PyGeoStudio 0.4.3
GeoStudio 2024.2
Python 3.11.7
Spyder 5.4.3
(1) 自动化数值模型处理过程
(2) 利用 Python 的强大功能来处理数据
(3) 使用matplotlib制作高质量绘图或导出到强大的后处理工具,如 Paraview 或 Tecplot
(4) 不使用图形界面修改分析
(5) 程序分析,例如相同的分析但具有不同的参数
(6) 执行自动校准、优化和敏感性分析
(7) 实时预测和预测模拟
import PyGeoStudio as pgs
src_file = "./Rapid drawdown.gsz"
geofile = pgs.GeoStudioFile(src_file)
geofile.showGeometries()
geofile.showReinforcements()
geofile.showAnalysisTree()
analysis = geofile.getAnalysisByID(4)
timestepping = analysis["TimeIncrements"]
instant_drawdown = geofile.getAnalysisByName("3-SEEP Instantaneous drawdown")
results = instant_drawdown["Results"]
(2) 参数化计算不同的渗透系数
mat = geofile.getMaterialByName("Embankment")
Kfunction = mat["Hydraulic"]["KFn"]
Ksats = [1e-7, 1e-6, 1e-5, 1e-4]
Ts = []
PWPs = []
for new_Ksat in Ksats:
actual_relK = Kfunction.getYData()
actual_Ksat = Kfunction.getYData()[0]
Kfunction.setYData(new_Ksat/actual_Ksat * actual_relK)
geofile.save()
instant_drawdown = geofile.getAnalysisByName("3-SEEP Instantaneous drawdown")
pgs.run(geofile, analyses_to_solve=[instant_drawdown])
T,PWP = instant_drawdown["Results"].getVariablesVsTime(
"PoreWaterPressure",
locations=[[25,2]]
)
Ts.append(T)
PWPs.append(PWP)