24小时热门版块排行榜    

查看: 2157  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhuqw

金虫 (小有名气)

[求助] 初应力和初始缺陷分析后没有位移结果

我用ANSYS求解经过激光器经过回流过程后,在回流后的应力和变形基础上,在工作状态下的应力,应变和变形情况。思路是分析一个含有初始缺陷(变形)及应力的分析。操作流程是:1.先求出回流后激光器的应力场并写入file.ist文件,并求出工作后的温度分布file2.rth  2.重新建模用UPGEOM更新模型的节点坐标,使模型的形状变成回流后的形状,再对所有节点的所有自由度约束,读入IST也就是回流后的应力文件,然后写出,获得节点反力(file1.rst)文件  3.再重新建模,读入回流后的应力场(file.ist)以及支点反力(file1.rst)  4.读入工作温度分布file2.rth,设置参考温度,求解热应力

现在问题是,最后求解完成以后节点结果只能看到应力应变的,没有位移结果,也就是说我希望获得在有初始变形和应力的情况下,由于温度分布不均匀器件内部应力分布和导致的变形(displacement)。错误文件在我进行第3步后,提示there are no active degrees of freedom。我没一步都是重新写一个文件的,第三部除了施加个别面的约束外并没有对所有节点施加约束,不知道为什么会没有位移结果?

我最后一步求解工作应力的求解代码:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建模画网格!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!更新节点坐标,读入回流应力场及支反力!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/PREP7
UPGEOM,,,,file,rst
/SOLU
INISTATE,READ,file,ist
ASEL,S,AREA,,3
D,ALL,UY,0
ASEL,S,AREA,,2
D,ALL,UZ,0
LDREAD,REAC,,,,,file1,rst
SOLVE
FINISH

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!读入温度场求解工作应力!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOLU
D,ALL,ALL
ANTYPE,STATIC
EQSLV,SPARSE
BCSOPTION,,INCOERE
TREF,25
LDREAD,TEMP,,,,,file2,RTH
ASEL,S,AREA,,3
D,ALL,UY,0
SOLVE
FINISH
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuqw

金虫 (小有名气)

引用回帖:
2楼: Originally posted by 小左儿 at 2016-03-21 16:31:28
看你的命令流:
ASEL,S,AREA,,3
D,ALL,UY,0
ASEL,S,AREA,,2
D,ALL,UZ,0
  !D表示对节点施加命令流,而你并没选择节点,只选择面,所以你施加约束施加在全部节点上而不是面的节点上。
如果你想施加在面的节点上 ...

哦哦,谢谢!之前在帮助手册里看到这个命令一直查不出来是什么意思,现在知道了。
但是我现在用这个命令去约束以后就一直程序崩溃,我觉得是因为不收敛,我把四周的面都约束了,很多不该约束的也约束了,还是一直崩溃。。。。不知道有什么办法能解决?
3楼2016-03-22 09:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

小左儿

银虫 (正式写手)

看你的命令流:
ASEL,S,AREA,,3
D,ALL,UY,0
ASEL,S,AREA,,2
D,ALL,UZ,0
  !D表示对节点施加命令流,而你并没选择节点,只选择面,所以你施加约束施加在全部节点上而不是面的节点上。
如果你想施加在面的节点上,这样:
ASEL,S,AREA,,3
nsla      
D,ALL,UY,0

不知道是不是这个问题
微信公众号:CAE技术分享
2楼2016-03-21 16:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小左儿

银虫 (正式写手)

引用回帖:
3楼: Originally posted by zhuqw at 2016-03-22 09:52:51
哦哦,谢谢!之前在帮助手册里看到这个命令一直查不出来是什么意思,现在知道了。
但是我现在用这个命令去约束以后就一直程序崩溃,我觉得是因为不收敛,我把四周的面都约束了,很多不该约束的也约束了,还是一直 ...

应该不是不收敛,不收敛程序也不会直接崩溃,而是算着算着就停了。
有可能是过约束了。
微信公众号:CAE技术分享
4楼2016-03-22 10:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuqw

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 小左儿 at 2016-03-22 10:10:49
应该不是不收敛,不收敛程序也不会直接崩溃,而是算着算着就停了。
有可能是过约束了。...

还真是。。。我把所有约束都取消了,算出来了,不过这样跟真实情况又不一样了。。因为真实情况底面是被约束的。。。

还有这样算下来还是没有位移结果。。。只有应力结果。。。
5楼2016-03-22 10:36:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见