24小时热门版块排行榜

返回列表

mzgnjust

金虫 (小有名气)

应助: 20 (小学生)
金币: 776.4
散金: 11
红花: 2
帖子: 180
在线: 90.8小时
虫号: 718140
注册: 2009-03-09
性别: GG
专业: 金属材料的力学行为

[求助] ansys焊接模拟过程中出现的错误

模拟过程中出现了如下错误：
The value of TEMP at node 434 is 1.594205862E+25.  It is greater than
the current limit of 1.E+25.  This generally indicates rigid body
motion as a result of an unconstrained model.  Verify that your model
is properly constrained.

查看节点温度，时该节点的温度超高，且该区域的温度都很高。以下是我的命令流，求高手帮忙看看，哪里出了问题，谢了~

/CLEAR !清除数据库，新建文件
/FILNAME, SY_FANG_02 !改变工作名
!开始热分析
!建立有限元模型
/NERR,,10000000
/PREP7 !进入前处理器
ET,1,PLANE77 !平面热单元
ET,2,SOLID5 !空间热单元
V=0.007 !焊接线速度
LENGTH=0.035 !钢板沿焊缝方向长度
L_NUMBER=20 !沿焊缝方向单元段数
TIME=0.01
TIME_INC=LENGTH/(L_NUMBER*V) !第一层焊缝时间增量
!
!
!
!
MPTEMP,1,0,300,600,1000,1300,2000 !定义材料模型，温度表
MPDATA,DENS,1,1,7859,7770,7659,7600,6750,6750 !密度
MPDATA,DENS,2,1,7859,7770,7659,7600,6750,6750
MPDATA,C,1,1,450,514,967,1370,2800,2000 !比热
MPDATA,C,2,1,450,514,967,1370,2800,2000
MPDATA,KXX,1,1,67,53,39,31,20,20 !传导系数
MPDATA,KXX,2,1,67,53,39,31,20,20
MPTEMP,1,0,300,600,1000,1300,1400,2000 !温度表
MPDATA,EX,1,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7 !弹性模量
MPDATA,EX,2,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7
MPDATA,PRXY,1,1,0.33,0.33,0.35,0.36,0.4,0.4 !泊松比
MPDATA,PRXY,2,1,0.33,0.33,0.35,0.36,0.4,0.4
MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6 !线膨胀系数
MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6
TB,BKIN,1,5
TBTEMP,15
TBDATA,1,280E6,2E9 !屈服极限，切向模量
TBTEMP,800
TBDATA,1,10E6,1E7
TBTEMP,2500
TBDATA,1,10E6,1E6
MPTEMP,,,,,,,,                   !初始化随温度变化的材料参数表格
MPTEMP,1,20                      !定义表格的温度数值
MPTEMP,2,250
MPTEMP,3,400
MPTEMP,4,500
MPTEMP,5,750
MPTEMP,6,800
MPTEMP,7,1000
MPTEMP,8,1200
MPTEMP,9,1500
MPTEMP,10,1700
MPTEMP,11,2500
MPDATA,HF,1,,100                   !定义对应温度下的对流换热系数数值
MPDATA,HF,1,,350
MPDATA,HF,1,,480
MPDATA,HF,1,,520
MPDATA,HF,1,,1000
MPDATA,HF,1,,1050
MPDATA,HF,1,,1500
MPDATA,HF,1,,2000
MPDATA,HF,1,,3000
MPDATA,HF,1,,3100
MPDATA,HF,1,,3500
MPTEMP,,,,,,,,                   !初始化随温度变化的材料参数表格
MPTEMP,1,20                      !定义表格的温度数值
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,1000
MPTEMP,7,1100
MPTEMP,8,1500
MPTEMP,9,1700
MPTEMP,10,2500
MPDATA,HF,2,,120                !定义对应温度下的对流换热系数值
MPDATA,HF,2,,210
MPDATA,HF,2,,300
MPDATA,HF,2,,460
MPDATA,HF,2,,540
MPDATA,HF,2,,1200
MPDATA,HF,2,,1200
MPDATA,HF,2,,2800
MPDATA,HF,2,,3000
MPDATA,HF,2,,3100
!
TBCOPY,BKIN,1,2 !复制材料BKIN数据
MP,REFT,1,20 !参考温度
MP,REFT,2,20
!
!
!
K,1,0,0
K,2,0.015,0
K,3,0.015,0.004
K,4,0,0.004
K,5,0.002,0.004
K,6,0.002,0
A,1,6,5,4
A,6,2,3,5
AGLUE, ALL
LESIZE,1,,,3
LESIZE,3,,,3
LESIZE,4,,,4
LESIZE,2,,,4
LESIZE,6,,,4
LESIZE,5,,,10,4
LESIZE,7,,,10,0.25
VEXT,1,2,1,0,0,0.035
VGLUE,ALL
LESIZE,14,,,30
TYPE,2
MAT,1    !2号材料为母材，1号材料为焊缝
REAL,
VSWEEP,1
TYPE,2
MAT,2
REAL,
VSWEEP,2
!
FINISH
/SOLU !进入求解器
   DA,7,SYMM    !对称面A12施加对称载荷!已改
   DL,5,,UZ
   DL,1,,UZ
   DL,4,,UX
   DL,12,,UX
