Try This one
finish
/clear
/VIEW,1,0,0,1
/RGB,INDEX,100,100,100,0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX,0,0,0,15
/UNITS,SI
V1=10.000000
Name=47.000000
alfp=0.400000
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! PARAMETERS !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
! GEOMETRY
*set, b, 0.005 ! [m]
*set, S, 20.
*set, a, 0.4 ! [m]
*set, l, 0.4 ! [m]
!
! 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
!
*set, E_b, 320.24e9 ! [Pa]
*set, nu_b, 0.26
*set,dens_b,3750
!
! MATERIAL BOT
*set, E_t, 105.70e9 ! [Pa]
*set, nu_t, 0.29
*set,dens_t,4429
!
! LOADING
*set, p, 1.
!
! FGM gradation law
*set, a_z, 1./b ! [m^-1]
*set, n_z, V1
!
! MESH
*set, v_numb_z, b/8 ! 20, 40
!transient analysis
*set,Tstep,0.01
*set,Tstop,2
*set,Tres,0.01
*set,ALPHAD,0
*set,BETAD,0.00015
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! 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-(alfp/2)*(z_coord+0.5)*(E_t+E_b)
nu_fgm = 0.3
den_fgm = (dens_t-dens_b)*grad_common + dens_b-(alfp/2)*(z_coord+0.5)*(dens_t+dens_b)
/gopr
*end !mat_fgm_K
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/prep7
!
!Damping
ALPHAD,0
BETAD,0.00015
! ELEMENT TYPE SELECTION
! ELEMENT TYPE SELECTION
et, 1, SOLID45
*do, j, -b/2,+b/2-v_numb_z,v_numb_z
block, 0, l, 0, a, j, j+v_numb_z
*enddo
/PREP7 ! Enter Preprocessor
*DIM, MyTable, TABLE, 1, 10, 1 ! Define a 4x3 table (4 rows, 3 columns, 1 var)
counter = 1
vsl=0
*do, j, -b/2,+b/2-v_numb_z,v_numb_z
z_coord = j
counter = counter+1
vsl=vsl+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,,vsl
VATT,counter,,1,0
vsel,all
*SET, MyTable(1,vsl), den_fgm
*enddo
vsel, all
ebig=2e-2
ESIZE,ebig
vmesh, all
nsel,all
nummrg,node,1e-5
ALLSEL
nsel,s,loc,x,0 ! Clamp left end of bimorph
d,all,ux,0,,,,uy,uz ! Defines DOF constraints, clamp the end
nsel,all
NSEL,ALL
finish
/SOLU
ANTYPE,MODAL ! Modal analysis
nmodes = 3
modopt,LANB,nmodes ! Block Lanczos solver
mxpand,nmodes,,,yes ! Calculate element results and
SOLVE
FINISH
SAVE
Finish
- Badges
- Science topic