有没有ANSYS温度应力分析的例题啊？
菜鸟求救~~~~~~~~~~~~~~~~~~~~~~~:hug: /filname,Q235/title,Q235WELDING 2-D FEA
/units,si
/prep7
et,1,plane55
et,2,surf151
keyopt,2,8,4
mptemp,1,20,250,500,750,1000,1500
mptemp,7,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35
mpdata,kxx,1,7,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350
mpdata,dens,1,7,7300,7090
mpdata,c,1,1,460,480,530,675,670,660
mpdata,c,1,7,780,820
mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45
mpdata,prxy,1,7,0.48,0.5
mpdata,alpx,1,1,1.10e-5,1.22e-5,1.39e-5,1.48e-5,1.34e-5,1.33e-5
mpdata,alpx,1,7,1.32e-5,1.31e-5
mpdata,ex,1,1,2.05e11,1.87e11,1.50e11,7.0e10,2.0e10,1.9e10
mpdata,ex,1,7,1.8e10,1.2e10
tb,biso,1,6,2
tbtemp,20
tbdata,,2.20e8,2.05e10
tbtemp,500
tbdata,,1.3e8,1.50e10
tbtemp,750
tbdata,,4.0e7,7.0e9
tbtemp,1000
tbdata,,2.50e7,2.0e9,
tbtemp,1700
tbdata,,1.0e6,1.8e9
tbtemp,2500
tbdata,,1.0e5,1.2e9,
pi=3.1415926
I=100 !焊接电流
U=23 !焊接电弧电压
effcient=0.85 !焊接热效率η
r=0.005 !电弧有效加热半径
Qm=3*I*U*effcient/(pi*r**2) !最大热流高斯密度
kk=0.85 !能量集中系数
Qmax=Qm*kk*0.25
ly=0.06
lx=0.3
lsize=0.005 !划分单元网格长度
v=0.004 !焊接速度
tm=lx/v !焊接总时间
tinc=lsize/v !单步时间长度
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
rectng,0,lx,0,ly
lesize,all,LSIZE
mshape,0,2d
mshkey,1
amesh,all
nsll,s
type,2
esurf
allsel,all
save
finish
!以下为输入高斯移动热源
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=Qmax/exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
finish
/solu
antype,trans
tintp,0.005,,,1,,, !设置向后积分
kbc,1
timint,off
autots,on
ESEL,S,TYPE,,2
SFE,ALL,,CONV,,100
SFE,ALL,,CONV,2,20 !表面效应单元施加对流系数
ALLSEL,ALL
tref,20
time,1e-6
d,all,temp,20
NSUBST,1
solve !进行稳态求解，得到初始条件
ddele,all,temp
timint,on
outres,all,all
tsres,erase
kbc,0
esel,s,type,,1
nsle,s
sf,all,hflux,%flux2%
allsel,all
time,73.75
autots,off
deltim,tinc,tinc,tinc
solve
sfdele,all,hflux
time,75
solve
time,725 !冷却过程到725s
autots,0
deltim,25,25,25
solve
finish
/prep7
etchg,tts !转换单元，进行结构分析
nsel,s,loc,x,0
nsel,r,loc,y,0
d,all,all
allsel,all
dl,1,,symm
finish
/solu
eresx,no !不选择积分点外推
tref,20
time,1e-6
LDREAD,TEMP,1,1,, ,'Q235','rth',' '
solve
*do,t,1.25,73.75,1.25
time,t
LDREAD,TEMP,,,t, ,'Q235','rth',' '
solve
*enddo
*do,i,1,26,1
time,50+i*25
LDREAD,TEMP,4,i,, ,'Q235','rth',' '
solve
*enddo
finish 这个例子不是太好 这例题就算了吧~ 楼上的两位想必是看过这个例子，说说它的优缺点呗，让大家知道下 这个例子虽然实现了应力场模拟的循环，但没有解决如何处理处于熔化状态的单元参与计算的问题，让处于液态的单元按照弹塑性理论计算肯定是不合理的。 我想问问,为什么要
time,1e-6
LDREAD,TEMP,1,1,, ,'Q235','rth',' '
solve time,1e-6
是设定的计算时间吧！ 先学习一下再说了，呵呵
页:
[1]