A New I 1 -Based Hyperelastic Model for Rubber Elastic Materials

Similar documents
New Approach to Predict the Micromechanical. Behavior of a Multiphase Composite Material

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE

Computational Methods of Solid Mechanics. Project report

Orbital Angular Momentum Eigenfunctions

Adsorption and Desorption Kinetics for Diffusion Controlled Systems with a Strongly Concentration Dependent Diffusivity

FE FORMULATIONS FOR PLASTICITY

1 Random Variable. Why Random Variable? Discrete Random Variable. Discrete Random Variable. Discrete Distributions - 1 DD1-1

8-3 Magnetic Materials

Numerical simulation of combustion frozen nanosized aluminum suspension with water

ATMO 551a Fall 08. Diffusion

30 The Electric Field Due to a Continuous Distribution of Charge on a Line

Application of Poisson Integral Formula on Solving Some Definite Integrals

Lecture 23: Central Force Motion

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007

Right-handed screw dislocation in an isotropic solid

Numerical Integration

A question of Gol dberg concerning entire functions with prescribed zeros

Some Remarks on the Boundary Behaviors of the Hardy Spaces

Probability Distribution (Probability Model) Chapter 2 Discrete Distributions. Discrete Random Variable. Random Variable. Why Random Variable?

Effect of axial prestretch on inflation instability in finitely extensible thin-walled tube

Some Ideal Convergent Sequence Spaces Defined by a Sequence of Modulus Functions Over n-normed Spaces

Study on GPS Common-view Observation Data with Multiscale Kalman Filter. based on correlation Structure of the Discrete Wavelet Coefficients

Vortex Initialization in HWRF/HMON Models

LECTURE 15. Phase-amplitude variables. Non-linear transverse motion

This is the published version of a paper published in Mechanics of composite materials. Citation for the original published paper (version of record):

JORDAN CANONICAL FORM AND ITS APPLICATIONS

ALOIS PANHOLZER AND HELMUT PRODINGER

Robust Spectrum Decision Protocol against Primary User Emulation Attacks in Dynamic Spectrum Access Networks

Finite Element Modeling and Free Vibration Analysis of Functionally Graded Nanocomposite Beams Reinforced by Randomly Oriented Carbon Nanotubes

Central limit theorem for functions of weakly dependent variables

Maximum Torque Control of Induction Traction Motor Based on DQ Axis Voltage Regulation

The Concept of the Effective Mass Tensor in GR. Clocks and Rods

BImpact of Supply Chain Coordination for Deteriorating Goods with Stock-Dependent Demand Rate

Stress, Cauchy s equation and the Navier-Stokes equations

τ TOT = r F Tang = r F sin φ Chapter 13 notes: Key issues for exam: The explicit formulas

2. Plane Elasticity Problems

F O R SOCI AL WORK RESE ARCH

Equations to Calculate Characteristic Frequencies of Multiple Chamber Aligned in Parallel Cavity Resonator (MCAP-CR)

Optimum Settings of Process Mean, Economic Order Quantity, and Commission Fee

LINEAR MOMENTUM Physical quantities that we have been using to characterize the motion of a particle

On the Correct Formulation of the Starting Point of Classical Mechanics

7.2.1 Basic relations for Torsion of Circular Members

Chapter 6 Differential Analysis of Fluid Flow

A new model for calculating impact force and energy dissipation based on the CR-factor and impact velocity

AMACHINE-to-machine (M2M) communication system

Simulation of Containment Atmosphere Stratification Experiment Using Local Instantaneous Description

KANTOROVICH TYPE INEQUALITIES FOR THE DIFFERENCE WITH TWO NEGATIVE PARAMETERS. Received April 13, 2010; revised August 18, 2010

Examine 3D Energy Balance Utility

Class 6 - Circular Motion and Gravitation

2 Governing Equations

On thermal boundary layers on a flat plate subjected to a variable heat flux

Chapters 1 8: Overview

Induction Motor Identification Using Elman Neural Network

Queuing Network Approximation Technique for Evaluating Performance of Computer Systems with Hybrid Input Source

Geometry and statistics in turbulence

Handout: IS/LM Model

Center of Mass and Linear

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

On the velocity autocorrelation function of a Brownian particle

Newton s E = mc 2 Two Hundred Years Before Einstein? Newton = Einstein at the Quantum Scale

Torque, Angular Momentum and Rotational Kinetic Energy

VLSI IMPLEMENTATION OF PARALLEL- SERIAL LMS ADAPTIVE FILTERS

Stress Intensity Factor

Distributed Adaptive Networks: A Graphical Evolutionary Game-Theoretic View

On the elastic properties of arteries

This paper is dedicated to the memory of Donna L. Wright.

