地应力平衡的一个简单例子.
地应力平衡总结
地应力平衡1、地应力平衡好坏评判标准1)地应力平衡后，位移云图中最大位移达到10-6量级或更低（接近于0）。
（主要判别条件）2)地应力平衡后，应力云图中应力有一定的数值。
（也就是应力不为0，但变形接近于0）2、进行地应力平衡的原因总的来说，如果不进行地应力平衡，而只施加重力，模型会在重力作用下产生变形，而实际工程中，我们施加荷载时，重力产生的而变形已经产生，实际上得到的是附加应力产生的变形。
1)我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。
2)地应力平衡中的外力和内力的问题。
地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
二力平衡在生活中的应用例子
二力平衡在生活中的应用例子二力平衡是物体在力的作用下保持静止的状态，是力学中的一个基本原理。
在生活中，二力平衡原理有着广泛的应用。
下面将列举十个符合标题要求的例子。
1. 吊车起重在建筑工地上，吊车经常用于起重重物。
吊车的吊臂和起重物之间存在着二力平衡，吊车的重力和起重物的重力相互平衡，保证了吊车和起重物的稳定。
2. 游乐设施安全游乐设施如过山车、摩天轮等在运行时需要保持平衡。
游乐设施的底座和车厢之间存在二力平衡，保证了游乐设施在高速旋转时的稳定性和安全性。
3. 骑自行车骑自行车时，车手和自行车之间存在二力平衡。
车手通过脚蹬和把手施加力，与自行车的重力和阻力相平衡，保持自行车的平衡和前进。
4. 桥梁结构桥梁是承载车辆和行人的重要交通设施，其结构需要保持稳定。
桥梁的桥墩和桥面之间存在二力平衡，保证了桥梁的整体稳定性。
5. 摆钟摆钟是利用钟摆的摆动实现时间测量的装置。
摆钟的摆杆和钟摆之间存在二力平衡，保证了钟摆的稳定摆动，从而实现准确的时间测量。
6. 网球运动打网球时，球拍和网球之间存在二力平衡。
球员通过挥拍和球的反弹，与球的重力和空气阻力相平衡，实现球的准确击打和运动。
7. 静止的书本当我们把书本放在桌子上时，书本与桌子之间存在二力平衡。
书本的重力和桌子对书本施加的支持力相互平衡，使书本保持静止状态。
8. 走在平衡木上当我们走在平衡木上时，我们需要保持平衡。
我们的脚踩在平衡木上，通过调整身体的重心，与重力相平衡，保持身体的稳定。
9. 插电器的稳定当我们插电器插头时，插头和插座之间存在二力平衡。
插头受到重力和插座对插头的支持力相互平衡，保证插头的稳定插入和电器的正常使用。
10. 乘坐电梯乘坐电梯时，电梯和乘客之间存在二力平衡。
电梯的重力和乘客的重力相互平衡，保证了电梯的稳定运行和乘客的安全。
通过以上十个例子，我们可以看到二力平衡在生活中的广泛应用。
无论是建筑工程、运动设施还是日常生活中的插电器、电梯等，都需要保持力的平衡，以确保物体的稳定和安全。
[测试成功]地应力平衡的一个简单例子(1)
![[测试成功]地应力平衡的一个简单例子(1)](https://img.taocdn.com/s3/m/d7fdd43e83c4bb4cf7ecd165.png)
地应力平衡方法熊志勇陈功奇第一部分地应力平衡方法简介地应力平衡有三种方法：（1）*initial conditions,type=stress,input=FileName.csv(或inp)该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。
（2）*initial conditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快,收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。
（3）*initial conditions,type=stress,geostatic,user该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。
第二部分 地应力平衡方法实例详解地应力平衡是岩土工程数值模拟分析的重要的内容，为了让师弟师妹们快点上手，我利用第一种方法做一个较简单的模型，希望对大家有用。
地应力平衡的一个简单例子
abaqus地应力平衡先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。
2. 地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
3. 地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场的方向是不是向上呢，这同样是我开始接触abaqus 的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪力之分，其方向也按是拉是压是顺时针或逆时针而分，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，则压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？aba中初始地应力场平衡一般在表面水平的情况下仅仅和密度相关，密度一样的话平衡的结果很好，别的参数改变之后经过计算，差别很小。
地应力平衡汇总
1、“地应力平衡”的含义、目的、作用我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
2、地应力如何平衡地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
abaqus的part模块绘图功能不是很强，因此常用AutoCAD绘出平面图后导入abaqus。
在abaqus6.10中，主导入要分为以下几个步骤：1，在AutoCAD中建好模型的平面图形，并且另存为dxf格式。
完整的地应力平衡方法
地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC）第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8 第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON) (请严格按照这个格式，注意大小写的字母)第四步：提交J0B，完成后第五步：按以下步骤,Roport---Report Field Output---选中S11,S22,S33,S12,S13,S23---Name:XX.INP---Write中选择Field Output-------------ok!!!第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP，保存格式内容单元号 S11 S22 S33 S12 S13 S23 （请注意，在保存内容中没有这一行的）1 , . , . , . , . , . , .2 , . , . , . , . , . , .. , . , . , . , . , . , .. , . , . , . , . , . , .这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些数字。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上：*initial conditions,type=stress,input=xx.inp 完成第八步：重新提交JOB，OK第九步：如果你还没有成功的话，那我只能说----------------我无语了。
：）ABAQUS的这项功能确实很不错。
：）这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要ABAQUS的这项功能，反正很多都用得着。
二力平衡的适用范围
二力平衡的适用范围嘿，朋友们！咱今天来聊聊二力平衡这档子事儿。
你说这二力平衡啊，就好像一场拔河比赛。
两边的力量势均力敌，谁也没法把对方拉过去，就这么僵持在那儿。
咱生活中到处都是二力平衡的影子呢！你想想看，咱站在地上稳稳当当的，不就是重力和地面给咱的支持力在那平衡着嘛！要是这俩力不平衡了，那咱不就得要么飘起来要么陷地里去啦，那多搞笑呀！再比如说，你挂个灯笼在那儿。
灯笼受到重力往下拉，可绳子拉住它往上拽，这不就是二力平衡嘛！要是不平衡了，灯笼要么掉地上摔个稀巴烂，要么嗖地一下飞上天，找不到影喽！二力平衡适用的范围可广啦！大到天上飞的飞机，小到咱手里拿的一支笔。
飞机在天上飞，那空气的升力和飞机的重力就得平衡好，不然咋能稳稳地飞呢？笔在咱手里，重力和手给它的摩擦力也得平衡呀，不然笔不就滑下去啦。
你说要是没有二力平衡，这世界得乱成啥样啊？车开着开着突然就飞起来了，房子立着立着就歪了。
哎呀呀，不敢想不敢想！咱平时骑自行车也是啊，你蹬车往前的力和地面的阻力得差不多平衡，这样速度才能稳定呀。
要是不平衡了，要么嗖地一下冲出去，要么就原地不动干着急。
还有啊，咱走路的时候，脚给地面一个向后的力，地面就得给咱一个向前的反作用力，这不也是一种平衡嘛。
要是不平衡了，咱还怎么走路呀，不就跟喝醉酒似的东倒西歪啦。
你看那起重机吊起重重的东西，不也是靠起重臂的拉力和重物的重力平衡嘛。
要是不平衡，重物掉下来那可不是开玩笑的！二力平衡就像是生活中的一个小魔法，让一切都变得井井有条。
它让我们能稳稳地站在地上，让各种东西能正常地工作和运行。
咱可得好好珍惜这个小魔法呀，不然这世界还不知道会变成啥样呢！所以说呀，二力平衡真的是太重要啦，无处不在，不可或缺呀！。
平衡力的条件与示例
平衡力的条件与示例平衡力是物体在静止或匀速直线运动时所受到的力之间的平衡状态。
要使一个物体处于平衡状态，必须满足一定的条件。
本文将探讨平衡力的条件，并通过一些示例来加深理解。
平衡力的条件1. 力的合力为零在物体处于平衡状态时，所有作用于它的力的合力应为零。
也就是说，物体所受到的所有力的矢量和为零。
只有当力的合力为零时，物体才能保持静止或匀速直线运动。
2. 力的合力和力矩为零平衡力的第二个条件是力的合力和力矩为零。
力矩是力绕某个点旋转的趋势，它取决于力的大小和作用点到旋转中心的距离。
当物体处于平衡状态时，力的合力和力矩都应为零。
平衡力的示例1. 钢轨上的火车在铁路上行驶的火车可以被视为平衡力的一个示例。
火车所受到的重力、摩擦力、风阻力等力都必须平衡，否则火车将无法保持匀速直线行驶。
当火车处于平衡状态时，所有作用于它的力的合力和力矩都为零。
2. 悬挂的钟摆钟摆悬挂在支点处，并且可以自由摆动，使其达到静止或匀速运动。
当钟摆静止时，它的平衡状态要求重力和张力的合力为零，并且力矩也为零。
这种平衡条件使得钟摆保持稳定。
3. 木桌上的书本将一本书放在水平放置的木桌上也可以看作是平衡力的示例。
在此情况下，书本所受到的重力和桌面对书本的支持力必须平衡。
只有当这两个力的合力为零时，书本才能保持静止。
4. 游泳运动员的平衡游泳运动员在水中比较容易保持平衡。
当游泳运动员在水中踢腿和划手时，他们必须通过调整身体的姿势来保持平衡。
只有当游泳运动员身体的各个部分受力平衡，才能保持稳定的游泳姿态。
总结平衡力的条件是力的合力和力矩为零。
只有当物体所受力的合力为零，并且力矩也为零时，物体才能保持平衡状态。
通过一些示例，我们可以更好地理解平衡力的概念。
无论是火车在铁轨上行驶，还是钟摆摆动，平衡力都起到了关键的作用。
在日常生活中，我们可以观察和应用平衡力的条件，以更好地理解物体的稳定性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
地应力平衡方法熊志勇陈功奇第一部分地应力平衡方法简介地应力平衡有三种方法:(1 *initial conditions,type=stress,input=FileName.csv(或 inp该方法中的文件 FILENAME.INP 获取方法为 :首先将已知边界条件施加到模型上进行正演计算 , 然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量 (也可以导出积分点处的应力分量 , 视要求平衡的精确程度而定。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况 , 适用范围广。
但由于外插的应力有一定误差 , 因此采用弹塑性本构模型时 , 可能会导致某些点的高斯点应力位于屈服面以外 , 当大面积的高斯点上的应力超出屈服面之后 , 应力转移要通过大量的迭代才能完成 , 而且有可能出现解不收敛的情况。
在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应 , 因此平衡后的效果不一定很理想 , 但无疑其适用性很强。
(2 *initial conditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则 , 表面大致水平 , 地应力平衡的好坏一般只受岩体密度的影响 , 无论采用弹性或弹塑性本构模型都能很好的达到平衡 , 可以不必局限于仅受泊松比的影响 , 能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快 , 收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型 , 需知道平整后模型的上覆岩体自重。
(3 *initial conditions,type=stress,geostatic,user该方法采用用户子程序 SIGINI 来定义初始应力场 , 可以定义其为应力分量为坐标、单元号、积分点号等变量的函数 , 要达到精确平衡需已知具体边界条件 , 在实际中应用较少。
第二部分地应力平衡方法实例详解地应力平衡是岩土工程数值模拟分析的重要的内容,为了让师弟师妹们快点上手, 我利用第一种方法做一个较简单的模型,希望对大家有用。
一、模型描述:3/2080m kg =ρ7101⨯=E35. 0=ν二、地应力平衡过程1. 启动 ABAQUS, 单击 Create ModelDatebase2.创建部件(Part在 Part 模块,单击创建部件按钮,弹出如右图的对话框,按图输入部件名:Part-soil ; 采用二维模型选择 2D Planar; Type 选择可变型(Deformable ; 基本特征选择壳体(Shell ;Approximate size输入 70,这个数值的大小,应根据模型的最大尺寸来确定:稍大于最大尺寸的 2倍。
比如本模型最大尺寸是 30那么我输入了 70, 但也不是绝对, 你当然也可以输入 65, 或 75等等。
最后单击 Continue, 继续下一步。
按照模型尺寸 (如图 , 建立模型部件, 双击鼠标中键,完成部件的建立。
60kpa203. 建立材料属性(P roperty在 Module 中切换到 P roperty模块,单击 ,输入材料名称(name : Material-soil, 单击 Density, 在弹出对话框中输入:密度 2080;然后单击Mechaniacal — Elasticity — Elastic, 在弹出的对话框输入图,单击 OK 完成材料的定义。
单击按钮 ,输入名字:S ection-soil,选择 Soild , Homogeneous,单击 Continue , OK, 完成截面的创立。
单击按钮 ,选取部件(单击或框选,选择后成粉红色,表示选中,单击 Done 或单击鼠标中键来确定。
在弹出的对话框中选中 Section-Soil ,单击 OK.4. 装配部件(Assembly在 Module 选择 Assembly 模块 , 单击 ,弹出对话框,采用默认值,单击OK.注意:1. 本模型只有一个部件所以自动选中,如有多个可按 Shift 键全选中;2.Instance Type 本例都无所谓,只有一个部件,但是如有多个部件,我比较喜欢选择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 。
单击 Edit Region ,选择整个模型,在 Component 2中输入 -9.8(重力加速度,单击OK7. 划分网格(Mesh (划分网格是一门艺术,本例只是简单的划分在 Module 选择 Mesh 模块,单击设定网格的种子,将全局种子大小设为 1,其余默认,单击OK. 单击 ,采用默认,单击 OK 。
单击 , Family 选择 Plane Strain,其余采用默认,单击OK单击 ,单击 Yes, 完成网格的划分,如下图:8. 在命令行中输入 mdb.models['Model'].setValues(noPartsInputFile=ON,按回车键9. 在 job 模块中创建名为 Job-NoInitialCondition 的分析步,提交分析。
应力图水平位移图10. 将分析得到的应力场保存为一个文本文件。
具体做法 :打开分析得到的ODB 文件, 选择菜单Report → Fiel d Output,在下图所示的对话框中,选中积分点上的各个应力分量(对于二维问题,应力分量 S11、 S22、 S33和 S12;对于三维问题,还应选中 S13和 S23。
图表 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 所要求的初始应力场文件格式,修改上述文件 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 → Edit keywords,在 *STEP语句之前添加以下语句:*initial conditions,type=stress,input=b.dat图表 5 修改前图表 6 修改之后13. 在 job 功能模块中将分析作业名 Job-WithInitialCondition , 重新提交分析。
注意, 初始应力场文件 b.csv 应该和 INP 文件 Job-WithInitialCondition.inp 位于同一个路径下,否则将会出现下列错误信息:The following file(s could not be located :b.csv(无法找到文件 b.csv14. 查看地应力平衡的结果。
打开 Job-WithInitialCondition.odb 初始状态下（0 时刻），模型就具有了一个初始应力场，这个应力场与上一个结果文件 Job-NoitialCondition.odb 中分析步结束时刻的应力场完全相同。
图表 7 平衡后的地应力场 15. 上面就已经完成了初始地应力平衡，接下来可以添加其他分析步（例如普通的静力分析步 Static，General），定义接触和实际的荷载，并去掉前面第一步中临时边界条件图表 8 施加静力的模型图表 9 应力云图图表 10 位移等值线图。