ooze7683/w s3.00+ .a0 Pqaluon Press Lid.
sdidr Slruefvrcr Vol. 20. No. WIO.PD. 897409. 1w Rioled inlhc U.S.A.
hf. 1.
DYNAMIC NECKING OF AN ELASTICPLASTIC CYLINDER UNDER UNIAXIAL STRETCHING 2. Q. HUANGt Departmentof Naval Architecture,Huazhong University of Science and Technology, Wuhan, Peoples Republic of China
and L. H. N. LaE Departmentof Aerospace and Mahlical E&n* University of Notre Dame, Notre tic,
IN 46556, U.S.A.
(Received7 June 1983;in m&d form I Nocrmbn 1983) AbrtnctThe axieymmctricnaking processof an cla&plastic cylinderunderuniaxialstrctchiig is analyzed by treatingthe problemas an initialtitigcnvduc problem.One end of the cylinder ir subjectedto a prescribedconstant axial velocity relativeto the other end and both ends arc shear free. The cigcnmodca of deviated motion departing from a uniform stretch beyond that of quasistatic bihcation by Hutch&m and Miles arc determined. Of all uniformly excited eigcnmdes of motion, only certainmodes will grow.The rateof nakhg depcadson the geometry, material propertiesand sketching velocity of the cylinder. 1. INTRODUCTION
The phenomenon of localized necking in an elasticplastic body under tensile loading has been observed and studied by many investigators [ 11.Most of the necking problems analyxed are based on Hill’s quasistatic theory of bifurcation and uniqueness[2]. In a quasistatic sense, the problem of the onset of necking in an elasticplastic cylinder under uniaxial tension has been analyxed[35]. Hutchinson and Miles[fl have shown that the state of uniform uniaxial tension is unique prior to the maximum support load of the specimen and that the true stress at bifurcation is greater than that at the maximum load by an amount which depends on the geometry and material properties. Dynamically, the necking process of the cylinder is influenced by an additional factor of stretching rate. There have been relatively few attempts to solve the dynamic problem. An approach based on a onedimensional formulation[q has been suggested for determining the dynamic axisymmetric necking process. In this paper, the axisymmetric necking process of an incompressible elasticplastic cylinder under uniaxial stretching is further analyzed by a twodimensional formulation. The ends of the cylinder are stretched in such a way that the ends remain free of tangential tractions and the lateral surface traction free. One end of the cylinder is subjected to a prescribed constant axial velocity relative to the other end. The axial velocity is relatively small in comparing with the elastic axial wave speed of the cylinder. For continuity, the state of quasistatic bifurcation as determined by Hutchinson and Miles[Sl is employed as a reference state. The dynamic necking problem is solved by employing a finite strain, dynamic, quasibifurcation theory[7,8] which is briefly described in the paper. Numerical results for a number of cylinders are presented and discussed. 2. QUASIBIFURCATION A certain motion of a system of n degrees of freedom is considered stable if, after a sufkiently small disturbance, the system remains to follow the undisturbed motion. In other words, the undisturbed motion is stable if the deviated motion, C,(r), I = 1,. . . , n, the tVisiting at the Departmentof Aerospace and Mahnical Engineering,University of Notre Dame, Notre Dame, Indiana, U.S.A. 897
898
Z. Q. HUNG and L. H. N. LEE
difference between the disturbed and undisturbed motions in n generalized coordinates, remains small. For any deviated motion, stable or not, there is the following identity[8]:
in which dots indicate differentiation with respect to time t and y, and y, are the initial values of the deviated motion. When the generalized forces, P,, of the system are configurationdependent and when the deviated motion is relatively small, the deviated motion is governed by the following variational equations[7]: &A,C,=AP,,
rl
,...,
n;
s=l,...,
n
(2)
where the coefficientsA,, are known functions of time depending on the undisturbed motion only. In eqn (I), without loss of generality, the initial condition may be taken as y, = 0 at a chosen time I = 0. Consequently, eqn (1) shows clearly that the undisturbed motion is unstable or the deviated motion grows monotonously when the quadratic term
is positive in a time period for any possible deviated motion. Furthermore, a variation of Q yields the following relationship: SQ = 24&5C,= 2t;Xr
(4)
which may be employed, by appropriate coordinate transformation, to obtain equations of motion in other generalized coordinates. The quadratic form Q at a given time can always be reduced to a linear combination of squares[9] such as: Q
amrIm
’ 9 m= l,...,n
(3
where a,,,(t) are the timedependent eigenvalues, PI,,,are related to C, by an orthogonal transformation &=lm’l)m,
r=
l,...,v
(6)
where I,,,,are the directional numbers of n mutually orthogonal vectors, q,,,,in the c, space known as eigenvectors or as eigenmode of motion. Such transformation is called a transformation to principal axes. In general, the directional numbers 1, may be time dependent[8]. For cases where I,,,,are always constants (or nearly constant), a variation of Q leads to the uncoupled ordinary differential equations &=a,q,,
m=l,...,n
(7)
where the underscores are placed under the indices to suspend the summation convention, If c, and q, are normalized, Q assumes the greatest value equal to the largest eigenvalue, a,, corresponding to the configuration qk in the c, space by a theorem of Weierstrasse[9]. Therefore, ak and & may also be determined by the extremum condition. If or,> 0 and a, > a, for r # k, then the eigenmode of motion, &, grows at the highest rate, provided that the initial disturbance has such a component. Such phenomenon is called a quasibifurcation phenomenon. The above concept and approach are applied to the solution of the present problem. Of an elasticplastic solid of a volume, V, let the undisturbed motion be described by
I)ynmic
necking of im &sicplastic
cylinderunderuniaxiulstretching
899
the displacement VAX,, 1) in the X, Cartesian coordinate system. The deviated (or additional) displacement is denoted by udXM, I). The corresponding Lagrangian strains are e, and efi and the corresponding ViolaKirchhoff stresses are S, and s,, respectively. When rc,,is relatively small, the additional strain, Ed, is given by
and the deviated motion is governed by the following equation of motion:
where p is the initial mass density, The deviated motion satisfies the following boundary conditions: ux(X,, 1) = 0 on that part of the boundary with prescribed kinematic conditions and
on that part of the boundary with prescribed surface traction TM, where NK is the outward unit normal to the surface. The quadratic term for the deviated motion of the elasticplastic continuum is given by
QW
=
IV
APMuMd V =
piMuM d V. sV
(11)
For a body which has no change in surface traction and no interior discontinuities of the variables, the quadratic term may be simplified as
Q(Q)= 
s
v (&LUM,L + SKL& dV.
(12)
The motion and its stability of such solid may be determined by analyzing the functional
QW By employing an appropriate constitutive relation, the function Q can be expressed in terms of the deviated displacement field uM.When the deviated motion is relatively small, the disturbed and undisturbed motions have nearly identical histories such that s, appear as small deviations from the stress path of S, in either loading or unloading condition. Thus, the stress deviations may be expressed in terms of the strain deviations by a general expression such as (13)
Hew CuN is the tensor of instantaneous elastic moduli, Bm is the unit tensor normal to the elastic domain in the strainincrement space, and g is a positive scalar which determines the strain hardening of the material and depends upon the strain history of the undisturbed motion. It is noted that the quasistatic bifurcation is characterixed by the occurrence of a nonzero solution, u,, to the variational equation SQ = 0 as shown by Hill [2]. 3. AXIAL
MOTION
Consider a solid circular cylinder having homogeneous elasticplastic properties. One end of the cylinder is subjected to a constant axial velocity u relative to the other end. It is assumed that the axial velocity is low enough such that the axial stress waves have
900
L.H.N.LEE
Z. Q. HUANGand
reflected many times throughout the length of the rod and that its dynamic state is the nearly uniform state of uniaxial tension at a time prior to the beginning of the necking process. It is also assumed that the ends are stretched in such a way that the ends remain free of tangential traction and the lateral surface traction free. Consider that there exist initial axisymmetric imperfections in the cylinder which are infinitesimal in magnitudes such that they produce negligible disturbances to the nearly uniform state of uniaxial tension prior to the bifurcation process. The imperfections may lead to an initial state of the deviated motion. The onset of quasistatic bifurcation in the cylinder, which has been determined by Hutchinson and Miles[Sj, occurs at a true axial stress cc. This state of stress is employed as a reference state which occurs at a time t = 0 when the length and radius of the cylinder are L and R, resp. If initial imperfections or disturbances are absent, the dynamic state of the cylinder at r > 0 is the nearly uniform state of axial tension with a true axial stress u and a length of (L + or). This state is considered as the undisturbed state. Assuming the material is incompressible, the second PiolaKirchoff axial stress, S,, may be expressed in terms of the true axial stress as (14)
c=
J E P
and
r=Cf
(15)
L’
Here, E is the Young’s modulus of elasticity and c the elastic wave speed. The true axial stress may be expressed in terms of s by employing a RambergOsgood relation[lO] between the true stress u and natural strain e such as ea =eY
3a* QY
+?
O(
<
(16) >
where ey and gy = Ee, are the effective yield strain and yield stress and n’is the hardening exponent. Let E, be the tangent modulus for an uniaxial increment of true stress according to E 1=x.da
(17)
The derivative of the tangent modulus at the reference state is found to be
where E: is the value of E, at Q,. Retaining only the first term of the Taylor’s series expansion of E, in the neighborhood of cc, it is found that (u  UC).
E, = E,’  (ii  1) 2
Following the definitions of E, and the additional natural strain e=
VT
vr
’
( >
c l+;
(20)
Dynamic neckingof an elasticplqtic cylinderunderuniaxialstretching
901
the integration of (19) yields
=
==c+(+
$&)[1 (I
+~r)+‘Ixl“)l
for ‘520
(21)
where
4. DEVIATED
MOTION
When the undisturbed motion is stable, all the disturbed motions will remain in the neighborhood of the undisturbed motion. When unstable, only certain deviated motions, however small in magnitude and whatever way excited initially, will grow and depart from the neighborhood. The axisymmetric, deviated motions, of the cylinder may be described in terms of the deviated displacements (u,, u,) in a cylindrical coordinate system (r, 8, z) with 0 Q r d R and 0 4 z < L associated with the refcrencc state. When the deviated motions are‘ relatively small, the constitutive relations of an incompressible elasticplastic solid having isotropic elastic properties and under continued plastic loading may be written as [Sj ~~=2cr,fr,+p s,=2G~,+&+p 2g
(22)
where G is the elastic shear modulus and p = f(s, + se + 8,). Here and in the remainder of the paper, the physical components of the deviated PiolaKirchhoff stresses (s,, s&s,, s,) and deviated Lqrangian strains (6, ~,e~e,J are used. For an incompressible material, 6, + tg + d = 0 and 3G = E. From the ddinition of E, and (22), it follows that[q 12[,,;(,E,)].
(23)
g
The equations of deviated motion and associated boundary conditions may be written as
(24)
s, = 0
sn = 0
for r = R
and
24,= 0
s, = 0
for2 ==O,L.
(25)
The equations of motion may be expressed in terms of the deviated diaplaamsnt components by using the constitutive equations (22) and the following relationships between strains and displacements:
(26)
5302
Z. Q. HUANG and L. H. N. LEE
Following the approach by Hutchinson and Miles[q, a function @(r,z, I), which ensures the condition of incompressibility, is introduced with
u*= 
[email protected] az
uZ=iy.
(27)
Furthermore, the function @ can be written in a separated form
kxz 9 = R+f~(r,i) sinL,
k = 1,2,3, . . ,
(28)
such that the boundary conditions on z = 0, L are satisfied. Using eqns (22X28), the tractionfree boundary conditions (25) on the lateral surface can be expressed as
(29)
where
(30)
The operator in (29) is defined by
(31)
It may be shown by further manipulation that eqns (29) may be reduced to a single constraint equation.
a4
(E,o+G)Ix+(E,aG)&=Oonc=l.
(32)
The equations of deviated motion (24) may be further reduced to a single hyperbolic equation in terms of the function 4 subject to the constraint by eqn (32). Instead of solving this equation directly, the following variational approach is employed.
5. EIGENVECTORS The quadratic functional Q, by (12), for this case of axisymmetric deformation, may be specialiml to the form
Dynamic
necking of an elasticplastic
cylinder
under uniaxial stretching
903
Employing qns (22), (26)(28) and integrating (33) with respect to z, we find
I#
+y2
1x1
( )I[ ( >
yz(p+
?+2
[
c
A,)
1+;7
I
@1)(1l,)
al
I ! 0,
2
+l
* *+Vr (
2 CdC*
c )
(34)
For small values of y, the solution of 4 may be represented by the following power series[5]
tw.7)=
i b”(7)
(n
_
.
nl
The dimensionless amplitudes (32), or
k’ 0 l)fn, 1
(1)l
p‘*
(35)
.
are subjected to the traction free boundary condition
b,(r)
7nI (l)+’ f,
(
b,(r)
;Y
>
(n  l)!n !
[(E,  a)n + G(n 
l)] = 0.
The eigenvalues and eigenvectors of the functional Q may be determined by the extremum condition. Substituting (35) into (34) and integrating with respect to [, the functional Q is reduced to a quadratic form in terms of b,(t). By introducing Lagrange’s multiplier bN+), the extremum condition and the condition of constraint may be combined as
ME,  0) + G(x  111 = 0.
For a chosen value of
7,
(37)
the variation of (37) with respect to each parameter b,(r) yields N+I
c
r1
m=l,...,
C&(7)=0,
N+l.
(38)
The eigenvalues, aP, may be determined by the characteristic equation
Icm apaM= 0.
(39)
The corresponding eigenvectors may be obtained by solving the simultaneous equations (Cm  a&J&=0,
ml
,...,
N+l,p=l,...,N+l.
(40)
The eigenvectors may be expressed by the orthogonal transformation bm(7 ) = Lvm(7
).
(41)
2. Q. HUNG and L. H. N. Lee
904
For a case of deviated motions from a uniform unperturbed axial motion, the directional numbers /,,,,,of the eigenvectors q, are practically constants. Substitution of (41) into (35) yields
(42) Substituting (42) into (34) and observing the orthogonal properties of l,,,,,the functional Q may be expressed in terms of the eigenmodes of motion as Q =  c nR2LEa,qi
(43)
P
where
ap=
(lY’
N+l
I
5( ,i 24
Y2
I
0
0
c
i
c [i n1
(n

l)!n!
*I
( (! N+,(l)~
+ @;,
[
f
c
n1
where S=g,
E1
2
(2n 
1
l)CV,
2
0 (n 
*=g,
EL
p3(yzp+4nz4n)fp
l)!n !
1 2
; znl
C=.=c E
6. NECKING The
)
czn21p )I
(n _ l)!nl
N+,(l)+’
+A,
2nr2”21pm
k’
0 ;
N+,(l)rl + 2hY2
(n  I)!n!
II1
2
f “I
PROCESS
functional Q may also be repre$ented by the integral at the right side of (1I), or
Dynamic
necking of an elasticplastic
cylinder
905
under uniaxial stretching
By substituting qns (27), (28), (35) and (41) into (45) and integrating, it is found that Q = npR’L‘c’x
j? %q,, p ‘dr’
where
(47) A variation of the functional Q with respect to the pth eigenmode of motion leads to the following uncoupled equations of motion $/I,(r)lr, 
_
p=l,...,N+l
where /I,(r) = 2
L ‘ap 0 R %’
(49)
Equations (48) describe an axisymmetric necking process accompanying with the axial stretching of a cylinder. The development of the necking process depends on the function A,(r). If A, < 0, the p th mode of motion is oscillatory. If /1, > 0 in a time interval, the p th mode of motion grows monotonously in amplitude in the time interval. When the value of A,(r) of a certain mode of motion becomes positive at the earliest time (T = 0), it marks the beginning of the necking process. For that particular mode of deformation, the condition A,(O)= 0 or s(O) = 0, which corresponds to the quasistatic bifurcation criterion, leads to the determination of the reference state having a true axial stress of 0,. For an elasti*plastic solid having isotropic elastic properties, the value of uC as given by Hutchinson and hGles[S]may be reduced to that by the following expression
(50) It has been shown by Hutchinson and Miles[q that the quasistatic bifurcation cannot take place prior to the maximum tensile load. Denoting the tangent modulus at the maximum load point by E,m,the maximum load condition yields he corresponding true stress urn= E,m. Furthermore, uC> urn,Et < E,“’and E,’ < u, for a material having properties described by qns (16X 19). In a dynamic process of stretching the cylinder to a state having a true stress u greater than u. as shown by qn (21), the corresponding E, is such that the value of (u  E,) increases as the stretching progresses beyond a,. The value of A,,of a certain eigenmode of motion increases positively as (a  E,) increases. Therefore, thrs is a principal mechanism pushing the necking process. In the necking process, all modes of deviated motion may be excited by the presence of unpreventable, Wnitesimal geometric or material imperfections or initial disturbances. However, only certain modes of motion having algebraically large values of AP,which depend on the geometry, material properties and rate of stretching, may grow in amplitude. To determine quantitatively the necking process of an elasticplastic cylinder, a numerical procedure and a computer program in double precision based on the foregoing approach have been developed. Eigenvaluw and eigenmodes of motion are
906
Z. Q. HUANGand L. H. N. LLS
determined by a direct method and the integration of eqn (Ss) is accomplished by employing the sixth order RungeKutta method. At time increment ATof less than l/50 of the shortest period of any oscillatory motion present in the deviated motion is employed. The power series by eqn (35) converges rapidly. Therefore, a limited number of terms may be sufficient to represent a solution. Accordingly, the number of distinct eigenmodes of motion will be dictated by the number of terms employed in the series representation. However, the eigenmodes of motion are determined by the extremum conditions of the functional Q and the respective values of A,,. Thus, the results of a limited series representation may cover sparsely a whole spectrum of deviated motion. Here, only first six terms of the series are employed in the computations. A number of cases with various combinations of other parameters have been considered and are described as follows. 7. NUMERICAL
RESULTS
Numerical results have been obtained by the foregoing procedure for a number of cylinders. Each of the cylinders has initial imperfections leading to an initial state of the deviated motion respresented by a displacement profile which has the components of all eigenmodes of deformation of an equal amplitude with q,(O)f 0.00001 and r&,(O) = 0, p= 1,. . . ,6. The typical results are shown in Figs. l6. All the cylinders referred to in Figs. 16 have the following material properties: q = 5, Et/E = 0.1 and various values of e, such that, for each cylinder, the quasistatic bifurcation stress 6, agrees with respective E,’ and y by eQn (50). Figure 1 shows the histories of the amplitudes of a number of eigenmodes of motion of Case 1 having the parameters: y = 0.1 and u/c = 0.001. The amplitudes of the 6rst and sixth eigenmodes of motion grow with time as soon as the axial true stress exceeds the quasistatic bifurcation stress eCor r > 0. The first and sixth eigenmodes of motion predominate in the necking process while all other eigenmodes of motion are oscillatory and remain to be of small amplitudes in an order of that excited initially. The fifth eigenmode of motion is not shown as it has a much shorter period of oscillation than that of q3 or g4_ Figure 2 shows the amplitudes of a number of eigenmodes of motion of Case 2 having the parameters: y = 0.2 and u/c = 0.001. Similar to that of Case 1, the first and sixth eigenmode of motion of this case also predominate, except that they grow at somewhat slower rates. It is noted that the product of W/Cgives the value of the average axial strain in addition to that of the reference state. Thus, Cylinder 2 developed a neck of a certain amplitude at a slightly larger average axial strain than that of Cylinder 1. Figure 3 shows
I 0
2
4
6
f,
IO
DIMENSIONLESS
12 TIME
Fig. 1. Eigenmodcsof motion of Case 1.
14
16
Dynamic necking of an elasticplastic cylinder under uniaxial stretching
I
0
2
4 T,
6 8 i0 12 DlMENSlDNtESS TIME
14
907
I6
Fig. 2. ET&emnodes of motion of Case 2.
r!
G 0;
1.20
z z E 1 2 p:
1.10
z: p! a’ 1
1.00 L
Id” 0.93
0
0.2
0.4
0.6
0.6
1.0
c , DIMENSIONLESS RADIAL POSITION
Fig. 3. Normalized radial strain profiles of Case 2.
the normalized radial strain profiles of the sii eigenmodes of motion of Case 2, where Z,= c,(~)/c,(O).For each of the cases considered, the eigenvahres and corresponding eigenmodes of motion are determined at each time increment, An eigemnodc shape (or &,,)does change with time but at an i~i~nt rate. The ins~~n~~ eigenvahtes are employed in integrating equation (48). The normaliz& radial strain profiles of the eigemnodes of motion of other cases have patterns and relative orders similar to this case. ,Figure 4 shows the histories of the amplitudes of a number of eigenmodes of motion of Case 3 having the parameters: y = 0.4 and u/c = 0.001. The second eigenmode of motion grows in amplitude continuously after r > 0. However, the first and sixth eigenmodes of motion grow at a later time. Cylinder 3 develops a neck of a certain amplitude at an average axial strain achy huger than that of Cylinder 2. Figure 5 shows the histories of the ampiitudes of a number of eigenmodes of motion of Case 4 having the parametes: y ~0.1 and+ ~O.~l.~cs~~~~~of~s~ is one tenth of that of Case t . The fmt and sixth eigenmodes of motion aiso predominate and grow in amplitude soon after r > 0 at rates also much slower than that of Case 1. The
Z. Q. HUANGand L. H. N. LFE
3
0
2
4 6 8 T, DIMENSIONLESS
Fig. 4. &gcnmodes
IQ
I2
I4
16
TIME
of motion of Case 3.
7
I ’ 0
I 5 r,
I I IO I5 DIMENSIONLESS
I 20 TIME
I 25
I 30
Fig. 5. Eigcmnodesof motion of Case 4.
effects of stretching velocity on the rate of necking are illustrated by Fig. 6. Which shows a neck of a certain amplitude of a cylinder is developed at a shorter time by applying a greater stetching velocity. However, a neck of a certain amplitude such as that having (0.00603 vs 0.0002), G(l) = 6 x lo’ occurs at a (e vs u/c) ratio of (0.00394 vs O.OOOl), (0.0173vs 0.001) or (0.078vs 0.01) for each of the four cases shown, respectitiy. It is noted that the application of the current approach to the case of relatively high stretching velocity of o/c = 0.01 may be questionable. It is shown here to indicate a trend. The trend shows that a necking failure is developed in a cylinder at the quasistatic bifurcation true stress a, and the corresponding average axial strain when the stretching velocity approaches to xcro. A number of other cases having material properties: q  5 and EJE = 0.01 have also been considered. Their necking patterns and processes are similar to the series presented here except they have slower necking rates.
Dynamic necking of an elasticplastic cylinder under uniaxial stretching
0
5
IO I5 20 25 r, DIMENSIONLESSTIME
30
35
909
40
Fig. 6. Rates of necking at various stretching velocities.
8. CONCLUDING
REMARKS
The computer program leading to the numerical results illustrated in paper has certain limitations. The program has been prepared for cases where the deviated motion is relatively small, the radius to length ratio y is relatively small and the wave effects are negligible. With additional modifications, the concept and approach presented herein may still be applicable to cases beyond the limitations. The results presented herein do indicate that the problem of dynamic necking of an elasticF_plasticcylinder may be treated as an initialvalueeigenvalue problem. The eigenmodes of motion can be obtained from a functional Q. The eigenmodes of motion depend on the geometry, material properties and stretching velocity of the cylinder. The initiation of an eigenmode of motion does require the presence of that particular mode of initial imperfection or disturbance, no matter how small it is. When all eigenmodes of motion are initially excited, only certain eigenmodes of motion predominate at a time depending on the characteristics of the functional Q which describes the intrinsic interactions between the undisturbed and disturbed motions of the cylinder. It is noted that the current approach is applicable to systems subjected to generalized forces which are configurationdependent. For stability problems involving velocitydependent forces, additional considerations are required. REFERENCES 1. V. Tvcrgaard. Bifurcation and imperfectionsensitivity at necking instabilities. DCAMM Rep. 157, Technical university of Denmark (1979). 2. R Hill, A general theory of uniqueness and stability in clastioplastic solids. L Mech. Phys. Soli& 6,23&249 (1958). 3. S. Y. Cbeng, S. T. Ariaratnam and N. Y. Dubey, Axisyrnrnetric bifurcation in an eluticplutic cylinder under axial load and lateral hydrostatic pressure. Quw. Appl. Marh. 29,41S] (1971). 4. A. Needlernan, A numerical study of necking in circular cylindrical bars. J. Mech. Phys. So/i& 20, 11ll 17 (1972). 5. J. W. Hutcbinson and J. P. Miles, B&c&ion analysis of the onset of necking in an clastieplutic cylinder under uniaxipl taision. 1. Mech. Phys. Solti 22, 6171 (1974). 6. N. K. Gupta and B. Karunes, A dynamic analysis of wking pbenomcna. Inr. /. Mech. Sd. 21, 387397 (1979). 7. L. H. N. Lee, Qursibiturcrrtionin dynamics of elsstieplastic continua. L Appl. Mech. 44.413418 (1977). 8. L. H. N. Lee, Flexural waves in rods witbin an tial plastic comprcasivc wave. Waue Marion 3.243255 (1981). 9. R. Courant and D. Hilbcrt, Mrllrodr of hiarhetna&al Physics, New York (1953). 10. W. Runberg md W. R. Osgood, Description of 8n curves by three pupmetnr. NACA TN902 (1943).