In the field of computational biomechanics, investigators have primarily used commercial software that is neither geared toward biological applications nor sufficiently flexible to follow the latest developments in the field. This lack of a tailored software environment has hampered research progress, as well as dissemination of models and results. To address these issues, we developed the FEBio software suite (http://mrl.sci.utah.edu/software/febio), a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. This paper provides an overview of the theoretical basis of FEBio and its main features. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures. Software verification is a large part of the development and maintenance of FEBio, and to demonstrate the general approach, the description and results of several problems from the FEBio Verification Suite are presented and compared to analytical solutions or results from other established and verified FE codes. An additional simulation is described that illustrates the application of FEBio to a research problem in biomechanics. Together with the pre- and postprocessing software PREVIEW and POSTVIEW, FEBio provides a tailored solution for research and development in computational biomechanics.
Issue Section:
Research Papers
References
1.
Anderson
, A. E.
, Peters
, C. L.
, Tuttle
, B. D.
, and Weiss
, J. A.
, 2005, “Subject-Specific Finite Element Model of the Pelvis: Development, Validation and Sensitivity Studies
,” J. Biomech. Eng.
, 127
, pp. 364
–373
.2.
Spilker
, R. L.
, Feinstein
, J. A.
, Parker
, D. W.
, Reddy
, V. M.
, and Taylor
, C. A.
, 2007, “Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries
,” Ann. Biomed. Eng.
, 35
, pp. 546
–559
.3.
Wang
, L.
, Zhang
, H.
, Shi
, P.
, and Liu
, H.
, 2006, “Imaging of 3D Cardiac Electrical Activity: A Model-Based Recovery Framework
,” Med. Image Comput. Comput. Assist. Interv.
, 9
, pp. 792
–799
.4.
Speelman
, L.
, Bohra
, A.
, Bosboom
, E. M.
, Schurink
, G. W.
, van de Vosse
, F. N.
, Makaorun
, M. S.
, and Vorp
, D. A.
, 2007, “Effects of Wall Calcifications in Patient-Specific Wall Stress Analyses of Abdominal Aortic Aneurysms
,” J. Biomech. Eng.
, 129
, pp. 105
–109
.5.
Portnoy
, S.
, Yarnitzky
, G.
, Yizhar
, Z.
, Kristal
, A.
, Oppenheim
, U.
, Siev-Ner
, I.
, and Gefen
, A.
, 2007, “Real-Time Patient-Specific Finite Element Analysis of Internal Stresses in the Soft Tissues of a Residual Limb: A New Tool for Prosthetic Fitting
,” Ann. Biomed. Eng.
, 35
, pp. 120
–135
.6.
Schmid-Schonbein
, G. W.
, and Diller
, K. R.
, 2005, “Transport Processes in Biomedical Systems: A Roadmap for Future Research Directions
,” Ann. Biomed. Eng.
, 33
, pp. 1136
–1141
.7.
Crawford
, R. P.
, Rosenberg
, W. S.
, and Keaveny
, T. M.
, 2003, “Quantitative Computed Tomography-Based Finite Element Models of the Human Lumbar Vertebral Body: Effect of Element Size on Stiffness, Damage, and Fracture Strength Predictions
,” J. Biomech. Eng.
, 125
, pp. 434
–438
.8.
Li
, L. P.
, and Herzog
, W.
, 2006, “Arthroscopic Evaluation of Cartilage Degeneration Using Indentation Testing—Influence of Indenter Geometry
,” Clin. Biomech. (Bristol, Avon)
, 21
, pp. 420
–426
.9.
Li
, L. P.
, and Herzog
, W.
, 2005, “Electromechanical Response of Articular Cartilage in Indentation—Considerations on the Determination of Cartilage Properties During Arthroscopy
,” Comput. Methods Biomech. Biomed. Eng.
, 8
, pp. 83
–91
.10.
Davids
, N.
, and Mani
, M. K.
, 1974, “A Finite Element Analysis of Endothelial Shear Stress for Pulsatile Blood Flow
,” Biorheology
, 11
, pp. 137
–147
. Available at http://www.ncbi.nlm.nih.gov/pubmed/4441640http://www.ncbi.nlm.nih.gov/pubmed/4441640.11.
Doyle
, J. M.
, and Dobrin
, P. B.
, 1971, “Finite Deformation Analysis of the Relaxed and Contracted Dog Carotid Artery
,” Microvasc. Res.
, 3
, pp. 400
–415
.12.
Janz
, R. F.
, and Grimm
, A. F.
, 1972, “Finite-Element Model for the Mechanical Behavior of the Left Ventricle. Prediction of Deformation in the Potassium-Arrested Rat Heart
,” Circ. Res.
, 30
, pp. 244
–252
.13.
Matthews
, F. L.
, and West
, J. B.
, 1972, “Finite Element Displacement Analysis of a Lung
,” J. Biomech.
, 5
, pp. 591
–600
.14.
Farah
, J. W.
, Craig
, R. G.
, and Sikarskie
, D. L.
, 1973, “Photoelastic and Finite Element Stress Analysis of a Restored Axisymmetric First Molar
,” J. Biomech.
, 6
, pp. 511
–520
.15.
Belytschko
, T.
, Kulak
, R. F.
, Schultz
, A. B.
, and Galante
, J. O.
, 1974, “Finite Element Stress Analysis of an Intervertebral Disc
,” J. Biomech.
, 7
, pp. 277
–285
.16.
Anderson
, A. E.
, Ellis
, B. J.
, and Weiss
, J. A.
, 2007, “Verification, Validation and Sensitivity Studies in Computational Biomechanics
,” Comput. Methods Biomech. Biomed. Eng.
, 10
, pp. 171
–184
.17.
Henninger
, H. B.
, Reese
, S. P.
, Anderson
, A. E.
, and Weiss
, J. A.
, 2010, “Validation of Computational Models in Biomechanics
,” Proc. Inst. Mech. Eng. Part H, J. Eng. Med.
, 224
, pp. 801
–812
.18.
Maas
, S. A.
, and Weiss
, J. A.
, 2007, FEBio User’s Manual, http://mrl.sci.utah.edu/software/febiohttp://mrl.sci.utah.edu/software/febio.19.
ASME, 2006, “
Guide for Verification and Validation in Computational Solid Mechanics
,” ASME Committee (PT60) on Verification and Validation in Computational Solid Mechanics.20.
Babuska
, I.
, and Oden
, J. T.
, 2004, “Verification and Validation in Computational Engineering and Science: Basic Concepts
,” Comput. Methods Appl. Mech. Eng.
, 193
, pp. 4057
–4066
.21.
Maas
, S. A.
, and Weiss
, J. A.
, 2007, FEBio Theory Manual, http://mrl.sci.utah.edu/software/febiohttp://mrl.sci.utah.edu/software/febio.22.
Bonet
, J.
, and Wood
, R. D.
, 1997, Nonlinear Continuum Mechanics for Finite Element Analysis
, Cambridge University Press
, Cambridge, NY
.23.
Matthies
, H.
, and Strang
, G.
, 1979, “The Solution of Nonlinear Finite Element Equations
,” Int. J. Num. Methods Eng.
, 14
, pp. 1613
–1626
.24.
Bonet
, J.
, and Wood
, R. D.
, 1997, Nonlinear Continuum Mechanics for Finite Element Analysis
, Cambridge University Press
, Cambridge, NY
.25.
Simo
, J. C.
, and Taylor
, R. L.
, 1991, “Quasi-Incompressible Finite Elasticity in Principal Stretches: Continuum Basis and Numerical Algorithms
,” Comput. Methods Appl. Mech. Eng.
, 85
, pp. 273
–310
.26.
Puso
, M. A.
, and Solberg
, J.
, 2006, “A Stabilized Nodally Integrated Tetrahedral
,” Int. J. Num. Methods Eng.
, 67
, pp. 841
–867
.27.
Betsch
, P.
, and Gruttmann
, F.
, 1996, “A 4-Node Finite Shell Element for the Implementation of General Hyperelastic 3D-Elasticity at Finite Strains
,” Comput. Methods Appl. Mech. Eng.
, 130
, pp. 57
–79
.28.
Mow
, V. C.
, Kuei
, S. C.
, Lai
, W. M.
, and Armstrong
, C. G.
, 1980, “Biphasic Creep and Stress Relaxation of Articular Cartilage in Compression: Theory and Experiments
,” J. Biomech. Eng.
, 102
, pp. 73
–84
.29.
Ateshian
, G. A.
, Ellis
, B. J.
, and Weiss
, J. A.
, 2007, “Equivalence Between Short-Time Biphasic and Incompressible Elastic Material Response
,” J. Biomech. Eng.
129
(3
), pp. 405
–412
.30.
Weiss
, J. A.
, Maker
, B. N.
, and Govindjee
, S.
, 1996, “Finite Element Implementation of Incompressible, Transversely Isotropic Hyperelasticity
,” Comput. Methods Appl. Mech. Eng.
, 135
, pp. 107
–128
.31.
Ateshian
, G. A.
, Maas
, S.
, and Weiss
, J. A.
, 2010, “Finite Element Algorithm for Frictionless Contact of Porous Permeable Media Under Finite Deformation and Sliding
,” J. Biomech. Eng.
, 132
, p. 061006
.32.
Laursen
, T. A.
, and Maker
, B. N.
, 1995, “Augmented Lagrangian Quasi-Newton Solver for Constrained Nonlinear Finite Element Applications
,” Int. J. Numer. Methods Eng.
, 38
, pp. 3571
–3590
.33.
Maker
, B. N.
, 1995, “Rigid Bodies for Metal Forming Analysis With NIKE3D
,” University of California, Lawrence Livermore Laboratory, Report No. UCRL-JC-119862, pp. 1
–8
.34.
Schenk
, O.
, and Gartner
, K.
, 2004, “Solving Unsymmetric Sparse Systems of Linear Equations With PARDISO
,” J. Future Gen. Comput. Syst.
, 20
, pp. 475
–487
.35.
Schenk
, O.
, and Gartner
, K.
, 2006, “On Fast Factorization Pivoting Methods for Symmetric Indefinite Systems
,” Elect. Trans. Numer. Anal.
, 23
, pp. 158
–179
. Available at http://etna.mcs.kent.edu/http://etna.mcs.kent.edu/.36.
Demmel
, J. W.
, Eisenstat
, S. C.
, Gilbert
, J. R.
, Li
, X. S.
, and Liu
, J. W. H.
, 1999, “A Supernodal Approach to Sparse Partial Pivoting
,” SIAM J. Matrix Anal. Appl.
, 20
, pp. 720
–755
.37.
Gupta
, A.
, “WSMP: Watson Sparse Matrix Package
,” IBM Research Report No. RC 21888 (98472), 2000.38.
Veldhuizen
, T.
, Jernigan
, M.
, Ishikawa
, Y.
, Oldehoeft
, R.
, Reynders
, J.
, and Tholburn
, M.
, 1997, “Will C++ be Faster than FORTRAN? Scientific Computing in Object-Oriented Parallel Environments
,” Lect. Notes Comput. Sci.
, 1343
, pp. 49
–56
.39.
Maker
, B. N.
, 1995, “NIKE3D: A Nonlinear, Implicit, Three-Dimensional Finite Element Code for Solid and Structural Mechanics
,” Lawrence Livermore Laboratory Technical Report No. UCRL-MA-105268.40.
Park
, S.
, Krishnan
, R.
, Nicoll
, S. B.
, and Ateshian
, G. A.
, 2003, “Cartilage Interstitial Fluid Load Support in Unconfined Compression
,” J. Biomech.
, 36
, pp. 1785
–1796
.41.
Ogden
, R. W.
, 1984, Non-Linear Elastic Deformations
, Ellis Horwood
, Chichester, UK
.42.
Erdemir
, A.
, Viveiros
, M. L.
, Ulbrecht
, J. S.
, and Cavanagh
, P. R.
, 2006, “An Inverse Finite-Element Model of Heel-Pad Indentation
,” J. Biomech.
, 39
, pp. 1279
–1286
.43.
Ateshian
, G. A.
, Warden
, W. H.
, Kim
, J. J.
, Grelsamer
, R. P.
, and Mow
, V. C.
, 1997, “Finite Deformation Biphasic Material Properties of Bovine Articular Cartilage From Confined Compression Experiments
,” J. Biomech.
, 30
, pp. 1157
–1164
.44.
Holmes
, M. H.
, 1986, “Finite Deformation of Soft Tissue: Analysis of a Mixture Model in Uni-Axial Compression
,” J. Biomech. Eng.
, 108
, pp. 372
–381
.45.
Armstrong
, C. G.
, Lai
, W. M.
, and Mow
, V. C.
, 1984, “An Analysis of the Unconfined Compression of Articular Cartilage
,” J. Biomech. Eng.
, 106
, pp. 165
–173
.46.
Oden
, J. T.
, 1972, Finite Elements of Nonlinear Continua
, McGraw-Hill
, New York
.47.
Treloar
, L. R. G.
, 1944, “Stress-Strain Data for Vulcanised Rubber Under Various Types of Deformation
,” Trans. Faraday Soc.
, 40
, pp. 59
–70
.48.
Bisshopp
, R. E.
, and D. C.
Drucker
, 1945, “Large Deflection of Cantilever Beams
,” Q. Appl. Math.
, 3
, pp. 272
–275
.49.
Batoz
, J. L.
, 1982, “An Explicit Formulation for an Efficient Triangular Plate Bending Element
,” Int. J. Numer. Methods Eng.
, 18
, pp. 1077
–1089
.50.
Moore
, S. M.
, Ellis
, B. J.
, Weiss
, J. A.
, McMahon
, P. J.
, and Debski
, R. E.
, 2010, “The Glenohumeral Capsule Should be Evaluated as a Sheet of Fibrous Tissue: A Validated Finite Element Model
,” Ann. Biomed. Eng.
, 38
(1
), pp. 66
–76
.51.
Drury
, N. J.
, Ellis
, B. J.
, Weiss
, J. A.
, McMahon
, P. J.
, and Debski
, R. E.
, 2011, “Finding Consistent Strain Distributions in the Glenohumeral Capsule Between Two Subjects: Implications for Physical Examinations
,” J. Biomech.
, 44
, pp. 607
–613
.Copyright © 2012
by American Society of Mechanical Engineers
You do not currently have access to this content.
