找回密码
 注册
Simdroid-非首页
查看: 6019|回复: 63

[土木] 保证通过的地应力平衡方法

[复制链接]
发表于 2007-3-31 15:54:23 | 显示全部楼层 |阅读模式 来自 上海
本人是个初学者,在论坛上看到地应力平衡法的贴子,发现有很多地方没讲清楚,现将地应力平衡法进行完善,调试通过,与大家分享。SC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent

地应力平衡方法:
第一步:建立模型,材料,分析步(GEOSTATIC

第二步:施加荷载,LOAD,选择施加重力GRAVITY,在你想施加重力的方向输入数值9.8SimWe
第三步:在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)  (请严格按照这个格式,注意大小写的字母)仿
第四步:提交J0B,完成后
第五步:按以下步骤,[Module]Visualization-Report---Report Field Output---下拉菜单里面选择
centriod,然后依次把s11s22s33s12s13s23点选上,setup页面选择报告文件的名字---***.txt
---Write中选择Field Output-------------ok!!!
第六步:
先打开excel,在excel中打开***.txt,选择“分隔符号选择“导入的起始行---下一步-----在出来的界面上点完成------保存为CSV逗号分隔文件***.csv(文件内容格式如下)

单元号     S11    S22   S33   S12   S13   S23 请注意,在保存内容中没有这一行的
   1        ,     .    ,   .    ,   .   ,    .   ,   .    ,
   2        ,     .    ,   .    ,   .   ,    .   ,   .    ,  
   .         ,     .    ,   .    ,   .   ,    .   ,   .    ,
   .         ,     .    ,   .    ,   .   ,    .   ,   .    ,  
这个结果文件是最重要的,在所保存的文件中只有数值部分,没有英文字母,没有上面那个单元号这一行,而且单元号前面也没有什么PART名字什么的,就是12........这些这些数字些数字。在CSV文件中,数字本身是用逗号隔开的,不需要另加逗号。


第七步:在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上:*initial conditions,type=stress,input=***.csv
完成仿



第八步:重新提交JOBOK



[ 本帖最后由 wmhappy2000123 于 2007-3-31 15:58 编辑 ]

评分

1

查看全部评分

发表于 2007-4-15 16:36:39 | 显示全部楼层 来自 上海虹口区
Simdroid开发平台
An initial condition has been specified on element set initialconditions - this element set is not active in the model是怎么回事?
发表于 2007-4-15 16:56:38 | 显示全部楼层 来自 德国
原帖由 iisltyf 于 2007-4-15 16:36 发表
An initial condition has been specified on element set initialconditions - this element set is not active in the model是怎么回事?




http://www.simwe.com/forum/viewthread.php?tid=746467
ABAQUS常见问题汇总(2.0版)

1.
论坛提问要诀
发表于 2007-4-15 18:06:59 | 显示全部楼层 来自 江苏徐州
第三步:在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)  (请严格按照这个格式,注意大小写的字母)仿SimWe
请问这个是在python命令行里面输入吗?菜鸟问
发表于 2007-4-15 18:07:59 | 显示全部楼层 来自 江苏徐州
请问这个是在python命令行里面输入吗?菜鸟问

第三步:在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)  (请严格按照这个格式,注意大小写的字母)仿SimWe

发表于 2007-4-15 22:33:46 | 显示全部楼层 来自 河南郑州
  回去试试
 楼主| 发表于 2007-4-16 11:33:07 | 显示全部楼层 来自 上海

回复 #2 iisltyf 的帖子

