论坛 › C04：ANSYS--流体软件（Fluen\CFX\Icepark) › 间接法求解焊接热应力

liujin2000161 发表于 2006-6-18 16:34:47

间接法求解焊接热应力

本人在作间接法求解焊接热应力，按错误提示修改过，但还是会出现下面错误！不知如何解决？求高手告知！
*** ERROR ***                           CP =      78.333   TIME= 16:08:53
One or more elements have become highly distorted.Excessive         
distortion of elements is usually a symptom indicating the need for   
corrective action elsewhere.Try incrementing the load more slowly   
(increase the number of substeps or decrease the time step size).You
may need to improve your mesh to obtain elements with better aspect   
ratios.Also consider the behavior of materials, contact pairs,      
and/or constraint equations.If this message appears in the first      
iteration of first substep, be sure to perform element shape checking.
源程序
FINISH
/CLEAR,START
/BATCH
/CWD,'d:\Ansyswork\EBM FEM Simulation for Metals\3D Temperture distrbution'
/FILNAME,MatCha_db,0
/TITLE,Change material performance db Method
!选择分析问题的类型：应力
!结构计算，在计算前应删除表面载荷并转换单元类型
/UNITS,SI
/PREP7                                                   !进入前处理
!0参数定义
!模型参数
!几何模型尺寸me-3=mm
*SET,Block_x,7e-3   
*SET,Block_y,3e-3
*SET,Block_z,0.2e-3
!划分网格的大小
*SET,LMeshx,1e-4!line1
*SET,LMeshy,1e-4!line4
*SET,LMeshz,1e-4!line9
!边界条件参数
*SET,Tinit,30      !初始温度
!1定义单元类型
ET,1,SOLID45       !or solid185
!2实常数
!Nothing
!3材料性能
toffst,273    !从绝对温度改为摄氏度

*DO,MNUM,1,2
MPTEMP
MPTEMP,1,24,400,1200,1600,2000,2500
MPDATA,DENS,MNUM,1,3983,4688.4,5406.8,5948.8,6908,6331          !密度
MPDATA,KXX,MNUM,1,1.6278,2.3812,3.6396,29.8,30.6,32      !热导率
mpdata,c,MNUM,1,472,535,679,859,945,1034                !比热

MPTEMP,,,,,,,,
MPTEMP,1,0,200,400,800,1200
MPTEMP,6,1400,1454,1800,2000
MPDATA,ENTH,MNUM,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9
MPDATA,ENTH,MNUM,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10        !热焓

MPTEMP,1,20,300,600,750,900,1200,1500,1800,2000
MPDATA,HF,MNUM,1,6,50,120,180,200,250,378,700,850 !对流系数

mp,murx,MNUM,1                  !***输入磁场相对渗透系数
mp,perx,MNUM,1                  !介质常数(PERY,PERZ)

MPTEMP,,,,,,,,
MPTEMP,1,20,200,400,600,800
MPTEMP,6,900,1000,1100,1200,1420,1460,3500
MPDATA,ALPX,MNUM,1,15.24e-6,16.43e-6,17.44e-6,18.21e-6,18.83e-6         !热膨胀系数
MPDATA,alpx,MNUM,6,19.11e-6,19.38e-6,19.66e-6,19.95e-6,20.7,20.7,20.7

MPTEMP
MPTEMP,1,20,100,200,300,400
MPTEMP,6,500,600,700,800,1000   
MPTEMP,11,1200,1400,3500
MPDATA,EX,MNUM,1,2.04E+11,1.96E+11,1.87E+11,1.81E+11,1.71E+11!   
MPDATA,EX,MNUM,6,1.6E+11,1.52E+11,1.44E+11,1.29E+11,6E+10
MPDATA,EX,MNUM,11,2E+10,1E+10,1E+10         !杨氏（弹性）模量


MPTEMP
MPTEMP,1,20,200,400,600,800
MPTEMP,6,900,1000,1100,1200,1420,1460,3500
MPDATA,PRXY,MNUM,0.267,0.29,0.322,0.296,0.262,0.24,0.229,0.223,0.223,0.223,0.223,0.223 !泊松比

!**应力应变求解设置**
TB,MISO,MNUM,13,4                ! Activate a data table
TBTEMP,20                   ! Temperature = 20
TBPT,DEFI,9.5931e-4,1.957e8    ! Strain, stress at temperature =2 0
TBPT,DEFI,2E-3,        2.86e8                !原来是2.06
TBPT,DEFI,1E-2,        3.11e8
TBPT,DEFI,0.2,        5.2e8

TBTEMP,100                   ! Temperature =100
TBPT,DEFI,7.4643e-4,1.463e8    ! Strain, stress at temperature =100
TBPT,DEFI,2E-3,        2.04e8                !原来是1.54
TBPT,DEFI,1E-2,        2.52e8
TBPT,DEFI,0.2,        4.47e8

