| 查看: 1052 | 回复: 12 | |||
[交流]
一个单层多道的例子，但是总是在进行一道之后，热源就不再移动了。
|
|||
|
找了一个单层多道的例子，但是总是在进行一道之后，热源就不再移动了。 /UNITS,SI /CONFIG,NRES,10000 /PREP7 ET,1,SOLID70 !单元类型，1表示此单元编号为1 MP,DENS,1,7930 !定义材料密度 MP,C,1,502 !定义材料比热 mptemp,1,20,100,500 mpdata,kxx,1,1,12.1,16.3,21.4 !kxx热导系数 /VIEW,1,1,1,1 !*************定义基板尺寸*********** lx=0.1 lz=0.1 ly=0.006 !基板大小 !*************定义焊接参数*********** layer=2 !总的焊接道数 lap=1/3 !两道之间搭接量 weld=0.01 !每道焊道宽度 wweld=layer*weld-(layer-1)*lap*weld !焊道的总宽度 wheight=0.004 !焊道的高度 remain=0.03 LSIZE=0.002 V=0.008 !速度 pi=3.1415926 U=25 I=180 Q=U*I !电源有效功率（在此，可以对Q通过UI来调整） R=0.006 !电源有效热半径 effect=0.8 Qmax=effect*Q/(pi*R*R) !中心处最大热流 tm=(lx-2*remain)/V !一道的总时间 tinc=LSIZE/V !每小段的时间间隔 !*****************建立模型，生成网格************** block,0,lx,0,ly,0,lz !创建一个块体积，可以根据基板尺寸调整坐标 wpoff,remain,ly,(lz-wweld)/2 !随着坐标系轴的变化平移工作平面 block,0,lx-2*remain,0,wheight,0,wweld vsbw,all,,delete !对block布尔运算 wpoff,,,wweld vsbw,all,,delete wprota,90,90 !转动工作面 vsbw,all,,delete wpoff,,,lx-2*remain vsbw,all,,delete save vglue,all !粘合所有的几何体 lsel,s,,,9,12,1 lsel,a,,,25,26,1 lsel,a,,,30 lsel,a,,,32 lsel,a,,,38 lsel,a,,,45 lsel,a,,,59 lsel,a,,,73 lesize,all,0.005,,,0.2,,,,1 allsel,all lsel,s,,,13,16,1 lsel,a,,,19,21,1 lsel,a,,,23,24,1 lsel,a,,,34 lsel,a,,,41,44,1 lsel,a,,,63,64,1 lsel,a,,,77,78,1 lsel,a,,,81,82,1 lesize,all,0.001,,,1,,,,1 allsel,all lsel,s,,,1 lsel,a,,,3,6,1 lsel,a,,,8 lsel,a,,,27 lsel,a,,,29 lsel,a,,,31 lsel,a,,,33 lsel,a,,,35 lsel,a,,,37 lsel,a,,,40 lsel,a,,,50 lsel,a,,,53 lsel,a,,,55 lsel,a,,,62 lsel,a,,,67 lsel,a,,,69 lsel,a,,,85 lesize,all,0.002,,,1,,,,1 allsel,all lsel,s,,,46 lsel,a,,,86 lsel,a,,,88 lesize,all,0.001,,,5,,,,1 allsel,all vsweep,all /PNUM,DEFA EPLOT fini !*******************进入后处理器********************* /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+wheight 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-2*remain)/LSIZE MAX_X=1+lx/LSIZE MAX_Z=1+lz/LSIZE !*************定义table数组 !****************************************** *DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME *do,n,1,layer,1 m=mod(n,2) !返回n/2的余数部分 *if,m,eq,1,then !n=1 *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=remain+V*(K-1)*tinc ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/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-weld/2,zcenter+weld/2 nsel,r,loc,y,ly,ly+wheight esln,s,0 ealive,all allsel,all !**********施加热流载荷********** esel,s,live nsel,s,loc,Y,ly sf,all,hflux,%heatflux% allsel,all !************求解***************** time,(n-1)*(tm+tinc)+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 *else !n=2 *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=lx-remain-V*(K-1)*tinc ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/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.001,xcenter+0.002 nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2 nsel,r,loc,y,ly,ly+wheight esln,s,0 ealive,all allsel,all !**********施加热流载荷********** esel,s,live nsel,s,loc,Y,ly sf,all,hflux,%heatflux% allsel,all !************求解***************** time,(n-1)*(tm+tinc)+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 *endif *enddo FINISH |
回复此楼» 猜你喜欢
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有12人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有8人回复
-
我不理解！
已经有24人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
国自然面上复盘~欢迎讨论
已经有19人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
» 本主题相关商家推荐: (我也要在这里推广)
简单回复
ansonal2楼
2014-05-26 09:19
回复
521majie(金币+1): 谢谢参与
baoer4103楼
2014-05-26 09:34
回复
521majie(金币+1): 谢谢参与
neu2344楼
2014-05-26 09:36
回复
521majie(金币+1): 谢谢参与
kingki5楼
2014-05-26 09:36
回复
521majie(金币+1): 谢谢参与
2014-05-26 09:39
回复
521majie(金币+1): 谢谢参与
7
假大空7楼
2014-05-26 09:39
回复
521majie(金币+1): 谢谢参与
521majie8楼
2014-05-26 09:45
回复
2014-05-26 09:55
回复
521majie(金币+1): 谢谢参与
henry83510楼
2014-05-26 10:08
回复
521majie(金币+1): 谢谢参与
铁观音茶叶11楼
2014-05-26 10:11
回复
521majie(金币+1): 谢谢参与
gaomeisdu12楼
2014-05-26 11:12
回复
521majie(金币+1): 谢谢参与
521majie13楼
2014-05-26 15:18
回复
20