ansys求助——温度载荷及后处理
小弟想做一个温度应力的东西，但是温度载荷加不上，载荷函数为5+34sin（t*3.14/6）希望大神给予帮助，下面是我建模及划分网格的一些命令，在此先谢谢了！
/prep7
et,1,42
mp,dens,1,5
mp,ex,1,2.7e10
mp,prxy,1,0.2
mp,dens,2,2.7
mp,ex,2,2.5e10
mp,prxy,2,0.2
mp,dens,3,2.2
mp,ex,3,2.1e10
mp,prxy,3,0.2
材料的温度属性
/prep7
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,2.37
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,2.37
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,2.37
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,3,REFT,,,
MPDATA,ALPX,3,,1e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,
MPDATA,ALPX,2,,1e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,
MPDATA,ALPX,1,,1e-5
csys,1
xmin=4900
xmax=5400
ymin=84.5
ymax=95.5
!定义骨料个数
num=68
!存放骨料位置及半径的数组
*dim,cir,array,num,3
cum=0
*do,i,1,10000000000000000
*if,cum,eq,num,then
*exit
*endif
x=rand(xmin,xmax)
y=rand(ymin,ymax)
*if,cum,lt,10,then
r=60
*elseif,cum,gt,9,and,cum,lt,40,then
r=30
*else,cum,gt,39,then
r=16
*endif
*if,x-r,gt,xmin,and,x+r,lt,xmax,then
*if,cum,eq,0,then
cum=cum+1
cir(cum,1)=x
cir(cum,2)=y
cir(cum,3)=r
*else
sum=0
*do,j,1,cum
*afun,deg
dist=sqrt((cir(j,1))**2+x**2-2*cir(j,1)*x*cos(cir(j,2)-y))
*if,dist,lt,cir(j,3)+r+10,then
*exit
*else
sum=sum+1
*endif
*enddo
*if,sum,eq,cum,then
cum=cum+1
cir(cum,1)=x
cir(cum,2)=y
cir(cum,3)=r
*endif
*endif
*endif
*enddo
!定义三个数组分别保存x,y,r
*dim,ncirx,array,cum
*dim,nciry,array,cum
*dim,ncirr,array,cum
*do,i,1,cum
ncirx(i)=cir(i,1)
nciry(i)=cir(i,2)
ncirr(i)=cir(i,3)
i=i+1
*enddo
cyl4,,,4900,84,5400,96
lesize,1,,,150
lesize,3,,,150
lesize,2,,,110
lesize,4,,,110
amesh,all
*do,i,1,16500 !(i代表的是单元)
k=0 !（用于计数判断材料属性如何定义）
*do,j,1,4 !(j代表的是单元上的8个点)
*do,s,1,68
*if,s,eq,cum,then
*exit
*endif
ndist=sqrt(nx(nelem(i,j))**2+ncirx(s)**2-2*cos(ny(nelem(i,j))-nciry(s))*nx(nelem(i,j))*ncirx(s))
*if,ndist,lt,ncirr(s),then
k=k+1
*endif
*enddo
*enddo
*if,K,eq,4,then
emodif,i,mat,3
*elseif,k,ne,0,and,k,ne,8,then
emodif,i,mat,2
*endif
*enddo
!加约束
dl,2,,ux
dl,2,,uy
dl,4,,ux
dl,4,,uy
!载荷
acel,,-9.8
sfl,1,pres,8.6
sfl,2,pres,1.5
sfl,4,pres,1.5
希望大神能帮助一下，如果结果得到结果，小弟会给予一定报酬的
页:
[1]