finite element analysis of human bone modelsbindel/present/2008-04-courant.pdf · bone basics bone...
TRANSCRIPT
![Page 1: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/1.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Finite Element Analysis of Human BoneModels
D. Bindel
Courant Institute for Mathematical SciencesNew York University
Biomath Lunch Seminar, 1 Apr 2008
![Page 2: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/2.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Outline
1 Bone basics
2 Bone measurement and modeling
3 BoneFEA software
4 Conclusion
![Page 3: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/3.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Outline
1 Bone basics
2 Bone measurement and modeling
3 BoneFEA software
4 Conclusion
![Page 4: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/4.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Why study bones?
Osteoporosis: 44M Americans, $17B / year> 55% of over 50 have osteoporosis or low bone mass350K hip fractures / year; over $10B / yearA quarter of hip fracture patients die within a year... and we’re getting older
![Page 5: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/5.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: macrostructure
![Page 6: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/6.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: microstructure
![Page 7: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/7.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: microstructure
![Page 8: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/8.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: trabecular microstructure
![Page 9: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/9.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: trabecular microstructure
(Scans from 23 and 85 year old females)
![Page 10: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/10.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone basics: orientation and remodeling
![Page 11: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/11.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Why study bones?
... because bone is a fascinating material!
Structurally complicated across length scalesStructure adapts to loads and changes over timeinhomogeneous, anisotropic, asymmetric, oftennonlinear
![Page 12: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/12.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Outline
1 Bone basics
2 Bone measurement and modeling
3 BoneFEA software
4 Conclusion
![Page 13: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/13.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Bone measurement
Diagnostic for osteoporosis: T-scores from DXAOrdinary microscopy on extracted coresQCT software: density profile, about 3 mm scaleMicro-CT and micro-MRI: O(10 micron)
![Page 14: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/14.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Micro-FE bone modeling
One vertebrate = 57M+ elements at 40 microns
![Page 15: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/15.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Whole bone modeling
Density only weakly predicts strengthWanted: Good effective constitutive relation
![Page 16: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/16.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Difficulties
Bone is:Variable over time and between individualsInhomogeneous and anisotropicDifferent in tension and compression
![Page 17: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/17.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Yielding and nonlinearity
Example difficulty:Trabecular network has beam and plate elementsSmall macro strains yield much larger micro strainsSmall-scale geometric nonlinearity a significant effect
![Page 18: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/18.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Yielding and nonlinearity
![Page 19: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/19.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
An approach
Micro-CT structure scans for orientationUse orientation indices + density to approximatematerial parametersProceed phenomenologically
![Page 20: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/20.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Outline
1 Bone basics
2 Bone measurement and modeling
3 BoneFEA software
4 Conclusion
![Page 21: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/21.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Diagnostic toolchain
Micro-CT scan data from patientInference of material propertiesConstruction of coarse FE model (voxels)Simulation under loadingOutput of stress fields, displacements, etc.
![Page 22: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/22.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
BoneFEA
Standard displacement-based finite element codeElastic and plastic material models (includinganisotropy and asymmetric yield surfaces)High-level: incremental load control loop,Newton-Krylov solvers with line search for nonlinearsystemsLibrary of (fairly simple) preconditioners; default is atwo-level geometric multigrid preconditioner
![Page 23: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/23.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Example analysis loop
mesh:rigid(mesh:numnp()-1, {z=’min’},function()return ’uuuuuu’, 0, 0, bound_disp
end)
pc = simple_msm_pc(mesh,20)mesh:set_cg{M=pc, tol=1e-6, max_iter=1000}for j=1,n dobound_disp = 0.2*jmesh:step()mesh:newton{max_iter=6, Rtol=1e-4}
end
![Page 24: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/24.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Example analyses
![Page 25: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/25.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Outline
1 Bone basics
2 Bone measurement and modeling
3 BoneFEA software
4 Conclusion
![Page 26: Finite Element Analysis of Human Bone Modelsbindel/present/2008-04-courant.pdf · Bone basics Bone measurement and modeling BoneFEA software Conclusion Finite Element Analysis of](https://reader034.vdocument.in/reader034/viewer/2022042309/5ed590c7bb43ae0d0e021b2e/html5/thumbnails/26.jpg)
Courant 08
Bone basics
Bonemeasurementand modeling
BoneFEAsoftware
Conclusion
Conclusion
Bones are interesting as well as important!Initial BoneFEA work done, in use by ON DiagnosticsPossible follow-up work for diagnostic toolPlenty of interesting research directions