srun.inp
来自「vibration simulation using ansys and mat」· INP 代码 · 共 1,672 行 · 第 1/5 页
INP
1,672 行
! srun.inp - suspension model
! documentation
! act.inp
! act.dwg
! output files:
! *.dat : transfer function ascii dump
! *.eig : eigenvector file
! *.grp1: transfer function plots
! *.grp2: pldi plots
! *.grp3: part plots
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 = 20 ! number of modes, eigenvalue expansion
! input variables for parameter studies
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 = 'suspr' ! watch length because of overall file name length limit
/filname,suspr,1 ! assigns name to all ANSYS files, opens new log and err files
name_ti = 'Reduced' ! title name
*else
filename = 'suspl' ! watch length because of overall file name length limit
/filname,suspl,1 ! assigns name to all ANSYS files, opens new log and err files
name_ti = 'Lanczos' ! title name
*endif
! input variables
! dynamics run
shkrmass = 1.0 ! 1 kg
shkrnode = 4407 ! node for shaker mass attach, front-middle weld at swage
gapnode = 22 ! gap node
shkrglvl = 1 ! shaker g level for transfer functions
zeta = 0.005 ! damping ratio, relative to critical damping
! view parameters, miscellaneous
! ************************************* /prep7 ********************************************
/title,Suspension Model
/prep7
/show,term
/view,1,1,-1,1 ! iso view
/angle,1,-60 ! iso view
/pnum,real,1 ! color by real
/numb,1 ! color on, numbers off
/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
!*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.71e-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
! real
r,1, ! slider
r,2, ! flexure
r,3, ! loadbeam
r,4, ! swage plate
r,5, ! actuator
r,6, ! bearing rings
r,7,kspring ! stiffness for bearing balls
r,8, ! coil bonding adhesive
r,9, ! coil
r,10,shkrmass,shkrmass,shkrmass ! shaker mass, ux, uy and uz components
! geometry, all dimensions in mm units
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
! 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
! input slider nodes, csys 21
csys,21
n, 1 , -1.498600, -1.100000, 0.000000
n, 2 , -1.036331, -1.100000, 0.000000
n, 3 , -0.690891, -1.100000, 0.000000
n, 4 , -0.345451, -1.100000, 0.000000
n, 5 , -0.195516, -1.100000, 0.000000
n, 6 , -0.000011, -1.100000, 0.000000
n, 7 , 0.195495, -1.100000, 0.000000
n, 8 , 0.345429, -1.100000, 0.000000
n, 9 , 0.690869, -1.100000, 0.000000
n, 10 , 1.036309, -1.100000, 0.000000
n, 11 , 1.346200, -1.100000, 0.000000
n, 12 , -1.498600, -0.876300, 0.000000
n, 13 , -1.036331, -0.876300, 0.000000
n, 14 , -0.690891, -0.876300, 0.000000
n, 15 , -0.345451, -0.876300, 0.000000
n, 16 , -0.195516, -0.876300, 0.000000
n, 17 , -0.000011, -0.876300, 0.000000
n, 18 , 0.195495, -0.876300, 0.000000
n, 19 , 0.345429, -0.876300, 0.000000
n, 20 , 0.690869, -0.876300, 0.000000
n, 21 , 1.036309, -0.876300, 0.000000
n, 22 , 1.346200, -0.876300, 0.000000
n, 23 , -1.498600, -0.658304, 0.000000
n, 24 , -1.036331, -0.658304, 0.000000
n, 25 , -0.690891, -0.658304, 0.000000
n, 26 , -0.345451, -0.658304, 0.000000
n, 27 , -0.195516, -0.658304, 0.000000
n, 28 , -0.000011, -0.658304, 0.000000
n, 29 , 0.195495, -0.658304, 0.000000
n, 30 , 0.345429, -0.658304, 0.000000
n, 31 , 0.690869, -0.658304, 0.000000
n, 32 , 1.036309, -0.658304, 0.000000
n, 33 , 1.346200, -0.658304, 0.000000
n, 34 , -1.498600, -0.431800, 0.000000
n, 35 , -1.036331, -0.431800, 0.000000
n, 36 , -0.690891, -0.431800, 0.000000
n, 37 , -0.345451, -0.431800, 0.000000
n, 38 , -0.195516, -0.431800, 0.000000
n, 39 , -0.000011, -0.431800, 0.000000
n, 40 , 0.195495, -0.431800, 0.000000
n, 41 , 0.345429, -0.431800, 0.000000
n, 42 , 0.690869, -0.431800, 0.000000
n, 43 , 1.036309, -0.431800, 0.000000
n, 44 , 1.346200, -0.431800, 0.000000
n, 45 , -1.498600, -0.215900, 0.000000
n, 46 , -1.036331, -0.215900, 0.000000
n, 47 , -0.690891, -0.215900, 0.000000
n, 48 , -0.345451, -0.215900, 0.000000
n, 49 , -0.195516, -0.195506, 0.000000
n, 50 , -0.000011, -0.276487, 0.000000
n, 51 , 0.195495, -0.195506, 0.000000
n, 52 , 0.345429, -0.215900, 0.000000
n, 53 , 0.690869, -0.215900, 0.000000
n, 54 , 1.036309, -0.215900, 0.000000
n, 55 , 1.346200, -0.215900, 0.000000
n, 56 , -1.498600, -0.000000, 0.000000
n, 57 , -1.036331, -0.000000, 0.000000
n, 58 , -0.690891, -0.000000, 0.000000
n, 59 , -0.345451, -0.000000, 0.000000
n, 60 , -0.276498, -0.000000, 0.000000
n, 61 , -0.000011, -0.000000, 0.000000
n, 62 , 0.276476, -0.000000, 0.000000
n, 63 , 0.345429, -0.000000, 0.000000
n, 64 , 0.690869, -0.000000, 0.000000
n, 65 , 1.036309, -0.000000, 0.000000
n, 66 , 1.346200, 0.000000, 0.000000
n, 67 , -1.498600, 0.215900, 0.000000
n, 68 , -1.036331, 0.215900, 0.000000
n, 69 , -0.690891, 0.215900, 0.000000
n, 70 , -0.345451, 0.215900, 0.000000
n, 71 , -0.195516, 0.195506, 0.000000
n, 72 , -0.000011, 0.276487, 0.000000
n, 73 , 0.195495, 0.195506, 0.000000
n, 74 , 0.345429, 0.215900, 0.000000
n, 75 , 0.690869, 0.215900, 0.000000
n, 76 , 1.036309, 0.215900, 0.000000
n, 77 , 1.346200, 0.215900, 0.000000
n, 78 , -1.498600, 0.431800, 0.000000
n, 79 , -1.036331, 0.431800, 0.000000
n, 80 , -0.690891, 0.431800, 0.000000
n, 81 , -0.345451, 0.431800, 0.000000
n, 82 , -0.195516, 0.431800, 0.000000
n, 83 , -0.000011, 0.431800, 0.000000
n, 84 , 0.195495, 0.431800, 0.000000
n, 85 , 0.345429, 0.431800, 0.000000
n, 86 , 0.690869, 0.431800, 0.000000
n, 87 , 1.036309, 0.431800, 0.000000
n, 88 , 1.346200, 0.431800, 0.000000
n, 89 , -1.498600, 0.768350, 0.000000
n, 90 , -1.036331, 0.768350, 0.000000
n, 91 , -0.690891, 0.768350, 0.000000
n, 92 , -0.345451, 0.768350, 0.000000
n, 93 , -0.195516, 0.768350, 0.000000
n, 94 , -0.000011, 0.768350, 0.000000
n, 95 , 0.195495, 0.768350, 0.000000
n, 96 , 0.345429, 0.768350, 0.000000
n, 97 , 0.690869, 0.768350, 0.000000
n, 98 , 1.036309, 0.768350, 0.000000
n, 99 , 1.346200, 0.768350, 0.000000
n, 100 , -1.498600, 1.100000, 0.000000
n, 101 , -1.036331, 1.100000, 0.000000
n, 102 , -0.690891, 1.100000, 0.000000
n, 103 , -0.345451, 1.100000, 0.000000
n, 104 , -0.195516, 1.100000, 0.000000
n, 105 , -0.000011, 1.100000, 0.000000
n, 106 , 0.195495, 1.100000, 0.000000
n, 107 , 0.345429, 1.100000, 0.000000
n, 108 , 0.690869, 1.100000, 0.000000
n, 109 , 1.036309, 1.100000, 0.000000
n, 110 , 1.346200, 1.100000, 0.000000
nall
nplo
! generate slider top nodes
nsel,s,node,,1,110
nmod,all,,,-dimptk
ngen,2,150,1,110,1,,,slidertk
nall
nplo
! input flexure nodes, non-tongue elements, csys 20
csys,20
n, 300 , -6.416125, -0.952500, 0.000000
n, 301 , -6.083300, -0.952500, 0.000000
n, 302 , -5.612377, -0.952500, 0.000000
n, 303 , -5.333989, -0.952500, 0.000000
n, 304 , -5.055602, -0.952500, 0.000000
n, 305 , -4.597357, -0.952500, 0.000000
n, 306 , -4.140200, -0.952500, 0.000000
n, 307 , -3.666882, -0.952500, 0.000000
n, 308 , -3.247657, -0.952500, 0.000000
n, 309 , -2.828432, -0.952500, 0.000000
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?