TUNIF,20,
TOFFST,273,
!
!
!!!!!!!!求解设置!!!!!!!!!!!!!!!!!!!!!!
ANTYPE,4 !瞬态分析
TRNOPT,FULL
LUMPM,0
NLGEOM,1
PSTRES,1 !打开应力刚化
NSUBST,5,20,1
OUTRES,ERASE
OUTRES,ALL,ALL !输出控制
AUTOTS,1 !打开自动时间步长
PRED,ON,,ON
LNSRCH,1 !打开线性搜索
!
!开始激活焊缝单元，进行瞬态热分析
TIMINT,0,MAG
TIMINT,0,ELECT
TIMINT,ON !打开时间积分
!
!
!!!!!!!!!!!!!!定义数组!!!!!!!!!!!!!!!!!!!!!!!
*DIM,QQ,ARRAY,8          !定义数组，存储单元节点热流密度
*DIM,E_NODE,ARRAY,20 !定义数组，存储单元所属节点编号
*GET,ELEM_MAX,ELEM,0,NUM,MAX !获取最大单元号
*GET,ELEM_MIN,ELEM,0,NUM,MIN !获取最小单元号
*GET,NODE_MAX,NODE,0, NUM,MAX !获取最大节点号
*GET,NODE_MIN,NODE,0,NUM,MIN !获取最小节点号
*DIM,NODE_AVE,ARRAY,NODE_MAX !定义数组，存储节点相邻单元数量
!
!
!
!
ESEL,S,MAT,,1 !选择焊缝单元!!!!!!!!!!!!!!!!!
EKILL,ALL !杀死单元!!!!!!!!!!!!!!!
!
!
alls!!!!!!!!!!!!!!
!
!
R=0.007 !电弧有效加热半径
  Q=140*20*0.8 !输入热量
  QMAX=3*Q/(3.1415926*R*R) !变量
  XC= 0 !第一层焊缝质心X坐标
  YC= 0.007 !第一层焊缝质心Y坐标
