ABAQUS 地应力平衡.docVIP

?ABAQUS 地应力平衡 abaqus 学习记录 2011-03-30 14:42:12 阅读107 评论0 ??字号：大中小?订阅 关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了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年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。