| 查看: 2747 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
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 |
回复此楼» 收录本帖的淘帖专辑推荐
 焊接有限元分析 |
» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有6人回复
-
博士申请都是内定的吗？
已经有7人回复
-
读博
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
ANSYS启动出现问题，寻找解决方案
已经有9人回复- ANSYS二维电场模拟电阻变化，求指教
已经有4人回复
- 关于使用ANSYS进行模拟的问题
已经有8人回复
- 求ANSYS用于焊接分析的生死单元的例子
已经有4人回复
- ansys 加载求助
已经有3人回复
- ansys温度场的模拟出现了问题 求大侠赐教
已经有5人回复
- ansys出现单元过度扭曲的错误
已经有10人回复
- 【原创】ansys使用心得
已经有11人回复
- 【求助】做摩擦焊的热力耦合 是用ansys好 还是abaqus好啊
已经有10人回复
- 【求助】ANSYS分析应力过程中存在的问题
已经有3人回复
3楼2012-05-30 22:42:08
4楼2012-05-30 22:48:48