arunpz.inp

vibration simulation using ansys and matlab 一书中的程序

!/batch

!	arun.inp - actuator model input
	
!	documentation
!		act.inp
!		act.dwg

!	output files:
!		*.dat		transfer function ascii dump
!		*.eig		eigenvector file
!		*.grp1		transfer function plots
!		*.grp2 		pldi plots
!		*.amass  	mass properties from static run

fini							!	cleans up any leftover runs


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  input required  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! eigenvalue extraction technique, 1 = reduced, 2 = block lanczos	 

eigextrm = 2

! which analyses to run, 1 = run, 0 = don't run

eigprob 	= 1				! eigenvalue problem
freqresp 	= 1				! frequency response
massprop 	= 0				! mass properties

! eigenvalue extraction technique, 1 = reduced, 2 = block lanczos	 

nummodes = 50			   	! number of modes, eigenvalue expansion

! input variables for parameter studies

!zht = dimptk + lbtk		! flat suspension
zht = 1.200					! distance from disk to arm/suspension interface
bump = 0.025				! "bump" along spring
offset = 0.025				! "kink" in spring, measured at slider end of spring

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  input required  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*if,eigextrm,eq,1,then

	filename = 'actrrpz'				! watch length because of overall file name length limit
	/filname,actrrpz,1				! assigns name to all AnSYS files, opens new log and err files
	name_ti = 'Reduced'				! title name

*else

	filename = 'actrlpz'				! watch length because of overall file name length limit
	/filname,actrlpz,1				! assigns name to all AnSYS files, opens new log and err files
	name_ti = 'Lanczos'				! title name

*endif

/title,Piezo Actuator/Suspension Model

! input variables

! dynamics run

coilfor1 = 25032
coilfor2 = 25037
coilfor3 = 25011
coilfor4 = 25016
gapnode1 = 22					! upper head gap node
gapnode2 = 10022				! lower head gap node
kball = 5.184E6					! mn/mm for 4 balls, nMB L950 Bearing

zeta = 0.005		   			! damping ratio, relative to critical damping

! view parameters, miscellaneous

! *************************************  /prep7  ********************************************

/prep7
/show,term
/view,1,1,-1,1     			! iso view
/angle,1,-60       			! iso view
/pnum,real,1       			! color by real
/type,1,2          			! hidden plot, centroid sorting
/pbc,u,,1	   				! show displacement bc
/pbc,m,,1					! show moment bc
/pbc,cp,,1					! show couple equations
/pbc,f,,1					! show applied forces
/triad,off         			! turns triad off
/uis,msgprop,0				! all warnings and errors displayed
/graphics,off				! powergraphics off
/numb,1   		         	! 1 is color on, numbers off, -1 is no colors or numbers
!*afun,deg					! degrees, not radians
shpp,on						! turns on error checking

! element types

et,1,45						! solid45 - 3d structural solid, brick
et,2,14,	          		! combin14, spring-damper, used as spring for bearing balls
et,3,63,					! plate elements
et,4,21						! point mass for shaker

! material properties - mn/mm^2, 

! stainless steel 300
mp,ex,1,180.0e6
mp,dens,1,8.0718e-6
mp,nuxy,1,0.32

! aluminum
mp,ex,2,71.0e6
mp,dens,2,2.77e-6
mp,nuxy,2,0.345

! ceramic for slider
mp,ex,3,392e6
mp,dens,3,4.24e-6
mp,nuxy,3,0.23

! coil
	
mp,ex,4,75.8e6
mp,dens,4,6.00e-6
mp,nuxy,4,0.343

! coil bonding adhesive

mp,ex,5,2.0e6
mp,dens,5,1.0e-6
mp,nuxy,5,0.3

! pzt material

mp,ex,6,5.7e6
mp,dens,6,7.7e-6
mp,nuxy,6,0.3

! real

r,1,								! slider
r,2,         						! flexure
r,3,								! loadbeam
r,4,								! swage plate
r,5,								! actuator
r,6,								! bearing rings
r,7,kball							! stiffness for bearing balls
r,8,								! coil bonding adhesive
r,9,								! coil
r,10,shkrmass,shkrmass,shkrmass		! shaker mass, ux, uy and uz components
r,11								! pzt material

! geometry, all dimensions in mm units, gives xc = -0.65299mm

tparmtpz = 7.450					! top arm top z dimension
tparmbtz = 6.075					! top arm bottom z dimension
btarmbtz = 1.050					! bottom arm bottom z dimension
coiltpz = 7.275						! coil top z dimension
coilbtz = 5.750						! coil bottom z dimension
coilclz = 0.5*(coiltpz+coilbtz)		! coil c/l dimension

