geostudio-seep模块

中仿科技有限公司主要 内容 财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容 数值模拟步骤SEEP/W分析基本原理几何和网格材料类型和属性边界条件分析类型数值问题结果解译数值模拟步骤1、运行软件，选择分析模块；2、模型分析设定：分析类型、分析控制、收敛及时步设定；3、绘制或导入几何模型（几何模型修正）；4、定义材料类型，输入材料参数；5、将材料赋给几何模型，网格查看与调整；6、定义边界条件，将边界条件赋给几何模型；7、检查模型，求解；8、结果查看；9、工况变更，再次求解。总结：数值模拟分析要素：几何、材料、网格、边界条件SEEP/W分析基本原理有效的利用SEEP/W软件来模拟分析问题，需要对一些关键的基本原理有清楚的认识Darcy’slaw总水头定义土的储水能力基质吸力或负孔隙水压力对于土体渗透性（渗透系数或渗透率）的影响Darcy’slaw在非饱和区的应用同饱和区是一样的q=流量k=渗透系数(材料属性)i=总水头梯度*SEEP/W的公式用总水头来 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 达总水头(H)=压力水头(h)+标高(y)*地表：标高yh*等式右端项给出土体储水变化率瞬态分析控制方程等式右端项=0即稳态分析*几何和网格GeoStudio中的几何对象区域（Region）、点和线自由点自由线依附于线的接触面单元面层（SurfaceLayer）圆孔网格生成结构化网格非结构化的四边形和三角形网格生成网格的一般原则区域（Region）特性当一个区域被定义时，有如下特性：只能赋予一种材料类型；只能有一种网格模式一个区域内所有单元阶数相同，同为一阶或者二阶单元同样的积分阶次接触面（Interfaceelementsonlines）将材料本构赋给线：SEEP/W中接触面模型可以被用来代表如截流幕墙或渗透通道，给几何线赋予渗透性材料。面层（SurfaceLayer）在地面条件改变以响应气候和气候条件在短时间内的剧烈改变，比如，在极短的时间内土体由干燥变饱和。为了数值上能够处理这种边界条件快速和剧烈的变化，在地表必须有好的网格离散。只能在地表定义创建与编辑几何对象GeoStudio软件的图形用户化界面导入CAD格式（.dxf）的文件精确点的定义(keyinpoint):输入坐标精确定义点的位置区域的修改与编辑（Modify、keyinregion）DXF文件导入GeoStudio软件的过程及技巧导入前的准备：设定单位和比例，红框内是重点注意的，这个决定了建模的尺度，如下图示，x向最大48.34m，y向最大35.64m，最小均为-5m，首先根据自己的模型尺度设定相应的值。cad的要求：要求导入的图形须是闭合的pl线，通常我们在cad作图是不会注意这一点，有一个好的方法是，即使不是闭合的pl线，用cad绘图菜单的边界命令可以轻松生成闭合的pl线（闭合region）实例操作模型修改精确点的捕捉模型检验网格生成GeoStudio的网格生成是全自动的，当材料赋给区域（Region）时网格自动生成，用户可以对生成的网格做用户化的设定与修改。优点：用户无需担心区域连接处的网格兼容性，也无需担心改变网格造成材料属性或边界条件的丢失。用户可以对基于区域（面）、点、线对象进行网格设定和调整，比如指定网格密度用实际单元长度，或全局网格比率，或者沿一条线等分等方法总的来说：推荐用户用默认的全局网格密度，如果需要，对局部进行网格加密或者粗糙化。结构化网格：需要更多的工作即创建多个区域，但不会对计算结果产生显著的改善混合网格（默认网格），推荐的网格生成网格的一般原则同数值模拟一样，网格剖分也是锻炼才能获得的技能。理想的网格的生成亦需要不断的实践和积累经验。经验会让你认识到网格与解之间的联系。当你能够预见解的近似分布时，对于网格剖分就会更精通。尽管没有硬性规定和快速准则，但还是有些是常有益的广泛准则：在开始分析时用尽可能少的单元，用不超过1000个单元计算，以验证概念并获得第一个近似解。当缩放比例是100%时，所有单元都是可见的。材料类型和属性创建材料类型及其属性，将材料赋给几何SEEP/W有两种主要的材料属性：土水特征曲线（VolumetricWaterContent）函数(function)土体渗透系数（Conductivity）饱和渗透系数saturatedconductivity(Ksat)函数（function）非饱和区渗透性是变化的，渗透系数亦非常数Ksat仅在饱和区应用的渗透系数当越多的空气充填孔隙，渗流变得越困难*材料属性：非饱和区渗透性AirEntryValueVolumetricWaterContentPore-WaterPressure(pwp)(-pwp)ResidualWaterContent0*材料属性：土水特征曲线（体积含水量函数VWC）如果土体会经历负孔隙水压力，那么渗透系数就非常数.渗透系数是含水量的函数间接的渗透系数是孔隙水压力(负)的函数渗透系数函数的基本形状由体积含水量函数派生而出*材料属性：渗透系数函数KsatAirEntryValueConductivity(logscale)0Pore-WaterPressure(+pwp)(-pwp)*渗透系数函数Whatishappeningassuctionincreases?-pwp,(kPa)Kq-pwp(kPa)SaturatedS=1FullyDrainedSrUnsaturated1>S>SrAir-EntryValueVWC*材料属性利用VWC样本函数EstimateVWCSamplefunctions*在SEEP/W中获得渗透性函数的方法之一EstimateVWCSamplefunctionsKfunctionEstimate*在SEEP/W中获得渗透性函数的方法之一利用VWC样本函数Suction(kPa)1S1001021041060SandSiltTillClay*用粒径信息估计VWC函数EstimateVWCGrain-sizeKfunctionEstimate*在SEEP/W中获得渗透性函数的方法之二在野外或实验室实测VWC函数MeasureVWCDatapointKfunctionEstimate*在SEEP/W中获得渗透性函数的方法之三Closedformoption1(FredlundandXing,1994)在SEEP/W中获得渗透性函数的方法之四Closedformoption1(VanGenuchten)参数拟合 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 材料属性：体积压缩系数（）体积压缩系数是体积含水量函数在孔隙水压力为正的部分的曲线斜率，其物理意义是描述在给定孔压改变的膨胀或收缩性能在大多数非饱和土体渗流分析中体积压缩系数不需准确定义，用户可以定义一个小的值，比如1.0e-5/kPa。如果是流固耦合分析，则这一参数需仔细输入。接触面模型及参数可以定义线的切向和法向的渗透系数，比如模拟排水板排水或者截流幕墙等在非饱和区（负孔隙水压力）和饱和区（孔压为正）都会发生流动渗透系数是基质吸力或者负孔隙水压力的函数体积含水量是基质吸力或者负孔隙水压力的函数*材料模型总结：**边界条件基本原则边界条件作用对象水头边界条件指定边界流量源和汇浸润面自由排水远场边界条件边界条件函数边界条件的基本原则指定边界条件是数值分析非常关键的部分！！！数值分析的解直接响应边界条件，没有边界条件就不可能获得解。边界条件本质上是驱动力。要点一：SEEP/W的边界条件只能是二选一，要么是水头要么是流量。要点二：需要完全理解当指定水头时，有限单元的解为流量，反之亦然。要点三：既没有指定流量亦未指定水头的结点，计算流量为0.边界条件作用对象在GeoStudio2007中所有边界条件必须直接作用在几何上，比如区域面、区域线、自由线和自由点。不可以作用在单元边或者节点上。水头边界条件总水头坝体水头边界条件之总水头坝体水头边界条件之总水头、压力水头指定边界流量Appliedunitflux(q)q=3e-5m/s1m2mQ在SEEP/W中，流入系统为正，流出系统为负。但结点流量Q为标量源和汇（Sourcesandsinks）有另外一类边界条件叫做源或者汇。这类边界条件也被成为第三类边界条件。典型的汇代表网格内部某点排水，源则可能是个补给井。关于源和汇一个重要的概念是它们代表了水流入或流出系统。在本例中指定的水头边界（H）代表汇（sink）水从该结点流出系统在网格内部，指定结点水头（H）或者流量（Q）常常代表源或者汇这种边界条件*浸润面（Seepagefaces）边界条件：PotentialSeepageFaceReview意义：流量与水头转化边界，通过迭代求解得到溢出点位置自由排水（Unitgradient）自由排水在某些情况下是很有用的边界条件，但也容易被用错。Kisch问题：粘土覆盖层下覆粗糙的沙土（aclaybarrieroveracoarsesandysoil）Kischpressureprofilewithheadtopandheadbaseboundarycondition远场边界（Farfieldboundaryconditions）在很多情况下，远离感兴趣点的边界条件是难以定义的，比如：远离抽水井或者开挖区域的边界条件，这被称为远场边界条件。处理方法之一：延展问题的几何，缺点：需要跟多的单元和跟多的计算时间，且有些情况会造成网格剖分困难，结果难于解释。处理方法之二：使用无限单元（infiniteregions）该边界条件上定义水头边界条件将会获得好的计算结果，非零的流量边界条件不推荐使用。使用建议：先不使用无限单元，在获得合理结果的基础上再用该方法以评价是否提供更符合预期的解边界条件函数（Boundaryfunctions）在瞬态分析时边界条件经常是时间的函数水头-时间（Headversustime）水头-体积的函数（Headversusvolume）结点流量-时间（NodalfluxQversustime）在时间相关的函数中，需确认时间的单位同渗透系数的单位(比如seconds,hours,days等)相一致。对于长期的瞬态分析，推荐将渗透系数定义为m/day以便说有的时间单位都是天，这会避免用864000s来代表10天。分析类型（AnalysisTypes）稳态分析（SteadyState）瞬态分析（Transient）时步（Timestepping）分步/多步分析轴对称分析平面视图（俯视图Planview）稳态分析（SteadyState）假定边界条件不随时间变化。对给定的特定的边界条件，模拟得到分析对象的水压力分布和流动状态。稳态分析不考虑时间，渗流方程被简化，同时收敛将变得困难——这取决于土体属性（渗透性函数）的非线性程度。稳态分析可以不定义体积含水量函数VWC-F，因为该函数的物理意义是当压力改变时土体吸水或释水性能。只有唯一的收敛结果。思考：完全饱和的材料？瞬态分析（Transient）分析多长时间土体响应用户定义的边界条件。例如：预测蓄水过程坝体内水位随时间变化，预测强降水的土体内水压力变化，等等。关键一点：瞬态分析必须定义初始条件（初始的各结点的总水头），不真实的初始条件会导致难以理解的、不符合实际的计算结果。初始条件的获得：定义初始水位线（Drawinitialwatertable）稳态分析的计算结果瞬态分析指定时步的计算结果其他GeoStudio软件的计算结果（SIGMA/WQUAKE/W）空间孔隙水压力稳态分析的计算结果注意事项定义初始水位线注意事项当水位线时水平时，用水位线的方法会给出精确的初始孔隙水压力，如果水位线是曲线，则给出的是近似的初始孔压。近似能否被接受取决于问题本身。空间孔隙水压力定义空间压力水头分布，未定义的区域软件插值得到。时步（Timestepping）瞬态分析必须对时步进行设定，合理的时步是问题求解的关键，准确的计算结果依赖于时步设定。时步设定基本规则:有限单元形状是重要的，三角形单元不应该有大约90度的内角，四边形单元可以有双倍于三角形单元的时间步长；由于单元尺寸与时间步长成正比，双倍的单元尺寸一位置可以设定双倍的时间步长，由此可知：减小单元尺寸却不减小时间步长将不会改善计算结果；时步设定基本规则:如果有瞬间荷载作用于系统，比如由于外部荷载导致的突然的孔隙水压力增加，应当将初始时步设定成近似于荷载作用土体的饱和渗透系数除以Mv值大小的数量级；如果模拟水库快速蓄水，可以使用水头-时间边界函数来模拟短期快速蓄水；如果在地表应用了流量-时间边界条件，且应用了许多不同值（比如模拟自然降水开/关），此时需要确保计算时步足够小以便允许求解器跟得上函数的性状。时步设定提示:在开始分析时，取总的分析时间的分5个相等的时步，运行模型，如果收敛观察需要多长时间达到收敛。如果不收敛，应该怀疑时步有问题，重新设定10个小的时步。模拟求解正常后，可以试验减小时步看收敛结果是否改变。唯一的方法是不断的试验分步/多步分析GeoStudio2007的新功能便是执行多步分析什么是多步分析？轴对称分析平面视图（俯视图Planview）（confinedaquiferonly）[K]includesconductivity(k){H}istheunknownButk=f(h)=f(unknownh)ThereforeiterativeprocedurerequiredStartwithKsat,thenuselastcomputedvalue-repeatedsubstitution*数值问题数值问题收敛判定不收敛收敛数值问题WaterbalanceerrorinatransientanalysisImprovingconvergenceEquationsolvers(directorparalleldirect)Automaticadaptivetimestepping结果解译流网（Flownets）渗透力（seepageforces）出溢比降（exitgradients）流网（Flownets）稳态的均质土体的渗流用拉普拉斯方程来描述，如下：图解这一方程就是我们熟知的流网，流网由等势线与流线交叉组成，在计算机和有限元数值方法被广泛应用之前，流网常常是唯一可用的求解渗流问题的方法。能够用流网来分析的问题其实是少数特别的情况。等势线近似的流网截面流量结构底部水压分布渗透力在所有的Geostudio软件中，当分析使用的是土体的总容重和边界力来代替浮力和渗透力，渗透力是不需考虑的。出溢比降（exitgradients）出溢比降和临界比降必须要考虑模型的几何、有效应力和土体剪切强度以及流速。计算所得的出溢比降为1并不意味着会发生破坏或管涌等问题。我们建议取溢出点位置合理深度范围内的比降平均值做参考，而非出溢点位置的值。掌握计算原理，清楚物理概念，更多的则是实践！Morepractice!