!应力场分析:
/CONFIG,NRES,1000000
/COM, Structural !打开结构分析
/units,SI !单位设置为国标
/PREP7 !前处理器
ET,1,SOLID45 !单元属性
!材料属性
MPTEMP,1,20,200,600,800,1000,1200,1374,1400,1460 !绘制温度曲线
MPDATA,KXX,1,1,13.31,16.33,22.38,25.07,27.53,29.76,30.1,31.95,320 !导热系数
MPDATA,C,1,1,470,508,592,634,676,719,740,12000,765 !比热
MP,DENS,1,8930 !密度
MP,EX,1,1.80e11 !弹性模量
MP,ALPX,1,15e-6 !热膨胀系数
MP,PRXY,1,0.288 !泊松比
!建模
BLOCK,-0.2,0.2,-0.002,0.002,0,0.003,
BLOCK,-0.2,0.2,0.002,0.1,0,0.003,
BLOCK,-0.2,0.2,-0.1,-0.002,0,0.003,
VGLUE,ALL
!划分焊缝网格
LESIZE,2, , ,50, , , , ,1 !线2划分50份
LESIZE,4, , ,50, , , , ,1 !线4划分50份
LESIZE,5, , ,50, , , , ,1 !线5划分50份
LESIZE,7, , ,50, , , , ,1 !线7划分50份
LESIZE,1, , ,1, , , , ,1 !线1划分1份
LESIZE,3, , ,1, , , , ,1 !线3划分1份
LESIZE,6, , ,1, , , , ,1 !线6划分1份
LESIZE,8, , ,1, , , , ,1 !线8划分1份
LESIZE,9, , ,1, , , , ,1 !线9划分1份
LESIZE,10, , ,1, , , , ,1 !线10划分1份
LESIZE,11, , ,1, , , , ,1 !线11划分1份
LESIZE,12, , ,1, , , , ,1 !线12划分1份
LESIZE,37, , ,20, , , , ,1 !线37划分20份
LESIZE,38, , ,20, , , , ,1 !线38划分20份
LESIZE,39, , ,20, , , , ,1 !线39划分20份
LESIZE,40, , ,20, , , , ,1 !线40划分20份
LESIZE,41, , ,20, , , , ,1 !线41划分20份
LESIZE,42, , ,20, , , , ,1 !线42划分20份
LESIZE,43, , ,20, , , , ,1 !线43划分20份
LESIZE,44, , ,20, , , , ,1 !线44划分20份
VSWEEP,1 !体扫略方式划分网格
!划分母材网格
ALLSEL
VSWEEP,4
VSWEEP,5
!将所有单元选中
ii=0
*DIM,NE,,50
*do,i,1,50
*if,esel(i),eq,1,then
ii=ii+1
ne(ii)=i
*endif
*enddo
*DIM,NC1,,50,1
*DO,I,1,50,1
NC1(I)=NE(I)
*ENDDO
!杀死焊缝单元
allsel,all
nropt,full !使用完全的Newton-Raphson,没有这个命令会报错
ESTIF,1e-056, !以前采取默认值1E-6造成很大的误差,单元好像并没有杀掉!
*DO,I,1,50
ekill,NC1(i)
esel,s,live
*ENDDO
*get,nlive,elem,,num,max
/SOLU
TIMINT,1,STRUCT
TIMINT,0,THERM
TIMINT,0,MAG
TIMINT,0,ELECT
TIMINT,0,FLUID
!设置位移约束
NSEL,ALL
NSEL,S,LOC,Y,0.08,0.085 !选择Y坐标大于0.08小于0.085的节点
NSEL,A,LOC,Y,-0.08,-0.085 !选择Y坐标大于-0.08小于-0.085的节点
D,ALL,ALL,0 !把选中的节点约束
ALLSEL,ALL
ANTYPE,4 !瞬态分析
NROPT,FULL !使用完全的Newton-Raphson
tunif,20 !初始温度
tref,20 !参考温度
kbc,1 !阶跃载荷,必须考虑相应时间内产生的温度梯度效应
tim=0
!一个点加载后马上将其删掉,并移到下一个点进行加载
!每层载荷循环的个数为规则单元的个数
*do,j,1,50,1
esel,all
ealive,NC1(j)
esel,s,live
esel,all
*get,nlive1,elem,,num,max
eplot
ldread,temp,,,tim, ,'file','rth',' ' !读取温度载荷作为体载荷
allsel
outres,all,all
tim=tim+2
time,tim
nsubst,1
autots,on
solve
save
ddele,nc1(j),temp,
allsel
*enddo
ldread,temp,,,1000, ,'file','rth',' ' !读取温度载荷作为体载荷
time,1000
solve
finish |