finish /clear,start /Title,Einzelentladungskrater bei der Mikrofunkenerosion /prep7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Parameter we=60e-6 ![J] d=20e-6 ![m] f=0.1 ![Hz] T=1/f ![s] te=1 !Anfangstemp.Phase te1=0.4e-6 ![s] te2=T-te1 te3=10 Qq=0.12*we/(3.1416*(d/2)*(d/2)*te1) !Wärmestromdichte [W/m²] Tu=20 !Umgebungstemperatur n=50 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Definition des Elemtetyps ET,1,PLANE77 !55 !!!!!!!!!!!!!!!!! MPTEMP,1,20,90,200,315,425,500 mpdata,kxx,1,1,7.2,7.5,8.5,9.8,11.0,12.1 !Thermische Leitfähigkeit [W/(m²K)] mpdata,c,1,1,560,570,600,650,700,740 !Thermische Kapazität [J/(kg K)] mp,dens,1,4450 !Dichte [kg/m³] !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Geometrie rectng,-2*d,6*d,0,3*d rectng,,4*d,0,d AOVLAP,1,2 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing Parameter LESIZE,5,,,4*n LESIZE,8,,,n LESIZE,4,,,n/4 LESIZE,2,,,n/4 LESIZE,3,,,n/2 LESIZE,9,,,n/6 LESIZE,10,,,n/6 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing MSHAPE,0,2D MSHKEY,1 AMESH,2 MSHAPE,1,2D MSHKEY,0 AMESH,3 finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Lastschritte... /SOLU ANTYPE,TRANS,NEW !Berechnung wird transient durchgeführt TINTP,,,,1 !Zeit-Integrationsparameter (Methode nach Euler Rückwärtseinsetzen) OUTRES,all,all autots,on lsel,s,,,9 lsel,a,,,4 lsel,a,,,3 lsel,a,,,2 lsel,a,,,10 nsll,s,1 d,all,temp,TU alls !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Lastschritt 1 (Zimmertemp) TIMINT,OFF !Ausstellen der transienten Efffekte time,te NSUB,5 !Zeit am Ende dieses Lastfalls D,ALL,TEMP,TU !Zimmertemp alls solve !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!2.Lastschritt-Aufwärmen TIMINT,on ddel,all,temp *Do,I,1,4 !!!a)Aufwärmen NSEL,S,LOC,Y,0,d/10 NSEL,r,LOC,X,(i-1)*d,i*d SF,all,hflux,Qq TIME,te+te1+(i-1)*(te2+te1) NSUB,20 neqit,200 kbc,1 alls Solve !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!b)Abkühlen NSEL,S,LOC,Y,0 NSEL,r,LOC,X,(i-1)*d,i*d SFDEL,all,HFLUX alls !Kratermittellinie über X-Y-Koordinaten auswählen NSEL,S,LOC,Y,0,d/10 NSEL,r,LOC,X,(i-1)*d,i*d SF,all,conv,1500,20 TIME,te+i*(te1+te2) NSUB,20 neqit,200 kbc,1 alls Solve *ENDDO !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !letzter Lastschritt vollst. Abkühlen NSEL,S,LOC,Y,0,d/10 NSEL,r,LOC,X,(i-1)*d,i*d SFDEL,all,conv alls TIME,te3 NSUB,20 neqit,200 kbc,1 alls solve !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Strukturanalyse /prep7 ET,2,PLANE82 !42 etchg,tts !Elementtyp von tehrm. auf Strukturanalyse umschaltenplane55 etlist !!!!!!!!!!!!!Materialdaten für Strukturanalyse MPTEMP,1,20,90,200,315,425,500 mpdata,alpx,1,1,8.8e-6,9.0e-6,9.3e-6,9.6e-6,9.8e-6,10e-6 mpdata,ex,1,1,113e+9,107.35e+9,101.7e+9,92.66e+9,84.75e+9,70.06e+9 mp,prxy,1,0.361 finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /solu antype,trans outpr,all,all nlgeom,on autots,on !!Einspannung lsel,s,,,3 nsll,s,1 d,all,uy !d,all,ux alls !!!Daten der thermischen Berechnung auf das Modell transferieren! !!!für jeden einzelnen SubStep der th.Sim. eine Strukturlösung berechnen! !!!!!1.Lastschritt TU=20°C K=0 *Do,i,1,1 *Do,j,1,4 ldread,temp,1,1,,,,rth k=k+1 TIME,(te*k)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!1.Entladung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!für Aufwärmvorgang k=0 *Do,i,2,2 *Do,j,1,4 ldread,temp,i,j,,,,rth k=k+1 TIME,te+(te1*k)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!!!!!!!für Abkühlvorgang R=0 *Do,l,3,3 *Do,m,1,4 ldread,temp,l,m,,,,rth R=R+1 TIME,te+te1+(te2*R)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!2.Entladung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!für Aufwärmvorgang k=0 *Do,i,4,4 *Do,j,1,4 ldread,temp,i,j,,,,rth k=k+1 TIME,te+te1+te2+(te1*k)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!!!!!!!für Abkühlvorgang R=0 *Do,l,5,5 *Do,m,1,4 ldread,temp,l,m,,,,rth R=R+1 TIME,te+2*te1+te2+(te2*R)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!3.Entladung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!für Aufwärmvorgang k=0 *Do,i,6,6 *Do,j,1,4 ldread,temp,i,j,,,,rth k=k+1 TIME,te+2*(te1+te2)+(te1*k)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!!!!!!!für Abkühlvorgang R=0 *Do,l,7,7 *Do,m,1,4 ldread,temp,l,m,,,,rth R=R+1 TIME,te+te1+2*(te1+te2)+(te2*R)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!4.Entladung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!für Aufwärmvorgang k=0 *Do,i,8,8 *Do,j,1,4 ldread,temp,i,j,,,,rth k=k+1 TIME,te+3*(te1+te2)+(te1*k)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!!!!!!!für Abkühlvorgang R=0 *Do,l,9,9 *Do,m,1,4 ldread,temp,l,m,,,,rth R=R+1 TIME,te+te1+3*(te1+te2)+(te2*R)/4 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!10.Lastschritt-Endabkühlen k=0 *Do,i,10,10 *Do,j,1,5 ldread,temp,i,j,,,,rth k=k+1 TIME,4*(te1+te2)+te+(te3-4*(te1+te2)-te)*k/5 neqit,200 kbc,0 NSUB,20 solve *EndDo *EndDo !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /post26 nsol,2,103,u,y ! Verschiebung in Y-Richtung !nsol,2,103,temp lines,1000 ! Anzahl der max.Zeilen prvar,2 ! Auflisten der Werte