! TITLE: reaction.mac
!
! Find reaction loads due to applied force-couple at link arm.
!
clear_start_new

/prep7

constants

Pa=380              ! Applied load
Llink=.9
Rg=1.00             ! Sector gear pitch radius
Faxial=Pa*Llink/Rg  ! Axial thrust load through worm screw

! Properties that apply to all elements
E=10.5E6
I=.05
A=.385

! End plate (element #1) properties (left end)
w1=1.37
t1=.100
A1=w1*t1
I1yy=w1**3*t1/12
I1zz=w1*t1**3/12

! Round boss properties (right end)
OD=.46
ID=.2
A4=pi/4*(OD**2 - ID**2)
I4=pi/64*(OD**4 - ID**4)

et,1,beam4
keyopt,1,6,1
keyopt,1,9,9

mp,ex,1,E

n,1,0,0,0
n,2,.55,0,0
n,3,1.55,0,0
n,4,6.432-.42,0,0
n,5,6.432,0,0
n,6,1.55,0,-.787
n,7,1.55,.658,-.787

!AREA, IZZ, IYY, TKZ, TKY, THETA, ISTRN, IXX, SHEARZ, SHEARY, SPIN, ADDMAS
r,1,A1,I1zz,I1yy
e,1,2

r,2,A,I,I
e,2,3
e,3,4
e,3,6
e,6,7

r,3,A4,I4,I4
e,4,5

d,1,all,0
d,5,ux,0,,,,uy,uz,roty,rotz

f,7,fz,Pa
f,3,fx,Faxial

numbers_toggle,5
numbers_toggle,6
/view, 1 ,1,1,1
/ang, 1
/pbc,all, ,1
/rep,fast

fini

/sol
solve
fini

/post1

/out,reaction,out
prnsol,dof
prnld,,,all
presol,elem
/out