关于高斯移动热源的问题(急急急)
最近在做一个移动加热温度场分析的问题。由于我是个初学者,对于移动热源的加载还是搞不清楚，看到论坛里有很多关于高斯移动热源的问题，可是需要积分，本人初来乍到，身无分文，但是急需解决这个问题！<BR>我所遇到的问题是这样,在立方体钢块上做等离子熔覆(近似于焊接),热源选用等离子弧(类似于激光),立方体底部考虑绝热,其他面于空气接触.主要的问题是不清楚怎么施加荷载,以及边界条件等等.请各位高手指教,先谢了! <IMG src="http://www.simu-online.com/bbs/images/emot/em13.gif" border=0> <BR>回复：（beckshjh）关于高斯移动热源的问题(急急急)
数组方式加载 timend <BR>num_x=pf_en+1<BR>num_y=phl_en+1<BR>num_time=timend/tinc+1<BR> <BR>*dim,flux3,table,num_x,num_y,num_time,x,y,time<BR> <BR>/com: timstep=pf_es/v !每个加热区单元的尺寸<BR>/com: phl_en !加热区的单元数<BR> <BR>radius=pfw_hr/2 !高斯热源作用的半径<BR> <BR>*do,k,1,num_time,1<BR> <BR> *do,j,1,num_y,1 !列y赋值<BR> <BR> flux3(0,j,k)=(j-1)*pf_es<BR> <BR> *enddo<BR> <BR> *do,ii,1,num_x,1 ! !行x赋值<BR> <BR> flux3(ii,0,k)=-0.5*(pf_en*pf_es)+(ii-1)*pf_es<BR> <BR> *enddo<BR> <BR> ydist=v*(k-1)*tinc<BR> <BR> *do,j,0,timend,tinc<BR> *do,i,1,num_x,1 !x坐标<BR> <BR> ycoor=v*j<BR> <BR> xcoor=-0.5*(pf_en*pf_es)+(i-1)*pf_es<BR> <BR> r=sqrt((ycoor-ydist)**2+xcoor*xcoor)<BR> <BR> *if,r,le,radius,then<BR> <BR> gauss=qmax/exp(k*r*r/1000000)<BR> flux3(i,j/tinc+1,k)=gauss !*10e5<BR> <BR> *endif<BR> <BR> *enddo<BR> *enddo<BR> <BR> flux3(0,0,k)=(k-1)*tinc !页time赋值<BR> <BR>*enddo <BR><BR>在划分单元时要在加热区划的密一些，因为不管怎样加载，都是想模拟出一个圆形的区域，随着时间不断变化的。起初可能是半园，而后多半个，整园，当扫描稳定后都是整园。当然这个园是锯齿形的。<BR>另外就是表定义要对其ROL,COL,PAGE进行赋值。循环赋值，当然哪一项作PAGE等，取决于你的程序编制。如是而已。模拟园的逼近程度取决于网格划分的好坏。当网格不适当时，可能会出现是椭圆的情况。所以网格尽量做到是正方形。回复：（beckshjh）关于高斯移动热源的问题(急急急)
推荐看看哈工大 武传松的《焊接热过程数值分析》回复：（beckshjh）关于高斯移动热源的问题(急急急)
采用Table型数组,定义模型表面各个节点在不同时刻的载荷值,我是参考论坛上面以前的帖子做的.<BR>详细的命令流如下:<BR>/UNITS,SI<BR>/CONFIG,NRES,10000<BR>/PREP7<BR>ET,1,SOLID70<BR>MP,DENS,1,3200<BR>MP,C,1,840<BR>MP,KXX,1,16.8<BR>/VIEW,1,1,1,1<BR>BLOCK,0,0.04,0,0.04,0,0.004<BR>SAVE<BR>LESIZE,1,,,40<BR>LESIZE,4,,,40<BR>LESIZE,9,,,4<BR>MSHKEY,1<BR>VATT,1,1,1<BR>VMESH,1<BR>/PNUM,DEFA<BR>EPLOT<BR>FINISH<BR>/SOLU<BR>LSIZE=0.001<BR>V=0.05<BR>TINC=LSIZE/V<BR>rb=0.002<BR>pi=3.1415926<BR>Q=600<BR>Qmax=Q/(pi*rb*rb)<BR>R=0.002<BR>MAX_TIME=41<BR>MAX_X=41<BR>MAX_Y=41<BR>*DIM,laser,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME<BR>*DO,K,1,MAX_TIME,1<BR>*DO,I,1,MAX_X,1<BR>laser(I,0,K)=(I-1)*LSIZE<BR>*ENDDO<BR>*DO,J,1,MAX_Y,1<BR>laser(0,J,K)=(J-1)*LSIZE<BR>*ENDDO<BR>*DO,I,1,MAX_X,1<BR>*DO,J,1,MAX_Y,1<BR>XCENTER=V*(K-1)*TINC<BR>YCENTER=0.02<BR>DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)<BR>*IF,DISTANCE,LE,R,THEN<BR>laser(I,J,K)=Qmax/EXP(2*DISTANCE**2/(0.002**2))<BR>*ELSE<BR>laser(I,J,K)=0<BR>*ENDIF<BR>*ENDDO<BR>*ENDDO<BR>laser(0,0,K)=(K-1)*TINC<BR>*ENDDO<BR>ANTYPE,TRANS<BR>TUNIF,27<BR>KBC,1<BR>*SET,DELT,1<BR>TIMINT,ON<BR>AUTOTS,OFF<BR>TT=0.8<BR>*DO,TM,1E-6,TT,TINC<BR>/VIEW,1,1,1,1<BR>/ANG,1<BR>/REP,FAST<BR>FLST,2,1,5,ORDE,1<BR>FITEM,2,2<BR>/GO<BR>!*<BR>TIME,TM<BR>AUTOTS,0<BR>DELTIM, , , ,1<BR>KBC,0<BR>!*<BR>TSRES,ERASE<BR>FLST,2,1,5,ORDE,1<BR>FITEM,2,2<BR>/GO<BR>!*<BR>!*<BR>SFA,P51X,1,HFLUX,%laser%<BR>SOLVE<BR>OUTRES,ALL,ALL,<BR>/PSF,HFLUX,,2<BR>/REPLOT<BR>*ENDDO<BR>OUTRES,ALL,ALL,<BR>FINISH<BR> 哇，我们做的东西很象的，看到与我联系好吗？讨论一下<BR>170293641我的QQ 大家好，都做的怎么样了？我现在也在做这方面的，不太懂，可以指教下吗？
谢谢大家！ 回复 2 # AaronSpark 的帖子
你好，你是在圆形路径上加载移动热源的吗？具体是怎么实现的啊？可以指教下吗？谢谢！
页:
[1]