本人在作间接法求解焊接热应力,按错误提示修改过,但还是会出现下面错误!不知如何解决?求高手告知!
*** 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 编辑 ] |