# Nonlinear Waves in Rods and Beams of Power-Law Materials.

1. IntroductionFree vibrations of rods and beams of power-law materials are considered. Analytic traveling wave solutions to the wave equations for power-law materials (see [1, 2]) are obtained which represent kink waves of single elevation that propagates without change of forms in time. It is shown that, unlike the wave equations for linear materials, the nonlinear wave equations do not allow arbitrary traveling wave forms in an infinite rod or beam. The results demonstrate that the traveling fronts of the waves may sharpen or flatten as the wave speeds increase depending upon the power-law index n and the bulk modulus. For n >1, the wave fronts sharpen, whereas for 0 < n < 1, the fronts flatten as the wave speeds increase. The solutions also demonstrate that the speeds of the nonlinear traveling waves depend not only on the material properties but also on the initial energy-level. It is well known that the speeds of waves for the linear elastic materials (n = 1, Hooke's law) depend only on the material properties in contrast to that of the waves in nonlinear materials. As far as we know these solutions are not available in literature, even though there are numerous research papers and books devoted to the discovery and study of traveling waves in elastic and plastic solids (see [3-7] for details). In the case of rods and beams of finite length, we also present some special solutions satisfying certain boundary and initial conditions. The closed formula solutions are expressed in terms of nonEuclidean sine functions (cf. [8]), which differ from the Euclidean sine functions corresponding to the waves in rods and beams of linear elastic materials.

The note is organized as follows. In Section 2, the powerlaw constitutive stress-strain equation is introduced. In Section 3, the potential energy and derivations of the wave equations of power-law materials are outlined. In Sections 4 and 5, closed-form solutions are derived. And, finally the results are summarized in Section 6.

2. Hollomon's Equation

It is well known that, in uniaxial state, the following powerlaw stress and strain relation is used for certain elastoplastic materials:

[sigma] = K[ansolute value of ([[epsilon].sup.n-1] [epsilon], 0 < n < [infinity], (1)

where [sigma] is the axial stress, e is the axial strain, and K and n are engineering constants with values depending on the specific material. The materials satisfying (1) sometimes are also referred to as Ludwick or as Hollomon's materials in literature (cf. [1,2]). Many heat-treated metals are well-known power-law materials. For a given annealed metal or alloy, K and n depend on the heat treatment received by the metal or alloy. The values of n are typically between 0 and 1 for such metals. For a comprehensive list of experimental values of K and n of common annealed industrial metals, see, for example, [9]. For some geological materials, such as certain rocks or ice, however, the values of n are greater than 1. In some biological tissues, experiments also indicate that the power-law index n satisfies 0 < n < 1 for bones such as tibia and femur, while n > 1 for cartilages such as common carotid artery and abdominal aorta (see, e.g., [10,11]). For a given value of 0 < n < 1, the stress-strain curve defined by (1) can result in a rapid increase in the yield stress for small strains or strain hardening. However, it can be the opposite for values of n > 1, for which large strains produce small stress or softening. For these reasons, n is called the strain-hardening or strain-softening exponent. Study of the mechanical properties of these heat-treated metals is very important in industries (see, e.g., [12], for stress analysis of beam columns made of Ludwick materials). If we allow n = 1, then (1) reduces to Hooke's law for linear elastic material and the constant K, also called the bulk modulus, equals the corresponding Young's modulus E. Power-law materials are a special case of a more general class of materials called Hencky plastics [13]. Physically, the constitutive equation (1) describes the hardening or softening of materials showing an elastic-plastic transition. In the following, bold letters are used to denote vectors or matrices. A vector is considered as a single row matrix. The transpose of a matrix A is denoted by [A.sup.[tau]], and the inner product of two vectors u and v by [uv.sub.[tau]]. The time derivative [partial derivative]u/[partial derivative]t is denoted by [??]. Let u(x, y, z, t) = (u(x, y, z, t), v(x,y,z,t),w(x,y,z,t)) denote the displacement vector,

[mathematical expression not reproducible] (2)

the strain components, and [[sigma].sub.x], [[sigma].sub.y], [[sigma].sub.z], [[tau].sub.xy], [[tau].sub.yz], and [[tau].sub.zx] the corresponding stress components. The following generalized power law can be derived from the Hencky total deformation theory [13]:

[mathematical expression not reproducible], (3)

where [parallel]D(u)[parallel] = [square root of ([[epsilon].sup.2.sub.x] + [[epsilon].sup.2.sub.y] + + [[epsilon].sup.2.sub.z] + 2[[gamma].sup.2.sub.xy] + 2[[gamma].sup.2.sub.yz] + 2[[gamma].sup.2.sub.zx])], where n, K, and v are the material constants; see also Wei [14]. Note that (3) is the three-dimensional version of (1). In the following two sections, wave equations of bars and beams made of the power-law elastoplastic materials are derived by (3) and the assumption of the Euler-Bernoulli beam theory. There are similar versions of generalized powerlaw stress-strain relations for strain-hardening or strainsofteningmaterialin theliteratureandsimilarwaveequations can be derived (see, e.g., [15-20]).

3. The Nonlinear Wave Equations

The potential energy for a power-law elastoplastic body occupying a three-dimension body V can by defined by

[mathematical expression not reproducible], (4)

where [epsilon] = ([[epsilon].sub.x], [[epsilon].sub.y], [[epsilon].sub.z], [[gamma].sub.xy], [[gamma].sub.xz], [[gamma].sub.yz]) and [sigma] = ([[sigma].sub.x], [[sigma].sub.y], [[sigma].sub.z], [[tau].sub.xy], [[tau].sub.xz], [[tau].sub.yz]). The Lagrangian energy functional 7(u) equals the kinetic energy T minus the elastoplastic potential energy U plus the work W done by external force. It can be written as

[mathematical expression not reproducible] (5)

where p is the density, [mathematical expression not reproducible] the velocity, f = ([f.sub.x], [f.sub.y], [f.sub.z]) the body force, and t = ([t.sub.x], [t.sub.y], [t.sub.z]) the surface force. See, for example, [21], for a standard definition of I(u). For a uniaxial bar of infinite length with cross-sectional area A(x), subject to axial force and zero surface force, we have [sigma] = ([[sigma].sub.x], 0, 0, 0, 0, 0), u = (u(x, t), 0, 0), [[sigma].sub.x] = K[[absolute value of ([[epsilon].sub.x])].sup.n-1][[epsilon].sub.x], f = (f(x, t), 0, 0), and t = (0, 0, 0). For an Euler beam of infinite length, it is assumed that the components of the displacement satisfy u(x, y, t) = -y([partial derivative]v/[partial derivative]x), v = v(x, t), w = 0, f = (0, r(x, t), 0), and t = (0,0,0). Therefore = [[epsilon].sub.x] [partial derivative]u/[partial derivative]x = -y([[partial derivative].sup.2]v/[partial derivative][x.sup.2]), [[epsilon].sub.xy] = (1/2)([partial derivative]u/[partial derivative]y + [partial derivative]v/[partial derivative]x) = 0, and [[epsilon].sub.y] = [[epsilon].sub.xz] = [[epsilon].sub.yz] = [[epsilon].sub.z] = 0. The potential energies for the bar and the beam are given by

[mathematical expression not reproducible], (6)

[mathematical expression not reproducible], (7)

respectively, where [mathematical expression not reproducible] is the generalized second moment of inertia of the beam. The x-axis is taken to be the axial direction of the bar and the beam. For rods and beams of finite length L, the corresponding Lagrangian functions are given by replacing -[infinity] and +[infinity] in (6) and (7) by 0 and L, respectively. Note that the assumptions made in this section on the elastoplastic bars and the beams are standard assumptions frequently made for elastic bars and beams (see, e.g., [22, 23], for details). The corresponding linear wave equations of elastic bars and beams corresponding to n =1 have been studied extensively.

For completeness, the derivation of the wave equations of the power-law materials given in [14] is outlined here. It is well known that Hamilton's principle seeks an equilibrium state in time dependent mechanical systems (see, e.g., [21]).

Specifically, Hamilton's principle requires that we seek a displacement u so that, for any time interval [[t.sub.1], [t.sub.2]], u([t.sub.1]) = u([t.sub.2]) and [??]([t.sub.1]) = [??]([t.sub.2]), and for all displacement of the form u + tv, where t is any real number, the first variation of the energy functional I satisfies

[mathematical expression not reproducible] (8)

for all v satisfying v([t.sub.1]) = v([t.sub.2]) = 0 and [??]([t.sub.1]) = [??]([t.sub.2]) = 0. The combination u(t) + [tau]v(t) is referred to as an admissible displacement for the mechanical system since it is required to satisfy some boundary conditions. It can be shown that if the displacement u satisfies (8) of Hamilton's principle, then it must also satisfy a differential wave equation under certain conditions. In particular, suppose that the cross-sectional area, denoted by A, is a nonzero constant, and then for the rod, we have

[rho] [[partial derivative].sup.2]u/[partial derivative][t.sup.2] = K[partial derivative]/[partial derivative]x ([[absolute value of ([partial derivative]u/[partial derivative]x)].sup.n-1] [partial derivative]u/[partial derivative]x) + f, x [member of] R, t [member of] [R.sup.+] (9)

and for the corresponding Euler beam

[mathematical expression not reproducible]. (10)

When n = 1, (9) reduces to the standard wave equation for the elastic bar

[rho][[partial derivative].sup.2]u/[partial derivative][t.sup.2] = K[[partial derivative].sup.2]u/[partial derivative][x.sup.2] + f, x [member of] R, t [member of] [R.sup.+]. (11)

and (10) to the standard wave equation for the elastic Euler beam

[rho]A [[partial derivative].sup.2]v/[partial derivative][t.sup.2] = -[[partial derivative].sup.2]/[partial derivative][x.sup.2] (KI [[partial derivative].sup.2]v/[partial derivative][x.sup.2]) + Ar, x [member of] R, t [member of] [R.sup.+]. (12)

The quantity In reduces to the second moment of inertia, [mathematical expression not reproducible] reduces to I when n = 1 in the elastic beam theory, and the material constant K becomes Young's modulus E for linear elastic materials. In deriving the wave equations (9) and (10), we have made the assumption that the solutions u and v are continuously differentiable and their appropriate lower order derivatives are bounded or vanishing when [absolute value of (x)] [member of] [infinity]. By (8), we get

[mathematical expression not reproducible]. (13)

Using integration by parts and interchange of the order of integration, with v([t.sub.1]) = v([t.sub.2]) = 0, and assuming that [lim.sub.x[right arrow][+ or -][infinity]][[absolute value of ([partial derivative]u(x, t) [partial derivative])].sup.n-1] ([partial derivative]u(x, t)/[partial derivative]x) is bounded by a constant independent of t and [lim.sub.x[right arrow][+ or -][infinity]]v(x, t) = 0 uniformly in t, we get the following:

[mathematical expression not reproducible] (14)

from (13). Since v, [t.sub.1], and [t.sub.2] are arbitrary and A [not equal to] 0, we then get (9) from (14). The corresponding beam equation (10) can be derived similarly which was reported in [14].

4. Traveling Waves in Rods and Beams of Arbitrary Length

In the following we will derive some traveling wave solutions to (9) and (10) for 0< n < [infinity] and n [not equal to] 1. As far as we know, these solutions are not available in literature, even though there are numerous research papers and books devoted to the discovery and study of traveling waves in elastic and plastic solids. For the study of traveling waves in nonlinear beam equations based on Hooke's law (n = 1) for elastic materials, see, for example, [3-5]. Also, see [6, 7, 24, 25], for more results of traveling waves in solids. Assuming that f = r =0 in (9) and (10), we have

[[partial derivative].sup.2]u/[partial derivative][t.sup.2] = [c.sup.2][partial derivative]/[partial derivative]x ([[absolute value of ([partial derivative]u/[partial derivative]x)].sup.n-1] [partial derivative]u/[partial derivative]x) + x [member of] R, t [member of] [R.sup.+] (15)

for the bar and

[[partial derivative].sup.2]v/[partial derivative][t.sup.2] = [[??].sup.2] [[partial derivative].sup.2]/[partial derivative][x.sup.2] ([[absolute value of ([[partial derivative].sup.2]v/[partial derivative][x.sup.2])].sup.n-1] [[partial derivative].sup.2]v/[partial derivative][x.sup.2]), x [member of] R, t [member of] [R.sup.+] (16)

for the beam, where [c.sup.2] = K/[rho] and [[??].sup.2] = -K[I.sub.n]/[rho]A, respectively. We look for traveling wave solutions of the form g(x - [lambda]t) for both (15) and (16), where A denotes a constant and g is a function to be determined. Let [phi](t) = [[absolute value of (t)].sup.n-1]t, where n is the index in power-law (3). The inverse of [phi] is [[phi].sup.-1](t) = [[absolute value of (t)].sup.(1-n)/n]t since

[mathematical expression not reproducible]. (17)

First, let [xi] = x - [lambda]t and substitute u(x, t) = g([xi]) into (15), so [[lambda].sup.2]g" = [c.sup.2](([phi](g'))'. After integration we get [[lambda].sub.2]g' = [c.sup.2][phi](g') + [c.sub.1], where [c.sub.1] is an arbitrary constant. Suppose that [lim.sub.x[right arrow][infinity]]([partial derivative]u/[partial derivative]x)(x, 0) = A. Since u(x, 0) = g(x) and g'(x) = ([partial derivative]u/[partial derivative]x)(x,0), we have [c.sub.1] = [[lambda].sup.2]A - [c.sup.2][phi](A). Looking for vnontrivial solutions for n [not equal to] 1 and assuming that A = [(c/[lambda]).sup.2/(1-n)]), we get [c.sub.1] = 0 and g' = [[+ or -](c/[lambda]).sup.2/(1-n)] which gives the following traveling wave solutions:

u(x,t) = [+ or -][(c/[lambda]).sup.2/(1-n)](x-[lambda]t) + [c.sub.2] (18)

for the bar equation (15). Note that solution (18) includes some physically meaningful solutions. For example, let us consider a semi-infinite bar with initial displacement

[mathematical expression not reproducible], (19)

and initial velocity

[mathematical expression not reproducible] (20)

and boundary condition [lim.sub.x[right arrow]+[infinity]]([partial derivative]u/[partial derivative]x)(x, t) = [(c/[lambda]).sup.2/(-1-n)] A particular solution satisfying these conditions is given by

[mathematical expression not reproducible] (21)

which is obtained from (18). The physical interpretation of the initial condition (19) is that half of the bar is initially subject to a constant stress and the other half is free of stress and fixed in position, and the second initial condition (20) means that the bar is initially moving at a constant speed and half of it is instantaneously stopped. Solution (21) explains that if a prestressed semi-infinite axial power-law rod subject to initial conditions [sigma](x, 0) = K[(c/[lambda]).sup.2n/(1-n)] and [??](x, 0) > 0 and boundary conditions u(0, t) = 0 for x = 0 and lower order derivatives are bounded or vanishing when [absolute value of (x)] [right arrow] [infinity], then the displacement in the interval [0, x] will be zero at time t = [[x.sup.1+n][absolute value of ([x.sup.1+n][[absolute value of ([??](x,0))] .sup.1-n]/[c.sup.2]]m+n) = x/[lambda] and the restoration of the deformed bar in interval [0, x(t)] to its undeformed configuration has a moving boundary x(t) which is expanding like a kink wave at a velocity A = c[(K/[sigma](x, 0)).sup.(1-n)/2n].

In the more general situation, for any value of q, the equation P(t) = [[lambda].sup.2]t - [c.sup.2][phi](t) - [c.sub.1] has at least one solution since it is continuous, [lim.sub.x[right arrow]+[infinity]]P(t) = +[infinity] and [lim.sub.t[right arrow]-[infinity]] P(f) = -[infinity]. Let P([[epsilon].sub.0]) = 0, and then g = e0 satisfies [[lambda].sup.2]g' = [c.sup.2][phi](g') + [c.sub.1]. We have the following similar solutions:

[mathematical expression not reproducible] (22)

satisfying the initial and boundary conditions

[mathematical expression not reproducible] (23)

For the linear elastic bar, n =1, [phi](g') = g', and if [c.sub.1] = 0, the equation [lambda]2g' = [c.sup.2][phi](g') is satisfied for any g' and also makes [lambda] = c. This shows that the linear elastic bar equation allows arbitrary traveling wave forms g in u(x, t) = g(x - [lambda]t), and however the wave can travel only at a fixed velocity [lambda] = C. If [c.sub.1] [not equal to] 0, then equation [[lambda].sup.2]g' = [c.sup.2][phi](g') + [c.sub.1], gives g' = [[epsilon].sub.0] = ([[lambda].sup.2] - [c.sup.2])/[c.sub.1] and the corresponding solution is (22), which is similar to the solutions for n [not equal to] 1.

The above shows that the difference between the nonlinear solution (n =1) and the linear case (n = 1) is that all the nonlinear traveling waves have the same shape and the traveling velocity depends not only on the material property but also on the initial stress-level while the linear traveling waves can take any form while keeping a fixed traveling velocity c that depends only on the material property.

Similarly, by substituting v(x, t) = g(x-[lambda]t) intothebeam equation (16), we get [[lambda].sub.2]g = [[??].sup.2]([phi](g"))". After integration twice, we get [[lambda].sup.2]g = [[??].sup.2][phi](g") + [c.sub.1][xi] + [c.sub.2], which gives [[lambda].sup.2]g = [[??].sup.2][phi](g") by setting [c.sub.1] = [c.sub.2] = 0. Let w = g'; we get w(dw/dg) = g" and [[lambda].sup.2]g = [[??].sup.2][phi](w(dw/dg)). From the last equation, we get

[mathematical expression not reproducible] (24)

which gives [(g').sup.2] = C - (2n/(1 + n))[[absolute value of ([lambda]/[??])].sup.2/n][[absolute value of (g)].sup.(1+n)/n], where C is the integration constant.

We assume g'(0) = 0 and g(0) > 0. So, the traveling wave solutions for the corresponding elastoplastic Euler beam are given implicitly

[mathematical expression not reproducible] (25)

which results in the following formula in terms of generalized trigonometric function defined in [26]

g(x - [lambda]t) = g(0) [sin.sub.2,1+1/n] (B (x - [lambda]t)), (26)

where B = [[absolute value of ([??])].sup.1/n]/[[absolute value of (g(0))].sup.(n+1)/2n/ g(0)[[absolute value of ([lambda])].sup.1/n][square root of (2n/(n + 1))]. Notice that for we obtain the well-known Euclidean sine traveling wave solution for the elastic Euler beam equation. We observe that the amplitude of the wave is determined by the initial condition v(0, 0). The traveling waves for elastoplastic beams can be applied to study piezoelectric robots; see [27].

5. Special Waves in Rods and Beams of Finite Length

Let us consider some special waves in rods and beams of finite length L with fixed ends. Equations (15) and (16) are solved for x [member of] (0, L) and t [member of] [R.sup.+] with homogeneous Dirichlet boundary conditions u(0, t) = u(L, t) = 0 and special initial conditions. We present some special solutions by using the generalized trigonometric functions developed by Drabek and Manasevich [26]. By using the separation of variables

u(x, t) = X(x)T(t) in (15) and using boundary conditions u(0, t) = u(L, t) = 0, we have

[mathematical expression not reproducible]. (27)

From the first equation and the boundary conditions, we get

[mathematical expression not reproducible]. (28)

By Thm 3.1 in [26], a sequence of solutions to the nonlinear eigenvalue problem (28) are given by [X.sub.k] = [sin.sub.n+12](k([x[pi].sub.n+12]/ L)), where [[lambda].sub.k] = -((k+1)kl2)[([[pi].sub.n+1,2]).sup.n+1] and [mathematical expression not reproducible]. Let us consider the initial conditions u(x, 0) = [X.sub.1](x) and [u.sub.t](x, 0) = 0. In this case we have T(0) = 1 and we can solve the second equation in (27). A special solution of this initial value problem is given by

[mathematical expression not reproducible]. (29)

The time evolution of the special solutions for c =1, L =1, n = 0.2, and n= 0.5 is presented in Figure 1, respectively.

Similarly, by using the separation of variables u(x, t) = X(x)T(t) in (16), we have

[mathematical expression not reproducible]. (30)

From the first equation and the boundary condition, we get

[mathematical expression not reproducible]. (31)

An analytic solution to (31) is not available and is an open problem. This is a nonlinear and nonhomogeneous eigenvalue problem which belongs to an active area of research beyond the scope of this paper, and we post it here as an open problem. Since superposition principle can not be applied to nonlinear problems, the solutions to (15) and (16) with general initial and boundary conditions require further investigations.

6. Conclusions

Two nonlinear wave equations are derived: one is for the longitudinal vibrations of a power-law bar and the other is for vertical vibrations of the power-law Euler beam. Analytic traveling wave solutions are found for these two equations for free vibrations in terms of generalized sine functions of two parameters. We recovered the linear elastic waves as special cases of our solutions. The traditional ways of determining the vibrations of a structure made of untreated metals do not apply to the structures of heat-treated metals with hardening and softening mechanical properties. The obtained results can be useful in engineering applications of the power-law materials, such as heat-treated metals and polyimide plastics. Further study of wave propagation and vibrations in structures made of the power-law nonlinear bars and beams seems necessary.

https://doi.org/10.1155/2017/2095425

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] K. Lee, "Large deflections of cantilever beams of non-linear elastic material under a combined loading," International Journal of Non-Linear Mechanics, vol. 37, no. 3, pp. 439-443, 2001.

[2] G. Lewis and F. Monasa, "Large deflections of cantilever beams of non-linear materials of the Ludwick type subjected to an end moment," International Journal of Non-Linear Mechanics, vol. 17, no. 1, pp. 1-6,1982.

[3] S. S. Antman, Nonlinear Problems of Elasticity, Springer-Verlag, New York, NY, USA, 1995.

[4] A. R. Champneys, P. J. McKenna, and P. A. Zegeling, "Solitary waves in nonlinear beam equations: stability, fission and fusion," Nonlinear Dynamics, vol. 21, no. 1, pp. 31-53, 2000.

[5] A. V. Porubov and M. G. Velarde, "Dispersive-dissipative solitons in nonlinear solids," Wave Motion, vol. 31, no. 3, pp. 197-207, 2000.

[6] W.-S. Duan and J.-B. Zhao, "Solitary waves in a quartic nonlinear elastic bar," Chaos, Solitons and Fractals, vol. 11, no. 8, pp. 1265-1267, 2000.

[7] V. I. Erofeev and N. V Klyueva, "Solitons and nonlinear periodic strain waves in rods, plates, and shells (a review)," Acoustical Physics, vol. 48, no. 6, pp. 725-741, 2002.

[8] D. Wei, Y. Liu, and M. B. Elgindi, "Some generalized trigonometric sine functions and their applications," Applied Mathematical Sciences, vol. 6, no. 121-124, pp. 6053-6068, 2012.

[9] J. F. Shackelford, in Introduction to Materials Science for Engineers, Prentice-Hall, Englewood Cliffs, NJ, USA, 5th edition, 2000.

[10] Y. C. Fung, Biomechanics: Mechanical Properties of Living Tissue, Springer-Valag, New York, NY, USA, 2nd edition, 1993.

[11] A. J. Grodzinsky, R. D. Kamm, and D. A. Lauffenburger, "Quantitative aspects of tissue engineering: basic issues in kinetics, transport and mechanics," in Principles of Tissue Engineering, R. Lanza, R. Langer, and W. Chick, Eds., R.G. Landes Co, Austin, Tex, USA, 1997.

[12] X. Teng and T. Wierzbicki, "Crush response of an inclined beam-column," Thin-Walled Structures, vol. 41, no. 12, pp. 1129-1158, 2003.

[13] H. Hencky, "Zur Theorie plastischer Deformationen und der hierdurch im Material hervorgerufenen Nachspannungen," Journal of Applied Mathematics and Mechanics, vol. 4, no. 4, pp. 323-334, 1924.

[14] D. Wei, "Nonlinear wave equations arising in modeling of some power-law elasto-plastic structures," in Proceedings of ICCP6 and CCP, pp. 248-251, Beijing, China, 2005.

[15] C. Atkinson and C. Y. Chen, "On interface dislocations between dissimilar materials with non-linear stress-strain laws (plane strain)," International Journal of Engineering Science, vol. 37, no. 5, pp. 553-573, 1999.

[16] J. Y. Chen, Y. Huang, K. C. Hwang, and Z. C. Xia, "Plane-stress deformation in strain gradient plasticity," Journal of Applied Mechanics, Transactions ASME, vol. 67, no. 1, pp. 105-111, 2000.

[17] C. H. Chou, J. Pan, and S. C. Tang, "Nonproportional loading effects on elastic-plastic behavior based on stress resultants for thin plates of strain hardening materials," International Journal of Plasticity, vol. 10, no. 4, pp. 327-346,1994.

[18] X.-L. Gao, "A mathematical analysis of the elasto-plastic plane stress problem of a power-law material," IMA Journal of Applied Mathematics, vol. 60, no. 2, pp. 139-149, 1998.

[19] A. E. Giannakopoulos, "Total deformation, plane-strain contact analysis of macroscopically homogeneous, compositionally graded materials with constant power-law strain hardening," Journal of Applied Mechanics, Transactions ASME, vol. 64, no. 4, pp. 853-860,1997.

[20] Z. K. Wang, X. X. Wei, and X. L. Gao, "A general analytical solution of a strain-hardening elastoplastic plate containing a circular hole subject to biaxial loading-with applications in pressure vessels," International Journal of Pressure Vessels and Piping, vol. 47, no. 1, pp. 35-55,1991.

[21] M. Petyt, Introduction to Finite Element Vibration Analysis, Cambridge University Press, Cambridge, UK, 1990.

[22] A. C. Ugural, Stresses in Plates and Shells, McGraw-Hill, New York, NY, USA, 1981.

[23] S. Timoshenko, D. H. Young, and W. Weaver Jr., Vibration Problems in Engineering, John Wiley & Sons, Hoboken, NJ, USA, 1974.

[24] J. Lamb, Elements of Soliton Theory, John Wiley & Sons, Inc, New York, NY, USA, 1982.

[25] M. Remoissenet, Waves Called Solitons-Concepts and Experiments, Springer, Berlin, Germany, 2nd edition, 1996.

[26] P. Drabek and R. Manasevich, "On the closed solution to some nonhomogeneous eigenvalue problems with p-Laplacian," Differential and Integral Equations, vol. 12, no. 6, pp. 773-788, 1999.

[27] H. Hariri, Y. Bernard, and A. Razek, "A traveling wave piezoelectric beam robot," Smart Materials and Structures, vol. 23, no. 2, Article ID 025013, 2014.

Dongming Wei, (1) Piotr Skrzypacz, (1) and Xijun Yu (2)

(1) Department of Mathematics, School of Science and Technology, Nazarbayev University, Astana 010000, Kazakhstan

(2) Institute of Applied Physics and Computational Mathematics, P.O. Box 8009, Beijing 100088, China

Correspondence should be addressed to Dongming Wei; dongming.wei@nu.edu.kz

Received 16 March 2017; Revised 4 May 2017; Accepted 4 June 2017; Published 13 July 2017

Academic Editor: Xin-Lin Gao

Caption: FIGURE 1: Vibrating power-law strings: n = 0.2 (a) and n = 0.5 (b) at t = 0, 0.2, 0.4, 0.5, 0.7, 0.8.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wei, Dongming; Skrzypacz, Piotr; Yu, Xijun |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 4564 |

Previous Article: | Hybrid Algorithm of Particle Swarm Optimization and Grey Wolf Optimizer for Improving Convergence Performance. |

Next Article: | Simulation of Wellbore Stability during Underbalanced Drilling Operation. |

Topics: |