这是因为你没有按照第三步运行
发表于 2007-4-16 16:13:07 | 显示全部楼层 来自 大连理工大学
在地质情况比较复杂或模型比较复杂的时候,最好是要用input的形式进行地应力平衡,而且往往要考虑侧向压力系数,个人理解上述方法好像就是在重力或体力作用下的静水压力平衡,没有考虑侧压系数,比如侧压系数为2.0。
那么该如何考虑侧压系数呢? 如果理解有误请指正。
发表于 2007-4-16 16:45:10 | 显示全部楼层 来自 天津
我想请教一下,用这种方法平衡地应力之后,怎么考虑固结呢?
发表于 2007-4-16 17:12:07 | 显示全部楼层 来自 大连理工大学
我按楼主操作存成csv格式的,运行提示错误:Error in job j3: in keyword *INITIALCONDITIONS, file "j3.inp", line 599: Unknown part instance name 134
Error in job j3: in keyword *INITIALCONDITIONS, file "j3.inp", line 599: Unknown part instance name 2-1
Error in job j3: in keyword *INITIALCONDITIONS, file "j3.inp", line 599: Unknown part instance name 3-2

加上逗号,后存成inp或csv格式的都可以运行,难道是我在导入excel的txt时候操作少了什么????
发表于 2007-4-16 17:33:30 | 显示全部楼层 来自 大连理工大学
如果没有好的办法,只好把导出的水平应力结果在excel中进行公式运算了。呵呵
发表于 2007-4-18 00:08:50 | 显示全部楼层 来自 天津
我模拟的是桩土相互作用,我想请教的是,是否在桩土装配之后,用该方法平衡地应力?还是先做土的地应力平衡,再加桩呢?先谢谢啦!
 楼主| 发表于 2007-4-19 00:10:17 | 显示全部楼层 来自 上海

回复 #9 chengshics 的帖子

那要自己计算地应力,输入*initial condition,type=stress,geostatic
单元名,值1,坐标1,值2,坐标2。
注意考虑水时,要注意饱和度影响孔隙水压的对平衡时的有效应力贡献,请参见 abaqus在土木工程应用一书中的固结部分,和P129页、P160页

评分

1

查看全部评分

 楼主| 发表于 2007-4-19 00:14:35 | 显示全部楼层 来自 上海

回复 #12 chengshics 的帖子

还是先做土的地应力平衡,再加桩。不过还有其他方法,那就是直接对不同单元输入应力值
发表于 2007-4-19 08:19:16 | 显示全部楼层 来自 天津

回复 #14 wmhappy2000123 的帖子

感谢wmhappy2000123的帮助,但是我想问的是用您帖子中提到的.inp文件的方法平衡地应力的话,是否先做土和桩的地应力平衡,然后输出s值,然后在做固结呢?
发表于 2007-4-27 20:54:34 | 显示全部楼层 来自 浙江杭州
请问语句mdb.models['模型名字'].setValues(noPartsInputFile=ON)  是什么意思?谢谢。
发表于 2007-5-27 23:05:18 | 显示全部楼层 来自 天津
考虑固结时的地应力平衡,得到的土底应力值是否应该和不考虑固结时的一样呢?
回复 不支持

使用道具 举报

发表于 2007-8-3 12:23:53 | 显示全部楼层 来自 辽宁大连

请楼主帮忙看看!

我的地应力平衡几何关系很简单,地表水平,只承受土体重力。土体分为三层,第一层线弹性模型,第二层D-P模型,第三层剑桥模型,总共20m深,地下水位线位于地表下一米处。由于地表水平。我选择了较为简单的type=stress,Geostatic来平衡地应力,但是结果很奇怪,尽管位移达到了e-4级,但是最大位移却发生在第二层与第三层的交接处。大家帮忙看看为什么会出现这种情况?
位移云图如下图:

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

发表于 2007-8-4 23:36:55 | 显示全部楼层 来自 上海
如果考虑固结,在最开始输入的inp文件中设置出示孔隙比和孔压,然后提交,后续同楼主所说,最后在inp文件中地应力的初始条件用文件输入即可
回复 不支持

使用道具 举报

发表于 2007-8-4 23:43:32 | 显示全部楼层 来自 上海
零点几毫米的位移已经算平衡了
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-1 21:35 , Processed in 0.070496 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表