| 查看: 1050 | 回复: 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 |
回复此楼» 猜你喜欢
面上意见出来了
已经有21人回复
-
能否申诉？
已经有8人回复
-
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
开心就散金
+5/1720
-
云南大学丁军桥团队2027年博士招生
+1/945
-
八年奋战，面上中了，同时一个读了8年的博士（中途接手）答辩通过了，都是涉险过关
+2/342
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/88
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/87
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/87
-
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/85
-
深圳大学刘翼振课题组诚聘博士后（CRISPR分子诊断）
+1/36
-
同一份本子，四份夸技术，一份骂意义——国自然评审的坐标系差异
+1/34
-
山东大学-山东字节信息智能技术实验室2026年招聘博士后
+2/32
-
【科研诚信】武汉轻工大学何江玲3篇论文：涉嫌抄袭、成组数据重复、制造虚假显著
+1/15
-
中南大学化学化工学院喻桂朋教授团队诚聘博士后（复合材料方向）
+1/15
-
中国科学院 理化技术研究所 招 动力工程及工程热物理 相近的保送硕士 博士
+1/9
-
湖南大学张旺课题组招收研究助理
+1/7
-
多重蛋白检测定量技术赋能心血管疾病基础与转化研究
+1/4
-
中国科学院广州生物医药与健康研究院课题组招收博士研究生1名（2027年9月入学）
+1/4
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+1/3
-
项目助理招聘启事
+1/2
-
中科院院士成会明湖南大学招收分离膜、高分子、AI4M等方向推免直博生（10天内有效）
+1/2
-
香港大学智能制造实验室招收27年春秋季入学全奖博士和硕士研究生
+1/1
简单回复
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
回复
投票:
