- 积分
- 0
- 注册时间
- 2018-5-17
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2019-9-19 17:25:37
|
显示全部楼层
来自 天津
地应力精确平衡和粘弹性边界静动力转换
现将模型侧面和下底面用传统约束边界Displacement约束,加载Geostatic步分析,在结果中输出场变量(field output),将积分点的S量中的S11、S22、S33、S12(平面应变)输出为inp文件,方法参见那本书,生成的inp就在工作目录中,找出后用Excel或者Ultraedit修改成书上那个格式(具体的不说了),然后导出成dat文件,比如就叫s11-s12.dat。
还是上次运算的那个模型,在关键词中于*STEP前加入“*initial conditions,type=stress,input=s11-s12.dat”,再次加载Geostatic步进行计算,算出结果后在后处理窗口中继续场变量输出,选择“唯一节点(Unique Nodal)”下的“反力RF(Reaction Force)”下的RF1、RF2项,输出为inp;然后修改该inp文件,这里我一般用Ultraedit修改,改的方法和1中类似,多余的说明文字要全去掉,只留下数字,格式为“name of part instance(部件名称),node(节点号),degree of freedom(自由度,RF1就是1之类的),magnitude(值)”的形式,注意RF输出结果为RF1和RF2在同一行,要把RF2移动到RF1下面去。
将约束边界撤销,加上粘弹性边界
在关键词**LOAD后面
使用*Cload
*include,input=xx.inp语句可以写入上面修改后的反力文件(初始应力场也要加进去)
(写入语句后,write input之后如果你import model代入该input文件你会发现语句消失,而load里节点出现很多反力。)
我是按这个方法做的,但是最终结果不对。我想问一下,导出约束反力的时候我是把所有节点的反力导出,再修改关键字加进去,为什么最后的位移特别大,但是应力是对的。这个方法应该是对的啊,第一个模型有约束进行地应力平衡后,位移忽略不计,第二个模型相同只是把约束改成约束反力加进去,但是位移特别大,应力同模型一。这个过程中我还没有加粘弹性边界。麻烦您看一下我的帖子,可能表达更明白一些。帖子为 静-动力人工边界转换,粘弹性人工边界。非常感谢! |
|