lbtk = 0.0762				! loadbeam thickness
fltk = 0.0508				! flexure thickness
swagetk = 0.3048			! swage plate thickness
bosstk = 0.6477				! boss height
slidertk = 0.6				! slider thickness
dimptk = 0.7524				! height of dimple top above disk
disktk = 1.25				! disk thickness

lbtipx = 0.2540				! dimple to loadbeam tip
flgx = 12.9540				! dimple to loadbeam flange end
spgx = 15.494				! dimple to swage end of spring

csys22z = 4.250				! csys0 at bottom of actuator to disk c/l

! input variables for parameter studies

!zht = dimptk + lbtk		! flat suspension
zht = 1.200					! distance from disk to arm/suspension interface
bump = 0.05					! "bump" along spring
offset = 0.025				! "kink" in spring, measured at slider end of spring

! calculate bump/offset beam dimensions

! loadbeam bottom at dimple to loadbeam bottom at swage plate

dzlb = zht - lbtk - dimptk		

! angle of loadbeam without offset wrt csys 20 at dimple, degrees

lbang0r = atan(dzlb/spgx)

! loadbeam bottom at dimple to loadbeam bottom at offset measurement point, no offset

dzof0 = flgx*tan(lbang0r)

! loadbeam bottom at dimple to loadbeam bottom at offset measurement point, with offset

dzof1 = dzof0 + offset

! calculate new angle of loadbeam

lbang1r = atan(dzof1/flgx)
lbang1d = lbang1r*180/3.141596

! loadbeam bottom at offset measurement point to loadbeam bottom at swage plate

dzspng = dzlb - dzof1

! local 20 coordinate system at dimple/loadbeam contact

csys20z = csys22z + 0.5*disktk + dimptk

local,20,0,52.293794,0,csys20z,0,0,lbang1d

! local 21 coordinate system at dimple/loadbeam contact, parallel to disk

local,21,0,52.293794,0,csys20z,0

! local 22 coordinate system for reflecting upper suspension

local,22,0,52.293794,0,csys22z

! local 23 cylindrical coordinate system at x location of swage boss c/l at csys 0 z location

local,23,1,34.259804,0,0

! input actuator nodes

csys,0

