ANSYS工字梁焊缝初应力例子

发布日期:[08-11-30 21:56:13] 浏览人次:[]

!建立好模型
…..
/SOLU
!*
ANTYPE,0
/INPUT,''input'',''txt'',
/INPUT,''output1'',''txt'',
!*
ISFILE,READ,1,txt, ,0
SOLVE
=====================
!定义残余应力矩阵
I=1
fy=310e6
!定义单元中心坐标数组,初应力数组
*DIM,EleCenter,ARRAY,NumNode,3,1, , ,
*DIM,EleIS,ARRAY,NumNode,1,1, , ,
*DO,I, 1, NumNode,1
!得到单元中心坐标
*GET,EleCenter(I,1,1),ELEM,I,CENT,X
*GET,EleCenter(I,2,1),ELEM,I,CENT,Y
*GET,EleCenter(I,3,1),ELEM,I,CENT,Z
!焊缝位置在X=1.0 到1.05 之间
*IF,EleCenter(I,1,1),GT,1,THEN
*IF,EleCenter(I,1,1),LT,1.05,THEN
*IF,EleCenter(I,2,1),GT,0.5,THEN
!上翼缘初始应力
EleIS(I,1,1)=-20*fy*EleCenter(I,3,1)**2+0.3*fy
*endif
*IF,EleCenter(I,2,1),LT,0.0,THEN
!下翼缘初始应力
EleIS(I,1,1)=-20*fy*EleCenter(I,3,1)**2+0.3*fy
*endif
!腹板初始应力
*if,EleCenter(I,2,1),GE,0.0,THEN
*IF,ELECENTER(I,2,1),LE,0.5,THEN
EleIS(I,1,1)=4.076*fy*(EleCenter(I,2,1)-0.5)**2
EleIS(I,1,1)=EleIS(I,1,1)-0.719*fy
*IF,EleIS(I,1,1),LT,-0.3*fy,THEN
ELeIS(I,1,1)=-0.3*fy
*ENDIF

*ENDIF
*ENDIF
*ENDIF
*ENDIF
*ENDDO
===============
!输出应力数值
*CFOPEN,''1'',''txt'','' ''
*DO,I,1,3840
!*IF,EleIS(I,1,1),ne,0.,then
*VWRITE,''!''
(A)
*VWRITE,''! Stress for element'',
(A,F)
*VWRITE,''!''
(A)
*VWRITE,''eis,'',I
(A,F)
P=EleIS(I,1,1)
*VWRITE,P,0.,0.,0.,0.,0.
(E,E,E,E,E,E)
!*ENDIF
*ENDDO
*CFCLOS

www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:ANSYS  工字梁  焊缝  初应力  例子
本文所属专题:
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务