1 / 23

2D Deformation and Creep Response of Articular Cartilage

2D Deformation and Creep Response of Articular Cartilage. By: Mikhail Yakhnis & Robert Zhang. Motivation. Articular cartilage transfers load between bones enables smooth motion along joints Cartilage has limited capacity for self repair

gerry
Télécharger la présentation

2D Deformation and Creep Response of Articular Cartilage

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 2D Deformation and Creep Response of Articular Cartilage By: Mikhail Yakhnis & Robert Zhang

  2. Motivation • Articular cartilage • transfers load between bones • enables smooth motion along joints • Cartilage has limited capacity for self repair • Applications: biomaterials, prosthetics, biomedical devices http://nigelhartnett.onlinemedical.com.au/images/articular%20knee%20injury.jpg

  3. Problem Description • Consider cartilage in an unconfined compression under constant load F • Analyze the 2D elastic deformation over time F Compression plate Articular Cartilage Frictionless Supports

  4. Material Background • Cartilage often modeled as a viscoelastic material • Viscous and elastic by superposition • Elasticity and viscosity can be linear or nonlinear • Established models: Kelvin-Voigt, Maxwell, Standard-Linear Solid http://www.allsealsinc.com/allseals/Orings/maxwell.gif

  5. Mathematical Model for Cartilage • We chose the Kelvin-Voigt model to focus on the creep response • The constitutive equation is Mechanical Analogue of Kelvin-Voigt Model http://en.wikipedia.org/wiki/File:Kelvin_Voigt_diagram.svg

  6. Assumptions for Model • Conditions • Constant force Fnormal to boundary B3 • No gravity (body force) • 2D, plane stress* • Confined in y-direction along B1andB3 • Confined in x-direction along B4 • Properties • c = 0.1m; L = 0.125m • Constant cross-sectional area A • Isotropic elasticity* * L c F B3 B4 B2 y B1 x

  7. Experimental Data • (Aggregate Modulus) Data Book on Mechanical Properties of Living Cells, Tissues, and Organs /. Tokyo ; New York : Springer, 1996. Print.

  8. Derivation of Weak Form • By definition, stress • Strain can be rewritten as gradient of displacement u • Our constitutive equation (in strong form) becomes

  9. Derivation of Weak Form (1) Take the gradient of the force equation (which equals zero) (2) Multiply by an arbitrary displacement (3) Integrate by parts to induce symmetry of and

  10. Decoupling a Transient Problem • We can decouple the formulation and assume the time and spatial variations are separate where is a function of time only and basis function N is function of space • The weak differential equation rewritten in matrix form is Reddy, J. N.. "Time-Dependent Problems." An introduction to nonlinear finite element analysis. Oxford: Oxford University Press, 2004. . Print.

  11. Displacement Equation for Creep Response • At each time step n • The equation for becomes

  12. Modeling Creep in MATLAB • Changes in Preprocessor.m • Provide initial displacement • Define time step • Adjust boundary conditions • Changes in Assemble.m • Assemble the damping matrix [C] • Changes in NodalSoln.m • Add initial condition, damping, time inputs • Modify reaction force and displacement equations

  13. Modeling Creep in MATLAB • Discussion: • MATLAB result converges toward experimental data farther away from initial time • 10% error at 6 seconds • MATLAB model reaches equilibrium faster than experimental data

  14. Modeling Creep in MATLAB

  15. Modeling Creep in ANSYS • A variety of models are available • Differences include suitability for primary and secondary creep • Usually of the form • Examples • Strain Hardening: • Time Hardening: ANSYS Advanced Nonlinear Materials: Lecture 3 – Rate Dependent Creep http://www.ansys-blog.com/wp-content/uploads/2012/06/Three-Types-of-Creep.png

  16. Considerations for ANSYS Model • What experimental data is available to us? • Can we fit the experimental data to the model? • Can we use the built-in Mechanical APDL curve fitting procedure? • Is there more emphasis on primary creep or secondary creep? • Does the model satisfy our constitutive equation?

  17. Parameters in the ANSYS Model • Experimental data provides aggregate modulus and Poisson’s ratio • Young’s Modulus can be derived from • The solution for time-dependent strain in the K-V model is • We can use the Modified Exponential Function in ANSYS where ; we can solve for and ANSYS Advanced Nonlinear Materials: Lecture 3 – Rate Dependent Creep

  18. ANSYS Results – Creep Response Short Term Response – 30 Seconds Long Term Response – 3000 Seconds

  19. Animation of Deformation in ANSYS

  20. Comparison of ANSYS and Experiment • Result: • Theoretical Model-Based ANSYS data tends to overshoot experimental data • Error is between 30% to 40% per data point • Experimental-based model performs better • Discussion: • Results demonstrate the limitations of ANSYS models • A combined primary-secondary model is ideal • Long term response in ANSYS is not accurate • Function models primary response Primary + Secondary Time Hardening

  21. ANSYS Model: Mesh and Time Refinement

  22. Sensitivity Analysis • Recall the creep model: • We varied each non-zero model constant by 50%* to perform a rudimentary sensitivity analysis: *The simulation did not converge at C2 +50% so C2 +10% was used instead

  23. 2D Deformation and Creep Response of Articular Cartilage By: DJ Mikey Mike & Big Rob Zhang Thank you for listening. Questions?

More Related