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