地应力平衡的一个简单例子.

地应力平衡方法熊志勇陈功奇第一部分地应力平衡方法简介地应力平衡有三种方法:(1*initialconditions,type=stress,input=FileName.csv(或inp该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定。其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。(2*initialconditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。计算速度快,收敛性好。缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。(3*initialconditions,type=stress,geostatic,user该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。第二部分地应力平衡方法实例详解地应力平衡是岩土工程数值模拟分析的重要的内容,为了让师弟师妹们快点上手,我利用第一种方法做一个较简单的模型,希望对大家有用。一、模型描述:3/2080mkg=ρ7101⨯=E35.0=ν二、地应力平衡过程1.启动ABAQUS,单击CreateModelDatebase2.创建部件(Part在Part模块,单击创建部件按钮,弹出如右图的对话框,按图输入部件名:Part-soil;采用二维模型选择2DPlanar;Type选择可变型(Deformable;基本特征选择壳体(Shell;Approximatesize输入70,这个数值的大小,应根据模型的最大尺寸来确定:稍大于最大尺寸的2倍。比如本模型最大尺寸是30那么我输入了70,但也不是绝对,你当然也可以输入65,或75等等。最后单击Continue,继续下一步。按照模型尺寸(如图,建立模型部件,双击鼠标中键,完成部件的建立。60kpa203.建立材料属性(Property在Module中切换到Property模块,单击,输入材料名称(name:Material-soil,单击Density,在弹出对话框中输入:密度2080;然后单击Mechaniacal—Elasticity—Elastic,在弹出的对话框输入图,单击OK完成材料的定义。单击按钮,输入名字:Section-soil,选择Soild,Homogeneous,单击Continue,OK,完成截面的创立。单击按钮,选取部件(单击或框选,选择后成粉红色,表示选中,单击Done或单击鼠标中键来确定。在弹出的对话框中选中Section-Soil,单击OK.4.装配部件(Assembly在Module选择Assembly模块,单击,弹出对话框,采用默认值,单击OK.注意:1.本模型只有一个部件所以自动选中,如有多个可按Shift键全选中;2.InstanceType本例都无所谓,只有一个部件,但是如有多个部件,我比较喜欢选择Independent,因为在这种情况下,所有的部件会出现在一个窗口,划分网格更方便,并且直观的看到不同部件连接处网格划分的是否协调。5.创建分析步(Step在Module选择Step模块,单击,弹出下面的对话框,创建分析步Step-1,选择Geostatic,单击Continue,在分析步编辑框Basic中选择On(如果会发生大变形的情况下要选择On,其他默认,单击OK6.施加荷载和定义边界条件(Load在Module选择Load模块,单击定义边界条件,需要对模型的左、右、下底面定义边界条件,分别命名BC-1、BC-2、BC-3,选择初始步(Initial,Displacement/Rotation,单击Continue,选择左、右、下底面边界,单击Done,分别选择U1、U1、U2(U1是水平方向,U2是竖直方向,单击OK完成边界条件的定义。单击定义重力荷载,Name:Load-grv,选择Step-1,Mechanical,Gravity,单击Continue。单击EditRegion,选择整个模型,在Component2中输入-9.8(重力加速度,单击OK7.划分网格(Mesh(划分网格是一门艺术,本例只是简单的划分在Module选择Mesh模块,单击设定网格的种子,将全局种子大小设为1,其余默认,单击OK.单击,采用默认,单击OK。单击,Family选择PlaneStrain,其余采用默认,单击OK单击,单击Yes,完成网格的划分,如下图:8.在命令行中输入mdb.models['Model'].setValues(noPartsInputFile=ON,按回车键9.在job模块中创建名为Job-NoInitialCondition的分析步,提交分析。应力图水平位移图10.将分析得到的应力场保存为一个文本文件。具体做法:打开分析得到的ODB文件,选择菜单Report→FieldOutput,在下图所示的对话框中,选中积分点上的各个应力分量(对于二维问题,应力分量S11、S22、S33和S12;对于三维问题,还应选中S13和S23。图表1输入常变量S11、S22、S33和S12单击此对话框中的Setup标签页,在Name文本框中输入要保存的文件名b.inp,取消对Appendtofile的选择(即创建一个新文件,在Write后面只选中FieldOutput(如下图所示。图表2设置输出场变量注意,此处输出的是当前增量步结束时的应力结果,因此上述对话框顶部的Step必须是Geostatic分析步,Frame必须是1。如果Frame是0,会看到输出的应力都是0。11.按照ABAQUS所要求的初始应力场文件格式,修改上述文件b.inp中的内容。具体方法为:用Excel打开上述文件b.inp,在”文本文件导入向导”的步骤1中选择“分隔符号”,在步骤2中选择“Tab”键和“空格”键,这样b.inp中的各列数据就成为Excel表格中的各个列。删除表格中开始几行的模型信息,再删除积分点编号所在的第2列数据(都为数字1,只保留单元编号和各个应力分量列,并将各个应力分量的科学计数法格式改为显示小数点后5位数字。修改前和修改后的数据如下:图表3修改前图表4修改后下面将上述数据输出为以逗号分隔的文本文件b.csv,具体的方法是:在Excel中单击菜单“文件”→“另存为”,将文件类型设置为“CSV(逗号分隔”,对于出现的提示信息,单击“是”,即可。12.为模型中定义初始应力场。在ABAQUS/CAE中无法直接定义初始应力,只能手工添加关键词,具体方法为:将原来的CAE模型另存为a.cae,选择菜单Model→Editkeywords,在*STEP语句之前添加以下语句:*initialconditions,type=stress,input=b.dat图表5修改前图表6修改之后13.在job功能模块中将分析作业名Job-WithInitialCondition,重新提交分析。注意,初始应力场文件b.csv应该和INP文件Job-WithInitialCondition.inp位于同一个路径下,否则将会出现下列错误信息:Thefollowingfile(scouldnotbelocated:b.csv(无法找到文件b.csv14.查看地应力平衡的结果。打开Job-WithInitialCondition.odb初始状态下（0时刻），模型就具有了一个初始应力场，这个应力场与上一个结果文件Job-NoitialCondition.odb中分析步结束时刻的应力场完全相同。图表7平衡后的地应力场15.上面就已经完成了初始地应力平衡，接下来可以添加其他分析步（例如普通的静力分析步Static，General），定义接触和实际的荷载，并去掉前面第一步中临时边界条件图表8施加静力的模型图表9应力云图图表10位移等值线图