ansys求助——温度载荷及后处理

liangsflai

小弟想做一个温度应力的东西，但是温度载荷加不上，载荷函数为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

liangsflai

希望大神能帮助一下，如果结果得到结果，小弟会给予一定报酬的