finish
/clear,all

vtop = 200

/prep7



! - Bimorph
block,0,9e-3,0,2e-3,10e-3,115e-3
block,0,9e-3,-2e-3,0,0,125e-3
block,0,9e-3,-4e-3,-2e-3,10e-3,115e-3

! - Skin
block,-9e-3,18e-3,2e-3,3e-3,0,125e-3


vglue,1,2,3


!- Bimorph Material


/com
/com	Pz21	Y-polarized
/com


/com	Stiffness
TB,	ANEL,	1	,	1	,0
TBDATA,	1,	1.1419E+11	,	7.2357E+10	,	7.5740E+10
TBDATA,	7,	1.1077E+11	,	7.2357E+10
TBDATA,	12,	1.1419E+11
TBDATA,	16,	2.6291E+10
TBDATA,	19,	2.6291E+10
TBDATA,	21,	1.9227E+10


/com	Piezo matrix
TB,	PIEZ,1
TBDATA,	2,	-2.9179
TBDATA,	5,	33.3852
TBDATA,	8,	-2.9179
TBDATA,	10,	16.1936
TBDATA,	15,	16.1936


/com	Permittivity
EMUNIT,	EPZRO,	8.85E-12
MP,	PERX,	1	,	2118
MP,	PERY,	1	,	1401
MP,	PERZ,	1	,	2118


/com	Density
MP,	DENS,	1	,	7780

/com
/com  -- Material properties for chosen middle Material (Glass Epoxy Composite)
/com


mp,ex,2,25e9				!Young modulus
mp,prxy,2,0.136				!Poisson coefficient
mp,dens,2,1900				!Density

/com
/com  -- Material properties for Skin (Elastomer Silicone)
/com

mp,ex,3,2.5e6				 !Young modulus
mp,prxy,3,0.48				!Poisson coefficient
mp,dens,3,1250				!Density

et,1,solid227,1001  	! Hexaedral Coupled-Field Solid

et,2,solid186   	    ! 3-D 20-Node Structural

et,3,conta174
    keyopt,3,1,5				! DOF: x,y,z,volt
    keyopt,3,4,1
    keyopt,3,5,1				! close gaps
    keyopt,3,9,3				! only offset
    keyopt,3,11,1			! consider shellthickness
    keyopt,3,12,5			! always connected

et,4,targe170
		keyopt,5,5,1			! Art der RB: Solid-Solid (Rot FHG sind frei)

vsel,s,loc,y,2e-3,3e-3
type,2
mat,2
real,2
mshape,1,3d				!Volume mesh using tetrahedral shape
mshkey,0				!Mapped meshing off
vmesh,all
allsel,all

vsel,s,loc,y,0,2e-3
type,1
mat,1
real,1
mshape,1,3d				!Volume mesh using tetrahedral shape
mshkey,0				!Mapped meshing off
vmesh,all
allsel,all


vsel,s,loc,y,-2e-3,0
type,2
mat,2
real,2
mshape,1,3d				!Volume mesh using tetrahedral shape
mshkey,0				!Mapped meshing off
vmesh,all
allsel,all

vsel,s,loc,y,-4e-3,-2e-3
type,1
mat,1
mshape,1,3d				!Volume mesh using tetrahedral shape
mshkey,0				!Mapped meshing off
vmesh,all
allsel,all

vsel,s,loc,y,2e-3,3e-3
asel,r,loc,y,2e-3
type,3
real,3
esurf
allsel,all

vsel,s,loc,y,0,2e-3
asel,r,loc,y,2e-3
type,4
real,4
esurf
allsel,all

!vsel,s,loc,y,0,2e-3
!asel,r,loc,y,2e-3
!cm,piezoCont,area
!allsel
!
!vsel,s,loc,y,2e-3,3e-3
!asel,r,loc,y,2e-3
!cm,TargCont,area
!allsel

!edbx,add,1,-9e-3,18e-3,2e-3,3e-3,0,125e-3
!edbx,add,2,0,9e-3,-4e-3,2e-3,10e-3,115e-3
!edcgen,ag,,,,,,,,,,,,,,2,1

nsel,s,loc,y,2e-3
nsel,r,loc,x,0,9e-3
d,all,volt,vtop
allsel,all

nsel,s,loc,y,-4e-3
d,all,volt,vtop
allsel,all

nsel,s,loc,y,0
d,all,volt,0
allsel,all

nsel,s,loc,y,-2e-3
d,all,volt,0
allsel,all

nsel,s,loc,z,0
d,all,ux,0
d,all,uy,0
d,all,uz,0

nsel,s,loc,z,125e-3
d,all,ux,0
d,all,uy,0
d,all,uz,0

allsel,all
/EOF
/Solu

solve

/post1

pldisp,2