TBTEMP,200                   ! Temperature =200
TBPT,DEFI,6.5535e-4,1.2255e8    ! Strain, stress at temperature =200
TBPT,DEFI,2E-3,        1.54e8               
TBPT,DEFI,1E-2,        2.19e8
TBPT,DEFI,0.2,        4.00e8


TBTEMP,300                   ! Temperature =300
TBPT,DEFI,6.0359e-4,1.0925e8    ! Strain, stress at temperature =300
TBPT,DEFI,2E-3,        1.15e8
TBPT,DEFI,1E-2,        2.09e8
TBPT,DEFI,0.2,        3.97e8

TBTEMP,400                   ! Temperature =400
TBPT,DEFI,5.8889e-4,1.007e8   ! Strain, stress at temperature =400
TBPT,DEFI,2E-3,        1.06e8
TBPT,DEFI,1E-2,        2.03e8
TBPT,DEFI,0.2,        3.97e8

TBTEMP,500                   ! Temperature =500
TBPT,DEFI,5.9375e-4,9.5e7    ! Strain, stress at temperature =500
TBPT,DEFI,2E-3,        1.00e8
TBPT,DEFI,1E-2,        1.90e8
TBPT,DEFI,0.2,        3.71e8

TBTEMP,600                   ! Temperature =600
TBPT,DEFI,5.8125e-4,8.835e7    ! Strain, stress at temperature =600
TBPT,DEFI,2E-3,        9.30e7
TBPT,DEFI,1E-2,        1.66e8
TBPT,DEFI,0.2,        3.11e8

TBTEMP,700                   ! Temperature =700
TBPT,DEFI,5.3438e-4,7.695e7    ! Strain, stress at temperature =700
TBPT,DEFI,2E-3,        8.10e7
TBPT,DEFI,1E-2,        1.27e8
TBPT,DEFI,0.2,        2.18e8

TBTEMP,800                   ! Temperature =800
TBPT,DEFI,4.4186e-4,5.7e7    ! Strain, stress at temperature =800
TBPT,DEFI,2E-3,        6e7
TBPT,DEFI,1E-2,        7.63e7
TBPT,DEFI,0.2,        1.09e8


TBTEMP,1000                   ! Temperature =1000
TBPT,DEFI,4.75e-4,2.85e7    ! Strain, stress at temperature =1000
TBPT,DEFI,2E-3,        3e7
TBPT,DEFI,1E-2,        4e7
TBPT,DEFI,0.2,        6e7

TBTEMP,1200                   ! Temperature =1200
TBPT,DEFI,4.75e-4,9.5e6    ! Strain, stress at temperature =1200
TBPT,DEFI,2E-3,        10e6
TBPT,DEFI,1E-2,        13.3e6
TBPT,DEFI,0.2,        2e7

TBTEMP,1400                   ! Temperature =1400
TBPT,DEFI,4.75e-4,4.75e6    ! Strain, stress at temperature =1400
TBPT,DEFI,2E-3,        5e6
TBPT,DEFI,1E-2,        6.67e6
TBPT,DEFI,0.2,        1e7

TBTEMP,3500                   ! Temperature =2000
TBPT,DEFI,4.75e-4,4.75e6    ! Strain, stress at temperature =3000
TBPT,DEFI,2E-3,        4.55e6
TBPT,DEFI,1E-2,        6.e6
TBPT,DEFI,0.2,        0.9e7
*ENDDO       
mp,reft,1,24
mp,reft,2,1400
!4几何模型
BLOCK,0,Block_x,0,Block_y,0,Block_z,
!5网格划分控制
!ESIZE,1e-4,
lesize,1,LMeshx,
lesize,4,LMeshy,
lesize,9,LMeshz,

MSHKEY,1
VMESH,ALL


*set,F,10e-3            !10mm/s
*SET,dt,LMeshx/F          !时间步长 即每步加热的时间       

/prep7
etchg,stt

ldread,temp,,,,,,rth

/post1
*dim,TB,char,120
TB(1)='TB1','TB2','TB3','TB4','TB5','TB6','TB7','TB8','TB9','TB10'
!TB(11)='TB11','TB12','TB13','TB14','TB15','TB16','TB17','TB18','TB19','TB20'
!TB(21)='TB21','TB22','TB23','TB24','TB25','TB26','TB27','TB28','TB29','TB30'
!TB(31)='TB31','TB32','TB33','TB34','TB35','TB36','TB37','TB38','TB39','TB40'
!TB(41)='TB41','TB42','TB43','TB44','TB45','TB46','TB47','TB48','TB49','TB50'
!TB(51)='TB51','TB52','TB53','TB54','TB55','TB56','TB57','TB58','TB59','TB60'
!TB(61)='TB61','TB62','TB63','TB64','TB65','TB66','TB67','TB68','TB69','TB70'
!TB(71)='TB71','TB72','TB73','TB74','TB75','TB76','TB77','TB78','TB79','TB80'
!TB(81)='TB81','TB82','TB83','TB84','TB85','TB86','TB87','TB88','TB89','TB90'
!TB(91)='TB91','TB92','TB93','TB94','TB95','TB96','TB97','TB98','TB99','TB100'
!TB(101)='TB101','TB102','TB103','TB104','TB105','TB106','TB107','TB108','TB109','TB110'
!TB(111)='TB111','TB112','TB113','TB114','TB115','TB116','TB117','TB118','TB119','TB120'
*do,J,1,120
set,J,last,1
etable,TB(J),temp
*enddo
finish

