| | | Dassault PLM Integrationen mit der 3DViewStation, eine Pressemitteilung
|
Autor
|
Thema: Laser Erwärmung - DFLUX (261 / mal gelesen)
|
samson05 Mitglied Wissenschaftlicher Mitarbeiter
Beiträge: 104 Registriert: 20.05.2013 Xeon E3-1245 v5 @3.50 GHz 64GB Ram 512 Gb SSD Inventor 2019, Ansys 2020R2, Abaqus 2020 Windows 7 64bit
|
erstellt am: 16. Okt. 2020 11:36 <-- editieren / zitieren --> Unities abgeben:
Hallo zusammen, ich versuche grad die Aufheizung einer Blechplatine mit einem Laser zu simulieren. Der Laser fährt auf einer Kreisbahn. Genutzt wird die User Subroutine DFLUX. Das funktioniert auch alles soweit. Das einzige was nicht funktioniert sind einigermaßen sinnvolle wärte für die später eingebrachte Temperatur und Temperaturverteilung im Blech. Wenn ich mit einem 4kW Laser auf mein Blech strahle, dann sollte da doch eigentlich nicht nur ein Temperaturdelta von 0.5 Grad heraus kommen. Mit den konsistenten Einheiten bin ich mir eigentlich sicher, dass ich hier keinen Fehler habe. Vielleicht liegt es irgendwo anders - ich bin mittlerweile ratlos. Hat jemand eine Idee? .inp und die Subroutine sind angefügt. Hier nochmal der code: Code:
subroutine DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS, & JLTYP,TEMP,PRESS,SNAME) C include 'ABA_PARAM.INC' C dimension FLUX(2), TIME(2), COORDS(3) CHARACTER*80 SNAME X=COORDS(1) Y=COORDS(2) Z=COORDS(3) TotalTime = 16 !bei der aktuellen Bogenlaenge entspricht das 5mm/s Dist = 0.0 R = 100 Bogenlsenge = 1/4*3.1415*R Geschw = Bogenlsenge/TotalTime Omega = (2*3.1415/(8*TotalTime)) Y_center = 0- R*cos(1*Omega*TIME(2)) X_center = 0.0 - R*sin(1*Omega*TIME(2)) YT=Y_center - Y XT=X_center - X Dist = SQRT((YT*YT)+(XT*XT)) QL=40000.d0 !LASER POWER [mW/mm2] R0=1d0 !BEAM RADIUS [mm] q=QL*EXP((XT**2+YT**2)/R0**2) !Gauss FLUX [mW] ? ! V=5.d0 !LASER BEAM VELOCITY [mm/s] if (Dist .le. R0) then FLUX(1) = q else FLUX(1) = 0.0 endif RETURN END
Viele Grüße Samson
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mustaine Ehrenmitglied V.I.P. h.c.
Beiträge: 3585 Registriert: 04.08.2005 Abaqus
|
erstellt am: 16. Okt. 2020 15:17 <-- editieren / zitieren --> Unities abgeben: Nur für samson05
|
rombikus Mitglied
Beiträge: 29 Registriert: 24.01.2013
|
erstellt am: 16. Okt. 2020 16:37 <-- editieren / zitieren --> Unities abgeben: Nur für samson05
Ausgehend aus dem E-Modul verwendest du ein MPa System. In diesem Fall stimmt die Dichte nicht und muss 7.8e-9 groß sein, wenn es sich um das Stahlmaterial handelt. Die spezifische Wärmekapazität ist J/(kg*K) mit J=kg*m^2/s^2. Bei 443 J/(kg*K) ergibt sich eine in diesem Einheitensystem umgerechnete anfängliche Wärmekapazität von 4.43e9 mm^2/(K*s^2). Prüfe dies noch ein Mal. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
samson05 Mitglied Wissenschaftlicher Mitarbeiter
Beiträge: 104 Registriert: 20.05.2013 Xeon E3-1245 v5 @3.50 GHz 64GB Ram 512 Gb SSD Inventor 2019, Ansys 2020R2, Abaqus 2020 Windows 7 64bit
|
erstellt am: 16. Okt. 2020 19:59 <-- editieren / zitieren --> Unities abgeben:
|
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|