应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys焊接模拟过程中出现的错误
51/1返回列表
查看: 2716  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

mzgnjust

金虫 (小有名气)

[求助] 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
回复此楼

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

焊接有限元分析

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-05-28 20:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

马小英

新虫 (小有名气)
楼主你好！在你每道焊缝的加载热源的过程中，质心坐标是怎么计算出来的？谢谢！
赞一下 回复此楼
选择大于一切
7楼2016-03-21 15:33:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

pengqiong

新虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与，应助指数 +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

金虫 (小有名气)
引用回帖:
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的回帖

mzgnjust

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

不好意思，不知道怎么贴图，我跟你描述一下吧，在第二个循环以后激活的焊缝单元上面一直是杀死状态的温度（室温），然后母材有部分在焊接半径之内的单元还是有相应的温度，就是第二个循环以后的焊缝处没有温度，不知道是不是单元激活的方法或者选单元的方法上有问题？
赞一下 回复此楼
4楼2012-05-30 22:48:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