| 查看: 2998 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 焊接有限元分析 |
» 猜你喜欢
材料与化工300求调剂
已经有26人回复
-
考研二轮调剂
已经有6人回复
-
材料复试求调剂
已经有17人回复
-
085601初试330分找调剂
已经有9人回复
-
291 求调剂
已经有12人回复
-
化学工程与技术324调剂
已经有17人回复
-
0856专硕求调剂 希望是a区院校
已经有12人回复
-
求调剂
已经有6人回复
-
求调剂 材料与工程 324分 专硕
已经有4人回复
-
博士自荐
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- ANSYS启动出现问题，寻找解决方案
已经有9人回复
- ANSYS二维电场模拟电阻变化，求指教
已经有4人回复
- 关于使用ANSYS进行模拟的问题
已经有8人回复
- 求ANSYS用于焊接分析的生死单元的例子
已经有4人回复
- ansys 加载求助
已经有3人回复
- ansys温度场的模拟出现了问题 求大侠赐教
已经有5人回复
- ansys出现单元过度扭曲的错误
已经有10人回复
- 【原创】ansys使用心得
已经有11人回复
- 【求助】做摩擦焊的热力耦合 是用ansys好 还是abaqus好啊
已经有10人回复
- 【求助】ANSYS分析应力过程中存在的问题
已经有3人回复
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
臭水沟: 金币+2, 谢谢交流~~ 2012-05-30 12:28:03
感谢参与，应助指数 +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
3楼2012-05-30 22:42:08
4楼2012-05-30 22:48:48
hegongda123
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 213
- 散金: 10
- 帖子: 100
- 在线: 26.4小时
- 虫号: 2071745
- 注册: 2012-10-19
- 专业: 机械摩擦学与表面技术
5楼2013-04-16 11:05:57
李军帅082104
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 987
- 散金: 135
- 帖子: 177
- 在线: 372.2小时
- 虫号: 2131364
- 注册: 2012-11-17
- 性别: GG
- 专业: 金属结构材料
6楼2013-08-18 08:43:17
7楼2016-03-21 15:33:03
8楼2017-04-25 08:08:47
