我的初始地应力平衡方法!

lfycsu

看了版上很多帖子,受益菲浅!
关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!
1.先施加重力荷载的作用,可以在cae中实现;
2.在inp文件中的output request中写上
*el print
s,
这样就会将施加重力荷载后的应力输出到*.dat文件中了;

3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如
ELEMENT  PT 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,比较简单,希望能对大家有用!

[ 本帖最后由 aba_aba 于 2006-10-10 21:48 编辑 ]

heerduo

个人觉得CAE中
Visualization模块中Report的Field Output的应力输出会比
*el print
s,
的方法更方便输出，个人小建议
大家参考一下吧!!

土木锋

在运行STEP——geostatic时，我有时也会出现如下错误：The fixed time increacement is too large。 不知道为何，请问如何修改？有时候能算，很奇怪。。。

Allen-2000

原帖由 pibird 于 2007-4-1 14:20 发表
我地应力平衡后，位移达到了e-6，但是地表的应力有30kpa。
不知道大家做的怎么样。

俺做的也是这样，地表应力不为0。差不多也是30KPA左右，

sunmoon

好东东

lyj198171

为什么我在geostatic步总是算不了？ 说步长过大

firce

geostatic相比于static有什么好处？为什么不直接用后者？

lfycsu

我碰到的情况,步长过大的原因可能是因为土体参数的值取的小了,由于重力的作用很大,所以不收敛了!static有人用过,不过既然用geostatic就可以做了,那不是更好!

airtanks

1、我的geostatic没办法算啊，总是出现Fixed time  increment is too large!
2、不知道上述命令在哪里输入？难道是写字板中？
恳请高人解答

airtanks

4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance.part.,即单元编号变为:instance名称.part名称.序号;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!
这步不理解啊，能不能详细解释一下。
6.最后在inp文件的step之前写上 *initial conditions,type=stress,input=文件名.csv即可
出现的错误是An initial condition has been specified on element set sand - this element set is not active in the model，什么原因呢？

airtanks

地应力平衡每建一种不同的模型是不是一定要重新取数据平衡一次？？
重力荷载不要删吗？？

chenpeilei

这个方法好，这样能把一些坡面下初始应力去掉。但我觉得用*initial conditions, TYPE=stress, input=*.inp比较方便

tzwjc

原帖由 lfycsu 于 2006-5-9 15:57 发表
单独存为一个*.dat文件
4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为:  instance名称.序号  ;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!
接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),
...

请问怎么“其中的1所在的列”去掉？存了后，打开就没有带逗号分隔符的格式，
请朋友们帮助！谢谢！

wmhappy2000123

我的贴子上有完整的步骤，希望对你有帮助

yuminhust2005

原帖由 firce 于 2006-5-9 20:43 发表
geostatic相比于static有什么好处？为什么不直接用后者？

听说是因为用geostatic只施加了应力不会有变形，而static会产生变形。
不过好像也没有什么大的区别！
希望真正了解内情的高手赐教。

pibird

我地应力平衡后，位移达到了e-6，但是地表的应力有30kpa。
不知道大家做的怎么样。

heerduo

我做完平衡后，位移达到了仍有几个mm的差异，不知道怎么做能改善变小一点，
因为我的总位移只有几公分，几个mm会严重影响精度。
还有我的地表应力也不为零，约有几个KPa，有办法让他是零?

guoyoulin0708

好东西    强烈顶起啊    感谢楼主

rocksoul

如果是多个part组成的model，导出的应应力的单元号好像应该有重号码的，那该怎么办，
还有采用input地应力平衡方法的时候，怎样考虑施加侧压系数呢？

rocksoul

上述步骤对多个part的model也适用，步骤也是一样。

[ 本帖最后由 rocksoul 于 2007-4-18 09:07 编辑 ]