fini /clear /title, Composite Platte /FILNAME, Platte aus Composite /com Alles in SI Einheiten /com Solid5 und Shell63 /com 4 Node Shell and 8 node Solid /com use same nodes /PREP7 ! Material 1 - Piezo Solid5 ! == Material parameters for piezo ! -- Elastic coefficients ! Alle Werte in N/m^2=Pa c11 = 120e9 c12 = 75.1e9 c13 = 75.1e9 c16 = 0 c14 = 0 c15 = 0 c22 = c11 c23 = c13 c26 = 0 c24 = 0 c25 = 0 c33 = 111e9 c36 = 0 c34 = 0 c35 = 0 c44 = 21e9 c45 = 0 c55 = c44 c66 = c44 c64 = 0 c65 = 0 ! Nur Eingabe des oberen Dreiecks, da auch für anisotrope Materialien gillt: ! cij=cji. deswegen nur unteres oder oberes Dreieck notwendig tb,ANEL,0 ! ,1 bedeutet Eingabe als Nachgiebigkeitsmatrix, 0 ist Steifigkeitsmatrix tbdata, 1, c11, c12, c13, c16, c14, c15 tbdata, 7, c22, c23, c26, c24, c25 tbdata,12, c33, c36, c34, c35 tbdata,16, c66, c64, c65 tbdata,19, c44, c45 tbdata,21, c55 ! -- Piezoelectric coefficients, Coulombs/m^2 e11 = 0 e21 = 0 e31 = -5.31 e12 = 0 e22 = 0 e32 = e31 e13 = 0 e23 = 0 e33 = 0 e33 = 15.78 e16 = 0 e26 = 0 e36 = 0 e14 = 0 e24 = 12.3 e34 = 0 e15 = e24 e25 = 0 e35 = 0 ! 0 -- ! Piezoelectric stress matrix [e] (used as supplied) ! 1 -- ! Piezoelectric strain matrix [d] (converted to [e] form before use) tb,PIEZ,0 tbdata, 1, e11, e21, e31 tbdata, 4, e12, e22, e32 tbdata, 7, e13, e23, e33 tbdata,10, e16, e26, e36 tbdata,13, e14, e24, e34 tbdata,16, e15, e25, e35 ! solid 226: dielectricity matrix in F/m ! hier aber Solid5: ! für solid 5 Eingabe als Absolutwerte nötig ! Epsilon_1_S / Epsilon_Null mp,PERx,1,916 mp,PERy,1,916 mp,PERz,1,830 ! -- Density, kg/m3 mp,DENS,1,7500 ! Daten sind für Gewebe ! Materialparameter für orthotropes Material Nr. 1 eingeben ! Mechanik MP,EX,2,26E9 ! ORTHOTROPIC MATERIAL PROPERTIES MP,EY,2,26E9 MP,EZ,2,10E9 MP,GXY,2,3.53E9 MP,GYZ,2,3.7E9 MP,GXZ,2,3.7E9 MP,PRXY,2,0.112 ! MAJOR POISSONS RATIO MP,PRYZ,2,0.311 ! MAJOR POISSONS RATIO MP,PRXZ,2,0.311 ! MAJOR POISSONS RATIO ! -- Density, kg/m3 mp,DENS,2,2000 ET,1,SOLID5,3 ! 8 Node Multi-Physics-Solid; UX, UY, UZ, VOLT if KEYOPT(1) = 3 ET,2,SHELL63 ! 4 NODE LAYERED SHELL; 6 DOF R,2,0.002 et,3,94,0 ! Circuit elemente für piezoelectricity keyopt(1)=0, also Widerstand r,3,1000 ! Real constant 1000 Ohm für Widerstand KEYOPT,3,6,0 ! keyoption(6) vom elementtyp 3 ist=0, also piezoelektrisch ! Definition of Analysis parameters frequency_start = 1 ! Untere Analysefrequenz frequency_stop = 401! ! Obere Analysefrequenz frequencysteps = 3 ! (frequency_stop - frequency_start) ! Anzahl der Frequenzanalysen Widerstand = 1000 ! Wert für elektrischen Widerstand definieren ! Dimensions of Aluminum plate/ GFK-plate length_alu = 200e-3 width_alu = 60e-3 height_alu = 2e-3 ! == Dimensions of piezo length_piezo = 0.02 width_piezo = 0.02 height_piezo = 0.0002 offset_piezo_x = width_alu/2-width_piezo/2 offset_piezo_y = 0.1 offset_piezo_z = height_alu ! Generate Piezo Volume block,offset_piezo_x,offset_piezo_x+width_piezo,offset_piezo_y,offset_piezo_y+length_piezo,offset_piezo_z,offset_piezo_z+height_piezo ! Generate area for plate shell wplane,1,0,0,offset_piezo_z,1,0,offset_piezo_z,0,1,offset_piezo_z Rectng,0,width_alu,0,length_alu wplane,1,0,0,0,1,0,0,0,1,0 ! Material affiliation, mesh mshkey, 1 ! use mapped meshing mshape, 0,2D ! hexaeder-elemente in 3d nehmen esize,0.002 real,2 type,2 mat,2 allsel amesh,7 ! aluplatte allsel mshkey, 1 ! use mapped meshing mshape, 0,3D ! hexaeder-elemente in 3d nehmen ESIZE,0.002 ! Wenig Elemente erzeugen real,1 type,1 mat,1 vsel,all vsel,s,loc,z,offset_piezo_z,offset_piezo_z+height_piezo vmesh,all ! piezo solid allsel,all ! Gemeinsame Knoten von Shell und Solid mergen nummrg,node,e-4,e-4 !================================================================ ! Erzeugen von Knoten für Widerstand n,90000,-0.1,0,0 n,90001,-0.1,0.1,0 ! type, 3 ! Pointer auf 3 = Circuit elements setzen ! mat,0 ! Material pointer setzen ! real, 3 ! Real constant set pointer ! e,90000,90001 ! Erzeuge Element vom Typ 3 aus den beiden Knoten allsel,all !================================================================ ! Knoten auf Oberfläche für obere Elektrode selektieren nsel,s,loc,z,(offset_piezo_z+height_piezo),(offset_piezo_z+height_piezo) ! Alle Knoten auf der Oberfläche der Keramik selektieren nsel,a,node,,90000 cp,1,volt,all *get,n_supply,node,0,num,min ! Niedrigste Knotennummer des definierten knoten-sets (cp) suchen allsel ! Knoten Oberfläche für untere Elektrode selektieren nsel,s,loc,z,offset_piezo_z,offset_piezo_z nsel,r,loc,x,offset_piezo_x,offset_piezo_x+width_piezo nsel,r,loc,y,offset_piezo_y,offset_piezo_y+length_piezo nsel,a,node,,90001 cp,2,volt,all ! Alle Knoten zu einem set nr. 2 zusammenfassen *get,n_ground,node,0,num,min allsel eplot ! Mechnical BCs: einseitig eingespannt allsel nsel,s,loc,y,0,0 ! unterer Rand fest eingespannt nsel,r,loc,x,0,width_alu d,all,all,0 allsel ematwrite,yes ! forces the writing of all element matrices to file.emat finish /solu eqslv,sparse ! Sparse solver verwenden, bei Circuit analysis funktioniert kein anderer solver nsel,s,loc,z,height_alu,height_alu ! Define Node for acting force nsel,r,loc,x,0.9*width_alu/2,1.1*width_alu/2 nsel,r,loc,y,1.0*length_alu,1.0*length_alu *get,fnode,node,,num,min allsel f,fnode,fz,-1 ! An fnode eine sinuskraft von 1 Newton ! in negative z-Richtung aufbringen antype,harmic ! harmonische Analyse starten hropt, full ! Specifies harmonic analysis options, full=alles harfrq,frequency_start,frequency_stop ! Defines the frequency range in the harmonic response analysis nsubst,frequencysteps ! Specifies the number of substeps to be taken this load step outres,all,all ! Controls the solution data written to the database outpr,all,all ! Controls the solution printout kbc,1 ! Specifies stepped or ramped loading within a load step ! 1=step changed mp,dmprat,2,0.001 ! damping ratio allsel solve finish