!/batch,list finish /clear,start /Title,Mehrfachentladung bei der Mikrofunkenerosion /prep7 !shpp,off ! stellt ShapeWarning aus !!!Parameterliste!!!!!!!!!!!!!!!!!!!!!!!!!!! we=40e-6 ! [J] d=20e-6 ! [µm] ![m] te1=0.3e-6 ! [µs] Qq=0.12*we/(3.1416*(d/2)*(d/2)*te1) !Wärmestromdichte [W/m²] v=8.33e-6 ! [m/s] V=0.5 mm/min h=0.2 ! [m] Tu=20 ! [°C] !Umgebungstemperatur a=10 b=a*d ! [µm] n=50 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Berechnung der Frequenz = F(h) f=1000 ! [Hz] Entladefrequenz des Drahtes x=h/d ! Entladungsanzahl entlang einer Höhenlinie f1=f/x ! Entladungsfrequenz entlang einer Linie in Vorschubrichtung des Drahtes T=1/f1 ! T=Zeit für einen Entlade- und Abkühlvorgang te2=10e-6 ! 1.Abkühlzeit te3=T-te1-te2! 2.Abkühlzeit zw. 2 Entladungen te=1e-6 ! 20°C Aufwärmzeit y=2 ! Substepanzahl vom Aufwärmen auf 20° und Endabkühlvorgang z=5 ! Substepzahl vom Entlade- und Abschreckvorgang !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Berechnung der Iterationsschrittanzahl m1=(b)/(v*T) m=nint(m1) ! Anzahl der sim. Entladungen tg=m*T+te+10 ! Gesamtzeit eb=d/n ! Elementbreite u1=(b)/m ! Verschiebungsweg nach jeder Entladung k1=u1/eb ! Anzahl Elemente/Weg k=nint(k1) ! ganzzahlige Elementanzahl u=k*eb ! normierter Verschiebungsweg te4=tg ! Endabkühlzeit !!!!!!!!!!!!!!!!!!!!!!!!!!thermische Umgebung 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 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Geometrie rectng,,b,0,d rectng,b,2*b,-b/2+d/2,b/2+d/2 aglue,1,2 !!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing Parameter LESIZE,1,,,10*n LESIZE,4,,,n LESIZE,9,,,n LESIZE,10,,,n !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing MSHAPE,0,2D MSHKEY,1 AMESH,1 MSHAPE,1,2D MSHKEY,0 SMRTSIZE,2 AMESH,3 finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!therm.SOLU /SOLU ANTYPE,TRANS,NEW !Berechnung wird transient durchgeführt TINTP,,,,1 !Zeit-Integrationsparameter (Methode nach Euler Rückwärtseinsetzen) OUTRES,ERASE OUTRES,ALL,NONE OUTRES,NSOL,ALL autots,off lsel,s,,,6 lsel,a,,,5 lsel,a,,,7 nsll,s,1 d,all,temp,TU alls physics,write,thermal !physics,clear ! Clear all bc's and options finish !!!!!!!!!!!!!!!!!!!!!!!!!Strukturumgebung /prep7 lsclear,all ET,1,PLANE82 !42 !!!!!!!!!!!!!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 !!!!!plastische Verformungsoption!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! tb,BKIN,1,6,,1 ! Bilinear Kinematic Hardening für 6 Temp. TBTEMP,20 TBDATA,1,900e+6,0.24e+9 !11.0e+9 !berechnete Werte! TBTEMP,90 TBDATA,1,810e+6,0.19e+9 !10.2e+9 ! TBTEMP,200 TBDATA,1,630e+6,0.33e+9 !9.4e+9 ! TBTEMP,315 TBDATA,1,570e+6,0.245e+9 !9.0e+9 ! TBTEMP,425 TBDATA,1,540e+6,0.192e+9 !8.3e+9 ! TBTEMP,500 TBDATA,1,450e+6,0.11e+9 !7.8e+9 ! TBLIST,BKIN,1 /XRANGE,0,0.02 TBPLOT,BKIN,1 finish !!!!!!!!!!!!!!!!!!!!!!!Struktur Solu!!!!!!!! /solu !!!statische Berechnung antype,static !OUTRES,ERASE !OUTRES,ALL,NONE OUTRES,ALL,ALL nlgeom,on autots,on !off !rescontrol,,all,1,5 !!Einspannung lsel,s,,,6 lsel,a,,,7 lsel,a,,,5 nsll,s,1 d,all,uy d,all,ux alls physics,write,struct finish !!!!!!!!!!!!Schleife für die Folgeentladungen *DO,i,0,m !!!für i=0 bis Entladungsanzahl!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!thermischer Teil!!!!!!!!! physics,read,thermal !therm.Umgebung laden !/assign,esav,SR,esav !Resultsfile fortsetzen !/assign,emat,SR,emat !/assign,rth,SR,rth /solu *IF,i,lt,1,then TUNIF,20 time,te !Zeit am Ende dieses Lastfalls NSUB,1 alls solve *ELSE !antype,trans,rest !!!a)Aufwärmen NSEL,S,LOC,Y,0 NSEL,r,LOC,X,(i-1)*u,d+(i-1)*u SF,all,hflux,Qq TIME,te+te1+(te1+te2+te3)*(i-1) NSUB,z neqit,200 kbc,1 alls Solve !!!b)Abkühlen NSEL,S,LOC,Y,0 NSEL,r,LOC,X,(i-1)*u,d+(i-1)*u SFDEL,all,HFLUX alls NSEL,S,LOC,Y,0 NSEL,r,LOC,X,(i-1)*u,d+(i-1)*u SF,all,conv,1500,20 TIME,te+te1+te2+(te1+te2+te3)*(i-1) NSUB,z neqit,200 kbc,1 alls Solve !!!c)Pause !NSEL,S,LOC,Y,0 !NSEL,r,LOC,X,(i-1)*u,d+(i-1)*u !SFDEL,all,conv TIME,te+i*(te1+te2+te3) NSUB,y neqit,200 kbc,1 alls Solve *ENDIF finish !/assign,esav !/assign,emat !/assign,rth !!!!!!!!!!!!!!!!!!!!!!!!!!!Strukturteil!!!!!!!! physics,read,struct !Strukturumgebung laden /assign,esav,sr,esav !Resultsfile fortsetzen /assign,emat,sr,emat !/assign,rst,sr,rst /solu *IF,i,lt,1,then !!!!!1.Lastschritt TU=20°C k=0 *Do,l,1,1 *Do,m,1,1 ldread,temp,l,m,,,,rth k=k+1 TIME,te*k/1 neqit,200 kbc,0 NSUB,2 alls solve *EndDo *EndDo *ELSE q=(i*3)-3 antype,static,rest,i+q,y !!!1.Entladung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!für Aufwärmvorgang k=0 *Do,l,1+q,1+q *Do,m,1,z ldread,temp,l,m,,,,rth k=k+1 TIME,te+(te1+te2+te3)*(i-1)+(te1*k)/z neqit,200 kbc,0 NSUB,5 alls solve *EndDo *EndDo !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !geschmolzene Elemente entfernen nach Spannung! !/post1 !set,last !etable,molten,s,eqv !esel,s,etab,molten,0.4e+9 !finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !/solu !antype,,rest !ekill,all !finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !/post1 !set,last !esel,s,live !plnsol,s,eqv !finish !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /solu !!!!!!!!!!!für Abkühlvorgang R=0 *Do,l,2+q,2+q *Do,m,1,z ldread,temp,l,m,,,,rth R=R+1 TIME,te+(te1*i)+(te2+te3)*(i-1)+(te2*R)/z neqit,200 kbc,0 NSUB,5 alls solve *EndDo *EndDo !!!Pause!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! R=0 *Do,l,3+q,3+q *Do,m,1,y ldread,temp,l,m,,,,rth R=R+1 TIME,te+(te1+te2)*i+te3*(i-1)+(te3*R)/y neqit,200 kbc,0 NSUB,2 alls solve *EndDo *EndDo !!!!!!!!!!!!! *ENDIF finish /assign,esav /assign,emat !/assign,rst *ENDDO !Ende der Entladungsschleife!!!!!!!!!!!!!!!!! finish save !save Database