ABAQUS模拟渗流要点

ABAQUS模拟渗流要点(2011-05-2422:56:10)转载标签：杂谈从其他网站看到一些这样的介绍，希望对饱和土方面模拟的朋友有帮助：当只进行渗流计算时：1、由于Abaqus中缺乏非耦合的孔压单元，这时可采用耦合单元，但要约束住所有位移的自由度。2、渗流材料参数选择。在CAE中都是在(Material-creat-other-porefluid)选项中。(1)Gel：定义凝胶微粒吸湿膨胀的发育过程，这在一般的岩土分析中应用不多。(2)Moistureswelling：定义由于吸湿饱和所引起的固体骨架体积膨胀(或负吸力引起的骨架收缩)。(3)Permeability：定义饱和介质的渗透系数，该渗透系数可以在type选项中定义为各向同性、正交各向异性和各向异性，并且可以根据VoidRatio定义为孔隙比的函数。在Suboptions中选择SaturationDependent参数来指定与饱和度相关性系数ks(s)，缺省设置为ks＝s3，而非饱和介质渗透系数k’=ksk。选择Velocitydependence参数可以激活Forchheimer定律，缺省的是Darcy定律。(4)PoreFluidExpansion：定义固体颗粒与流体体积热变化效应。(5)PorousBulkModuli：定义固体颗粒与流体体积模量。(6)Sorption：定义负孔隙压力与饱和度之间的相关性。当type=Absorption时，定义吸湿曲线，type=Exsorption时定义排水曲线。3、载荷及边界条件(1)通过(Load-creat-step-fluid-surfaceporefluid)选项定义沿着单元表面的外法线方向的渗流速度vn，当考虑降雨影响时可采用此载荷(2)边界条件(Boundarycondition-creat-other-porepressure)选项定义孔压边界条件，此时要先假定浸润面的位置，然后定义浸润面上的孔压为零，Abaqus会在后续的分析计算中自动计算出浸润面的位置。Abaqus默认的是不透水边界。(3)当渗流自由面遇到临空的自由排水面时，需要定义一个特殊的边界条件。此时可以通过在inp文件中加入*Flow或*Sflow来定义(4)初始条件的定义。初始条件中一般要定义以下几种：*initialcondition,type=saturation初始饱和度*initialcondition,type=porepressure初始孔压*initialcondition,type=ratio初始孔隙比当进行耦合分析时，基本步骤同上，但要去掉除边界条件之外的约束，同时还要在边界上加上流体压力。***这些问题都是很有意义的问题和很关键的技术，一下子讨论和理解的人也不是很多，如能给出一些例子来进一步阐述和说明（越详细越好），我想可以作为精华帖。***likuiming编辑上午只看了一个damseep.inp文件，大坝渗流。只是了解其中每个命令的含义，没有运行看其结果。*MATERIAL,NAME=soil*DENSITY2000.0,0.0*ELASTIC,TYPE=ISOTROPIC1000.0E9,0.2,0.0*****E,r,温度**********PERMEABILITY，TYPE=ISOTROPIC（默认），ANISOTROPIC，SATURATION，（饱和）VELOCITY（速度）******PERMEABILITY,SPECIFIC=10000.0***流体比重specificweight（FL–3）**2.11E-4,***渗透系数K，空隙率e,温度****SORPTION**定义材料吸附作用TYPE=ABSORPTION(default)**-100000000.0,0.04**孔隙压力（负压），饱和度s（0.01=s=1.0）**-10000000.0,0.050.0,1.0*INITIALCONDITIONS,TYPE=SATURATION****定义初始饱和度，默认为1**ALLNODE,1.0*****节点集或节点号，饱和度（默认为1）*******定义初始孔隙水压力（initialporefluidpressures）foracoupledporefluiddiffusion/stressanalysis******INITIALCONDITIONS,TYPE=POREPRESSUREPOR,196000.,0.0,0.0,20.**节点集，流体第一孔隙压力，垂直向第一孔隙压力，流体第二孔隙压力，垂直向第二孔隙压力**POR11,0.0POR12,0.POR13,0POR14,0.POR15,0.*INITIALCONDITIONS,TYPE=RATIO***定义初始孔隙比**ALLNODE,1.0***节点，第一孔隙比，垂向，第二，垂向***BOUNDARYALLNODE,1ALLNODE,2*RESTART,WRITE*STEP,INC=100,UNSYMM=YES*SOILS,UTOL=0.2**调用自动时间步**0.2,1.,1.0E-5,0.3**初始时间增量，计算总时间，最小，最大时间增量，孔隙压力随时间变化率***DLOADDAM,GRAV,9.8,0.0,-1.,0.*BOUNDARY66,8,,0.067,8,,19600.68,8,,39200.69,8,,58800.70,8,,78400.71,8,,98000.72,8,,117600.73,8,,137200.74,8,,156800.75,8,,176400.6,8,,196000.26,8,,58800.027,8,,39200.28,8,,19600.29,8,,0.0*FLOW**定义渗透系数和孔隙压力？？**FSIDE,Q2D,0.01**单元集，渗流类型，起始孔隙压力，起始渗透系数***CONTROLS,ANALYSIS=DISCONTINUOUS*RESTART,WRITE,FREQUENCY=1000*ENDSTEP书中说了shift键选择多个项目，比如线段，我发现在无意中多选了一个项目后，可以用Ctrl键取消