Qualitative Analysis for Solutions of a Class of. Nonlinear Ordinary Differential Equations

Similar documents
Available online at ScienceDirect. Procedia Engineering 91 (2014 ) 32 36

Chapter 6 Thermoelasticity

About Some Inequalities for Isotonic Linear Functionals and Applications

4.2 Boussinesq s Theory. Contents

A COMPARISON OF MEMBRANE SHELL THEORIES OF HYBRID ANISOTROPIC MATERIALS ABSTRACT

s c s (b) Hence, show that the entropy for rubber-like materials must have the separable form

Micro-scale adhesive contact of a spherical rigid punch on a. piezoelectric half-space

arxiv: v1 [hep-th] 6 Jul 2016

Two dimensional polar coordinate system in airy stress functions

Friedmannien equations

Fourier-Bessel Expansions with Arbitrary Radial Boundaries

Modelling of Low Velocity Impact Damage in Laminated Composites

SPA7010U/SPA7010P: THE GALAXY. Solutions for Coursework 1. Questions distributed on: 25 January 2018.

On Natural Partial Orders of IC-Abundant Semigroups

International Journal of Pure and Applied Sciences and Technology

Radial geodesics in Schwarzschild spacetime

EXISTENCE OF THREE SOLUTIONS FOR A KIRCHHOFF-TYPE BOUNDARY-VALUE PROBLEM

On Polynomials Construction

9.4 The response of equilibrium to temperature (continued)

On the Eötvös effect

Right-indefinite half-linear Sturm Liouville problems

Michael Rotkowitz 1,2

RELATIVE KINEMATICS. q 2 R 12. u 1 O 2 S 2 S 1. r 1 O 1. Figure 1

( ) D x ( s) if r s (3) ( ) (6) ( r) = d dr D x

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING FLUID MECHANICS III Solutions to Problem Sheet 3

Investigations of Boundary Treatments in Incompressible Smoothed Particle Hydrodynamics for Fluid-Structural Interactions

igid nd non-leky two-comptment building. Yu et l [8] developed non-line govening equtions by consideing the effect of bckgound lekge. Howeve, thee e n

Research Article Hermite-Hadamard-Type Inequalities for r-preinvex Functions

A comparison principle for nonlinear heat Rockland operators on graded groups

Hadamard and Caputo-Hadamard FDE s with Three Point Integral Boundary Conditions

Effect of Heat Generation on Quasi- Static Thermal Stresses in a Solid Sphere

1 Using Integration to Find Arc Lengths and Surface Areas

A Parametric Study on the Centrifugal Force-Induced Stress and Displacements in Power-Law Graded Hyperbolic Discs

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

Electric Potential. and Equipotentials

Numerical Investigation of Flow in a New DC Pump MHD

Important design issues and engineering applications of SDOF system Frequency response Functions

Solution of fuzzy multi-objective nonlinear programming problem using interval arithmetic based alpha-cut

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Class Summary. be functions and f( D) , we define the composition of f with g, denoted g f by

Research Article Modeling of Thermal Distributions around a Barrier at the Interface of Coating and Substrate

This immediately suggests an inverse-square law for a "piece" of current along the line.

Discrete Model Parametrization

Integrals and Polygamma Representations for Binomial Sums

Elastic limit angular speed of solid and annular disks under thermomechanical

Answers to test yourself questions

SURFACE TENSION. e-edge Education Classes 1 of 7 website: , ,

General Physics II. number of field lines/area. for whole surface: for continuous surface is a whole surface

Simple analytical solutions for underground circular and elliptical openings

Poroelasticity of a covalently crosslinked alginate hydrogel under compression

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function

D-STABLE ROBUST RELIABLE CONTROL FOR UNCERTAIN DELTA OPERATOR SYSTEMS

New Expansion and Infinite Series

IMECE DRAFT

Topics for Review for Final Exam in Calculus 16A

Influence of Hydraulic Hysteresis on Effective Stress in Unsaturated Clay

Two Dimensional Inertial Flow of a Viscous Fluid in a Corner

STUDY OF THE UNIFORM MAGNETIC FIELD DOMAINS (3D) IN THE CASE OF THE HELMHOLTZ COILS

ANALYTICAL SOLUTION TO CONVECTION-RADIATION OF A CONTINUOUSLY MOVING FIN WITH TEMPERATURE-DEPENDENT THERMAL CONDUCTIVITY

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution

Tests for Correlation on Bivariate Non-Normal Data

