我的初始地应力平衡方法!
看了版上很多帖子,受益菲浅!关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了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 编辑 ] 个人觉得CAE中
Visualization模块中Report的Field Output的应力输出会比
*el print
s,
的方法更方便输出,个人小建议
大家参考一下吧!! 在运行STEP——geostatic时,我有时也会出现如下错误:The fixed time increacement is too large。 不知道为何,请问如何修改?有时候能算,很奇怪。。。 原帖由 pibird 于 2007-4-1 14:20 发表 http://www.simwe.com/forum/images/common/back.gif
我地应力平衡后,位移达到了e-6,但是地表的应力有30kpa。
不知道大家做的怎么样。
俺做的也是这样,地表应力不为0。差不多也是30KPA左右, 好东东 为什么我在geostatic步总是算不了? 说步长过大 geostatic相比于static有什么好处?为什么不直接用后者? 我碰到的情况,步长过大的原因可能是因为土体参数的值取的小了,由于重力的作用很大,所以不收敛了!static有人用过,不过既然用geostatic就可以做了,那不是更好!
具体做法?
1、我的geostatic没办法算啊,总是出现Fixed timeincrement is too large!2、不知道上述命令在哪里输入?难道是写字板中?
恳请高人解答 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,什么原因呢? 地应力平衡每建一种不同的模型是不是一定要重新取数据平衡一次??
重力荷载不要删吗?? 这个方法好,这样能把一些坡面下初始应力去掉。但我觉得用*initial conditions, TYPE=stress, input=*.inp比较方便 原帖由 lfycsu 于 2006-5-9 15:57 发表
单独存为一个*.dat文件
4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为:instance名称.序号;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!
接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),
...
请问怎么“其中的1所在的列”去掉?存了后,打开就没有带逗号分隔符的格式,
请朋友们帮助!谢谢! 我的贴子上有完整的步骤,希望对你有帮助 原帖由 firce 于 2006-5-9 20:43 发表 http://www.simwe.com/forum/images/common/back.gif
geostatic相比于static有什么好处?为什么不直接用后者?
听说是因为用geostatic只施加了应力不会有变形,而static会产生变形。
不过好像也没有什么大的区别!
希望真正了解内情的高手赐教。
问题
我地应力平衡后,位移达到了e-6,但是地表的应力有30kpa。不知道大家做的怎么样。 我做完平衡后,位移达到了仍有几个mm的差异,不知道怎么做能改善变小一点,
因为我的总位移只有几公分,几个mm会严重影响精度。
还有我的地表应力也不为零,约有几个KPa,有办法让他是零? 好东西 强烈顶起啊 感谢楼主 如果是多个part组成的model,导出的应应力的单元号好像应该有重号码的,那该怎么办,
还有采用input地应力平衡方法的时候,怎样考虑施加侧压系数呢? 上述步骤对多个part的model也适用,步骤也是一样。
[ 本帖最后由 rocksoul 于 2007-4-18 09:07 编辑 ]