n, 20001 ,-31.196109,  0.000000,  0.000000
n, 20002 ,-31.154875, -1.603435,  0.000000
n, 20003 ,-31.031281, -3.202632,  0.000000
n, 20004 ,-30.825654, -4.793362,  0.000000
n, 20005 ,-30.538537, -6.371420,  0.000000
n, 20006 ,-30.170690, -7.932636,  0.000000
n, 20007 ,-29.924915, -8.444614,  0.000000
n, 20008 ,-29.506420, -8.828528,  0.000000
n, 20009 ,-28.975194, -9.029343,  0.000000
n, 20010 ,-28.407387, -9.018275,  0.000000
n, 20011 ,-26.218205, -8.551795,  0.000000
n, 20012 ,-24.029023, -8.085315,  0.000000
n, 20013 ,-21.839841, -7.618834,  0.000000
n, 20014 ,-19.650659, -7.152354,  0.000000
n, 20015 ,-17.461477, -6.685873,  0.000000
n, 20016 ,-15.272295, -6.219393,  0.000000
n, 20017 ,-13.083113, -5.752913,  0.000000
n, 20018 ,-12.611490, -5.564512,  0.000000
n, 20019 ,-12.229765, -5.229535,  0.000000
n, 20020 ,-11.895722, -4.786381,  0.000000
n, 20021 ,-11.895722, -4.285849,  0.000000
n, 20022 ,-11.895722, -2.940254,  0.000000
n, 20023 ,-11.895722, -1.579227,  0.000000
n, 20024 ,-11.895722,  0.000000,  0.000000
n, 20025 ,-11.895722,  1.579227,  0.000000
n, 20026 ,-11.895722,  2.940254,  0.000000
n, 20027 ,-11.895722,  4.285849,  0.000000
n, 20028 ,-11.895722,  4.786381,  0.000000
n, 20029 ,-12.229765,  5.229535,  0.000000
n, 20030 ,-12.611490,  5.564512,  0.000000
n, 20031 ,-13.083113,  5.752913,  0.000000
n, 20032 ,-15.272295,  6.219393,  0.000000
n, 20033 ,-17.461477,  6.685873,  0.000000
n, 20034 ,-19.650659,  7.152354,  0.000000
n, 20035 ,-21.839841,  7.618834,  0.000000
n, 20036 ,-24.029023,  8.085315,  0.000000
n, 20037 ,-26.218205,  8.551795,  0.000000
n, 20038 ,-28.407387,  9.018275,  0.000000
n, 20039 ,-28.975194,  9.029343,  0.000000
n, 20040 ,-29.506420,  8.828528,  0.000000
n, 20041 ,-29.924915,  8.444614,  0.000000
n, 20042 ,-30.170690,  7.932636,  0.000000
n, 20043 ,-30.538537,  6.371420,  0.000000
n, 20044 ,-30.825654,  4.793362,  0.000000
n, 20045 ,-31.031281,  3.202632,  0.000000
n, 20046 ,-31.154875,  1.603435,  0.000000
n, 20047 ,-31.196109,  0.000000,  0.000000
n, 20051 ,-32.696109,  0.000000,  0.000000
n, 20052 ,-32.652892, -1.680533,  0.000000
n, 20053 ,-32.515429, -3.432560,  0.000000
n, 20054 ,-32.307841, -5.023841,  0.000000
n, 20055 ,-32.006919, -6.677777,  0.000000
n, 20056 ,-31.621385, -8.314060,  0.000000
n, 20057 ,-31.129834, -9.338017,  0.000000
n, 20058 ,-30.292844,-10.105844,  0.000000
n, 20059 ,-29.230393,-10.507475,  0.000000
n, 20060 ,-28.094778,-10.485339,  0.000000
n, 20061 ,-25.905597,-10.018859,  0.000000
n, 20062 ,-23.716415, -9.552378,  0.000000
n, 20063 ,-21.527233, -9.085898,  0.000000
n, 20064 ,-19.338051, -8.619418,  0.000000
n, 20065 ,-17.148869, -8.152937,  0.000000
n, 20066 ,-14.959687, -7.686457,  0.000000
n, 20067 ,-12.770505, -7.219976,  0.000000
n, 20068 ,-11.827258, -6.843176,  0.000000
n, 20069 ,-11.063807, -6.173222,  0.000000
n, 20070 ,-10.549975, -5.235440,  0.000000
n, 20071 ,-10.395722, -4.285849,  0.000000
n, 20072 ,-10.395722, -2.940254,  0.000000
n, 20073 ,-10.395722, -1.579227,  0.000000
n, 20074 ,-10.395722,  0.000000,  0.000000
n, 20075 ,-10.395722,  1.579227,  0.000000
n, 20076 ,-10.395722,  2.940254,  0.000000
n, 20077 ,-10.395722,  4.285849,  0.000000
n, 20078 ,-10.549975,  5.235440,  0.000000
n, 20079 ,-11.063807,  6.173222,  0.000000
n, 20080 ,-11.827258,  6.843176,  0.000000
n, 20081 ,-12.770505,  7.219976,  0.000000
n, 20082 ,-14.959687,  7.686457,  0.000000
n, 20083 ,-17.148869,  8.152937,  0.000000
n, 20084 ,-19.338051,  8.619418,  0.000000
n, 20085 ,-21.527233,  9.085898,  0.000000
n, 20086 ,-23.716415,  9.552378,  0.000000
n, 20087 ,-25.905597, 10.018859,  0.000000
n, 20088 ,-28.094778, 10.485339,  0.000000
n, 20089 ,-29.230393, 10.507475,  0.000000
n, 20090 ,-30.292844, 10.105844,  0.000000
n, 20091 ,-31.129834,  9.338017,  0.000000
n, 20092 ,-31.621385,  8.314060,  0.000000
n, 20093 ,-32.006919,  6.677777,  0.000000
n, 20094 ,-32.307841,  5.023841,  0.000000
n, 20095 ,-32.515429,  3.432560,  0.000000
n, 20096 ,-32.652892,  1.680533,  0.000000
n, 20097 ,-32.696109,  0.000000,  0.000000
n, 20101 ,-34.196109,  0.000000,  0.000000
n, 20102 ,-34.150910, -1.757631,  0.000000
n, 20103 ,-34.007139, -3.590036,  0.000000
n, 20104 ,-33.790028, -5.254319,  0.000000
n, 20105 ,-33.475301, -6.984133,  0.000000
n, 20106 ,-33.072080, -8.695484,  0.000000
n, 20107 ,-32.334754,-10.231420,  0.000000
n, 20108 ,-31.079269,-11.383160,  0.000000
n, 20109 ,-29.485592,-11.985607,  0.000000
n, 20110 ,-27.782170,-11.952403,  0.000000
n, 20111 ,-25.592988,-11.485922,  0.000000
n, 20112 ,-23.403806,-11.019442,  0.000000