On generalized Laguerre matrix polynomials

Assessment of the general non-linear case ..=(LK).", J'

Research Article Approximation of Signals (Functions) by Trigonometric Polynomials in L p -Norm

TO THE STATISTICS OF DOUBLE STARS

Fracture Mechanics: Linear Elastic Fracture Mechanics 2/2

To Feel a Force Chapter 7 Static equilibrium - torque and friction

A New Approach to General Relativity

VELOCITY PROFILE MODELING FOR NON-ISOTHERMAL FLOWS INSIDE A CIRCULAR TUBE

Long-range stress re-distribution resulting from damage in heterogeneous media

HAM for finding solutions to coupled system of variable coefficient equations arising in fluid dynamics

FARADAY'S LAW dt

Determining solar characteristics using planetary data

Modeling Fermi Level Effects in Atomistic Simulations

Chapter 31 Faraday s Law

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

Verified Solution for a Statically Determinate Truss Structure with Uncertain Node Locations

Psychometric Methods: Theory into Practice Larry R. Price

On Bounds for Harmonic Topological Index

Black Body Radiation and Radiometric Parameters:

rect_patch_cavity.doc Page 1 of 12 Microstrip Antennas- Rectangular Patch Chapter 14 in Antenna Theory, Analysis and Design (4th Edition) by Balanis

Chem 453/544 Fall /08/03. Exam #1 Solutions

5. Pressure Vessels and

Article : 8 Article : 8 Stress Field. and. Singularity Problem

Multi-relational Weighted Tensor Decomposition

Introduction to Money & Banking Lecture notes 3/2012. Matti Estola

, the tangent line is an approximation of the curve (and easier to deal with than the curve).

Electrostatics (Electric Charges and Field) #2 2010

Value Prediction with FA. Chapter 8: Generalization and Function Approximation. Adapt Supervised Learning Algorithms. Backups as Training Examples [ ]

Charged particle motion in magnetic field

SIMPLIFIED APPROACH FOR THE ANALYSIS OF A VISCOELASTIC PLATE IMPACT RESPONSE USING FRACTIONAL DERIVATIVE CONSTITUTIVE EQUATIONS

A GENERALIZATION OF A CONJECTURE OF MELHAM. 1. Introduction The Fibonomial coefficient is, for n m 1, defined by

r ˆr F = Section 2: Newton s Law of Gravitation m 2 m 1 Consider two masses and, separated by distance Gravitational force on due to is

Static equilibrium requires a balance of forces and a balance of moments.

Chapter 8: Generalization and Function Approximation

Transcription:

A New I -Based Hypeelastic Model fo Rubbe Elastic Mateials Osca Lopez-Paies SES Octobe -4, Evanston, IL

Neo-Hookean odel W ìï ï ( I - ) = ( if l + l + l - ) lll = = í ï + othewise ïî (*). Matheatical siplicity (aenable to analytical solutions fo fundaental BVP and ipleentation in FEM codes). Physical intepetation (Gaussian statistics ). Good ageeent with expeients (fo soe elastoes and up to easonably lage defoations) Eve since its deivation in the 94 s, thee have been nueous (oe than 4!!!) efineents of the Neo-Hookean odel (*). Many of the do not confo with featues,, and/o Vahapoglu & Kaadeniz ()

Poposed constitutive odel A oe appopiate basic easue of stain -a a a j( I ; a) = ( I - ), a Î a Note: i) it genealizes the linea NH easue ii) it lineaizes popely (unlike ( I - ) a ) iii) it neglects I fo atheatical siplicity Poposed odel WI M -a a a ( ) = å ( I - ) = a, a eal-valued ateial paaetes j ( I ; ) = ( I - ) M tes to be included in the suation Lopez-Paies ()

. Matheatical siplicity The esulting Cauchy stess is siply given by M W æ ö -a a- T T = - pi = I å FF - pi F ç è = ø with pincipal stesses æ M ö -a a- t = I i ç å l - p ( i =,, ) i çè = ø And the inceental tangent odulus by M M W æ ö æ ö -a a- -a a- L = = ( a - ) I FÄF+ I F çå è çå = ø è= ø These closed-fo quantities ae needed fo the solution of BVPs, and fo the ipleentation in nueical codes such as ABAQUS.

. Physical significance of ateial paaetes Non-Gaussian statistical echanics odel of Beatty () T T =Y( I )FF -pi Hee, Y( I ) contains infoation about the statistical distibution of the undelying polyeic chains Exaple: Auda-Boyce odel with enegy Y = AB nkt W I C I ( ) = ( ) AB å - = - ( I / h ) nkt I / h - Note: By choosing = nktc and a = WI educes to the AB odel M -a a a ( ) = å ( I - ) = a

