《地应力平衡地一个简单例子.docx》由会员分享，可在线阅读，更多相关《地应力平衡地一个简单例子.docx（20页珍藏版）》请在三一办公上搜索。
1、abaqus地应力平衡先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际 边坡尺寸一致，但我们可以夸X-点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似 目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果 能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进展计算，变形后形状和现状边坡形状一致，其内力也就是 初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模
2、型就是边坡现在 的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。如果我们知道现状边坡的内力， 将其提取出来作为几何模型的内力，再和外力重力平衡，如此我们建立的模型才能算和实际模型一致。真实地知道现状边坡的内力是很难的，我们采取的方法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进展计 算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何 模型，再施加重力当然边界条件也应根本一致以平衡，这样才算建立了与现状模型根本一致的模型，其下的计算才成为可能。这就 是所谓“地应力平衡的含
3、义、目的、作用。2.地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内 力，仅有外力重力，没有内力是不可能的，同样，仅有内力专指初始应力场而不受重力也是不可能的，否如此，整个体系的力不会 平衡。这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。为的是内力和外力平衡。3.地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场的方向是不是向上呢，这同样是我开始接触 abaqus的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪 力之分，其方向也按是拉是压是顺时针或逆时针而分
4、，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，如此 压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？aba中初始地应力场平衡一般在外表水平的情况下仅仅和密度相，藩度一样的话平衡的结果很好，别的参数改变之后经过计算， 差异很小。外表不水平的情况如此最好通过文件导入初应力的情况进展平衡。一般来讲，外表不平的时候有很多因素造成误差很大。cdstudio 版主的总结：=908486&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D470地应力平衡的正确时间点：所谓地应力平衡是指,当我们建任何东西或挖任何东
5、西之前，地表的位移都是零，但是土体的应 力却存在，这种无位移但有应力的时间点叫地应力平衡.那么正确的时间点应该是什么呢，就是不管土原来是什么样的，高山，河流， 丘陵，平原都可以，关键是在我们对它做任何扰动之前.有些人问桩土作用如何平衡？这个问题首先就是错的，因为地应力平衡的时候 是不能有桩的.一旦有桩就跨过了地应力平衡这个时间点，数值模拟要忠实于实际情况，不能想当然，所以正确的选择时间点是地应 力平衡计算中最最关键的.处理地应力平衡的三种方法:1.地表水平土层分层水平的情况下的地应力平衡地表水平而且土层 分层水平的情况下的地应力平衡是最简单的情况，事实上也是大多数计算用到的平衡方法，这个时候可
6、以用ABAQUS提供的*initial conditions， type=stress，geostatic方法来做.具体怎么写这个语句，请自行参考文档.2.其它的情况下的地应力平衡地表不水平或 土分层不规如此的情况下只能使用导入地应力数据的方式进展地应力平衡，其方法的核心是给土体加重力，不必理会其土体位移，只是 单纯的提取每个单元的平均应力，再将这个应力读回ABAQUS作为初始应力，以达到平衡土体位移的效果.具体操作参考这个贴 forum.simwe./viewthread.php?tid=762963&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2 3.不预
7、设地应力所谓的不进展地应力平衡是 指计算前不设定土体位移归零，而是记录下土体在重力作用下每个结点的位移值，在进展后继计算后，将位移结果减去重力作用下每个 结点的位移值就能得到当前计算步的净位移，这种方法只适用于小变形分析，不适用于几何非线性的情况，好处是收敛性好，精度高， 因为记录的是结点位移，缺点是麻烦，因为要提前记录，还要处理数据.1法是对单元集进展赋值，与结点集无关;常见的问题1.有结构物的情况如何平衡?有结构物的情况下就不是正确的平衡点，如果一定要这样平衡，比如原来有一个隧道，现在在上面建房子， 要估计这种情况下的追加位移是多少，这种情况下可以用上面的方法2和3来完成.2.多层土的情况
8、下的平衡?多层土的情况可以根据不 同的情况选择方法1或2进展.3.有地下水的情况下的平衡?有水的情况下最主要是要搞清楚土的三相，按照不同的孔隙率计算土的干密 度和湿密度，这种平衡目前最好是只用方法1来完成.4.边坡加重力无法平衡?边坡加重力无法平衡表示边坡本身是不稳的，请检查模型 是否合理，塑性参数是否合理.aba_aba大哥解答：word你的inp文件中有part和instance,而geo.dat中的节点编号前面没有instance 名称。有两种方法：1。建立没有没有instance 名称。见.simwe./forum/viewthread.php?tid=7089422。建立有加上 in
9、stance 名称。见先建立好你的模型，按照你的课题把边界条件、接触什么的都建好，然后使用命令mdb.models你的模型名字 .setValues（noPartsInputFile=ON），建立没有 part 的 inp 文件abaqus 不认带 part 的 istress 数据zhuyi:要写你的 model 的名字 在CAE窗口顶部，Model右面，默认名称为Model-1, Model-2等等，不是CAE或JOB的名称。你的mode l的名字前后都要加上单引 号地应力平衡方法熊志勇陈功奇第一局部地应力平衡方法简介地应力平衡有三种方法：1*initial conditions,type
10、=stress,input=FileName.csv（或 inp）该方法中的文件FILENAME.INP获取方法为:首先将边界条件施加到模型上进展正演 计算，然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量（也 可以导出积分点处的应力分量，视要求平衡的准确程度而定）。其所采用的几何模型可以 考虑地表起伏不平的情况以与岩土材料极其不均匀的情况,适用X围广。但由于外插的 应力有一定误差,因此采用弹塑性本构模型时，可能会导致某些点的高斯点应力位于屈 服面以外，当大面积的高斯点上的应力超出屈服面之后，应力转移要通过大量的迭代才 能完成,而且有可能出现解不收敛的情况。在仅考虑自重情况
11、下只能考虑受泊松比的影 响带来的侧压力系数效应，因此平衡后的效果不一定很理想，但无疑其适用性很强。2 *initial conditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力与其相应坐标。所采用的 几何模型一般较规如此,外表大致水平，地应力平衡的好坏一般只受岩体密度的影响，无 论采用弹性或弹塑性本构模型都能很好的达到平衡，可以不必局限于仅受泊松比的影响, 能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。计算速度快，收敛性好。 缺点就是不能够很好平衡具有起伏外表的几何模型,需知道平整后模型的上覆岩体自 重。3*initial c
12、onditions,type=stress,geostatic,user该方法采用用户子程序SIGINI来定义初始应力场，可以定义其为应力分量为坐标、 单元号、积分点号等变量的函数,要达到准确平衡需具体边界条件,在实际中应用较少。60kpa(0,20)(30,20)地基(0,0)(30,0)12川12m 一 |Model: Model-1 i PartModule: Part(XV2|:第二局部地应力平衡方法实例详解地应力平衡是岩土工程数值模拟分析的重要的内容，为了让师弟师妹们快点上手, 我利用第一种方法做一个较简单的模型，希望对大家有用。一、模型描述：P = 2080kg / m3E =1x
13、107v = 0.3530m二、地应力平衡过程1. 启动 ABAQUS,单击 Create Model Datebase2. 创建部件Part在Part模块，单击创建部件按钮，弹 出如右图的对话框，按图输入部件名：Part-soil； 采用二维模型选择2D Planar； Type 选择可变型Deformable； 根本特征选择壳体Shell； Approximate size输入70，这个数值 的大小，应根据模型的最大尺寸来确 定：稍大于最大尺寸的2倍。比如本 模型最大尺寸是30那么我输入了 70， 但也不是绝对，你当然也可以输入65, 或75等等。最后单击Continue,继续下一步。 按
14、照模型尺寸如图，建立模型部件, 双击鼠标中键，完成部件的建立。:机囹斑 RE1 U i B 1 i b ajJ： : S O : ，.，= i ，-:wE* -!日 ModeliB Pm P-n -solilej安 beetu cMsia.i：ftiI. BlfadvL1a hiix bwn jTMtsd在Module中切换到P roperty模块，单击 *3. 建立材料属性P roperty，输入材料名称name：Material-soil,单击Density，在弹出对话框中输入：密度2080 ；然后单击 Mechaniacal一Elasticity一Elastic,在弹出的对话框输入图，单
15、击OK完成材料的 定义。Edit MaterialName: Material-soiName Material-soilDescription:Material BehaviorsIDescriptinn:Materia I SehfiftiicnEDensityGeneral Mechanical Thermal OtherDensityGereral Medianical Thermal OtherDepvairRegularizationUser MaterialUser Defined FieldUser Output VariablesDenuityElasticityElasti
16、cE Use terNumber cDataPlasticityHyperelastirDamage for Ductilt MetalsHyperfc-dmDa m ag,e for Trad on Sep arati m Laws bHypoelasticDa m a叩 for E ber-Reinforced Compasites kpornus Ela rticDefonnatiori PlasticityViscoelastcDamping ExpansionBrittle CrackingElssti cType: |lsotropic Use tern perature de p
17、endernt dataGulboptici nzsNumber of field vsriabloT-Long-termModuli time scale foir visco el a EticityJ s;B No compiressiomI I No tern si o nOKDataVcung Mcdul UPoi&nac Ratio,le70.2SCancel单击按钮一，输入名字：S ection-soil,选择Soild , Homogeneous，单击 Continue, OK，完成截面的创立。单击按钮，选取部件单击或框选，选择后成粉红色，表示选中，单击）one或单击鼠标中键来
18、确定。在弹出的对话框中选中Section-Soil，单击OK.Edit Section AssignmentSectionSe-cti on; Sr-dti on -scilCircate.Note: Lkt contains cnly EBctions日 ppli 匚abJm to the selectedeg io ns.Type:Go-lidj HomogeneoueMateri a I: Materi a I -soilRegionR.egicn: (Picked)Shell OffsetI Offset: IOKCancel4.装配部件Assembly在Module选择Assembl
19、y模块，单击，弹出对话框，采用默认值，单击OK.ParisIns-tsince8 Typed- Dependent mesh o-m part) Independent (me-sh on instance)Note: To change a Dependent instances mesh you must edit its parts mesh.Auto-otiset from oth&r instancesOKApplyCancel注意：1.本模型只有一个部件所以自动选中，如有多个可按Shift键全选中；2.Instance Type本例都无所谓，只有一部件，但是如有多个部件，我比拟喜欢
20、 选择Independent,因为在这种情况下，所有的部件会出现在一个窗口，划分网 格更方便，并且直观的看到不同部件连接处网格划分的是否协调。5.创建分析步Step在Module选择Step模块，单击，弹出下面的对话框，创建分析步Step-1,选择Geostatic,单击Continue，在分析步编辑框Basic中选择On （如果会发生大变 形的情况下要选择On）,其他默认，单击OK6.施加荷载和定义边界条件Load在Module选择Load模块，单击定义边界条件，需要对模型的左、右、下底面定义边界条件，分别命名BC-1、BC-2、BC-3，选择初始步 Initial,Displacement
21、/Rotation,单击Continue，选择左、右、下底面边界, 单击Done，分别选择U1、U1、U2U1是水平方向，U2是竖直方向，单击OK完 成边界条件的定义。Create Bcjndsry Conditio单击一定义重力荷载，Name:Load-grv,选择 Step-1, Mechanical, Gravity, 单击Continue。单击Edit Region，选择整个模型，在ponent 2中输入-9.8重力 加速度，单击OKTypes for Sle-cted StepCategciry例| MechanicalThermal。Acaustic wtdb.modelsL Mo
22、del-1 .setValues(ncPsrtsInputFi 15=011) 9. 在job模块中创建名为Job-NoInitialCondition的分析步，提交分析。H-艮Hws|Wg: 7爽）IT T.&iLc!捋 由山皿-1. H14-bl T.NLe* -J.：=-Lc + DS R- 十w - .如十 w H.WM。山-l-d.Cd 14-bd应力图郴哼皿KDaDiDLDiIMDLIMgDL - - - a - B嘴器密归 _L Ik 1J Bu P_I4&.?! m I m2.d.c.B.L LLL Lz.12, -.1 T - - - - - a - d.d_一 1一一 一水
23、平位移图10. 将分析得到的应力场保存为一个文本文件。具体做法:打开分析得到的ODB文件， 选择菜单ReportTField Output，在如下图所示的对话框中，选中积分点上的各个 应力分量对于二维问题，应力分量S11、S22、S33和S12；对于三维问题，还应 选中S13和S23。|Sr Abaqus/CE Ve rsien &.7-1 - Model Debase: C:UsersxzyDocumentl?a q us.ca e 图 File Model Viewport View Result Hot Animate Report Ctptions TcField Output.Res
24、ultsModule: Visualization N ODB:eRe port Field Output图表1输入常变量S11、S22、S33和S12单击此对话框中的Setup标签页，在Name文本框中输入要保存的文件名b.inp,取 消对Append to file的选择即创建一新文件，在Write后面只选中Field Output 如如下图所示。图表2设置输出场变量注意，此处输出的是当前增量步完毕时的应力结果，因此上述对话框顶部的Step 必须是Geostatic分析步，Frame必须是1。如果Frame是0，会看到输出的应力都是0。11. 按照ABAQUS所要求的初始应力场文件格式，修
25、改上述文件b.inp中的内容。具体 方法为：用Excel打开上述文件b.inp，在文本文件导入向导的步骤1中选择“分 隔符号，在步骤2中选择“Tab 键和“空格键，这样b.inp中的各列数据就成 为Excel表格中的各个列。删除表格中开始几行的模型信息，再删除积分点编号所在的第2列数据都为 数字1，只保存单元编号和各个应力分量列，并将各个应力分量的科学计数法格 式改为显示小数点后5位数字。修改前和修改后的数据如下：1ln:era:ionpirrval-iesiromsourcesortedcoLumnEl cm ent Label.Outputatintegrationp-c-intsforp
26、art :FAIR Rm entF 一甲1S. S33F 一 f!1 MLabel(5Loc1Los1LocL 111-1. 54E+05-2. S7E+05-1. 54E+?E-7G6. 44B|1I21-7.B3E+04-1. 47E+05-7. BBE+?4-193. 5DZ31-7. 77E+04-1.46E+C5-7. B2E+D4377. 094111. 6)E+0S-3. OOE+05-1. 61E+D5288. 727b1-1. blH+Ub-2. B4E+05-1. 53E+D56Z0. 3Z3A1-1.43E+05-2. G2E+05-1 . 4-4F+1FRA4-. g7
27、1-1. 14EiOS-2. 13Ei C!3-1. 14Ei-4%. 60101-1.42E+05-2. 66E+C5-1. 43E+D5-202. 391-1.63E+05-3. 03E+05-1. 63E+D58. 14437_U1-S. 73E+U4-1. B0E+C5-9. 71E+047_9. luj11一1 一加+04-2. S6E+04-1.汛F+M-5i. 053B121-2. 03E+05-3. 78Ei C5-2. 03E+D5137. 831_31-5. 04E+03-9.12E+03-4. 95E+D3-14. iy141-1.7DE+05-3.1GE+05-1. 7
28、0E+?E-41. 75SL1-7.1DE+04-1. 32E+05-7. 09E+?4-17. 665B161-5.OBE+03-5. 50E+03-5. 10E+D3-17. 07B9171-9. 27E+04-1. 72E+05-9. 28E+D440. 6349101-5.13E+Q3-9. 56E+03-5. 1BE+D3-4b. OhioA1-F； 一 R 邛+GR-=i.Fi. *F+祚一舛一 A71 1图表3修改前11-1.54E+05-2.87E+05-1.54E+05-766. 44：82-7.83E+04-1. 4=7E+O5-7.88E+04-193.55233-7.7
29、7E+04-1.46E+05-7.S2E+04377.09444-1.60E+05-3.OOE+05-1.61E+05238. 727|155-1.54E+05-2.84E+05-1.53E+05620.32366-1.43E+05-2.68E+05-1.44E+05864.95977-1.14E+05-2.13E+05-1.14E+05-472. 601B8-1.42E+05-2.S6E+O5-1.43E+05-202. 399-1.63E+05-3.03E+05-1.63E+05B. 1443710r io-9.73E+04-1.80E+05-9.71E+04719. 10911u-1.5
30、3E+04-2.86E+04-1.54E+04-54- 05SS1212-2.03E+05-3.7SE+O5-2.03E+05137. B3113-5.04E+03-9.12E+03-4.95E+03-14. 19214r 14-1.70E+05-3.16E+05-1.70E+05-41. 75811515-7.10E+04-1.32E+05-7.09E+04-17. 66581616-5.06E+03-9.50E+03-5. 10E+03-17. 07891717-9.27E+04-1.72E+05-9.28E+0440. 634=918r is-5.13E+03-9.G6E+O3-5. 1
31、8E+03-45. 85161919-5.34E+03-9.90E+03-5.33E+03-34. 97112020-1.04E+05-1.93E+05-1.04E+05IB. 60&821-5.38E+03-9.90E+03-5.35E+03-9. 567522 Jr 22-1.37E+05-2.54E+05-1.37E+05-29. 34&G2323-1.47E+05-2.73E+05-1.47E+05-B5. 08&824-2.14E+05-3.98E+QS-2. 14E+0532. 46L4十耳明-9 1耳IT+门耳-? ClnF + ri-9 1 riF+i-i图表4修改后下面将上述
32、数据输出为以逗号分隔的文本文件b.csv,具体的方法是:在Excel 中单击菜单“文件一“另存为，将文件类型设置为“CSV逗号分隔，对于出 现的提示信息，单击“是，即可。12. 为模型中定义初始应力场。在ABAQUS/CAE中无法直接定义初始应力，只能手工添加 关键词，具体方法为：将原来的CAE模型另存为a.cae,选择菜单ModelTEdit keywords,在*，丁?语 句之前添加以下语句：。Abaqus/CAE Vension 6.7-1 - Model Database: C:Us图 FileModel ViewportManager.M&delSessionCopy Model C
33、ogy Objects.Ol0 n gnEdit Attributes k | Edit Keyvord& h BlsEluRenamekDeletekRainbowmView Result Plotlule: Visualization! 、图表5修改前图表6修改之后13. 在job功能模块中将分析作业名Job-WithInitialCondition,重新提交分析。注意, 初始应力场文件b.csv应该和INP文件Job-WithInitialCondition.inp位于同一个 路径下，否如此将会出现如下错误信息：The following file(s) could not be loc
34、ated :b.csv(无 法找到文件 b.csv)14, 查看地应力平衡的结果。打开Job-WithInitialCondition.odb初始状态下0时刻, 模型就具有了一个初始应力场，这个应力场与上一个结果文件Job-NoitialCondition.odb中分析步完毕时刻的应力场完全一样。图表7平衡后的地应力场15. 上面就已经完成了初始地应力平衡，接下来可以添加其他分析步例如普通的静力 分析步Static, General，定义接触和实际的荷载并去掉前面第一步中临时边 界条件图表8施加静力的模型E, H栖g; 75窗I+1.77i- + DS -+l.b224D& . - +1.4E+ 区 - +l.Dt-=4B& -+1.M3+ 谜 -+.ld-= + rd . - +.E：,J&= + W + b.7t-= + M +5.167 + M +.bll-= + M +2. DJE-4-M 1-十45匚叫:十匚门图表9应力云图Ll, M n iLudeiI 中曰.dEHliEDE4-1.7J t-02 4-2.E?-D2 4-Z.&15JS-DZ 十2,混4-1 E=tJ-D2 十Ld宁三土-DF 4-l.lllc-C.Z 4-7.-4 7JC-CQ_ 4- D. CiCiCl- 4- DD图表10位移等值线图
