finally with my colleagues G.Ferrari and I.Lelli I finished the simulation work about friction stir welding...We are very excited because the simulation gave us the expected results. The maximum temperature does not exceed the melting of the material, the pattern of isotermee is congruent with the results shown in the literature and various checks seem to justify the work done ...so now we are waiting the verdict of the professor with fingers crossed...
fabrizio leone
I am working on friction
I am working on friction stir welding process in that i have to move
tool shape heat source in plate to find out temperature distribution in
plate . I have make circular shape heat source on the surface of the
plate but i want to make cylindrical shape pin as a heat source in that.
presenyly i am using dflux subroutine pls help me.
that is my soubroutine
SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,
1 TEMP,PRESS,SNAME)
INCLUDE 'ABA_PARAM.INC'
DIMENSION FLUX(2),TIME(2),COORDS(3)
CHARACTER*80 SNAME
integer DZ
REAL V,Q1,Q2,Q3,CX,CY,DX,DY,r
V = 0.002
Q1 = 4640014.831
Q2 = 2000000
Q3 = 1202051
CX = 0.05
CY = V*TIME(2)
DX = COORDS(1)-CX
DY = COORDS(2)-CY
r = sqrt(DX**2 + DY**2)
!
IF(COORDS(3).LT.2e-04) THEN
IF ((r.ge.0.004).AND.(r.le.0.0125))THEN
flux(1)=Q1*r
ENDIF
ENDIF
IF((COORDS(3).GT.2e-04).AND.(COORDS(3).LT.41E-04)) THEN
IF ((r.ge.0.0039).AND.(r.le.0.0040))THEN
flux(1)=Q2*r
ENDIF
ENDIF
IF((COORDS(3).GT.41E-04).AND.(COORDS(3).LT.43E-04)) THEN
IF (r.le.0.0040)THEN
flux(1)=Q3*r
ENDIF
ENDIF
RETURN
END
i am able to restrict pin to some length