abaqus怎么计算土体的地应力？ 请问abaqus建模时，采用修正剑桥模型建立的非水平地面怎么...
abaqus\u600e\u4e48\u8ba1\u7b97\u571f\u4f53\u7684\u5730\u5e94\u529b\u5206\u6790\u90e8\u53ef\u4ee5\u5b9e\u73b0\u81ea\u52a8\u5e73\u8861
\u8fd8\u6709\u53ef\u4ee5\u5148\u8ba1\u7b97\u51fa\u81ea\u91cd\u5e94\u529b\uff0c\u5728\u5bfc\u5165aba
\u653e\u5761\u6bb5\u548c\u5730\u5e94\u529b\u5e73\u8861\u6ca1\u6709\u5173\u7cfb
\u4f60\u7684\u6a21\u62df\u987a\u5e8f\u5e94\u8be5\u662f\u5730\u5e94\u529b\u5e73\u8861\uff0c\u7136\u540e\u751f\u6b7b\u5355\u5143\u6740\u6b7b\u653e\u5761\u6bb5\u7684\u90e8\u5206\u6bb5\u5143\u5427
\u81f3\u4e8e\u4e0d\u80fd\u8ba1\u7b97\u60c5\u51b5\u53ef\u4ee5\u7528\u5f39\u6027\u8ba1\u7b97\u91cd\u529b\uff0c\u7136\u540e\u5012\u56de\u53bb\uff0c\u6a21\u578b\u518d\u52a0\u4e0a\u5851\u6027
初始地应力平衡应该满足：
（1） 竖向位移为零或者很小，一般至少10的负6次方，甚至更小
（2） 竖向应力s22要等于土体深度*密度*g
关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!
1.先施加重力荷载的作用,可以在cae中实现;
2.在inp文件中的output request中写上
*el print
s,
这样就会将施加重力荷载后的应力输出到*.dat文件中了;
3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如
ELEMENT T FOOT- S11 S22 S33 S12
NOTE
1 1 -1.2598E+05 -1434. -3.1852E+04 892.7
2 1 -1.2249E+05 -6287. -3.2194E+04 1223.
3 1 -1.1795E+05 -497.7 -2.9611E+04 1664.
4 1 -1.1210E+05 -7240. -2.9834E+04 1992.
5 1 -1.0485E+05 579.0 -2.6068E+04 2600.
6 1 -9.5803E+04 -8272. -2.6019E+04 3031.
7 1 -8.4709E+04 1915. -2.0698E+04 4083.
8 1 -7.0634E+04 -9746. -2.0095E+04 4339.
9 1 -5.1088E+04 5401. -1.1422E+04 8519.
10 1 -2.4353E+04 -1.1150E+04 -8876. 1.2126E+04
11 1 -1.2847E+05 268.1 -3.2050E+04 738.1
12 1 -1.2786E+05 -9868. -3.4433E+04 629.1
13 1 -1.2938E+05 -4224. -3.3402E+04 502.5
14 1 -1.3039E+05 -3458. -3.3461E+04 165.9
单独存为一个*.dat文件,
4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为: instance名称.序号 ;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!
5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),
6.最后在inp文件的step之前写上
*initial conditions,type=stress,input=文件名.csv即可!
这种方法不需要用python,比较简单,希望能对大家有用!
先说为什么要施加地应力：
1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。这就是所谓“地应力平衡”的含义、目的、作用。
2. 地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。为的是内力和外力平衡。) q0 F3 q6 H1 O#
3. 地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场的方向是不是向上呢，这同样是我开始接触abaqus的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪力之分，其方向也按是拉是压是顺时针或逆时针而分，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，则压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？
不怕各位笑话，以上几点在高手看来是很简单的问题，却是我经历了漫长而艰辛的摸索才得到的，今天也写给初学者，不要再走我的老弯路了。
aba中初始地应力场平衡一般在表面水平的情况下仅仅和密度相关，密度一样的话平衡的结果很好，别的参数改变之后经过计算，差别很小。表面不水平的情况则最好通过文件导入初应力的情况进行平衡。一般来讲，表面不平的时候有很多因素造成误差很大。
地应力平衡方法：
第一步：建立模型，材料，分析步（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 （请注意，在保存内容中没有这一行的）
, . , . , . , . , . ,
, . , . , . , . , . , .
, . , . , . , . , . , .
, . , . , . , . , . , .
这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有!
上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2.
.......这些数字。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面
加上：*initial conditions,type=stress,input=xx.inp 完成,
第八步：重新提交JOB，OK
第九步：如果你还没有成功的话，那我只能说----------------我无语了。
ABAQUS的这项功能确实很不错。：）这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要ABAQUS的这项功能，反正很多都用得着。
希望我的建议能帮助初学者（其实我也是初学者，呵呵）很快掌握这个技巧。
地应力平衡的正确时间点
所谓地应力平衡是指, 当我们建任何东西或挖任何东西之前, 地表的位移都是零, 但是土体的应力却存在, 这种无位移但有应力的时间点叫地应力平衡. 那么正确的时间点应该是什么呢, 就是不管土原来是什么样的, 高山, 河流, 丘陵, 平原都可以, 关键是在我们对它做任何扰动之前. 有些人问桩土作用如何平衡? 这个问题首先就是错的, 因为地应力平衡的时候是不能有桩的. 一旦有桩就跨过了地应力平衡这个时间点, 数值模拟要忠实于实际情况, 不能想当然, 所以正确的选择时间点是地应力平衡计算中最最关键的.
处理地应力平衡的三种方法
1.地表水平土层分层水平的情况下的地应力平衡
地表水平而且土层分层水平的情况下的地应力平衡是最简单的情况, 事实上也是大多数计算用到的平衡方法, 这个时候可以用ABAQUS提供的*initial conditions, type=stress, geostatic方法来做. 具体怎么写这个语句, 请自行参考文档.
2. 其它的情况下的地应力平衡 [' s4 x& ]% I8 A# x* c1 Y$ q
地表不水平或土分层不规则的情况下只能使用导入地应力数据的方式进行地应力平衡, 其方法的核心是给土体加重力, 不必理会其土体位移, 只是单纯的提取每个单元的平均应力, 再将这个应力读回ABAQUS作为初始应力, 以达到平衡土体位移的效果. 具体操作参考这个贴
http://forum.simwe.com/viewthread.php?tid=762963&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2,
3. 不预设地应力
所谓的不进行地应力平衡是指计算前不设定土体位移归零, 而是记录下土体在重力作用下每个结点的位移值, 在进行后继计算后, 将位移结果减去重力作用下每个结点的位移值就能得到当前计算步的净位移, 这种方法只适用于小变形分析, 不适用于几何非线性的情况, 好处是收敛性好, 精度高, 因为记录的是结点位移, 缺点是麻烦, 因为要提前记录, 还要处理数据.
1法是对单元集进行赋值, 与结点集无关;
常见的问题
1. 有结构物的情况如何平衡?
有结构物的情况下就不是正确的平衡点, 如果一定要这样平衡, 比如原来有一个隧道, 现在在上面建房子, 要估计这种情况下的追加位移是多少, 这种情况下可以用上面的方法2和3来完成.
2. 多层土的情况下的平衡?
多层土的情况可以根据不同的情况选择方法1或2进行.
3. 有地下水的情况下的平衡?
有水的情况下最主要是要搞清楚土的三相, 按照不同的孔隙率计算土的干密度和湿密度, 这种平衡目前最好是只用方法1来完成.
4. 边坡加重力无法平衡?
边坡加重力无法平衡表示边坡本身是不稳的, 请检查模型是否合理, 塑性参数是否合理./
我是新手、自学者，总结了一点经验，主要是给和我一样的新手
*initial condition,type=stress,geostatic 就不说了，啃啃手册就明白了，或者在论坛搜索下 geostatic
我主要说一下*initial conditions,type=stress,input=XXX.XXX的方法：
1、先建立好你的模型，按照你的课题把边界条件、接触什么的都建好，然后使用命令mdb.models['你的模型名字']setValues(noPartsInputFile=ON)，建立没有part的inp文件（abaqus不认带part的istress数据）
2、建立*geostatic的step，只施加重力荷载，然后求解
3、完成之后用viewer打开odb，选report->field report,下拉菜单里面选择 centriod，然后依次把s11、s22、s33、s12、s13、s23点选上，去setup页面选择报告文件的名字（例如istres.rpt)把下面的和值之类的额外输出点掉，要不还要手动删除
4、编辑你的报告文件，把数据行之外的东西都删除。（正确的格式：单元号，s11，s22，s33，s12，s13，s23），使用excel或者ultraedit在每个数据间加上一个逗号（ , ），数据之间的空格不要紧，但是要注意不要有空行或者逗号之间是空白，否则会出现element 0错误
5、重新建立你的step，写入新的inp文件之后，编辑它，在材料属性之后、step之前加入*initial conditions,type=stress,input=XXX.XXX
搞定，位移一般都被控制在10-6数量级，满足要求。
在此感谢论坛的各位前辈，也建议新人多搜索一下论坛的经典旧帖
---------------------------------------
补充一下：
使用mdb.models['你的model的名字'].setValues(noPartsInputFile=ON) 时注意:
- 要写你的model的名字（在CAE窗口顶部，Model右面，默认名称为Model-1, Model-2等等），不是CAE或JOB的名称。
- 你的model的名字前后都要加上单引号
- 这个命令要在CAE窗口底部输入（ 先点左下角的 >>> ）
- 输入此命令后，要在Job模块中用write job来生成INP文件（loveestboy1981 是不是没作这一步？）
- 生成的INP文件名和Job名一样，而不是CAE文件或Model的名称。
- 生成的INP文件在当前的工作目录下（不一定是TEMP文件夹）
你是要用地应力平衡吗？
绛旓細濡傛灉妯″瀷姣旇緝瑙勫垯锛屾病鏈夋々鍩烘垨鑰呭湴闈㈠缓绛戯紝鍙互鍦╨oad涓柦鍔犲姏body force锛岀劧鍚庡湪Key word閲岄潰杈撳叆鐩稿簲骞宠 鍛戒护銆傚鏋滄ā鍨嬫瘮杈冨鏉傦紝灏辩◢寰夯鐑︼紝浣嗘讳綋涓嶇敤瀵煎嚭鑺傜偣鍔涳紝鏈澶氳緭涓懡浠ゅ弽澶嶅钩琛″嚑娆°
绛旓細鐪熷疄鍦扮煡閬撶幇鐘惰竟鍧＄殑鍐呭姏鏄緢闅剧殑,鎴戜滑閲囧彇鐨勫姙娉曟槸,鐢ㄦ垜浠墍寤虹珛鐨勫嚑浣曟ā鍨嬫柦鍔犲拰瀹為檯妯″瀷涓鑷寸殑閲嶅姏鍜岃竟鐣屾潯浠惰繘琛璁＄畻,寰楀埌鍙樺舰鍚庢垨鍙樺緱鏇村皬鎴栦笌鐜扮姸杈瑰潯涓嶅畬鍏ㄤ竴鑷寸殑杈瑰潯鍐呭姏杩戜技鐨勪綔涓虹幇鐘惰竟鍧＄殑鍐呭姏,骞堕噸鏂板皢鍏舵柦鍔犱簬涓庣幇鐘惰竟鍧′竴鑷寸殑鍑犱綍妯″瀷,鍐嶆柦鍔犻噸鍔(褰撶劧杈圭晫鏉′欢涔熷簲鍩烘湰涓鑷)浠ュ钩琛,杩欐牱鎵嶇畻寤虹珛浜嗕笌鐜扮姸...
绛旓細閲囩敤鍓戞ˉ妯″瀷鏃,寮规ч儴鍒嗛渶閫夋嫨澶氬瓟寮规фā鍨,鍦ㄨ繘琛屽湴搴斿姏骞宠 鏃,鎴戠殑鍔炴硶鏄:绗竴姝,鐢ㄥ脊鎬фā鍨,鍔犻噸鍔涜嵎杞,绠楀嚭鍦熶綋搴斿姏;绗簩姝,鐢ㄥ墤妗ユā鍨,鎶婄涓姝ョ殑鍒濆搴斿姏鐢╥nitial conditions,type=stress,input=inistress.txt鏉ュ鍏,浣嗘暟鎹鏌ョ殑鏃跺欒繕鏄彁绀篨XXmust be positive,璇烽棶浣犲綋鏃舵槸鎬庝箞瑙ｅ喅鍓戞ˉ妯″瀷鐨勫湴搴斿姏骞宠 闂...
绛旓細Plaxis銆ABAQUS杞欢銆侾LAXIS2D/3D绋嬪簭鏄敱鑽峰叞PLAXISB.V.鍏徃鎺ㄥ嚭鐨勪竴绯诲垪鍔熻兘寮哄ぇ鐨勯氱敤宀╁湡鏈夐檺鍏璁＄畻杞欢锛屽凡骞挎硾搴旂敤浜庡悇绉嶅鏉傚博鍦熷伐绋嬮」鐩殑鏈夐檺鍏冨垎鏋愪腑锛屽:澶у瀷鍩哄潙涓庡懆杈圭幆澧冪浉浜掑奖鍝嶃佺浘鏋勯毀閬撴柦宸ヤ笌鍛ㄨ竟鏃㈡湁寤虹瓚鐗╃浉浜掍綔鐢ㄣ佸ぇ鍨嬫々绛忓熀纭(妗ユ々鍩虹)涓庨偦杩戝熀鍧戠殑鐩镐簰褰卞搷銆佹澘妗╃爜澶村簲鍔涘彉褰㈠垎鏋愩佸簱...
绛旓細杩欎釜 姣旇緝澶嶆潅 鏈変笁绉嶆濊矾 1 鏌ユ壘鐩稿叧鏂囩尞 杩涜鍙傛暟纭畾 2 鏍规嵁浜嬪疄瀹為獙缁撴灉 閫氳繃澶ч噺绠椾緥杩涜璋冭瘯 3 鐪熷疄鍋氭尟鍔ㄥ疄楠岋紝姹傜殑鍦熶綋鍙傛暟
绛旓細鍦ㄧ矘鍦熶腑鎵撴々涔熶細寮曡捣鍦拌〃闈㈤殕璧凤紝鏍规嵁涓嶅悓寤虹瓚鑽疯浇瑕佹眰鍙婂満鍦版潯浠讹紝妗╅暱10锝24m锛岀蹇冮敜閲25锝50kN锛岃惤璺3锝5m锛屽簲鍦ㄥ紑鎸栧墠鎵撴々锛堜絾闇閫佹々锛屼細闄嶄綆鎵撴々鏁堢巼锛夛紝鎵撴々浼氬紩璧峰潙搴曢殕璧枫傚姞澶ф々璺濆彲鍑忓皯妗╂暟銆傚湪娣卞熀鍧戝唴锛岃繕闇娉ㄦ剰鍒版墦妗╂椂鍦熶綋涔嬮棿鐩告尋鍘嬮犳垚闅嗚捣鍙婃柇妗╃瓑闂锛屾婚殕璧烽噺澶х害鐩稿綋浜庣兢妗...
绛旓細abaqus/cae娌℃湁鐩存帴鐨勬柟娉曢夋嫨鍐呴潰锛岄渶瑕佸涓繃婊ゆ柟寮忛厤鍚堣捣鏉ラ夋嫨锛岃瑙佷笅鍥撅紙绗簩鍥剧墖鏄夋嫨鐨勭粨鏋滐紝涓涓暱鏂逛綋鍐呯殑鍐呴儴绌虹悆闈級锛
绛旓細鍦熶綋鍗曞厓涓篊PE8RP瀛斿帇鍗曞厓锛岄瀹氫箟鍦轰腑鍒濆闃舵璁剧疆瀛旈殭姣旓紝鍏蜂綋鍙湅璐瑰悍鑰佸笀鐨勩ABAQUS宀╁湡宸ョ▼瀹炰緥璇﹁В銆7,8绔狅紝缃戜笂鏈夊搴旂殑cae鏂囦欢
绛旓細棣栧厛瑕佹槑纭竴浠朵簨鎯咃細涓涓崟鍏冨彧鑳芥湁涓缁勫脊鎬фā閲忋傞偅涔堣瀹炵幇绾挎у彉鍖栵紝鍙兘灏介噺澶氱殑鍒掑垎鍗曞厓锛岀劧鍚庢瘡涓崟鍏冩寚瀹氬悇鑷脊鎬фā閲忋傚綋鐒讹紝杩欐牱鍋氬緢楹荤儲锛屾墍浠ユ姌涓殑鍔炴硶鏄寜鐓ч珮搴﹀垝block锛屾瘡涓猙lock鎸囧畾涓绉嶆潗鏂欏睘鎬с傚鏋滆繕鏄珜楹荤儲锛屽彲浠ョ湅鐪abaqus鐨鍛戒护娴佹枃浠讹紝鍙互鎻愰珮鍋氭绫婚噸澶嶅伐浣滅殑鏁堢巼 ...
绛旓細绗竴姝ュ厛璁＄畻鍒濆鍦板簲鍔涳紝鍔犻噸鍔涘悗灏嗗簲鍔涘鍑猴紝鐒跺悗璇诲叆锛屽啀绠椾竴閬嶏紝杩欐椂鍙互灏嗗湪閲嶅姏浣滅敤涓嬬殑浣嶇Щ鎶垫秷銆傚湪abaqus涓涔熸槸鐢ㄨ繖涓柟娉曪紝浣嗗彲浠ョ敤*geostatic杩涜鑷姩骞宠 锛屽叾鍦板簲鍔涚殑绯绘暟骞朵笉鏄儚ansys鏄氳繃鍦熶綋娉婃澗姣璁＄畻鐨銆傜劧鍚巈kill涓姝ヤ竴姝ョ殑鎸栦簡銆傛瘡涓姝ラ兘瑕佺敓鎴愰噸鍚姩鏂囦欢锛屾敞鎰弐estart鏂囦欢涓嶄繚瀛樺崟鍏冪殑...