. Pedictive capabilities The poposed enegy adits the polynoial epesentation: -i é M æi- ö ù WI ( ) = å ( a j) - ( I- ) i= i êåç ú! = è ë j= ø û Thus, in the liit of sall defoations it educes to a - M M WI ( ) = å ( I- ) + å ( I- ) + O( I-) = = å = i ( ) In the sall defoation egie: it lineaizes popely and has no dependence on a with M = In the odeate defoation egie: polynoial dependence on ( I - ) In the lage defoation egie: stong dependence on a povides the eans to captue the typical lock-up

Copaisons with elastoes Fo deonstation puposes, conside the two-te odel -a -a a a a a WI ( ) = ( I - ) + ( I - ) a a Note: the ateial paaetes, a,, a will be deteined fo uniaxial data only by least-squae fitting and by enfocing the pope lineaization condition + = Fo copaison puposes, we also conside the esponse of the Gent odel J I W ( I ) ln é - ù =- - G ê J ú ë û whee and J ae ateial paaetes

Copaisons with a vulcanized ubbe S un (MPa) 6 5 4 Uniaxial Tension Two - Te Model Gent Model Data S (MPa) 4.5.5.5.5 Biaxial Tension Pue Shea Tension l = Two - Te Model Gent Model Shea Data Biaxial Data 4 5 6 7 8 l 4 5 6 l Teloa (944)

Copaisons with a silicone ubbe Uniaxial Tension & Copession Biaxial Tension.5.6.4. Two - Te Model Gent Model Data S un (MPa) -.5 S bi (MPa).8.6 - Two - Te Model Gent Model Data.5.5 l.4...4.6.8. l Meunie et al. (8)

Copaisons with a silicone ubbe Pue Shea Copession l = Pue Shea Tension l = -.4.8 Two - Te Model Gent Model Data S ps (MPa) -.8 -. S ps (MPa).6.4 -.6 Two - Te Model Gent Model Data.5.6.7.8.9 l...4.6.8. l Meunie et al. (8)

Copaisons with a Michelin elastoe Uniaxial Tension Siple Shea 4.5.5 S un (MPa).5.5 S ss (MPa).5 Two - Te Model Gent Model Data....4.5.6.7.8 l.5 Two - Te Model Gent Model Data....4.5.6.7.8 g Lahellec et al. (4)

Constitutive estiction on and Necessay and sufficient conditions fo stict polyconvexity a W '( I ) > and W '( I ) + IW ''( I ) > Necessay and sufficient conditions fo stong ellipticity W '( I ) > and ( l l - ) W '( I ) + I - - W ''( I ) > ( i =,, ) i i Sufficient conditions fo stict polyconvexity (and theefoe fo stong ellipticity) in tes of the ateial paaetes > and a > ( =,,..., M)

Cavitation instabilities Undefoed FE odel l Mesh nea cavity l l x 5 l l l = = Soe specifics: 64,8 8-node bick eleents f 5x 4 Initial volue faction of cavity: f = p/6».5-9 -9 f Cavitation ensues wheneve f Nakaua & Lopez-Paies () 5 = f.5.5.5 s/

FEM Results t > t t= > t t= t f / f 5 unde fixed f / f = 5 Nea Cavity

Onset-of-cavitation suface: Michelin Rubbe Axisyetic Loading.5..5 s.5..5. s t/ - - - - - - t /.5 t = t = t - - t / Hydostatic stess: Shea stesses: t s t + t + t = t -t t -t =, t = Nakaua & Lopez-Paies ()

Genealizations The poposed incopessible, isotopic odel constitutes a pactical platfo fo which to account fo oe levels of coplexity to odel elastoes Although weakly, W is expected to depend on the second invaiant I. At the expense of sacificing atheatical siplicity it ay be of inteest to conside M -a N -bs a a bs bs = å - + u å - s = a s= bs WI (, I) ( I ) ( I ) Copessibility effects ay be eadily added in a nube of diffeent ways, e.g., M -a M a a ' WI J I J J å å = a = (, ) = ( - )- ln + ( -) Othe effects that can be added include Mullins effect, hysteesis, as well as ate and theal effects

Final eaks We have poposed a new constitutive odel fo ubbe like ateials that: i) is atheatically siple, ii) contains paaetes that ay be given a physical significance, and iii) chaacteizes and pedicts accuately the esponse of a vaiety of elastoes. In view of its functional siplicity, it is staightfowad to ipleent it in coecial finite eleent packages (e.g., ABAQUS) fo the study of stuctual pobles. In addition, the poposed I -based odel peits to eadily check fo conditions of polyconvexity and stong ellipticity, needed to undestand the developent of instabilities