cantbeam.inp

来自「vibration simulation using ansys and mat」· INP 代码 · 共 139 行

/title, cantbeam.inp, 0.2 thick x 2 wide x 20mm long steel cantilever beam, 10 elements

/prep7

et,1,4            	! element type for beam

! steel

ex,1,190e6			! mN/mm^2
dens,1,7.83e-6		! kg/mm^3
nuxy,1,.293         

! real value to define beam characteristics

r,1,0.4,0.1333,0.0013333,0.2,2		! area, Izz, Iyy, TKz, TKy

! define plotting characteristics

/view,1,1,-1,1      ! iso view
/angle,1,-60        ! iso view
/pnum,mat,1         ! color by material
/num,1              ! numbers off
/type,1,0           ! hidden plot
/pbc,all,1          ! show all boundary conditions

csys,0				! define global coordinate system

! nodes

n,1,0,0,0			! left-hand node
n,11,20,0,0			! right-hand node

fill,1,11			! interior nodes

nplo

! elements

type,1
mat,1
real,1
e,1,2
egen,10,1,-1

! constrain left-hand end

d,1,all,0			! constrain node 1, all dof's

! constrain all but uz and roty for all other nodes to allow only those dof's
! this will give 10 nodes, node 2 through node 11, each with 2 dof, giving a total of 20 dof
! can calculate a maximum of 20 eigenvalues if don't use Guyan reduction to reduce size of
! eigenvalue problem, maximum of 10 eigenvalues if use Guyan reduction

nall
nsel,s,node,,2,11
d,all,ux
d,all,uy
d,all,rotx
d,all,rotz

allsel
nplo
eplo

! *******************************  eigenvalue run  ******************************************

fini						! fini just in case not in begin

/solu						! enters the solution processor, needs to be here to do editing below

allsel						! default selects all items of specified entity type, typically nodes, elements

nsel,s,node,,2,11
m,all,uz

antype,modal,new
modopt,reduc,10				! method - reduced Householder, number of modes to extract
expass,off					! key = off, no expansion pass, key = on, do expansion
mxpand,10,,,no				! nummodes to expand
total,10,1					! total masters, 10 to be used, exclude rotational dofs

allsel

solve						! starts the solution of one load step of a solution sequence, modal here

fini

! plot first mode

/post1

set,1,1

pldi,1


! *******************************  output frequencies  **************************************

/output,cantbeam,frq         ! write out frequency list to ascii file .frq

set,list

/output,term				! returns output to terminal
 
! *******************************  output eigenvectors  *************************************

! define nodes for output:  forces applied or output displacements

nall
!nsel,s,node,,11				! cantilever tip

/output,cantbeam,eig         ! write out eigenvectors to ascii file .eig

*do,i,1,10
	set,,i
	prdisp
*enddo

/output,term

! *******************************  plot modes  **********************************************

! pldi plots

/show,cantbeam,grp,0
allsel

/view,1,,-1,,			! side view for plotting
/angle,1,0
/auto

*do,i,1,10
	set,1,i
    pldi
*enddo

/show,term