/VIEW,1,1,2,3   
/ANG,1
/REP,FAST

/prep7
etchg,tts

/solu
ANTYPE,TRANS
nlgeom,on
tref,24
cnvtol,u,,0.5
CUTCONTROL,PLSLIMIT, 10,
lnsrch,on
nropt,full
neqit,2500
EQSLV,ICCG,,2
kbc,0       
TIMINT,ON,STRUC
autots,on
outres,all
outpr,all,last
NSUBST,10,25,5
!NSUBST,1,1,1

local,11,0,0,0,,,,,,
nsel,s,loc,z,0
nsel,r,loc,y,Block_y/2
nsel,r,loc,x,Block_x/2
d,all,all,0
allsel                                        !约束自由度

/solu
*set,tim,0
*set,i,0
!进入主循环
*DO,Y,1,3         !三道扫描线
*do,X,1,40       !扫描线长度5mm

tim=tim+dt
i=i+1
time,tim

allsel
ESEL,S,ETAB,TB(i),1395
mpchg,2,all
ekill,all                        !杀高温单元
esel,inve
ealive,all
eplot
allsel

ldread,temp,i,,,,,rth!读入时间tim时的温度分布 ,（i步时）
/pbf,temp,,1                !Shows magnitude of body force loads on displays.
solve

eplot
/auto,1
/replot

*enddo
*ENDDO

save

[ 本帖最后由 liujin2000161 于 2006-8-12 12:09 编辑 ]

半个苹果 发表于 2006-6-21 21:52:38

55555

我也遇到此问题    希望高手给个解决的步骤

半个苹果 发表于 2006-6-21 21:53:27

！！！

是不是应该首先检查一下网格形状呀

liujin2000161 发表于 2006-6-21 22:49:07

我后来把步长改为25就可以了！但是在第三道有出现这个警告，计算不下去了！

[ 本帖最后由 liujin2000161 于 2006-8-12 12:01 编辑 ]

liujin2000161 发表于 2006-8-12 11:59:32

自己顶一下！

xiguaxiaojie 发表于 2006-8-15 21:08:45

liujin2000161，半个苹果，
你们好，我在计算到第二步时也遇到同样的问题，不知你们现在解决的怎么样了，是不是只要增加子步数就可以了呢？谢谢告知。

liujin2000161 发表于 2006-8-16 12:40:24

就算增加到100，也还会出现这个问题阿！郁闷中

xiguaxiaojie 发表于 2006-8-16 20:57:52

呵呵，是啊，我现在子步数已经是60了，增大到120，不知能不能得到解决，那要计算多久啊，急！

liujin2000161 发表于 2006-8-17 13:37:39

看来得修改一下模型阿,重新划分网格

后起之秀 发表于 2006-8-21 09:36:27

不杀死高温单元，相应地修改高温时的材料参数

liujin2000161 发表于 2006-8-21 10:58:39

有液相产生阿！修改材料参数，更难收敛！

liujin2000161 发表于 2006-9-8 09:40:06

顶起来！

zxw19851015 发表于 2009-4-3 15:43:14

论坛里面高手很多,我只有学习的份

acailaogu 发表于 2009-4-3 16:15:56

增加子布、，加密网格！！！！

hanxue 发表于 2009-4-3 23:25:44

还可以试试修改一些非线性选项的默认值，还有就是网格质量会影响计算结果的收敛。

wangwenna218 发表于 2009-4-21 16:15:49

我也遇到过，不知大家解决了吗

guotao102 发表于 2009-4-22 14:06:50

楼主，不知道呢问题现在解决了没？

guotao102 发表于 2009-4-22 14:07:13

我现在也有这样的问题，加我QQ，我们讨论下

guotao102 发表于 2009-4-22 14:07:41

忘记说QQ号了，522029871

飘叶和流星 发表于 2009-5-18 11:07:42

10# 后起之秀
你好
想问您关于焊接残余应力方面的问题
http://forum.simwe.com/thread-880314-1-1.html
这是我遇到的一个问题
我的QQ401021563
欢迎指导交流

查看完整版本: 间接法求解焊接热应力