24小时热门版块排行榜    

查看: 1922  |  回复: 0

Liv_Guo

新虫 (初入文坛)

[求助] ANSYS焊接移动热源APDL

ANSYS焊接高斯移动热源,需要导入模型进行模拟,写的APDL结果都不行,求助~~~
附上我改的APDL,能找错也可以~最好帮忙提供一份能根据导入模型的APDL
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
MP,DENS,1,7930
MP,C,1,502
mptemp,1,20,100,500
mpdata,kxx,1,1,12.1,16.3,21.4
/VIEW,1,1,1,1
pi=3.1415926
i=160
u=15.7
efficient=0.8
r=0.03
Qm=3*i*u*efficient/(pi*r*r)
kk=0.85
Qmax=Qm*kk
ly=0.4
lx=3
lz=0.025
lsize=0.09
v=0.0018
tm=lx/v
tinc=lsize/v

/solu
antype,trans,new
!nlgeom,on
timint,0,struct
timint,1,therm
timint,0,mag
timint,0,elect
tref,25
nropt,auto
autos,on
kbc,0
pred,on
lnsrch,on

nsel,s,loc,y,ly,ly+0.004
esln,s,1
ekill,all
allsel,all

esel,s,live
eplot

esel,s,live
nsel,s,loc,z,0
nsel,a,loc,z,lz
nsel,a,loc,x,0
nsel,a,loc,x,lx
sf,all,conv,10,25

MAX_TIME=1+tm/tinc
MAX_X=1+lx/lsize
MAX_Z=1+lz/lsize

*DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
        *DO,K,1,MAX_TIME,1
        *DO,I,1,MAX_X,1
        heatflux(I,0,K)=(I-1)*lsize
        *ENDDO
        *DO,J,1,MAX_Z,1
        heatflux(0,J,K)=(J-1)*lsize
        *ENDDO
        *DO,I,1,MAX_X,1
        *DO,J,1,MAX_Z,1
        XCENTER=v*(K-1)*tinc
        ZCENTER=lz/2
        DISTANCE=SQRT(((I-1)*lsize-XCENTER)**2+((J-1)*lsize-ZCENTER)**2)
        *IF,DISTANCE,LE,r,THEN
        heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2))
        *ELSE
        heatflux(I,J,K)=0
        *ENDIF
        *ENDDO
        *ENDDO
        heatflux(0,0,K)=(K-1)*tinc
        esel,s,live
        nsel,s,loc,x,xcenter-0.002,xcenter+0.001
        nsel,r,loc,z,zcenter-0.005,zcenter+0.005
        nsel,r,loc,y,ly,ly+0.004
        esln,s,0
        ealive,all
        allsel,all
        esel,s,live
        nsel,s,loc,y,ly
        sf,all,hflux,%heatflux%
        allsel,all
        time,K*lsize/v
        nsubst,2
        SOLVE
        nsel,s,loc,y,ly
        sfdele,all,hflux
        allsel,all

        esel,s,live
        eplot

        OUTRES,ALL,ALL,
        /PSF,HFLUX,,2
        /REPLOT
        OUTRES,ALL,ALL,
        *ENDDO
FINISH
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Liv_Guo 的主题更新
信息提示
请填处理意见