finish /clear /prep7 ! !************************************** !********** Parametereingabe ********** !************************************** ! /prep7 laenge_naht=0.5 ! [m] nahtvorb=0.005 ![5mm] breite=0.2 ! [m] nahtspalt=0.002 dicke=0.011 nahtueberhoeh=0.0028 ! [m] nahtrad=0.011 ! [m] Krìmmungsradius der Naht !leistung_unten=3.07571172*10**10 ! 10000 Watt !leistung_oben=1.751072751*10**10 leistung_unten=0.125*3.07571172*10**10 ! 5000 watt leistung_oben=0.125*1.751072751*10**10 ! !************************************** !*********** Knoteneingabe ************ !************************************** ! /prep7 k,1,0,0,0 ! Ursprung k,10000,0,0,5 ! Hilfs-Achse zur Drehung ... k,10001,0,0,-5 ! ...im Zyl-Kos --> Geometriegenerierung k,2,0.5*nahtspalt,0,0 k,3,0.5*nahtspalt+nahtvorb,dicke,0 k,4,0.5*nahtspalt+nahtvorb,0,0 k,5,0.5*nahtspalt+nahtvorb+breite,0,0 k,6,0.5*nahtspalt+nahtvorb+breite,dicke,0 k,7,0,nahtueberhoeh+dicke,0 ! !************************************** !*********** Linieneingabe ************ !************************************** ! /prep7 larc,7,3,1,nahtrad l,1,2 l,1,7 l,2,3 l,2,4 l,3,4 l,3,6 l,5,6 l,4,5 ! !************************************** !*********** Unterteilung ************* !************************************** ! /prep7 lesize,1,,,1,,,,, lesize,2,,,1,,,,, lesize,3,,,2,,,,, lesize,4,,,2,,,,, lesize,5,,,2,,,,, lesize,6,,,2,,,,, lesize,7,,,20,4,,,, lesize,8,,,2,,,,, lesize,9,,,20,4,,,, ! !************************************** !************* Flaechen *************** !************************************** ! /prep7 al,1,2,3,4 al,4,5,6 al,6,7,8,9 ! !************************************** !******** Elementtypdefinition ******** !******* und Materialdefinition ******* !************************************** ! /prep7 et,1,solid278 mp,kxx,1,40 ! [W/m*K] mp,dens,1,7800 ! [kg/m3] mp,c,1,500 ! [J/kg*K] ! ! voffst,1,laenge_naht,0 voffst,2,laenge_naht,0 voffst,3,-laenge_naht,0 ! ! nummrg,all numcmp,all ! delete all unused numbers for lines ! !************************************** !************ Netzgroesse ************* !************************************** ! /prep7 lesize,14,,,40,,,,, lesize,15,,,40,,,,, lesize,16,,,40,,,,, lesize,17,,,40,,,,, lesize,20,,,40,,,,, lesize,24,,,40,,,,, lesize,25,,,40,,,,, ! ! vsymm,x,1,3,1,0,0 vglue,1,4 numcmp,all lesize,43,,,1,,,,, lesize,44,,,1,,,,, lesize,23,,,20,0.25,,,, lesize,40,,,20,0.25,,,, vglue,1,6 vmesh,all ! ! nummrg,all ! ! /view,,1,1,1 /ang,,30 /dist,all,auto, eplot ! !************************************** !************* Schleife *************** !******* untere Knotenschichten ******* !************************************** ! /prep7 allsel *do,i,1,41,1 ! vsel,s,volu,,1,6,5 nslv,s,1 nsel,r,loc,z,0.5/80-(0.5/40)*(i-1),-0.5/80-(0.5/40)*(i-1) ! Selektion des i-ten layer der Schweissnahtknoten (NDIV) nsel,u,loc,y,0.00691,0.015 cm,schweissnahtknoten_unten_%i-1%,node ! Speicherung der selektierten Knotenlayer als gruppe Schweissnahtknoten_%i-1% *enddo ! !************************************** !************* Schleife *************** !******* obere Knotenschichten ******** !************************************** ! /prep7 allsel *do,i,1,41,1 ! vsel,s,volu,,1,6,5 nslv,s,1 nsel,r,loc,z,0.5/80-(0.5/40)*(i-1),-0.5/80-(0.5/40)*(i-1) ! Selektion des i-ten layer der Schweissnahtknoten (NDIV) nsel,u,loc,y,-0.001,0.001 cm,schweissnahtknoten_oben_%i-1%,node ! Speicherung der selektierten Knotenlayer als gruppe Schweissnahtknoten_%i-1% *enddo ! !************************************** !************* Schleife *************** !******* untere Elementschicht ******** !************************************** ! /prep7 allsel *do,i,1,40,1 cmsel,s,schweissnahtknoten_unten_%i-1% cmsel,a,schweissnahtknoten_unten_%i% esln,s,1,active ! Abw½hlen der vorherigen Reihe an Elementen (inelegante l‚sung) cm,schweissnahtelemente_unten_%i-1%,elem ! Speichern der Ausgew½hlten Elementschicht als gruppierung *enddo ! ! *do,i,1,40,1 *if,i,eq,1,then, cmsel,s,schweissnahtelemente_unten_%0% *else cmsel,a,schweissnahtelemente_unten_%i-1% *endif *enddo cm,schweissnahtelemente_unten,elem ! !************************************** !************* Schleife *************** !******** obere Elementschicht ******** !************************************** ! /prep7 allsel *do,i,1,40,1 cmsel,s,schweissnahtknoten_oben_%i-1% cmsel,a,schweissnahtknoten_oben_%i% esln,s,1,active ! Abw½hlen der vorherigen Reihe an Elementen (inelegante l‚sung) cm,schweissnahtelemente_oben_%i-1%,elem ! Speichern der Ausgew½hlten Elementschicht als gruppierung *enddo ! ! *do,i,1,40,1 *if,i,eq,1,then, cmsel,s,schweissnahtelemente_oben_%0% *else cmsel,a,schweissnahtelemente_oben_%i-1% *endif *enddo cm,schweissnahtelemente_oben,elem ! ! ! cmsel,s,schweissnahtelemente_unten cmsel,a,schweissnahtelemente_oben cm,schweissnaht,elem ! ! ! allsel nsel,s,loc,x,-0.207, cm,lagerung_links,node allsel nsel,s,loc,x,0.207, cm,lagerung_rechts,node ! !************************************** !********** Schreiben der ************* !******* thermischen Umgebung ********* !************************************** ! allsel physics,write,thermal physics,clear ! !************************************** !********** Schreiben der ************* !****** strukturellen Umgebung ******** !************************************** ! allsel etchg,tts mp,ex,1,210e9 ! [N/m3] mp,prxy,1,0.3 ! [/] mp,dens,1,7800 ! [kg/m3] mp,alpx,1,12e-6 ! [m/m*K] allsel physics,write,struct physics,clear finish ! !************************************** !************* Schleife *************** !******** Wärmeeintrag unten ********** !********** und Berechnung ************ !************************************** ! /solu allsel antype,trans ! Transiente Analyse physics,read,thermal ! Lesen der thermischen Physik nropt,full outres,all,all ! Alle Zischenresultate fØr alle Lastschritte speichern time,1e-6 ! Anfangszeit der Anfangsbedingung nsubst,2 ! Substeps timint,on ! Transiente Effekte aus tunif,20 ! Alle Knoten 20 grad Celsius Anfangsbedingung allsel ! W½hle alles aus (alive und gekillte Elemente, gekillte elemente sind aber unwirksam) solve ! L‚sung finish ! Schliessen des L‚sers /solu antype,,restart timint,on ! Transiente Effekte durch Zeitintegration ein cmsel,s,schweissnahtelemente_oben cmsel,a,schweissnahtelemente_unten bfe,all,temp,,20 !nsle,s,active, !ice,all,temp,20 ekill,all esel,s,live allsel time,1e-4 solve ! ! ! *do,i,1,40,1 ! Schleife ìber alle Elementlayer /solu ! Im Loesungsmodus outres,all,all esel,s,live ! Alle live Elemente auswaehlen time,8*i ! Zeitinkrement mit Schleife gekoppelt nsubst,1,5,1 ! Substep der Zeit !cmsel,s,schweissnahtelemente_%i-1% ! Selektiere Gespeicherte Elementlayer_%i-1% !ealive,all ! Aktiviere die neue Elementschicht !eplot *if,i,eq,1,then ! Falls erte Elementschicht cmsel,s,schweissnahtelemente_unten_0 ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_unten_0 bfedele,schweissnahtelemente_unten_0,temp nsle,s,active, ddele,all,temp ! Mittlere Knotenschicht wird von oberen und unteren Elementen geteilt ! das heisst dass nur in der obersten Knotenreihe T=20 bleibt, di mittleren Knoten sind ! in den unteren (zu dem Zeitpunkt aliven) Elementen aktiv und erfahren die volle Leistung ! folglich erfahren die oberen gekillten Elemente ein Temperaturanstieg an der Unterseite ! abhilfe: vielleicht RB als element-volumenlast oder estif matrix-vorfaktor stärker bfe,schweissnahtelemente_unten_0,hgen,,leistung_unten ! Fìr alle anderen Knotenschichten 350 watt *elseif,i,eq,40 ! Falls letzte Elementschicht cmsel,s,schweissnahtelemente_unten_39 ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_unten_39 nsle,s,active ddele,all,temp bfe,schweissnahtelemente_unten_39,hgen,,leistung_unten ! Fìr alle anderen Knotenschichten 350 watt *else cmsel,s,schweissnahtelemente_unten_%i-1% ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_unten_%i-1% nsle,s,active ddele,all,temp bfe,schweissnahtelemente_unten_%i-1%,hgen,,leistung_unten ! Fìr alle anderen Knotenschichten 350 watt *endif !esel,s,live allsel eplot solve /solu allsel bfedele,all,all !esel,s,live !eplot *enddo ! !************************************** !************* Schleife *************** !********* Wärmeeintrag oben ********** !********** und Berechnung ************ !************************************** ! /solu *do,i,1,40,1 ! Schleife ìber alle Elementlayer /solu ! Im Loesungsmodus outres,all,all esel,s,live ! Alles ausw½hlen time,320+8*i ! Zeitinkrement mit Schleife gekoppelt nsubst,1,5,1 ! Substep der Zeit !cmsel,s,schweissnahtelemente_%i-1% ! Selektiere Gespeicherte Elementlayer_%i-1% !ealive,all ! Aktiviere die neue Elementschicht !eplot *if,i,eq,1,then ! Falls erte Elementschicht cmsel,s,schweissnahtelemente_oben_0 ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_oben_0 nsle,s,active ddele,all,temp bfe,schweissnahtelemente_oben_0,hgen,,leistung_oben ! Fìr alle anderen Knotenschichten 350 watt *elseif,i,eq,40 ! Falls letzte Elementschicht cmsel,s,schweissnahtelemente_oben_39 ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_oben_39 nsle,s,active ddele,all,temp bfe,schweissnahtelemente_oben_39,hgen,,leistung_oben ! Fìr alle anderen Knotenschichten 350 watt *else cmsel,s,schweissnahtelemente_oben_%i-1% ! Selektiere Erstes Gespeichertes Elementlayer ealive,schweissnahtelemente_oben_%i-1% nsle,s,active ddele,all,temp bfe,schweissnahtelemente_oben_%i-1%,hgen,,leistung_oben ! Fìr alle anderen Knotenschichten 350 watt *endif !esel,s,live allsel eplot solve /solu allsel bfedele,all,all !esel,s,live !eplot *enddo bfedele,all,all time,660 nsubst,20,500,20 solve