01 January 1970 0 3K Report

I have completed the APDL as follows, but the moving heat source do not any work. Does anyone have idea about that? !%%%%�������ݽ���%%%%%%%%%%%%%%%%%%%%%%%%%%%% /clear,start /filename,Gaussian moving heat source /title,Gaussian moving heat source1_2D /PREP7 !�����ƶ���Դ��Gauss heat source !V=2mm/s,power=5w,d=4um !time=25s,tinc=esize/V=0.5s !the whole computational time whole_t=25 tinc=0.5 esize=1 V=2 ! Welding speed d=4 ! Effective heat arc radius R=d/2 LL=100 ! Weldment length WW=50 ! Weldment width HH=1 ! Weldment height MAX_TIME=51 MAX_X=101 MAX_Y=51 Q=5 FL=0.001 !the initial temperature Int_tem=20 dt=0.1 !%the number of steps n_step=whole_t/dt+1 !%%the material properties%% KK=0.025 CC=658 rho=7.6e-6 ET,1,PLANE55 MP,KXX,1,KK MP,C,1,CC MP,DENS,1,rho RECTNG,0,LL,0,WW !% mesh%%%% LSEL,S,LOC,y,0 LSEL,A,LOC,y,WW LESIZE,all,,,100,1, , , ,1 LSEL,S,LOC,x,0 LSEL,A,LOC,x,LL LESIZE,all,,,50,1, , , ,1 /PSYMB,LDIR,1 ASEL,ALL MSHKEY,1 AMESH,ALL ALLSEL,ALL *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)*ESIZE *ENDDO *DO,J,1,MAX_Y,1 LASER(0,J,K)=(J-1)*ESIZE *ENDDO *DO,I,1,MAX_X,1 *DO,J,1,MAX_Y,1 XCENTER=100-V*(K-1)*TINC YCENTER=25 DISTANCE=SQRT(((I-1)*ESIZE-XCENTER)**2+((J-1)*ESIZE-YCENTER)**2) *IF,DISTANCE,LE,R,THEN LASER(I,J,K)=Q*EXP(-DISTANCE**2/(R**2)) ! ����ı������� *ELSE LASER(I,J,K)=0 *ENDIF *ENDDO *ENDDO LASER(0,0,K)=(K-1)*TINC *ENDDO !!!!definite the initial temperature������������ TUNIF,Int_tem !!!!!the Down boundary condition!!!!!!!!!!!! LSEL,S,LOC,Y,0 SFL,ALL,HFLUX,FL ALLSEL,ALL !!!!the Down boundary condition���������������������������� !!!!!the Upper boundary condition!!!!!!!!!!!! LSEL,S,LOC,Y,WW SFL,ALL,HFLUX,FL ALLSEL,ALL !!!!the Upper boundary condition���������������������������� !!!!!the Right boundary condition!!!!!!!!!!!! LSEL,S,LOC,X,LL SFL,ALL,HFLUX,FL ALLSEL,ALL !!!!the Right boundary condition���������������������������� !!!!!the Left boundary condition!!!!!!!!!!!! NSEL,S,LOC,X,0 D,ALL,temp,Int_tem ALLSEL,ALL !!!!the Left boundary condition���������������������������� /SOL ANTYPE,4 TRNOPT,FULL ! Complete set method LUMPM,0 *do,t,1e-6,10,tinc time,t autots,off deltim,,,,1 !Use final time step from previous load step kbc,1 tsres,erase sf,ALL,hflux,%laser% solve outres,all,all *enddo outres,all,all finish

Similar questions and discussions