关于高斯移动热源的问题(急急急)

beckshjh · 发表于 2006-4-30 17:21

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

最近在做一个移动加热温度场分析的问题。由于我是个初学者,对于移动热源的加载还是搞不清楚，看到论坛里有很多关于高斯移动热源的问题，可是需要积分，本人初来乍到，身无分文，但是急需解决这个问题！ 我所遇到的问题是这样,在立方体钢块上做等离子熔覆(近似于焊接),热源选用等离子弧(类似于激光),立方体底部考虑绝热,其他面于空气接触.主要的问题是不清楚怎么施加荷载,以及边界条件等等.请各位高手指教,先谢了! <IMG src="http://www.simu-online.com/bbs/images/emot/em13.gif" border=0>

AaronSpark · 发表于 2006-5-1 07:40

数组方式加载 timend num_x=pf_en+1 num_y=phl_en+1 num_time=timend/tinc+1 *dim,flux3,table,num_x,num_y,num_time,x,y,time /com: timstep=pf_es/v !每个加热区单元的尺寸 /com: phl_en !加热区的单元数 radius=pfw_hr/2 !高斯热源作用的半径 *do,k,1,num_time,1 *do,j,1,num_y,1 !列y赋值 flux3(0,j,k)=(j-1)*pf_es *enddo *do,ii,1,num_x,1 ! !行x赋值 flux3(ii,0,k)=-0.5*(pf_en*pf_es)+(ii-1)*pf_es *enddo ydist=v*(k-1)*tinc *do,j,0,timend,tinc *do,i,1,num_x,1 !x坐标 ycoor=v*j xcoor=-0.5*(pf_en*pf_es)+(i-1)*pf_es r=sqrt((ycoor-ydist)**2+xcoor*xcoor) *if,r,le,radius,then gauss=qmax/exp(k*r*r/1000000) flux3(i,j/tinc+1,k)=gauss !*10e5 *endif *enddo *enddo flux3(0,0,k)=(k-1)*tinc !页time赋值 *enddo 在划分单元时要在加热区划的密一些，因为不管怎样加载，都是想模拟出一个圆形的区域，随着时间不断变化的。起初可能是半园，而后多半个，整园，当扫描稳定后都是整园。当然这个园是锯齿形的。 另外就是表定义要对其ROL,COL,PAGE进行赋值。循环赋值，当然哪一项作PAGE等，取决于你的程序编制。如是而已。模拟园的逼近程度取决于网格划分的好坏。当网格不适当时，可能会出现是椭圆的情况。所以网格尽量做到是正方形。

AaronSpark · 发表于 2006-5-1 07:45

推荐看看哈工大武传松的《焊接热过程数值分析》

AaronSpark · 发表于 2006-5-1 07:46

采用Table型数组,定义模型表面各个节点在不同时刻的载荷值,我是参考论坛上面以前的帖子做的. 详细的命令流如下: /UNITS,SI /CONFIG,NRES,10000 /PREP7 ET,1,SOLID70 MP,DENS,1,3200 MP,C,1,840 MP,KXX,1,16.8 /VIEW,1,1,1,1 BLOCK,0,0.04,0,0.04,0,0.004 SAVE LESIZE,1,,,40 LESIZE,4,,,40 LESIZE,9,,,4 MSHKEY,1 VATT,1,1,1 VMESH,1 /PNUM,DEFA EPLOT FINISH /SOLU LSIZE=0.001 V=0.05 TINC=LSIZE/V rb=0.002 pi=3.1415926 Q=600 Qmax=Q/(pi*rb*rb) R=0.002 MAX_TIME=41 MAX_X=41 MAX_Y=41 *DIM,laser,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME *DO,K,1,MAX_TIME,1 *DO,I,1,MAX_X,1 laser(I,0,K)=(I-1)*LSIZE *ENDDO *DO,J,1,MAX_Y,1 laser(0,J,K)=(J-1)*LSIZE *ENDDO *DO,I,1,MAX_X,1 *DO,J,1,MAX_Y,1 XCENTER=V*(K-1)*TINC YCENTER=0.02 DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2) *IF,DISTANCE,LE,R,THEN laser(I,J,K)=Qmax/EXP(2*DISTANCE**2/(0.002**2)) *ELSE laser(I,J,K)=0 *ENDIF *ENDDO *ENDDO laser(0,0,K)=(K-1)*TINC *ENDDO ANTYPE,TRANS TUNIF,27 KBC,1 *SET,DELT,1 TIMINT,ON AUTOTS,OFF TT=0.8 *DO,TM,1E-6,TT,TINC /VIEW,1,1,1,1 /ANG,1 /REP,FAST FLST,2,1,5,ORDE,1 FITEM,2,2 /GO !* TIME,TM AUTOTS,0 DELTIM, , , ,1 KBC,0 !* TSRES,ERASE FLST,2,1,5,ORDE,1 FITEM,2,2 /GO !* !* SFA,P51X,1,HFLUX,%laser% SOLVE OUTRES,ALL,ALL, /PSF,HFLUX,,2 /REPLOT *ENDDO OUTRES,ALL,ALL, FINISH

andy_yier · 发表于 2006-5-1 09:26

哇，我们做的东西很象的，看到与我联系好吗？讨论一下 170293641我的QQ

守护蓝天 · 发表于 2011-1-3 15:10

大家好，都做的怎么样了？
我现在也在做这方面的，不太懂，可以指教下吗？
谢谢大家！

守护蓝天 · 发表于 2011-1-3 15:11

回复 2 # AaronSpark 的帖子

你好，你是在圆形路径上加载移动热源的吗？具体是怎么实现的啊？可以指教下吗？谢谢！

账号		自动登录	找回密码
密码			我要加入

关于高斯移动热源的问题(急急急)

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（beckshjh）关于高斯移动热源的问题(急急急)

回复：（beckshjh）关于高斯移动热源的问题(急急急)

回复：（beckshjh）关于高斯移动热源的问题(急急急)