lfycsu 发表于 2006-5-9 15:57:37

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

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

3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如
ELEMENTPT 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 发表于 2007-4-9 00:10:03

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

土木锋 发表于 2009-1-13 17:48:43

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

Allen-2000 发表于 2007-4-4 11:43:05

原帖由 pibird 于 2007-4-1 14:20 发表 http://www.simwe.com/forum/images/common/back.gif
我地应力平衡后，位移达到了e-6，但是地表的应力有30kpa。
不知道大家做的怎么样。
俺做的也是这样，地表应力不为0。差不多也是30KPA左右，

sunmoon 发表于 2006-5-9 16:33:26

好东东

lyj198171 发表于 2006-5-9 18:40:47

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

firce 发表于 2006-5-9 20:43:19

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

lfycsu 发表于 2006-5-9 21:46:30

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

airtanks 发表于 2006-7-13 16:22:29

具体做法？

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

airtanks 发表于 2006-7-14 11:21:16

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 发表于 2006-7-16 22:08:44

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

msnsky 发表于 2006-11-13 17:30:06

chenpeilei 发表于 2006-11-13 19:57:52

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

tzwjc 发表于 2007-1-25 00:03:33

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

wmhappy2000123 发表于 2007-3-31 19:07:42

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

yuminhust2005 发表于 2007-4-1 11:42:34

原帖由 firce 于 2006-5-9 20:43 发表 http://www.simwe.com/forum/images/common/back.gif
geostatic相比于static有什么好处？为什么不直接用后者？
听说是因为用geostatic只施加了应力不会有变形，而static会产生变形。
不过好像也没有什么大的区别！
希望真正了解内情的高手赐教。

pibird 发表于 2007-4-1 14:20:44

问题

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

heerduo 发表于 2007-4-5 11:42:40

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

guoyoulin0708 发表于 2007-4-8 15:17:51

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

rocksoul 发表于 2007-4-16 17:14:45

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

rocksoul 发表于 2007-4-16 17:19:58

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

[ 本帖最后由 rocksoul 于 2007-4-18 09:07 编辑 ]
页: [1] 2 3 4
查看完整版本: 我的初始地应力平衡方法!