save,pre_mod,db
/prep7
/com,       Units SI! (m, kg, s, N ...)
/COM,
/COM,       MATERTIAL PROPERTIES FOR THE FINITE ELEMENT SOLUTION:
/COM,       CONSTITUTIVE MATRICES FOR PZT-4 (POLAR AXIS ALONG Y)
/COM,
/COM,   [c11 c13 c12  0  0  0 ]        [ 0 e31  0 ]      [ep11  0    0  ]
/COM,   [c13 c33 c13  0  0  0 ]        [ 0 e33  0 ]      [  0  ep33  0  ]
/COM,   [c12 c13 c11  0  0  0 ]        [ 0 e31  0 ]      [  0   0   ep11]
/COM,   [ 0   0   0  c44 0  0 ]        [e15 0   0 ]
/COM,   [ 0   0   0   0 c44 0 ]        [ 0  0  e15]
/COM,   [ 0   0   0   0  0 c66]        [ 0  0   0 ]
/COM,

MP,PERX,10,728.5                 ! PERMITTIVITY AT CONSTANT STRAIN
MP,PERY,10,634.7
MP,PERZ,10,728.5

TB,ANEL,10                       ! ANISOTROPIC ELASTIC STIFFNESS
TBDA,1,13.9E10,7.43E10,7.78E10  ! c11,c13,c12
TBDA,7,11.5E10,7.43E10          ! c33,c13
TBDA,12,13.9E10                 ! c11
TBDA,16,2.56E10                 ! c44
TBDA,19,2.56E10                 ! c44
TBDA,21,3.06E10                 ! c66

TB,PIEZ,10                       ! PIEZOELECTRIC STRESS COEFFICIENTS
TBDA,2,-5.2                     ! e31
TBDA,5,15.1                     ! e33
TBDA,8,-5.2                     ! e31
TBDA,10,12.7                    ! e15
TBDA,15,12.7                    ! e15

et,10,226,1001
cmsel,s,piezo
emodif,all,type,10
emodif,all,mat,10

/com, Elektroden generieren 

allsel,all 

cmsel,s,electrode1 		! Elektrode 1 auswählen 
cp,1,volt,all 		! Potential aller Knoten koppeln 
*get,n_phase,node,,num,min 	! Knotenbezeichnung: Phase 
allsel,all 

allsel,all 
cmsel,s,electrode2 		! Elektrode 2 auswählen 
cp,2,volt,all 		! Potential aller Knoten koppeln 
*get,n_gnd,node,,num,min 	! Knotenbezeichnung: Masse/ Ground 

allsel,all 
d,n_gnd,volt,0 		! Masse auf 0 Volt legen 
!d,n_phase,volt,0 
allsel,all 


R=100000000	
R,1,R 			
ET,11,CIRCU94,0 	
TYPE,11 
E,n_phase,n_gnd 	
alls

allsel,all 
outres,all,all 
finish
/solu