| 查看: 673 | 回复: 0 | |||
[求助]
ansys 焊接移动热源的施加
|
|
以下命令流中施加的热源无法移动，不知是哪里的错误！！望各位给指导指导！！！！ finish /clear,nostart /filname,1.11thermal_analysis /title,3D_weld_thermal_analysis /units,si /prep7 len=100e-3 wid=50e-3 thic=6e-3 b=0.002 h=0 wweld=15e-3 lnum=10 I=200 U=30 v=0.01 arf=60 pi=3.1415926 eta=0.7 K=1.2e4 themp=25 qmax=0.24*k*eta*I*u/pi percent=0.99 radius=abs(sqrt(1/k*log(1/(1-percent)))) et,1,solid70 mptemp,1,0,100,200,300,400,500 mptemp,,1200,1250,1350,1400,1450,1500 mptemp,,1505,1705,1905,2105,2305,2500 mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 mpdata,kxx,1,,30,30,30,30,30,30 mpdata,kxx,1,,30,30,30,30,30,30 mpdata,c,1,1,450,469,481,508.5,536,569 mpdata,c,1,,700,2172,5116,6589,8061,9533 mpdata,c,1,,9533,7757,5982,4206,2431,700 *afun,deg K,1,0,0,0 K,2,wweld,0,0 K,3,wid,0,0 K,4,wid,thic,0 K,5,thic*tan(arf/2)+wweld,thic,0 K,6,thic*tan(arf/2),thic,0 K,7,0,0,len A,1,3,4,6 A,1,2,5,6 cyl4,0,0,thic/cos(arf/2),arf,0,90 cyl4,0,0,thic/cos(arf/2)*sqrt(2)/2,arf,0,90 L,1,7 Vdrag,all,,,,,,14 Vptn,all Vglue,all Nummrg,all aadd,24,32 lccat,40,49 lccat,42,52 lsel,s,line,,17,21,2 lsel,a,line,,24,26,2 lsel,a,line,,35,36 lsel,a,line,,14,30,16 lesize,all,2e-3 lsel,s,line,,6,25,19 lsel,a,line,,2,18,16 lesize,all,,,6 lsel,s,line,,41,43,2 lsel,a,line,,50,53,3 lesize,all,1.5e-3 lsel,s,line,,45,47,2 lsel,a,line,,44,46,2 lesize,all,,,18,0.8 lsel,s,line,,7,27,20 lesize,all,,,12,0.7 lsel,a,line,,5,23,18 lesize,all,,,12,10/7 lsel,s,line,,11,34,23 lsel,a,line,,48,51,3 lesize,all,,,4 alls mshkey,1 mshape,0,3d vsel,all vatt,,,1 aslv,s aatt,,,1 type,1 vmesh,all alls /view,,-3,2,-5 v1=5 v2=7 a1=18 a2=31 esel,s,type,,1 *get,nemax,elem,,num,max *get,nemin,elem,,num,min vsel,s,,,v2 eslv ekill,all esel,s,live eplot nsel,s,loc,x,0 dsym,symm,x alls /solu antype,trans trnopt,full nropt,full,,on pred,on timint,on tref,25 t=0 dt0=0.001 dt1=0.501 tinc=1 time,dt0 timint,off kbc,1 nsubst,1 esel,s,type,,1 nsle,s ic,all,temp,25 alls outres,all,all solve finish esel,s,type,,1 esel,r,live eplot /prep7 *do,im,0,lnum,1 esel,s,type,,1 esel,r,live vsel,s,,,v1 eslv,r nsle asel,s,,,a1 nsla,r,1 esln,r,0 c=V*tinc*im b=thic/cos(arf/2)*sqrt(2)/2 tm=dt1+im*tinc /solu time,tm antype,4,rest timint,on autots,on kbc,1 deltim,0.1,0.01,0.1 *do,i,nemin,nemax,1 *if,i,le,nemax,then xsy=centrx(i) ysy=centry(i) zsy=centrz(i) rr=abs(sqrt(xsy*xsy+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c))) qr=qmax*exp(-k*rr*rr) *if,rr,le,radius,then sfe,i,nmface(i),hflux,,qr *else sfe,i,nmface(i),hflux,,0 *endif *endif *enddo allsel outres,all,all solve alls esel,s,type,,1 esel,r,live vsel,s,,,v1 eslv,r nsle asel,s,,,a1 nsla,r,1 esln,r,0 *do,i,1,6,1 sfedele,all,i,hflux *enddo *enddo alls esel,s,live eplot@wuming524 |
回复此楼» 收录本帖的淘帖专辑推荐
 ANSYS命令流 |
» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复