!
!
!
!!!!!!!!!!!!从头至尾激活焊缝，开始
*DO,J,1,L_NUMBER
alls
SFEDELE,ALL,ALL,HFLUX !删除以前施加的热流密度
VSEL,S,,,1 !选择第KK层焊缝所对应的体!!!!!!!!!!!!!!!!!!!!!!
ESLV,S,1 !选择体上单元即第KK层焊缝单元!!!!!!!!!!!!!!!!!!!!
NSEL,R,LOC,Z,(J-1)*LENGTH/L_NUMBER,(J+1)*LENGTH/L_NUMBER!!!!!!!!!!!!!!!!!!!!!
!选择第J行单元的中间节点!!!!!!!!!!!!!!!!!!!!!
ESLN,R !选择第KK层第J行单元!!!!!!!!!!!!!!!!
EALIVE,ALL !激活单元!!!!!!!!!!!!!!!!!
ALLS !!!!!!!!!!!!!!!
ZC=J*LENGTH/L_NUMBER !第KK层第J行单元质心Z坐标
*DO,I,ELEM_MIN,ELEM_MAX !计算施加单元热流密度，循环开始
   *IF,ESEL(I),EQ,1,THEN
   ESEL,S,,,I !选择第I个单元
         *DO,K,1,8 !计算第I单元前8个节点热流密度
         *GET,E_NODE(K),ELEM,I,NODE,K !获取第I个单元第K个节点的编号
         SS=E_NODE(K)
         DISP=SQRT((NX(SS)-XC)*(NX(SS)-XC)+(NY(SS)-YC)*(NY(SS)-YC)+(NZ(SS)-ZC)*(NZ(SS)-ZC))
         EEE=3*DISP*DISP/(R*R) !变量
         *IF,DISP,LT,R,THEN
         QQ(K)=QMAX*EXP(-EEE)/8 !第K个节点的热流密度
         *ELSE
         QQ(K)=0
         *ENDIF
         *ENDDO
      *IF,QQ(1)+QQ(2)+QQ(3)+QQ(4),NE,0,THEN
      SFE,I,1,HFLUX,,QQ(1), QQ(2), QQ(3),QQ(4) !在第I个单元第一个面上加热流密度
      *ENDIF
      *IF,QQ(1)+QQ(2)+QQ(5)+QQ(6),NE,0,THEN
      SFE,I,2,HFLUX,,QQ(1), QQ(2), QQ(6),QQ(5) !在第I个单元第二个面上加热流密度
      *ENDIF
      *IF,QQ(6)+QQ(2)+QQ(3)+QQ(7),NE,0,THEN
      SFE,I,3,HFLUX,,QQ(2), QQ(3), QQ(7),QQ(6) !在第I个单元第三个面上加热流密度
      *ENDIF
      *IF,QQ(1)+QQ(4)+QQ(5)+QQ(7),NE,0,THEN
      SFE,I,5,HFLUX,,QQ(1), QQ(4), QQ(7),QQ(5) !在第I个单元第五个面上加热流密度
      *ENDIF
      *IF,QQ(5)+QQ(6)+QQ(7)+QQ(8),NE,0,THEN
      SFE,I,6,HFLUX,,QQ(5), QQ(6), QQ(7),QQ(8) !在第I个单元第六个面上加热流密度
      *ENDIF
      ALLS!!!!!!!!!!!!!!!
      *ENDIF
   *ENDDO !计算施加单元热流密度，循环结束
   ALLS
!
!
!
TIME=TIME+TIME_INC !时间变量
!
TIME,TIME !设定分析时间
SOLVE
*ENDDO
SAVE !保存
FINI

回复此楼

» 收录本帖的淘帖专辑推荐

焊接有限元分析

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2012-05-28 20:16:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

pengqiong

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 34.3
帖子: 8
在线: 2.8小时
虫号: 1771239
注册: 2012-04-22
专业: 机械结构强度学

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
臭水沟: 金币+2, 谢谢交流~~ 2012-05-30 12:28:03

不知道具体哪里的问题，仅把我认为不妥的地方说下。
1、solid5是8节点的六面体单元，是没有中间节点的，你用77号单元分的面网格，sweep后的Z方向应该没有中间节点，这个你自己试下，所以激活焊缝时候可能出错。
2、施加热流密度时候为啥只加了5个面？有点不理解，第四个面呢？
3、没有细看你的命令流，好像你没有施加热分析的边界条件，这样，热只会进不会出。
4、焊接这个东西，三维的话，就用顺序耦合吧，看你的命令流应该是借鉴了别人的，但又用了不同的耦合方法。

赞一下

回复此楼

2楼2012-05-30 01:03:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mzgnjust

金虫 (小有名气)

应助: 20 (小学生)
金币: 776.4
散金: 11
红花: 2
帖子: 180
在线: 90.8小时
虫号: 718140
注册: 2009-03-09
性别: GG
专业: 金属材料的力学行为

引用回帖:

2楼: Originally posted by pengqiong at 2012-05-30 01:03:48
不知道具体哪里的问题，仅把我认为不妥的地方说下。
1、solid5是8节点的六面体单元，是没有中间节点的，你用77号单元分的面网格，sweep后的Z方向应该没有中间节点，这个你自己试下，所以激活焊缝时候可能出错。
2 ...

1.用的单元式solid5耦合单元，77号单元当初是为了体扫略网格划分的时候的一个过渡单元，我把77号单元换成不带中间节点的面单元重新计算以后，好像还是不行。
2.热分析边界条件加在了材料的属性里（具体是给材料赋予了对流换热属性）
3.这个命令流的确是借鉴了别人根据自己的情况改的，用顺序耦合的方法的我去试试

我后来在solve之前加了
ESEL,S,LIVE
NSLE,S
只对激活单元和节点进行计算，上面的错误倒是没有了，但是在第二个循环以后激活的焊缝单元上面温度总是上不去，看着就像是跟其他单元隔绝的一样，一会儿我把图贴出来

赞一下(1人)

回复此楼

3楼2012-05-30 22:42:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