Wave Generation by Oscillating Wall in Static Media

Math 4318 : Real Analysis II Mid-Term Exam 1 14 February 2013

Optimization. x = 22 corresponds to local maximum by second derivative test

Comparative Studies of Law of Gravity and General Relativity. No.1 of Comparative Physics Series Papers

DIRECT AND INVERSE PROBLEMS FOR THE EVERSION OF A SPHERICAL SHELL

U>, and is negative. Electric Potential Energy

PX3008 Problem Sheet 1

Theory of Elasticity Exam Problems and Answers Lecture CT5141 (Previously B16)

dx was area under f ( x ) if ( ) 0

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

Data Structures. Element Uniqueness Problem. Hash Tables. Example. Hash Tables. Dana Shapira. 19 x 1. ) h(x 4. ) h(x 2. ) h(x 3. h(x 1. x 4. x 2.

Experimental evaluation of the process of decohesion of adhesive joints with polymer films

Electricity & Magnetism Lecture 6: Electric Potential

On the integration of the equations of hydrodynamics

International Journal of Technical Research and Applications e-issn: , Special Issue 19 (June, 2015), PP.

Received 2 August 2014; revised 2 September 2014; accepted 10 September 2014

Lecture 11: Potential Gradient and Capacitor Review:

3.1 Magnetic Fields. Oersted and Ampere

7.2.1 Basic relations for Torsion of Circular Members

FI 2201 Electromagnetism

Collection of Formulas

Perturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion

defined on a domain can be expanded into the Taylor series around a point a except a singular point. Also, f( z)

(a) Counter-Clockwise (b) Clockwise ()N (c) No rotation (d) Not enough information

Creep modeling in functionally graded rotating disc of variable thickness

Week 8. Topic 2 Properties of Logarithms

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

10 m, so the distance from the Sun to the Moon during a solar eclipse is. The mass of the Sun, Earth, and Moon are = =

Geometry of the homogeneous and isotropic spaces

Fluids & Bernoulli s Equation. Group Problems 9

Supplementary material for " Coherent and Tunable Terahertz Radiation from Graphene Surface Plasmon Polarirons Excited by Cyclotron Electron Beam "

A CYLINDRICAL CONTACT MODEL FOR TWO DIMENSIONAL MULTIASPERITY PROFILES

Mathematical formulation of the F 0 motor model

THE STUDY OF THE SYSTEMS USED FOR THE POWER SOURCES COUPLING

arxiv: v2 [hep-th] 28 Jan 2016

1 INTRODUCTION. Dao Van Dung a Hoang Thi Thiem a, *

A Permanent Magnet Device for Measuring the Coercive Force

Electronic Supplementary Material

Transcription:

Adv. Theo. Appl. Mech., Vol. 7, 2014, no. 1, 1-7 HIKARI Ltd, www.m-hiki.com http://dx.doi.og/10.12988/tm.2014.458 Qulittive Anlysis fo Solutions of Clss of Nonline Odiny Diffeentil Equtions Juxin Li *, Wei Zho *, **, Cong Ren * nd Jing Zhng * * School of Science, Dlin Ntionlities Univesity, Dlin 116600, Chin ** Stte Key Lbotoy of Stuctul Anlysis fo Industil Equipment, Deptment of Engineeing Mechnics, Dlin Univesity of Technology, Dlin 116024, Chin Copyight 2014 Juxin Li, Wei Zho, Cong Ren nd Jing Zhng. This is n open ccess ticle distibuted unde the Cetive Commons Attibution License, which pemits unesticted use, distibution, nd epoduction in ny medium, povided the oiginl wok is popely cited. Abstct The boundy vlue poblem (BVP) fo clss of nonline odiny diffeentil equtions is exmined. It cn be used to descibe the evesion poblem fo spheicl shell composed of clss of tnsvesely isotopic incompessible Mooney-Rivlin mteils. The solutions of the nonline eqution descibing the eltion mong defomed thickness, initil thickness nd mteil pmetes of the spheicl shell e obtined. The effects of stuctue pmete nd mteil pmetes on the thickness of the eveted spheicl shell e discussed by numeicl exmples. Keywods : incompessible hypeelstic mteil, spheicl shell, evesion, defomtion, implicit solution 1. Intoduction Recently, the evesion poblems of bodies composed of hypeelstic mteils (such s ubbe nd ubbe-like mteils) widely ise fom engineeing design nd eospce fields, nd so on. These poblems cn be descibed s BVP fo clss of

2 Juxin Li, Wei Zho, Cong Ren nd Jing Zhng nonline odiny diffeentil equtions. Mny eseches studied the eveted spheicl shell theoeticlly nd expeimentlly. In 191, Amnni [1] fistly consideed the evesion poblem of spheicl shell composed of clss of compessible hypeelstic mteils. Eicksen [2] studied the existence nd uniqueness of the eveted solutions fo spheicl shell composed of incompessible Mooney-Rivlin mteil. Antmn [], Szei [4] gve futhe nlysis fo the existence of the eveted solutions fo compessible spheicl shells. In 199, Jeemih [5] mde compison between the qulittive fetues of the egul nd eveted shells composed of specil compessible mteil unde intenl pessue. It is found tht the stess distibution in both cses is diffeent nd lge pessue cn be sustined by the eveted shell. Fo some incompessible hypeelstic spheicl shells, Chen nd Hughton [6] poved tht, if the mteil stisfied the Bke-Eicksen inequlities, no mtte how the thickness ws, thee ws unique spheiclly symmetic eveted solution, nd they gve sufficient condition of cvity fomtion fo the eveted spheicl shell. Othewise, they lso found tht thicke spheicl shells could undego bifuction on evesion. In this ppe, the evesion poblem of spheicl shell composed of clss of tnsvesely isotopic incompessible Mooney-Rivlin mteils is consideed. The coesponding mthemticl model cn be teted s BVP. The implicit solution is obtined by using the incompessible condition nd the semi-invese method. The effects of stuctue pmete nd mteil pmetes on the thickness of the eveted spheicl shells e discussed by numeicl exmples. 2. Mthemticl model nd Solutions Hee we conside the solution of the evesion poblem fo n incompessible hypeelstic thin-wlled spheicl shell. Since the defomtion is spheicl symmety, the defomtion configution is given by = (R), 0 A R B, 0 < b, θ = π Θ, ϕ = Φ, (1) whee R nd e the dii of undefomed nd defomed shell, nd b e the inne nd oute dii, espectively. Note tht b = (A), = (B). The pincipl stetches e given by d λ = =, λθ = λϕ = = λ, (2) dr R whee λ is pmete descibing the xil stetch te. Since the mteil is

Qulittive nlysis fo solutions of clss of nonline ODE incompessible, we get λλλ θ ϕ = 1, then = B R. Moeove, b = B A. Fo diffeent mteils, thee e diffeent defomtion modes. So fo solving this poblem, we need to give specific stin-enegy function. In this ppe, suppose tht the spheicl shell is consist of dil tnsvesely isotopic incompessible Mooney-Rivlin mteil, whose stin-enegy function is [8] μ 1 W = + β 1 2 1 2 2 2 1 2 ( I ) + β ( I ) + α( λ ), () whee μ > 0 is the she modulus in infinitesiml defomtion, β ( 12 β 12) is dimensionless mteil constnt, α is dimensionless pmete mesuing the degee of nisotopy of the mteil, I nd 2 1 = λ + λ θ + λϕ I = λ λ + λ λ + λ λ 2 θ θ ϕ ϕ e the pincipl invints of the ight Cuchy defomtion tenso. The coesponding pincipl components of the Cuchy stess tenso e given by 2 1 1 2 2 1 1 σ = μ( ) ( + β + 2( β) λ + α) p, σθθ = μλ ( + β + ( β)(( ) + λ )) p, (4) whee p is the hydosttic pessue. In the bsence of body foce, the equilibium diffeentil equtions cn be educed to dσ 2 + ( σ σ θθ ) = 0. d (5) Integting Eq. (5) with espect to, we hve () () 2 σθθ σ σ = σ + d. (6) Assume tht the inne nd oute sufces of the eveted spheicl shell e tction-fee, the boundy conditions e s follows Then fom Eq. (6), we get σ ( ) = σ ( b) = 0. (7) b σθθ σ d = 0. (8)

4 Juxin Li, Wei Zho, Cong Ren nd Jing Zhng Substituting Eq. (4) into Eq. (8), we hve b 1 1 2 2 d ( ( ) ) + β + β λ α( ) = 0 λ. (9) Fo convenience, we intoduce the following nottions Then we obtin A δ =, B m =. (10) B b 1 δ + m M A δ 1 = = ( ), m η = = ( ) b 1 δ + m 1, B + 1 d dλ B + dλ λ = = ( 1), = =. (11) R λ λ 1+ λ λ Using the bove nottions, we ewite Eq. (9) s 1 1 1 1 1 1 1 1 1 α ( β) λ ( + α + β) + ( + α + β) + ( β) + ln( λ+ 1) 4 2 λ 4 2 λ λ α 1 1 1 + + = 4 2 2 ( ln(( λ ) ) ctn(2 ( λ ))) 0 M m. (12) Obviously, Eq. (12) descibing the finite defomtion of the eveted spheicl shell is nonline eqution with espect to δ nd m. Fom Eqs. (11), (12), we cn get the eltions mong the stuctue pmete δ, the mteil pmetes α, β nd the thickness η of the eveted spheicl shell.. Numeicl simultions Figs. 1-2 show the effects of the initil thickness δ, the mteil pmetes α nd β on the finite defomtion of the eveted spheicl shell composed of the incompessible mteil ().

Qulittive nlysis fo solutions of clss of nonline ODE O 5 () Fig. 1 Reltion cuves of δ ~ (b) m, η fo diffeent vlues of α ( β = 0.2 ) () (b)( Fig. 2 Reltion cuves of δ ~ m, η fo diffeent vlues v of β ( α = 2 ) Fom Fig. 1(), it is seen tht, fo the given vlues of β, the inne dius off the eveted spheicl shell inceses with the incesing initil thickness δ, nd the incesing eltions e line l ppoximtely. Menwhile, the inne dius off the eveted spheicl shell inceses withh the incesing the mteil pmete α, tht is to sy, the influence of the mteil pmete α on the inne dius nd thickness of the eveted thickness is significnt. Fom Fig. 1( (b), it is seen tht, fo the given vlues of β, the thickness of the eveted spheicl shell inceses with the incesing initil

6 Juxin Li, Wei Zho, Cong Ren nd Jing Zhng thickness δ. Howeve, the influence of the mteil pmete α on the thickness of the eveted thickness is not obvious. Fom Fig. 2(), it is seen tht, fo the given vlues of α, the inne dius of the eveted spheicl shell inceses with the incesing initil thickness δ. Moeove, the effect of the mteil pmete β on the inne dius of the eveted spheicl shell is not obvious. Fom Fig. 2(b), it is seen tht, fo the thickness of the eveted spheicl shell, thee e some simil popeties to the inne dius in Fig. 2(). 4. Conclusions The nonline equtions descibing the finite defomtion fo n eveted spheicl shell composed of dil tnsvesely isotopy Mooney-Rivlin mteil is consideed. The implicit solution is obtined by using the incompessible condition nd the semi-invese method. The conclusions show tht (1) The thinne the initil spheicl shell is, the bigge the defomed inne dius is. (2) The thinne the initil spheicl shell is, the thinne the eveted spheicl shell is. () The influence of the pmete α mesuing the degee of nisotopy of the mteil on the inne dius nd thickness of the eveted thickness is significnt. (4) The influence of the mteil pmete β on the inne dius nd thickness of the eveted thickness is not obvious. Acknowledgement This wok ws suppoted by the Ntionl Ntul Science Foundtion of Chin (No. 112200), the Pogm fo Lioning Excellent Tlents in Univesity (No. LR2012044) nd the Fundmentl Resech Funds fo Centl Univesities (No. DC120101124). Refeences [1] G.S. Amnni. Defomzioni finite dei solidi elstici isotopi, II. Nuovo Cimento, 10(191), 424-427. [2] J.L. Eicksen. Invesion of pefectly elstic spheicl shell. A. Angew Mth Mech, 5(1955), 82-85.

Qulittive nlysis fo solutions of clss of nonline ODE 7 [] S.S. Antmn. The evesion of thick spheicl shells. Ach Rt Mech Anlysis, 70(1979), 11-12. [4] A.J. Szei. On the eveted stte of spheicl nd cylindicl shells. Q Appl Mth, 48(1990), 49-58. [5] G.M. Jeemih. Infltion nd evesion of spheicl shells of specil compessible mteil. Jounl of Elsticity, 0(199), 251-276. [6] D.M. Hughton, Y.C. Chen. On the evesion of incompessible elstic spheicl shells. ZAMP, 50(1999), 12-26. [7] D.M. Hughton, Y.C. Chen. Asymptotic bifuction esults fo the evesion of elstic shells. ZAMP, 54(200), 191-211. [8] P. Chdwick. The existence nd uniqueness of solutions of two poblems in the Mooney-Rivlin theoy fo ubbe. Jounl of Elsticity, 2(1972), 12 128. Received: My 1, 2014