---

Hallo zusammen,

ich versuche zurzeit einen Laser zu simulieren, was mir allerdings noch nicht ganz gelingt.
Wenn ich in meiner Subroutine die Effizienz von 1 ändere bekomm ich den Fehler zero heat flux. Ich weiß leider überhaupt nicht wo da das Problem liegt, da ich blutige Anfängerin bin. Die Subroutine  und den Inputfile teilweise hier drunter hereinkopiert. Ich hoffe sehr mir kann jemand helfen.

Lg

Subroutine:
      SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
    1 TEMP,PRESS,SNAME)
C   
      INCLUDE 'ABA_PARAM.INC'
C   
      DIMENSION FLUX(2), TIME(2), COORDS(3)
      CHARACTER*80 SNAME
     
      d=0.00010d0 !Laserdiameter [m]
     
      x0=0.0006d0     !STARTING COORDINATE OF THE LASER BEAM
      y0=0.00055d0        !STARTING COORDINATE OF THE LASER BEAM
      z0=0.0d0          !STARTING COORDINATE OF THE LASER BEAM
     
   

   
      h=0.0001d0 !penetrationdepth [m]
      n=1d0
      W=120d0 !Laserpower [W]

      !components of vector of laser velocity Geschwindigkeit [m/s]
      vx=0.00d0
      vy=0.00d0
      vz=0.7d0
     
     
      xarc=vx*time(1)+x0
      yarc=vy*time(1)+y0
      zarc=vz*time(1)+z0
     
      xf=coords(1)-xarc
      yf=coords(2)-yarc
      zf=coords(3)-zarc
     
     
      f=(2/h)*(1+(yf/h))
      Q=(8*n*W)/(3.1415926d0*(d**2))
     
      heat=Q*f*exp(-((8*(((xf)**2)+((zf)**2)))/d**2))
 
     
      IF(SQRT(xf*xf+yf*yf+zf*zf)<d) THEN
 
            FLUX(1)=heat
            ELSE
            FLUX(1)=0.
      End IF
     
     
      RETURN
      END

Inputfile:
*Heading
** Job name: X30Mn22z Model name: Model-1
** Generated by: Abaqus/CAE 2017
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
      1,    0.000665, 0.000549999997,          0.
      2,    0.000665, 0.000549999997, 0.00400000019
      3,    0.000665, 0.000449999992, 0.00400000019
      4,    0.000665, 0.000449999992,          0.
      5, 0.00079999998, 0.000549999997, 0.00400000019
      6, 0.00079999998, 0.000449999992, 0.00400000019
      7, 0.00079999998, 0.000449999992,          0.
      8, 0.00079999998, 0.000549999997,          0.
      9, 0.00039999999, 0.000549999997, 0.00400000019
    10, 0.00039999999, 0.000449999992, 0.00400000

                    .......

*Surface, type=ELEMENT, name=Surf-7
_Surf-7_S2, S2
_Surf-7_S5, S5
_Surf-7_S6, S6
_Surf-7_S1, S1
*End Assembly
**
** MATERIALS
**
*Material, name=X30Mn22
*Conductivity
  14.,  0.
17.5, 300.
  21., 600.
  23., 900.
  27.,1200.
  33.,1350.
  42.,1500.
*Density
6906.,
*Latent Heat
333.7,1650.,1750.
*Specific Heat
493.31,  20.
  500.1, 100.
508.73, 300.
535.51, 500.
554.14, 700.
573.42, 900.
593.37,1100.
  603.6,1200.
**
** PHYSICAL CONSTANTS
**
*Physical Constants, absolute zero=-275., stefan boltzmann=5.67e-08
** ----------------------------------------------------------------
**
** STEP: Welding
**
*Step, name=Welding, nlgeom=NO, inc=100000
*Heat Transfer, end=PERIOD
4e-05, 0.00571, , ,
**
** LOADS
**
** Name: bodyheatflux  Type: Body heat flux
*Dflux
Set-8, BFNU, 1.
**
** INTERACTIONS
**
** Interaction: Konvektion
*Sfilm
Surf-5, F, 23., 0.65
** Interaction: Strahlung
*Sradiate
Surf-6, R, 23., 0.7
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1_termocouple01
**
*Output, history
*Node Output, nset=tempcheckpoint
NT,
*End Step

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

zero heat flux heißt hier, dass keine Wärme eingebracht wird, also:
- die Routine DFLUX wird nicht aufgerufen
- oder die Routine DFLUX liefert 0 zurück.
Ich vermute letzteres.

Dein Laserdurchmesser hat 0.1 mm und die Lasertiefe hat ebenfalls 0.1 mm. Für die Laser, die ich kenne liegt beides eher in der Größe 1-10 mm. Gleichzeitig liegt Dein 1. Knoten bei (0.665 mm, 0.55 mm, 0 mm) und der 2. bei (0.665 mm, 0.55 mm, 4.0 mm). Damit scheint der Knotenabstand deutlich größer als Dein Laserspot zu sein. Die Vernetzung (auch in Tiefenrichtung) unter dem Laserspot muß aber deutlich feiner als der Laserspot ausfallen, da die Wärmeeinkopplung nur an den Integrationspunkten(?) erfolgt. So aber strahlt Dein Laser ins numerische Nirwana und es wir keine Wärme eingebracht.

Anmerkung:
absolute zero liegt bei -273.15 °C

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP