CC BY
9
УДК 517.958:539.3
Exact Analytical Solution of One Problem on Planar Deformation of Nonlinear-Elastic Media
Georgiy M. Sevastyanov*
Institute of Machine Engineering and Metallurgy FEB RAS Metallurgov, 1, Komsomolsk-on-Amur, 681005
Russia
*
Received 10.05.2015, received in revised form 15.06.2015, accepted 10.07.2015
An analytical solution of a problem on planar deformation of isotropic incompressible nonlinear-elastic (rubber-like) media with a cylindrical cavity is constructed in quasi-static approximation. A contour of the cavity is a smooth symmetrical curve. The special kind of follower load provides purely radial movement of material. Mass forces are neglected. A physical model of medium is given by elastic potential, which is analogous of Mooney-Rivlin strain energy potential (with a difference in the used finite strain tensor). The obtained solution is exact since in equations connecting Cauchy stress tensor and Almansi finite strain tensor all nonlinear terms are kept (for accepted medium model the maximum is the fourth power of components of displacement gradient tensor).
Keywords: planar deformation, nonlinear elasticity, incompressibility, Mooney-Rivlin solid, finite strain, Almansi strain tensor, exact analytical solution. DOI: 10.17516/1997-1397-2015-8-3-352-355
Let incompressible nonlinear-elastic medium occupies unlimited space and has a cylindrical cutout of infinite length (Fig. 1), on the surface of which the special kind of an uneven follower load (Fig. 2) is given. We assume the contour R0 of the cutout border in undeformed state a smooth curve, symmetric about two perpendicular lines. Consider the problem of elastostatic equilibrium of such media on the assumption of plane strain.
Fig. 1. Computational domain R(0) > R0(0) Fig. 2. On the determination of boundary load
* [email protected] © Siberian Federal University. All rights reserved
0 = Я /
Ь;го(ф)
RoW
We introduce a cylindrical coordinate system (R, 6, Z) centered at the point of intersection of axes of symmetry of the cutout (an axis Z is perpendicular to the plane of the Fig. 1), counting a polar angle 6 will lead from one of the lines of symmetry. Assume that points of continuum are displaced only radially during a loading. We define a motion of continuous medium in the Euler form as R = r — ur (r,4>), 6 = Z = z, where (R, 6, Z) are coordinates of points of medium in the initial state, (r, z) - same in the deformed state.
Almansi finite srain tensor is A = [(Vu) + (Vu)T — (Vu)T (Vu)] /2, where the displacement gradient tensor in cylindrical coordinates in view of u = ur (r, 4>) has the form
/ dur 1 dur o dr r
ur
0 ^ 0
r
Vu =
\ 0 0 0
Its nonzero components are expressed through displacement as follows:
1 dur ( 8ur\ 1 ur ( ur \ 1(1 dur x 2
2 ^ dr V2 dr)' a00 2 •r A2 r) 2 V r ^ 8$ y (1)
1/1 dur\( dur x
&r& = a.0r = X - -• - 1 -
r0 - "0r =2 Ar • \1 - ~dT
Incompressibility condition is (1 — 2- A1) • (1 — 2^ A2) — 1, where A1, A2 are the principal values of the Almansi strain tensor (A3 — 0). Integration of the incompressibility equation relative to displacement (assuming ur/r < 1, dur/dr < 1) gives a view of the solution:
ur (r, 4>) — r — \Jr2 + (2)
thus the problem of determining the kinematics of media is reduced to identification of
We define physical relations by elastic potential of isotropic incompressible medium (which
is analogous to the Mooney-Rivlin potential [1]) in the form W(I1,I2) — —b (I1 — a I2), where
a, b is elastic constants, b > 0 is in pressure units, a G (0, 2) is dimensionless quantity; Ii, I2 are
linear and quadratic invariants of the Almansi tensor, I1 — tr(A), 2-12 — (tr(A))2 — tr (A2).
The components of the Cauchy stress tensor a are determined by dependencies [2] dW
aij — — P• Sij + —--(Skj — 2- akj), sum by index k, P is unknown function of additional hy-
daik
drostatic pressure, Sij is Kronecker delta.
In expanded form, in view of Ii = arr + a00, I2 = arr • a00 — ar0- a^r, dIi/daij = S.
3I2/da.ij = Iy Stj — atj and hence dW/datj = —b- (Stj — a- (Ii- Stj — atj)), we have:
2
K0
—P — b- (1 — a- a.00\ (1 — 2- arr ) + 2- a- b- a'20, aT.0 = b- (2 — a)- a.r0,
(3)
r 0 '
The equilibrium equations (without mass forces):
&00 = —P — b- (1 — a• arr)• (1 — 2- a.00) + 2- a• b- aI
Barr + Orr — &00 + 1 dor0 = o (4)
dr r r
dar0 , 1 da00 2
+----TuT +---Or0 = 0. (5)
dr r r
We write the auxiliary relations resulting from (2):
dur ur dur 1
dr r — ur'
a
2 r — ur
then from (1)
— _ 1 dar0 _ r (6)
_ 4 ^ r2 + dr _ 2 ^ (r2 + ' (6)
the prime denotes a derivative with respect to
Integrating (5) in view of (6) in relation to a^, we have:
^=^ (c <r)+i° <r2+*) — ),
where C(r) is some arbitrary function.
Now the function of additional hydrostatic pressure (and, consequently, the stress component arr) can be expressed from (3) by C(r) and and the difference (arr — a^) and dartp/d^ — by After the necessary transformations (4) becomes
1 dC 2-r2 • + +4-_ 8-
.4
0. (7)
2 dr 8^ r3• (r2 +
Differentiating (7) with respect to we have:
2- r4 ' + 4^ ) + 2- r2• ^ (^ ''' + 4^ ^') + • (2^ ^ ^ '' - ^'2 + 4^ = 0.
The last equation is divided in the case of = 0 into two autonomous differential equations with respect to
U '" + 4- = 0,
^ - '2 + 4^ = 0.
Note that d (2- ^ i>" - ^'2 + 4^ /d^ = 2- ' + 4^ ), hence, any non-trivial solution of the second equation of the system identically satisfies the first equation. Thus, equation (7) has an exact solution of the form:
= k cos2(^ + ^0), C'(r) = (k - 2^ r2) /r3,
where k and are constants of integration.
Symmetry conditions ar0 = ar0 = 0 are satisfied by ¿0 = 0 (or ¿>0 = n/2, solutions
0=0 0=n/2
are identical up to rotation of the coordinate axes through n/2). In the absence of stress at infinity, we have C(r) = - ln (r2) - k/ (2-r2). The constant of integration k can be considered as time-similar parameter of loading.
Obviously, for compliance with adopted kinematics (u = ur (r, the load on the surface of the cylindrical cutout can not contain functional arbitrariness.
Normal aL and tangent tl components of boundary load have the form aL = n- nT, tl = s^ a nT, where the value of the stress tensor a should be taken at the boundary L : r0(^>) in the actual configuration in view of r0(^) = \JR2(^) - n = (sin w, cos w, 0),
s = (-cosw, sinw, 0) are normal and tangential to L unit vectors, respectively. The function w(^) is the angle between the extension of the radius-vector and the tangent to L, it is defined
r0(^)
by the following relationship [3, Ch. VII, §2.232]: w(^>) = arctan
r0(W
Consider the case when the contour of the cavity in the undeformed state has the form of an ellipse with semi-axes equal to 1 and 1/2: R0(^>) = 1/\/sin2 ^ + 4- cos2 ^. Let the elastic modulus to be a = 1. The Fig. 3 shows graphs aL/b and TL/b for the loading parameter
Fig. 3. Components of boundary load Fig. 4. Stress intensity (von Mises)
value k = -3/4 (it corresponds to a deformed state with r0(0) = 1). The Fig. 4 shows isolines of tangential stress intensity (according to von Mises), which is also divided by the elastic modulus b.
Thus, the maximum stress intensity falls on the point at the boundary of computational domain with polar angle ф = 0. From this point the propagation of elasto-plastic boundary begins, when a load reaches some critical value.
Obtained solution can be useful for testing algorithms and programs for numerical integration of equations of the nonlinear theory of elasticity.
References
[1] M.Mooney, A theory of large elastic deformation, J. of App. Phys., 9(1940), no. 11, 582-592.
[2] A.I.Lurie, Nonlinear theory of elasticity, Amsterdam, North-Holland, 1990.
[3] G.M.Fikhtengolts, Calculus, V.1, M., Fizmatlit, 1962 (in Russian).
Точное аналитическое решение одной задачи о плоской деформации нелинейно-упругой среды
Георгий М. Севастьянов
Построено аналитическое решение задачи о плоской деформации изотропной нелинейно-упругой несжимаемой (резиноподобной) среды с цилиндрической полостью в квазистатическом приближении. Контур полости представляет собой гладкую симметричную кривую. Определенный частный вид следящего нагружения обеспечивает чисто радиальное перемещение материала. Массовые силы не учитываются. Физическая модель среды задана упругим потенциалом, аналогичным потенциалу Муни-Ривлина (с различием в используемом тензоре конечных деформаций). Полученное решение является точным: в уравнениях связи тензоров напряжений Коши и конечных деформаций Альманси сохранены все нелинейные члены (максимальная степень при выбранной модели среды — четвертая по компонентам градиента вектора перемещений).
Ключевые слова: плоская деформация, нелинейная упругость, несжимаемость, среда Муни-Ривлина, конечные деформации, тензор Альманси, точное аналитическое решение.
