Technical Report of Mathematical Models and Engineering Analysis of Howitzer
Technical Report of Mathematical Models and Engineering Analysis of Howitzer
Technical Report of Mathematical Models and Engineering Analysis of Howitzer
FOR
4 j ENGINERIMLiSIS NOESIGN
2~SQF ROIL-)M2TJ
HWI 19LIGHT.#,TOWED-'19
TECHNI,,CAL. REPORT
~
- ~MiQ'-IAEL C. NERDH
JERRYW FRANTZ
MAY1973
~,07 009OO
Distribution, 1 d- t S.
S Governnt
agencies only; Ta vauaion of itay
Hardware pres ;Hay
requsts fo is documont imst be o ered
e
_to Arti Weapons Sys toew-Directora
U. Weaponsf Commnd.
DISPOSITIO1K INSTMUCTIONS
TECHNICAL REPORT
R-RRA-S-2-2 8-73
MATHEMATICAL MODELS
FOR
ENGINEERING ANILYSIS AND DESIGN
OF
HOWITZER, LIGHT, TOWED; 105MM SOFT RECOTT,, XM204 (
By Vq,
Michael C. Nerdahl
Advanced Systems Division ii
Artillery Weapon Systems Directorate
andI
Jerry W. Frantz
Systems Performance Division "
Research Directorate* •
DISTRIBUTION:
D ibution limited to U. . Go nt
agenc Test and EvSw on of
military har ed; May 1973.
Other requests his t must be
referred rtillevy Weapon Sys
Dir ate, Rock Island Weapons Laboratory.
- *- ' -
[ i
".
-.
..,
..
.cx
..
,..
ti
....
€,
.,
#;,
.
..
..
.-
i.x
.
..*
'"<"
"<"'2,
. ."":,
:,.,..
....
)zt.
,.:,,
•.:..
..2.
2.•
ABSTRACT
FOREWORD
r.AnFor
h#. ,- - ,f
iiIi ::i!j
CONTENTS
FUNCTIONS ..................................... 26
REFERENCES ....................................... . 79
APPENDIX ..................
*..... .......... 88
Siii I
k
LIST OF ILLUSTRATIONS
iv
-. d .7 - a, .W
INTRODUCTION
TIME
saC
.ATT.RY
POSITION
CONVENTIONAL RECOIL
0.0 , PULL LANYARD TO FIRE
RECOIL
. MAXIMUM RECOIL
I .COUNTER RECOIL
,.1- END OF CYCLE
.,,"r.SOFT RECOIL
B
POS:TION
- COUNTER RECOIL
1.4 END OF CYCLE
FIGURE 1
L1
Developoent of a reliable soft-recoil mechanism for modern
artillery began in 1957 with the modification of an MI01 Howitzer,
105mm. While evaluation tests of this mouified weapon demonstrated
the feasibility and value of a soft-recoil system, mechanism
reliability was unacceptable. (Reference 1). In March 1964,
design and fabrication of an experimental firing fixture having an
improved recoil mechanism was initiated. Extensive firing tests of
this fixture were conducted to confirm feasibility of this weapon
concept, to determine accuracy and durability characteristics, and
to idantify and examine functional problems. (References 2 & 3).
___k,_V_
...
Design of the Engineering Development Prototypes that
will be used for DT /OT = has been initiated. As a result
of the changes in the cannon and the increase in ammunition impulse
levels, new recoil and counterrecoil orifice designs were required
and the effect of carriage flexibility had to be reaqse3sed.
Consequently, the original mathematical models were modified
to more closely represent the current design configuration. At
the same time, an attempt was made to generalize the form so that
minor alterations would allow for their use in the analysis of
future soft-recoil weapon concepts.
|4
"i!
II
CLi
U6F
00
0'0
.61.
I z 3
I
7A1
4-
ILA
I C6
yf z
a. F
'V - t.. C4
P3
FiguA
5a
Reloti2
Poito n irnDd~WVUie no
xO
As how, Fr AX 6 ,a 1 ~ 2 Iforxx.a: 2
aUII ONTRICI
Figure
ReltivPoitin flbude
F andDigram iigIgaaig
(As how , XO
Fo A~X,,ts,<,,;Forxjx,, a=A2
RIGID BODY MODEL OF A SOFT RECOIL MECHANISM
Lix
recoiling parts (as in the case of a cook-off). This spear buffer
restricts fluid flow between the pressures P and P2 during part
of the cycle. As illustrated schematically in Figure 5, this
restriction exists thile
>X IIx
-with "x " being the value of "x" for which a becomes equal to A2
For any "x", the flow diagram is shown in Figure 6 with flow
from P to P4 when "x" is increasing (defined as recoil). The direction
of flow is reversed when "x" is decreasing (during run-up and during
counterrecoil). Values for the orifice areas (a.'s)
1 will vary with
direction of fluid flow, position of the spear buffer relative to
the piston, and functioning of the velocity sensor.
Then,
(2) P1 - P2 - g(vi)
(3) F2 - P3 g(v 2 )
(4) P - P
3 4 g(v 3 )
17 i
.............
..... ,.
.......
..
.,
wie re
NA i Na v
(5) v =
1 a1
Nagz 1 + NA2 a Na 2 v 2
Na 2 v 2 =NAx + a 3v 3
N
(7) v3 = a x *
.,42
4
Now,
a3 v3 - A4 Cx- ,)
NARX Ax
4- 4y
A4
(8) " A4 N
A4 NA R
A4_-
A4
(10) A 4 - NAR
A4
(11) P3 P +
&(v 3 4 3
= + + +
(12) P2 P g(v) - P g(v) g(v3)
( 2 P2 3 2 4*4 2 + 3
From the free body diagrams (Figures 2 and 3), with "x zx (spear
buffer rescricting fluid flow)
4
~ANPN+ 4P 4 +NA 3 P 3 NA 2 2 11P
t' lI,
9
" II I i i I I i i I i I ii I
and
"
(13) Mp V Wp sin + APN + Ff. sgn(- A4P 4
A PN + A4 P 4 + NA 3P 3 - NA2 P 2 - NAP I
A
(15.1) M 4 R "x W sin N + sgn() A .+
SA p PN Ffp A4p4
r + A4 - NAR 1
(16) A MJ x =B(t) + (W +W sin~
- (NFp + F ) sgn(l)
(NF + F ) sgn(k) + NA 3 P4
10
r
Collecting terms
[M1 A NM
)4 x B(t) + (WR + Wp) sinY
- (Np + F agn(i)
p g
-N(A-
1 + A- A3 ) p4
- N(Al I A2 -A 3) g(v3)
(15.2) P
4
iW
A4 p
sin X+ --A4 F
fp
sgn(;)
AN - M NAR
A N AH
44 4
R
_714 M B(t) + (WR + Wp) gin
-(NF
p
+ F)
9
sgn(;) - A4~w~ sin
NA R NA
A4 fp A PN
4 4
NAAN A4 - NA R
AA
44
- - (A+
(g(v) A2) g(v 2 ) - N g(v3 )
N g3
L1
' : .- : ..
, ! - 7 ' - . :.-. ,.. '-
.. .. ,' * .,
or, rearranging tetva
NJ M NAR)
If the initial value is VV the gas volume for any displacement "x"
is written as
VN = vo - N (x -Y) -N A4
Then, assuming adiabatic gas laws, the gas pressure for any displacement
"x" is determined from
NN
k
0 0
k
IvI. I
~(~)k
Voo
(18) PN o° _ _ _ __
0 x
A
"4
RN
From the free body diagram (Figure 4), the force cn the recoil rod is
given by
(19) R - NA P+ NA 2 - NAP + % gn()
33 p
N(Ai+A 2) - Nlv
AA2
)= a1 x
1'
12 p
I'
Therefore, for "x <~ x e", Equation 14.1 should be written as
+
-NPN A4p4 + N A3p3
- NAp - UlPl
]
- (NF + F) sgn(i)
p g
- (NFp + ) ( )9
+ NA P +NA
S 4 A3 g(v3 )
Al(.
+ A2) P4
- N(A+) P- ".%+%
N(A + A) {1l)
g(V
13
1.
Ij
nn.i n u - . . . . - - -, -
Therefore, collecting terms and noting that AR -A, I A2 -A 3
-(NF + F )sgn(i)
p g
-N P - N(A1 + A2 ) g(v,)
(17.1) L+ A4 NA)MJ
+( B(t) +
+ WR +4 A4t )Wpjsn
sinI
NARAR
-
A N
4
NA )[gv+gv)
- '1 + 9( 12 g~2)
- AR g(v3 )
By defining
H 1i-if x x
(20'
Hii 0 if x xe
14
Equations 5 and 5.1 may be written as
(5.2) v1 a X
ai'
+ I N
- (NFP + Fg + Ffp) s 1
L N(A HA~g(v)
- g(v 3 )
Then, by using Equations 5.2 and 18.2 along with Equation 20,
a single system of equations defines the entire cycle as long as
he eanism remins filled wit fluid.
=
TheP'3 P2
' 'l -o0
4 (v 3 RP+-
15
and Equations 14 and 15 must be rewritten as
+ A4 g(V3 )
with
or
(23) 3 a4 _
and
P
(24) PN P
(25) R- Y sg(i)
p
16
I
The governing differential equations of motion are solved
by standard numerical methods, with the system controls of the
actual mechanism being simulated by logic decisions of the computer.
These logic decisions are defined in the following manner.
• 17
The orifice area may be defined as
+al..k
+" af Xe
1
a"+l f x> x" and "x < 0"
ax +al..k if WX> xe i n i
By making "a A2 " for "x < Xe" the definition siuplifies to
a +a of
x leak unless "x > x 7 x"and " > 0"
a ,
7. Define
H -
>
P3 0 or I3 !: 0
18
f.i
'I
TABLE I
A. EQUATIONS OF MOTION
-(N? + F
pg
+ NRF
A4 fp sn)
- NAR g(v3 )
1 if "x c x
(20) H
,,,
0 if fix >__ x e
(18 ) PN pP .. . 0
V I
o
(:vl ) 2
A4
(v)
(1) g(vi) " "gn(v
(5.2) va A +a
19
~...
TABLE I (cont'd)
A + HA
1 2
(6) v2 " 2
- 8
(7)
A4 - NAR
(i0) y - A4 Z
4 4
A N
mA 4
(11) P3 -P 4 + g(v3 )
- fp ( - - - A 4 S(v 3 )
+ g(v3 )
A4
(23) v a3
20
TABLE I (cont'd)
(24) PN o ( V ANxY
f
a~b whenever x> 0
a,
21
. . . . . . ..
A.
DESIGN OF RECOIL AND COUNTERRECOIL CONTROL ORIFICE
(17) [ (A 2
4: P mB(t) +R + (A 4 -R]
- (NF + Fg + F)
-RA sgP)
p p A4 N
L(26)
This equation has the form
where
M~ x A(t) - D(t)
(27) H (AM
(28) W -( 4 A ) p sin
(5) v1 -
a1
(6) V a2 x
(7) v.3 -
i 22
3
~i.
NA"A
+
A4 PN -N(Al A2) g(v2 )
- NA 3
(32)
v sgnvg)
and (32) ~0 (KVnTC
gvvl (v7) sngv
and
Ax
1
(33) aI
v1
NARAN
A4 N NAI g(vl)
23
41 d.
with "x < 0" (during counterrecoil), this may be written
(after assuming c :; C 2 0 c3 " ) as
(35) M(NF+ F + _ F
Meffh p 9 A4 fp
NARAM
-- P
N
A4
3 (NA) 2
[. NAI N(A 1+2 A)
2 a a 2
29c 2 al
3
NA N(A, + ANA):' .2
2 a12 + 2 a 2
2g[ 2
(6 23af2
NAR N"RA
+NF+F - Ff) - pN
p g A4 fp A4 N
(38) NA jP):2(A)
1 T23a f
24
Actually, the minimilm value for a1 is determined by the
necessary clearance between the spear buffer and the recoil piston
head. Then, at latch position,
a aleak
a ad
25
io
. D(t) Aft)
x |
FIGURE 7
Free Body Diagram of Recoiling Mass
Since the firing cycle begins with the recoiling parts at rest, the
initial conditions are:
x 0 and - 0 when t - 0
26
Le . . . ... ,', , , . ... . ,-. .. ,.-. . .. 1,, ,,.' ... o . . -- , ,,! :'.z'. , . - : : 2 *, C -- M'
Mosent Area Equations
ti
0
t- Sw~) D(t) dt
'JJ 0
where
ti
(43) A(t) dt - Area under curve of A(t) from t - 0 to t t
J
t
(44) D(t) dt - Area under curve of D(t) from t - 0 to t - ti
0
i. 0
27
- - --- - --. -
D er a Value of D(t) at end of recoil (t tr)
Ar
r - Specified rise and fall time for D(t)
B(t) t -t r A
eAccelerating forces
Decelerating forces
FIGURE 8
Assumed Force System for a
28
I.
k !
Ir, I
Equations 41 and 42 say be solved situltaneously for values
of the umknown DR and t .
Meff j+ + De- Do
(47) t 6 r
r 21 + (De -D O) ar
I + Wef t r -5 (D° + D e )A
I ff rt
(48) D =
a FAccelerating forces
(t) '
tf
eff
-t
D
t t
r
Decelerating forces
FIGURE 9
Assumed Force System for an
Ideal Soft Recoil Cycle
29
with
x - 0 and 0 at t - 0 (Battery)
z- Z and a x at t - tf (Firing)
The momnt area equa ione (Equations 41 and 42) can be solved
for
t - Cycle time
Xf Firing velocity
(49) D - 4( f I + i) 'f
eff
(50) t - D
--
1f ++
(51) f 'N ff %ef f
30
Villi i i i -ail
(3) returned to latch by the gas force.
31
"°-": .-. . . . . . . -
i
2 "
4A ~ ~ ~ p 4f )
- P - MA (P-P 22)
A4 N1U
PN =: P4 P3 -2 P2
2
(53) + H (N~P + F + 2i i
NA A P - NA (Pl- PN)
A4 N
N"
A4 N + AlPN
. N %i
D- (N + F+ +"- F)
p g A4 fp
32
With counterrecoil control occurring near the latch position where
PN - P , a liniting design value for "D" is given by
)
FN(A - A-A - P-
(54) D- (NF + F + N>
p S A4 fp +N A 4j- ~0
The shape oL' the control function r.-y be chosei,n some essentially
arbitrary fashion so long as adequate flow control can be maintained.
One such choice is shown (Figure 10).
TT
Beginning of c_
counterrecoil
control by spear
buffer
b t
FIGURE 10
Counterrecoil Control Function
1
ii
3- .
Then, from "t 0 to t T "1 otion i defined by
21
2CR
Me f f z " T ' - t
where
H +,
DCI N~ $ A4 LA" 4 o
M eff 2 A4
NOiE: .'t' M
s aseewmed zero when counterrecoil control begins.
TheR,
DA R
+ + K1
DCR t 2
Ieff x - + K1 t+K 2
D3
T 3
Since
X x v and x x. at t 0
K1 eff immn
K2 "Meff Xv
34
II
Therefore, at t aandx-X
D 2
CR T
T. 4 a ft =in
4 ef f M
and
M X DCR
eff~~
- T38 + ef M1ff xv
DCRmin + Tx +x
24 Me 2
2 DCR
eff CR TC
with
x - Kand x - X at t
NOTE: Here, "t" has been assumed zero when the control function
begins
to decrease.
2DD 2
DCR t3
Meft x DCR t 2-2 -T -3 1" Kl t+ K2
2. ? eff
35
M
As originally specified, when "D" retutes to zero
ZI and x
=1- -
Thcn, ac t T
T
2 I
- D T DCR T2
Meff 2 CR f -4
T
DCRT
+* eff 414ef f + Meff
m ;i n
DT
(55) -T DCR T
aff
CR T CR T
eff XT 2 3T -
+
+ 4 T
DCR T Mef2f Xmin T
PC T
2
T T
eff T ' 24 e T4o1
36
-=
.... -
D CR T 2
2 T
CR DCR
Me
f T - + Mef T T -2 + Meff xv
(57) DCR 2
4- Meff '- T - Meff (x v - ) -0
-A7
Z
o
wz 0
WU
x
lAM
AI
XEz
0.38
CARRIAGE RESPONSE TO FIRING LOADS
T-n; :A coordinate system rotating with the mass N and having its
origin at the elevation trunnion. At the time, t - 0, this
system has been rotated through the elevation angle (Y)
from its original position parallel to the 0-YZ coordinate
system.
YA a YT + nA cOs (Y + ) -A sin (Y + 6)
ZA - ZT + n s n Y
(Y + ) +A Coo (Y + 8)
39
I
.. -... ,.-.
.. .....
....-
. ... . . ... .. ". -.. -' -;
0 'A
400
BA -YT COs -ZT i 0 + nA cos (y + 6+ ) "CA sin (y + 8 + *)
CA W YT sin 0 + ZT Cos 0 + nA sin (y + 8 + 4) + CA cos (y +48+ )
Similarily:
YB - YT + nB cos (y + 0) -B sin (y + 8)
ZB - ZT + "B (Y + 6) + B cos (y + )
B = YT Cos 0 - ZT sin 0 + nB cos (y + B + ) -C. sin (Y + B + *)
C - Y sin * + ZT cos + nB sin (Y + + 0) + B (Y +8+ *)
And,
3 C(15-Zy
D -YD Bil
C YD sin + ZD cos 0
BA - - s in 0 + k Cos
+ C; + I nA Cos ( + 0 + 0) - A sin (y + 8 + )]
+ A sin (y + a + 0)
JY sinl * + 7T coo 0]
- (8 + ) [ sin C, + 8 + ) + CBcoo (y + 6 + ,)
LYT Co 0- 7 sin 0
+ (0 + coo (y + + 0) S ( + 8 +*)]
41
Differentiating again wiUh respect to time gives:
- (8 + $) [2; A sAin (Y + + ]
S
;2 [YT coS 0 - ZT sin *]
-(8 + $)2[InA cos (y + 8 + 0)-CAsin (y + + )]
* sin 0]
CA" 0 [YT cos - ZT
+sA Sin (y + + )
+(t + $) [2; COA C + + )J
-2[Y T sin + ZT Cos
BB 0 IT sin - + ZT CosB0
- (8 + * 1
r
,B sin (y + e + 9) + CB cos (y +8 + ,)]
- T c
coS - ZT sin 0]
IA
42
BD " [YD sin 0 + ZD coo 4]
.2
- c[D coo - ZDs
F'D = *[D Co - ZD n ]
- .2
[D sn + ZD cos
By defining
(8) 3T]- F
1 •2 .2
- MA(BA +CA]
1 • 2 .2
+ Y MB EBB +CdBI
1 2 .2
2 KDLBD +CDI
1 2
+ 2 IA+ IB3 [; +
1 .2
+ ID
43
*
After terms are expanded and collected, the following expression
for kinetic energy is obtained:
2
(60) T A [2 T2 + ZT2 + n 2 A2+ (nAYT + CAZT) con (y + e)
2
+ (nAZT - CAYT) sin (y + 01J
.2 2 2 .2
A+n +
* A ]+ nA
+
AZT
~Ysin
AYT
(Y + e)]
2 7 2 (Y + cos (y
- Ae A - A A~;- T si 8) + zT + 8))
N [2 2 2 2
+ - 'T +ZT + nB + C + 2 (n BYT + tsrB ) Cos (Y + 8
.2 2 2
2
M 2 2 1_+_B_2_D_
Y +~
.4zI + 2 + -n.;
D 2 2
441
Differentiating Equa~tion 60 with respect to
MAPYT2 +T A2 A T A T
+ 2(Ti ~Y
"T
AT
sin (y +e)
AZT
2 2 .
+ 8 EnI + A+~AT ~Z)Cos (y + 8
E; y'A-~sin (y + 8) +- ZT Cos (T + 8) j
+ . 5 T 2 +C 2 (J
~ T+ + B B + (YT + 413ZT) Cos (y + 0
. 2 2 +2*
+H[Y *ZD+ +V CTAcs y+6
B B5
And then with respect to time
(6)d r 2 + z2~+
[Mtily 2+ 2 2(n + ~Z)cr Y+
(6) dt TI A T A A AT AY o y+e
S [YD 2 + ZD2 ] + [A + IB +
ID
+'M [IIA +
CA +
(n AYT +
^ZT) cos (y
+ ZnAzAT -
+8)
(Y
sYTsin + 0)]
I
+ MB [B 2.+ 4B2 + (nBYT + ;BZT) COS (Y + %)
Y )
+ (I"Z-; - sin ( + 8)]
1B ]
[(+
+
A A (AYT Si + 0) + Z Coo (y )
46 .
Similarily, differentiating Equation 60 first vith respect to 0
A7 A A
+
- qA j
'A
A ; A%2
+ M riB 2+ B2
+ B + 2
[B2 + (nBY T + BZT) cos (y + O + (nBZT By ) sin (y + 8)]j
+ +'A
+ I I [ + )
oL d r4, 1
J 2
(62)a IT .i 'r[ +Z A T AT c (
+ [IA + IB]]
2 2 + 2 +2+1A + n"A
+ 8 M n A+ A B
I nB+ B + "MA Y
47
U
[:
II
+ A 10 1
VG - Potenttal energy stored in the effective spring at the front support point.
Thac is,
48L
A ~- !
:;i
iL i
e0
VW CA A + CS WB + CD W D
I
or
(66)
LIi
VW [YT sin + ZT cooB + nA sin (Y +) + A (Y + + A
+ [YT Cos .
sin + n B cos
T + c sin (Y + e + 0)) WB
+ TB Cos (Y + + 0) %sin (y + 8+ )] WB
49
4
By defining static equilibrium as the refertnce position, the static
load in the ground spring is given by
B W + b W BD W
static A A B B D D
L
A -- L
The minus sign decreases the spring deflection as 0 increases when
L > 0 (i.e. - when the pivot point is at the rear), and i-creeres
the sprtng deflection as * increases when L < 0 "i.e. - when the pivot
point is in front of the spring). This allows cie total spring
deflection to be written as:
G Wt G~~
static
1 A 1 2
VG m G
(KA) KG G
+ A 2
V KG [A
22* CP
I + static
2 G KG
2 KGLJJ
-KG 2Pstatic 12
2 Kc L
50 1
LJ ---
Letting
-static
KL
static L
or
3 W + B W + B W I
$
static ~A A B B
Y'GL2 D D
(69) VG 1 - staLi
2
(70) " K - static
and
(71) --. =0
51
Z1 - ZT +n sLin (y + 0) + c Cs
BEL '
EL E2
2 EL
LXI
LEL Y + 0
a (Y
(Y yYL 2
2 +z-(ZT _ zE) )2++ n2+
+
2
2 ;EL 2
y+O FL T aL a ELI
.'.
52
"
...
.. . - . . f. . !.
: t_: .? . ' t
:i :: "' /
,- ... .,4 ,, l ,'¢ --. ,, , ,,
By defining the constants:
2 + +
(72) EL (YT - YEL 1- (ZlT - EL)2 EL EL
1 /2
[ELI + 2EL2 cos (y + 8) - 2EL3 sin (Y +
(75) LEL Y + -
53
A spring equilibrator (compression type) is shown
schematically in Figure 13.
FIGURE 13 L!
Spring Equilibrator - Coiupression Type
Let
54[
Then, the total force for two parallel equilibrators
is given by
With
energythe reference
stored in theposition definedisbygiven
equilibrator 8 -by
0, the potential k
LEL +
KEQ 2
2 [(PEq + KEQLEL,) (LEL y +e)
LEL Y
2 (LELy + - LEL y
Then
i - K Q LELEL
LEL + 0
KEQ y+ 8 a
ae ]
I
2!ES 2[LP EQ - K~ (LEL+ - LEL 0 ) 3E
(77)
(78) L=q 0
LEL y+e
FIGURE 14
Scl~eaatic Diagram of Elevation Strut
56
Treating an increase in length as a positive deflection of the
spring assembly and defining a tension load as positive enables
one to depict graphically the force in one strut (Figure 15).
FEL I-
2 K
H Pr
2 s I
s
FEL- "O
2 I
-=======---Deflection
II
FIGURE 15
Load Deflection Curve for one Elevation Strut 1
Note that the value of FEL is taken as zero at time 0. Then the
spring deflection is given by
=LL + LEL i
Ii
-. j'-- ..-
Since the load in the elevation strut must be defined in a piecewise
fashion, the potential energy will also be defined in this manner.
With the equilibrators approximately balancing the tipping parts when
e - 0, the potential energy in the elevation struts (both) may be
defined as twice the area between the load deflection curve and the
A axis from A - 0 to a The curve chosen will be dependent on
whether the spring assembly is being compressed or releasing its energy.
+
(79) VEL 1A - r]2 K + 2[A - r] Ps H
Since A is independent of 0
3VEL
(80) E -o
However,
VEL A a
S2 [A - r] K + 2p
-2 K [A - r) + Ps
and since
A- L~ 8 - LELy
The equation
-LELy - r] + P] y+
(81) -- -2 K [LELy +
is obtained
58
iI
- T .1
~I1
V m VW + VG + VEQ + VEL
- - VEL
3V 3VV V 5V 3
- G EQ +-
3e Be 36 386 e
+6
+ 2 jPEQ-KF (LELy + "LEL )] 3
3LVL
+ 2[K(LEL -LEL -r)+P] * +
59'
1
[ELI 2EL2 cos (v+ ) 2EL3 sin (y + 0)] /2
60
By the resolution of the breech force [B(Ct)] into a force along
the n - axis and a couple around the point T and by including the
dissipation functions representing dmping in the ground spring
and in the elevation trunnion, the generalized forces may be
.aken as the torques around the appropriate pivot points (0 and T).
Then,
or
=
(85) Fe B(t) q C0
Then, since the Lagrange Equation for the pitch motion (defined by 0) is
d ITT + .- F
-" I -
and for the relative rotation of masses MB and MD (defined by f) is
d [Ia T 3V F
61
in obtaining the equaLion for linear motion of the mas center, Newton's
Equations may be used. From Figure 11
BA. + CA k
A sinl
B~ (Y + 8 + 4) + CA Cos y +8-i+ 0)
+ n
2 $ nA 6 nA +
-"A - A [hA +T cos (Y + ) + ZT sin (Y + 6)]
and
'- [~T c*
.. rv.Y ( +
oo (Y e)
+)Z + zTri
sin (Y'+o
+8) + nAAl
+' CnA]
-2$8 A+2 8
AA+ 2 A
62
.ii
. .:. .. . . ,_, -,, , -' . ,. ,. . - - ,. . . *1..
,.. . . - ,*. .
I
+
(87) XF MA Q, Sl S2 - WA cos (Y + +)
L,
+ R(t) [
A
(%-
-.
2'
) A +
A2 "2) cG)G]
Then, the equation for linear motion of the recoiling parts is (by
substituting for Q in Equation 86 1
Equation 91 (See Page 64)
Equations 89, 90, and 91 are the governing motion equations which define
the variables
* * a and nA
A
63
*.
i;
EQUATION 89 E
MA2 2 7AT A
2
+Zf
I'+~ IA+ 41+ (AYT 4 CA Z)COS V4-A T) sin(T4.+e8))( MAI ~
A+E rC~
{
+ M.[
+ [01I'A
12iL,' C.2
18~
2(jY,+4CZ,XCos(74-8J-#+ 2(1 ZI-
10LMAA
8 Y) Slr1
+
+ 0)] +C2+%Y
2'± 4I
MAI'AA"AC nT 7+9)(AZ7TCS
[n'2
8Y + C2Z )sn(Y, Cer -y
4-+{s+[,) -
)
YT 0(CA
ns-
9)s +82K[
COTSn 7'40
-E~o }) MAI
-[L
I
os'+C
+M1)1
-2~A, fA Y: i r(-+*) 4- Zr Cs 7+ 0) +
+~[O ++~A
2 MA 'Ar+CZIS"- + O('?ACA ]
O( + 0)] *2 KL [L r2A~o
20 YT
MA A
R-) CA
EQU4TION 90
~Z~r
ZAe CAYT)in(7+e) A + C+ C ,Yr+VAl)COS(7+ ) + (IAZ CAY)SIP(? 9)
7+0)+ 2
(q 8 2 + t12
+~ Xfasi(+6] ~?oy+ )f ~ 8 ~- STSY 7
8
+xi
+1
+81
+Y (zTCr)sn(7 4 9)] + {4- [A~~ + M
4[8 n + +[A#]
++ +fA[iMAaA
{
-CASi"
2 r f 2L2cosj,,+ 0 -2EL_
+ [LI os(+ )- E~snj
0)-C~~sin(T-#-L E0). [EL1Esie +ELcr+O
+ - Stin1 7 +'
T CoT)C )
'Jy jEL2smT(7+,) 4- EL3cosq+ e)
M (~YT+
MI1 ["A CAZT)S i "(7 + 0) -(11A ZT CATIC S(?+8
EQUATiON 91
0
[gll+U1l+Fsni)-B Ai(r+G MA +0 iMA'A + eAA + Y 1COS(7+ 9)1+ Z7 Sinqy~9
c:4
I
Note that the breech force [B(t)] is the only forcing function
included if Equations 89 and 90. 'owever, the following functions
are included in Equation 91:
(92" x nAo- nA
(93) -
(9-,)) - *;A
65
From Equation 87
or, S 2 = WA cos (y + e + ) + MA
U
SI + S 2 - W A cos (Y +
+
+ 0)
aA-2 A +2 A +2hA
From Equation 88
]
S () - S ("-a) B(t) [B - F
1 2 1 A ~B(A 5)
+ R(t) [A-
+ [f +
CL f2 (;A"
A-
'2)2 +
(CA "G )
FG
66
Let
+ R(t) k -r ] - +
A 2 A
Then,
I + S -KYZ
~1 2
Solving simultaneously,
XYZ () - XYZl - RAIL
(98) S2 =
b
(99) S1 a XYZ - S
67
The force in the two equilibrators has been defined 8 (page 55)
and the force in the two elevation struts is given by (See Figure 15)
FEL + 21(A - r) K+ P ]
or
/ 2
(75) LELy + 0 - [El + 2EL2 cos (Y + 8) - 2EL3 sin (Y + E))L
in the n direction.
0
R -B sin (Y+ + ) + Cos (Y+ + 0)
B B
68
Substituting for BB and CB and simplifying allows writing these components
in the following form
-
"
- YT ns (Y + 0) + ZT cos (Y + 8)
- B - n
•.2 .2
- [B +
YT cos (Y + e) + ZT sin (Y+O)] - b B
and
Then, from the free body diagram of the cradle [Figure 12(b)]
=
FF 1K =TH FEQ+scos T - R(t) + BF
+ [fI + f2 + F -
and
F R Tv + FEQ+eslT -S-S 2
69
: 1
- !I
tan- y
Zl- + (y + 8
Then, I-..
(102) T - " B - YT sin (y + 8) + ZT cos (y + a)] + O€ + 2;6n B
.2 *2 rig]
And
- ¥ sin (y + 0) + Z cos \ 1 -
70
The free body diagram of the mass M9 is shown in Figure 16. Since this
mass is also assumed to simply rotate aroumd point 0, the constants,
(104)
I
TV tan- I
zj
Y-D
,o;
and
2
(105) R3 Y
[y + zD21/
Static
Then, since
and
Pstatic
(106) KG L
Ptatic
1 71
Ii
-- i i V~~i .... j W -ii -
.....
r -~ Ii
I:
0
44
00
~ 4
morn
~ou
2.4
- U
.i.O
01
- I I.
i~.
[ p.
-
b.
a
I
I
I>
72
I,
77i
Therefore,
G - L s tatic)
by
The reaction at the effective spring is given
Vspring - KG G [5
tB " D YD s in U~ + 0)
T cos (Y + 8 + 0) - Tv sin (Y - 8 + )
F L cos (T Y
-F .+ Cos (T-y- 8- ) -1
y + !
73
...... ...
'P MD t& o ~~0
i
V sp ring FEL Ty+ 6 sin (r-y- e- )
-w + V pivot
D
Therefore,
+ +
(108) vpivo t ,TH sin (Y + )+ TvCos (y + )
- FEL+Y+ 6 sin (T - y - 6 -
(109) H - TN cos (y + e + *) - Ts in (y + 6 + )
74
EVALUATION OF STATIC REACTIONS
ZF - 0 -S 1 + S 2 - WA cos (Y)
S I + S2 - WA Cos (Y)
"Hass Center =
L - 2 (' 1 - 2
Solving simultaneously
(ill) Si WA Cos Y - S2
75
To define the static luad In the clevatin meehanism (which was consid'red
approximately zero whcn developing Models A aud B), an examination of
Figures 11 and 12 shows that, for the tipping parts (MA and B)
ZM 0= W (n cos Y- sin Y)
T A~ AO A
+ W (n cos Y - B sin Y)
B B B
Z + I sin Y + cos Y - Z
(113) T - tan - T EL EL EL + y
YEL - YT- nEL Cos Y + EL sin Y
76
In a vximilar manner, the following equations allow for determination
of the trunniou reactions.
- W sin y -W sin y
A B
- GA cos y- WB cos Y
Therefore,
]
-WD [YD
EF; 0= H
77
I
Therefore, under static conditions
(117) H 0-
WA[YT + nAO cos y -A sin yl + WB[YT + B cos - CB sin y]
+ W Y
D D
(118) V n " L
When the recoiling mass returns to the latch position at the end
oi any firing cycle or when it reaches the limit of forward
travel in the case of a misfire, additional relative Motion between
the meases MA and ME is prevented by mechanical stops. Consequently,
the preceding system of equations must be modified by eliminating
Equation 91 completely and holding the value of n (and x)
constant from this time on. At the same tize, in he remaining
equations,
A 0
'I"
A 0
sgn
78
FIT
REFERENCES
79
- a4,
ri
SYMBOL TABLE
in& 'es
pounds
seconds
radians
80
aleak Flow area due to clearance between Fig. 5
De Value of D(t) at t - t R 8
Fig.
Dr Constant level for D(t) during recoil Fig. 8
81
EQ Subscript relating variable to equilibrator
82
I Impulse - Area under B(t)
LELa a!
Length of elevating struts at elevation - 0*
LEL
Y
Length of elevating struts at initial
elevation (y) i
LELY+ Length of elevating struts at rotation 0 J
MR Mass of recoiling parts without floating Fig. 2
piston
ii
83
N Number of recoil cylinders
chamber
84
'1
T Elevation trunnion - origin of T - Fig. 11
coordinate system
T Kinetic energy
T Trunnion reaction
Vi
- f to T - 4 axis Fig. 12
t Time variable
t Time of recoil
r
tf Time of firing
V Potential energy
V Vertical reaction at pivot point 0 Fin.
C 11l
pivot
Vspring Vertical reaction at ground spring Fig. 11
85
f Firing velocity
x Recoil displacement at which terminal
coumerrecoil velocity is reached
a C coordinate of fI Fig. 12
coordinate of f2 Fig. 12
Angle locating mass center of non-elevating Pig. 16
parts
86
I
I
AAA.
riA Defines position, velocity, and acceleration Fig. 12
of mass M A with respect to tho elevation
trunnion
Coefficient of friction
87
k'I.
APPENDIX
881
I
2
BF - ABPB - 2 K 3 (x - x,) + 2S'
where
89
.4
and, since
2 2
C"
Al 2g 2
a
Therefore,
2 2
2 (x
BF" AB 299 2 -2 2K3 xB)+ 2S'
BF--3 aa 2
2
BF sgnc ~2 K (x xB) + S'
90
i-
00
fA,
-i .
L16.
W6~k~
XI
a 91
DISTRIBUTION
Comander
Headquarters, U. S. Army Materiel Command
Washington, D. C. 20315
2 ATTN: AMCRD-WC
Commander
Headquarters, U. S. Army Combat Developments Comand
Fort Belvoir, Virginia 22060
1 ATTN: CDCMR-W
1 ATTN: Technical Libra-
1 Commander
U. S. Army Combat Developments Comand
Combat Systems Group
Fort Leavenworth, Kansas 66027
1 Comander
U. S. Army Combat Developments Coumand
Institute of Combined Arms & Support
Fort Leavenworth, Kansas 66027
Comander
Headquaters, U. S. Army Weapons Command
Rock Island, Illinois 61201
1 ATTN: AMSW-E-RD
1 ATTN: AMSWE-RDP
1 ATTN: AMSWE-LCA
1 ATTN: AMSWE-LCD
1 ATTN: AMSWE-LMC
1 ATTN: AHCPM-CAWS
Commander
Headquarters, U. S. Army Test & Evaluation Command
Aberdeen Proving Ground, Maryland 21005
1 ATTN: AY.STE-FA
_ ~. -.--.
seclit Ctasuafcatiee
DOCUMENT CONTROL DATA.- 2S D
4swwUt It"0., It.* #1 .b.,,o
0,11C010 #DeglMar i....M eft sou~~u
vftop -*w sow~w
o.bm 0040 lee
to..S
tof
a None
IvIIi
NoneU. S. AryWeapons Command
Mathematical models used in design of the XK204 Howitzer (AD and ED Prototypes)
are described in this report. The physical basis for the mathematical
representation is presented along with the derivation of the required equations.
While these models have been generalized to allow their use in other weapon
design situations, some modification will be necessary to include features not
specifically consideredl. Systems of equations which will provide for the
definition of required control functions as well an the prediction of recoil
mechanism functioning and weapon motions are sumarized.
DDNwVi.1473 I UNCLASSIFIEDSNN
Id. LINK " LINKB LINK C
Kay WORO-
1OLE ROL
01. WT 10 LE FT
Artillery
D04 Howitzer Design
Mathematical Model
Soft Recoil Mech&nasm Design
UNICLASIFIED
I
t. N
I
- - a -
- . ,A:......
,i &A: