DEAR ALL
could any one tel me whats wrong with this code
its an fgm beam with piezopatch
FINISH
/CLEAR
/title, Controlled smart structure - K.Bendine
/nopr ! Supresses extended interpreted data input
/com, ! Places comments in the output
!/FILNAME,fgmbeam_control,0 ! This sets the jobname to 'beam_control'
!/CONFIG,NRES,3000
! GEOMETRY
*set,Lb,700e-3 ! Support Length
*set,Hb,3e-3 ! Support Thickness
*set,Wb,20e-3
! [m]
*set,bpe,0.0001
*set, ele_numb_x, 12 ! 30, 70
*set, ele_numb_y, 12 ! 20, 40
*set, ele_numb_n, 100 ! 30, 70
*set, ele_numb_z, 5 ! 20, 40
*set, length_e_x, Lb/ele_numb_x
*set, length_e_y, Wb/ele_numb_y
*set, length_e_z, Hb/ele_numb_z
!
! MATERAL DUMMY
*set, nu_dummy, 0.3
*set, E_dummy, 71.700E9 ! [Pa]
*set, alpha_dummy, 23E-6 ! [1/K]
*set, kxx_dummy, 121. ! [W/mK]
*set, dens_dummy, 1210
!
! MATERIAL TOP
! MATERIAL TOP
*set, E_t, 15.1e10 ! [Pa]
*set, nu_t, 0.3
*set,dens_t,3000
!
! MATERIAL BOT
*set, E_b, 7.0e10 ! [Pa]
*set, nu_b, 0.3
*set,dens_b,2707
!
! LOADING
*set, p, 1.
!
! FGM gradation law
*set, a_z, 1./Hb ! [m^-1]
*set, n_z, 1
!
! MESH
! 20, 40
!transient analysis
*set,Tstep,0.01
*set,Tstop,2
*set,Tres,0.01
*set,ALPHAD,0
*set,BETAD,0.00015
*set,ks,5
*set,kc,10
*set,kv,5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! mat_fgm_K MACRO !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*create, mat_fgm_K
/nopr
z_coord = arg2
grad_common = (z_coord+0.5)**n_z
E_fgm = (E_t-E_b)*grad_common + E_b
nu_fgm = (nu_t-nu_b)*grad_common + nu_b
den_fgm = (dens_t-dens_b)*grad_common + dens_b
/gopr
*end !mat_fgm_K
!
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! PREP 7 !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/PREP7
SHPP,SILENT,ON !Shape checking in silent mode
SHPP,MODIFY,1,40
et, 1, solid5
et,2,solid45
! ====================================================
! Define problem parameters - PZT
! ====================================================
! - PZT Geometry (All PZTs QP16n)
L=60e-3 ! Length, m
H=1e-3 ! Thickness, m
W=20e-3 ! Width,m
! - Material properties for PVDF
s11=16.4e-12
s22=s11
s12=-5.74e-12
s13=-7.22e-12
s23=s13
s33=18.8e-12
s44=47.5e-12
s55=47.5e-12
s66=44.3e-12
pidns=7700 ! Density of PZT
d31=-179e-12 ! Piezoelectric strain coefficients, C/N
d32=-179e-12 ! Manufacturers data combined with PZT5A mat props
d33=350e-12
d24=584e-12
d16=584e-12
ept11=1730 ! Relative permittivity at constant stress
ept22=1730
ept33=1700
! ====================================================
! Materials and element types
! ====================================================
! Defining elements and materials:
! PZT-------------------------------------------------
et,1,SOLID5,3 ! 3-D piezoelectric element
MP,DENS,1,pidns
MP,PERX,1,ept11 ! Permittivity at constant stress
MP,PERY,1,ept22
MP,PERZ,1,ept33
tb,ANEL,1,,,1 ! Assembling elastic compliance matrix
tbda,1,s11,s12,s13
tbda,7,s22,s23
tbda,12,s33
tbda,16,s66
tbda,19,s55
tbda,21,s44
tb,PIEZ,1,,,1 ! Piezoelectric strain matrix
tbda,3,d31
tbda,6,d33
tbda,9,d32
tbda,14,d24
tbda,16,d16
! ====================================================
! Define problem parameters - BEAM
! ====================================================
! - Support Beam Thickness and material props
Lb=700e-3 ! Support Length
Hb=3e-3 ! Support Thickness
Wb=20e-3
! Support Width
! ====================================================
! Define problem parameters - PZT Placement
! ====================================================
! - Top PZT 1 (Starting from the clamped end)
L1=14e-3 ! Starting point of PZT from the left edge
W1=10e-3 ! ((in reality 9)Starting point of piezo (Measured from the lower side, when clamped at the left end)
! - Top PZT 2 (Starting from the clamped end)
L2=90e-3 ! Starting point of PZT from the left edge
W2=10e-3 ! Starting point of piezo (Measured from the lower side, when clamped at the left end)
! - Top PZT 3 (Starting from the clamped end)
L3=322e-3 ! Starting point of PZT from the left edge
W3=10e-3 ! Starting point of piezo (Measured from the lower side, when clamped at the left end)
! - Bottom PZT 1 (Starting from the clamped end)
L4=14e-3 ! Starting point of PZT from the left edge
W4=10e-3 !(in reality 9)! Starting point of piezo (Measured from the lower side, when clamped at the left end)
! ====================================================
! Solid model
! ====================================================
! Beam
*set, v_numb_z, Hb/8
*do, j, -Hb/2,+Hb/2-v_numb_z,v_numb_z
BLOCK,0,Lb,0,Wb,j, j+v_numb_z
vsel,all
*enddo
BLOCK,L1,L1+L,0,W,Hb/2,Hb/2+H ! Create the active QP16 piezo wafer P1
BLOCK,L2,L2+L,0,W,Hb/2,Hb/2+H ! Create the active QP16 piezo wafer P2
BLOCK,L3,L3+L,0,W,Hb/2,Hb/2+H ! Create the active QP16 piezo wafer P3
VSEL,All
!vglue,all
! Divinding the base aluminum beam for better mesh distribution
!WPROTA,,-90,
!WPOFFS,,,W1
!VSBW,ALL
!WPOFFS,,,W
!VSBW,ALL
!WPAVE,0,0,0,
!WPROTA,,,90
!WPOFFS,,,L1
!VSBW,ALL
!WPOFFS,,,L
!VSBW,ALL
!WPOFFS,,,L2-L1-L
!VSBW,ALL
!WPOFFS,,,L
!VSBW,ALL
!WPOFFS,,,L3-L2-L
!VSBW,ALL
!WPOFFS,,,L
!VSBW,ALL
!WPCSYS,,0
! Top piezos
NUMCMP, volu
allsel,all
!---------
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
! FGM MATERIAL PROPERTIES
counter = 0
*do, j, -Hb/2,+Hb/2-v_numb_z,v_numb_z
z_coord = j
counter = counter+1
*use, mat_fgm_K,, z_coord
mp, ex, counter, E_fgm
mp, prxy, counter, nu_fgm
mp, dens, counter, den_fgm
VSEL,S,volu,,counter
VATT,counter,,1,0
vsel,all
*enddo
VSEL,S,LOC,Z,Hb/2,Hb/2+H !Selecting piezos
VATT,1,1,2,0 !Assigning piezo material props
vsel,all
lsel, s, tan1, x, -1
lsel, a, tan1, x, +1
lesize, all, , , ele_numb_x
lsel, s, tan1, y, -1
lsel, a, tan1, y, +1
lesize, all, , , ele_numb_y
lsel, s, tan1, z, -1
lsel, a, tan1, z, +1
lesize, all, , , ele_numb_z
lsel, all
vmesh, all
NUMMRG,KP ! Merges duplicate solid model
! entities
numstr,node
!-------------
!Go and merge some nodes (eg. use glue layers of materials together)
!----------------------
nsel,all
nummrg,node,1e-5
!----------------------
!Clamp left end of beam
!----------------------
nsel,s,loc,x,0 ! Clamp left end of bimorph
d,all,ux,0,,,,uy,uz ! Defines DOF constraints, clamp the end
nsel,all
!----------------------
!Define upper electrodes on piezos
!----------------------
!P1
asel,s,loc,z,Hb/2+H
asel,r,loc,x,L1,L1+L
nsla,s,1
*get,p1,node,0,num,min
cp,1,volt,all
d,p1,VOLT,0,0
allsel
!P2
asel,s,loc,z,Hb/2+H
asel,r,loc,x,L2,L2+L
nsla,s,1
*get,p2,node,0,num,min
cp,2,volt,all
d,p2,VOLT,0,0
allsel
!P3
asel,s,loc,z,Hb/2+H
asel,r,loc,x,L3,L3+L
nsla,s,1
*get,p3,node,0,num,min
cp,3,volt,all
d,p3,VOLT,0,0
allsel
!P4
asel,s,loc,z,Hb/2
asel,r,loc,x,L1,L1+L
nsla,s,1
*get,p3,node,0,num,min
cp,4,volt,all
d,p4,VOLT,0,0
allsel
!----------------------
!Define lower electrodes on piezos (bonded side)
!----------------------
!P1
asel,s,loc,z,Hb/2
asel,r,loc,x,L1,L1+L
nsla,s,1
*get,p5,node,0,num,min
cp,5,volt,all
d,all,volt,0,0 ! Ground bottom electrode,
allsel
!P2
asel,s,loc,z,Hb/2
asel,r,loc,x,L2,L2+L
nsla,s,1
*get,p6,node,0,num,min
cp,6,volt,all
d,all,volt,0,0 ! Ground bottom electrode,
allsel
!P3
asel,s,loc,z,Hb/2
asel,r,loc,x,L3,L3+L
nsla,s,1
*get,p7,node,0,num,min
cp,7,volt,all
d,all,volt,0,0 ! Ground bottom electrode,
allsel
!----------------------
!Define gravitational field
!----------------------
!acel,0,-9.81,0
fini
! ====================================================
!!!!!!!!!!!!!!!!!!!!!!!modal analysis!!!!!!!!!!!!!!!!!!!!!!
! ====================================================
/SOLU
ANTYPE,MODAL ! Modal analysis
nmodes = 10
modopt,LANB,nmodes ! Block Lanczos solver
mxpand,nmodes,,,yes ! Calculate element results and
SOLVE
*get,f1,mode,1,freq
FINISH
SAVE
/COM, **** REVIEW THE RESULTS OF THE MODAL SOLUTION ****
!
/POST1
SET,LIST