Rib fractures under anterior–posterior dynamic loads: Experimental and finite-element study

The purpose of this study was to investigate whether using a finite-element (FE) mesh composed entirely of hexahedral elements to model cortical and trabecular bone (all-hex model) would provide more accurate simulations than those with variable thickness shell elements for cortical bone and hexahedral elements for trabecular bone (hex–shell model) in the modeling human ribs. First, quasi-static non-injurious and dynamic injurious experiments were performed using the second, fourth, and tenth human thoracic ribs to record the structural behavior and fracture tolerance of individual ribs under anterior–posterior bending loads. Then, all-hex and hex–shell FE models for the three ribs were developed using an octree-based and multi-block hex meshing approach, respectively. Material properties of cortical bone were optimized using dynamic experimental data and the hex–shell model of the fourth rib and trabecular bone properties were taken from the literature. Overall, the reaction force–displacement relationship predicted by both all-hex and hex–shell models with nodes in the offset middle-cortical surfaces compared well with those measured experimentally for all the three ribs. With the exception of fracture locations, the predictions from all-hex and offset hex–shell models of the second and fourth ribs agreed better with experimental data than those from the tenth rib models in terms of reaction force at fracture (difference <15.4%), ultimate failure displacement and time (difference <7.3%), and cortical bone strains. The hex–shell models with shell nodes in outer cortical surfaces increased static reaction forces up to 16.6%, compared to offset hex–shell models. These results indicated that both all-hex and hex–shell modeling strategies were applicable for simulating rib responses and bone fractures for the loading conditions considered, but coarse hex–shell models with constant or variable shell thickness were more computationally efficient and therefore preferred.     

A New Computational Efficient Approach for Trabecular Bone Analysis using Beam Models Generated with Skeletonized Graph Technique

Micro-finite element (FE) analysis is a well established technique for the evaluation of the elastic properties of trabecular bone, but is limited in its application due to the large number of elements that it requires to represent the complex internal structure of the bone. In this paper, we present an alternative FE approach that makes use of a recently developed 3D-Line Skeleton Graph Analysis (LSGA) technique to represent the complex internal structure of trabecular bone as a network of simple straight beam elements in which the beams are assigned geometrical properties of the trabeculae that they represent. Since an enormous reduction of cputime can be obtained with this beam modeling approach, ranging from approximately 1,200 to 3,600 for the problems investigated here, we think that the FE modeling technique that we introduced could potentially constitute an interesting alternative for the evaluation of the elastic mechanical properties of trabecular bone.     

A new computational efficient approach for trabecular bone analysis using beam models generated with skeletonized graph technique

Micro-finite element (FE) analysis is a well established technique for the evaluation of the elastic properties of trabecular bone, but is limited in its application due to the large number of elements that it requires to represent the complex internal structure of the bone. In this paper, we present an alternative FE approach that makes use of a recently developed 3D-Line Skeleton Graph Analysis (LSGA) technique to represent the complex internal structure of trabecular bone as a network of simple straight beam elements in which the beams are assigned geometrical properties of the trabeculae that they represent. Since an enormous reduction of cputime can be obtained with this beam modeling approach, ranging from approximately 1,200 to 3,600 for the problems investigated here, we think that the FE modeling technique that we introduced could potentially constitute an interesting alternative for the evaluation of the elastic mechanical properties of trabecular bone.     

Adaptive glenoid bone remodeling simulation

Glenoid prosthesis loosening is the most common cause for revision total shoulder arthroplasty. Stress-induced bone remodeling may compromise long-term prosthesis fixation and significantly contribute to loosening. Realistic, robust analysis of bone-prosthesis constructs need to look beyond initial post-implantation mechanics provided by static finite element (FE) simulation. Adaptive bone remodeling simulations based on Wolff's law are needed for evaluating long-term glenoid prostheses fixation. The purpose of this study was to take a first step towards this goal and create and validate two-dimensional FE simulations, using the intact glenoid, for computing subject-specific adaptive glenoid remodeling. Two-dimensional glenoid FE models were created from scapulae computed tomography images. Two distinct processes, “element” and “node” simulations, used the forward-Euler method to compute bone remodeling. Initial bone density was homogeneous. Center and offset load combinations were iteratively applied. To validate the simulations we performed location-specific statistical comparisons between predicted and actual bone density, load combinations, and “element” and “node” processes. Visually and quantitatively “element” simulations produced better results (p>0.22), and correlation coefficients ranged 0.51–0.69 (p<0.001). Having met this initial work's goals, we expect subject-specific FE glenoid bone remodeling simulations together with static FE stress analyses to be effective tools for designing and evaluating glenoid prostheses.     

Simplified boundary conditions alter cortical-trabecular load sharing at the distal radius; A multiscale finite element analysis

High-resolution peripheral quantitative computed tomography (HR-pQCT) derived micro-finite element (FE) modeling is used to evaluate mechanical behavior at the distal radius microstructure. However, these analyses typically simulate non-physiologic simplified platen-compression boundary conditions on a small section of the distal radius. Cortical and trabecular regions contribute uniquely to distal radius mechanical behavior, and various factors affect these regions distinctly. Generalized strength predictions from standardized platen-compression analyses may not adequately capture region specific responses in bone. Our goal was to compare load sharing within the cortical-trabecular compartments between the standardized platen-compression BC simulations, and physiologic BC simulations using a validated multiscale approach. Clinical- and high-resolution images were acquired from nine cadaveric forearm specimens using an HR-pQCT scanner. Multiscale FE models simulating physiologic BCs, and micro-FE only models simulating platen-compression BCs were created for each specimen. Cortical and trabecular loads (N) along the length of the distal radius micro-FE section were compared between BCs using correlations. Principal strain distributions were also compared quantitatively. Cortical and trabecular loads from the platen-compression BC simulations were strongly correlated to the physiologic BC simulations. However, a 30% difference in cortical loads distally, and a 53% difference in trabecular loads proximally was observed under platen BC simulations. Also, distribution of principal strains was clearly different. Our data indicated that platen-compression BC simulations alter cortical-trabecular load sharing. Therefore, results from these analyses should be interpreted in the appropriate mechanical context for clinical evaluations of normal and pathologic mechanical behavior at the distal radius.     

Concept and development of an orthotropic FE model of the proximal femur

PURPOSE: In contrast to many isotropic finite-element (FE) models of the femur in literature, it was the object of our study to develop an orthotropic FE "model femur" to realistically simulate three-dimensional bone remodelling. METHODS: The three-dimensional geometry of the proximal femur was reconstructed by CT scans of a pair of cadaveric femurs at equal distances of 2mm. These three-dimensional CT models were implemented into an FE simulation tool. Well-known "density-determined" bony material properties (Young's modulus; Poisson's ratio; ultimate strength in pressure, tension and torsion; shear modulus) were assigned to each FE of the same "CT-density-characterized" volumetric group.In order to fix the principal directions of stiffness in FE areas with the same "density characterization", the cadaveric femurs were cut in 2mm slices in frontal (left femur) and sagittal plane (right femur). Each femoral slice was scanned into a computer-based image processing system. On these images, the principal directions of stiffness of cancellous and cortical bone were determined manually using the orientation of the trabecular structures and the Haversian system. Finally, these geometric data were matched with the "CT-density characterized" three-dimensional femur model. In addition, the time and density-dependent adaptive behaviour of bone remodelling was taken into account by implementation of Carter's criterion. RESULTS: In the constructed "model femur", each FE is characterized by the principal directions of the stiffness and the "CT-density-determined" material properties of cortical and cancellous bone. Thus, on the basis of anatomic data a three-dimensional FE simulation reference model of the proximal femur was realized considering orthotropic conditions of bone behaviour. CONCLUSIONS: With the orthotropic "model femur", the fundamental basis has been formed to realize realistic simulations of the dynamical processes of bone remodelling under different loading conditions or operative procedures (osteotomies, total hip replacements, etc).     

Experimentally tested computer modeling of stress fractures in rats

The objective of this study was to develop a finite-element (FE) modeling methodology for studying the etiology of a stress fracture (SF). Several variants of three-dimensional FE models of a rat hindlimb, which differed in length or stiffness of tissues, enabling the analyses of mechanical strains and stress in the tibia, were created. We compared the occurrence of SFs in an animal model to validate locations of peak strains/stresses in the FE models. Four Sprague-Dawley male rats, age ~7 wk, were subjected to mechanical cyclic loads of 1.2 Hz and ~6 N, which were delivered to their hindlimb for 30 min, 3 times/wk, up to 12 wk, by using a specially designed apparatus. The results showed that 1) FE modeling predicted the maximal strains/stresses (~220,0 με and ~29 MPa, respectively) between the mid- and proximal thirds of the tibia; 2) in a longer shin, greater and more inhomogeneous tensile strains/stresses were evident, at the same location; 3) anatomical variants in shin length influenced the strain/stress distributions to a greater extent with respect to changes in mechanical properties of tissues; and 4) bone stiffness was more dominant than muscle stiffness in affecting the strain/stress distributions. In the animal study, 35,000 loading cycles were associated with the formation of a SF. The location of the identified SF in the rat limb verified the FE model. We find the suggested model a valuable tool in studying various aspects of SFs.     

Finite element modeling of α-helices and tropocollagen molecules referring to spike of SARS-CoV-2

The newly developed finite element (FE) modeling at the atomic scale was used to predict the static and dynamic response of the α-helix (AH) and tropocollagen (TC) protein fragments, the main building blocks of the spike of the SARS-CoV-2. The geometry and morphology of the spike’s stalk and its connection to the viral envelope were determined from the combination of most recent molecular dynamics (MD) simulation and images of cryoelectron microscopy. The stiffness parameters of the covalent bonds in the main chain of the helix were taken from the literature. The AH and TC were modeled using both beam elements (wire model) and shell elements (ribbon model) in FE analysis to predict their mechanical properties under tension. The asymptotic stiffening features of AH and TC under tensile loading were revealed and compared with a new analytical solution. The mechanical stiffnesses under other loading conditions, including compression, torsion, and bending, were also predicted numerically and correlated with the results of the existing MD simulations and tests. The mode shapes and natural frequencies of the spike were predicted using the built FE model. The frequencies were shown to be within the safe range of 1–20 MHz routinely used for medical imaging and diagnosis by means of ultrasound. These results provide a solid theoretical basis for using ultrasound to study damaging coronavirus through transient and resonant vibration at large deformations.     

Differences in 3D vs. 2D analysis in lumbar spinal fusion simulations

Lumbar interbody fusion is currently the gold standard in treating patients with disc degeneration or segmental instability. Despite it having been used for several decades, the non-union rate remains high. A failed fusion is frequently attributed to an inadequate mechanical environment after instrumentation. Finite element (FE) models can provide insights into the mechanics of the fusion process. Previous fusion simulations using FE models showed that the geometries and material of the cage can greatly influence the fusion outcome. However, these studies used axisymmetric models which lacked realistic spinal geometries. Therefore, different modeling approaches were evaluated to understand the bone-formation process.Three FE models of the lumbar motion segment (L4–L5) were developed: 2D, Sym-3D and Nonsym-3D. The fusion process based on existing mechano-regulation algorithms using the FE simulations to evaluate the mechanical environment was then integrated into these models. In addition, the influence of different lordotic angles (5, 10 and 15°) was investigated. The volume of newly formed bone, the axial stiffness of the whole segment and bone distribution inside and surrounding the cage were evaluated.In contrast to the Nonsym-3D, the 2D and Sym-3D models predicted excessive bone formation prior to bridging (peak values with 36 and 9% higher than in equilibrium, respectively). The 3D models predicted a more uniform bone distribution compared to the 2D model.The current results demonstrate the crucial role of the realistic 3D geometry of the lumbar motion segment in predicting bone formation after lumbar spinal fusion.     

A phase-contrast microscopy-based method for modeling the mechanical behavior of mesenchymal stem cells

We present three-dimensional (3D) finite element (FE) models of single, mesenchymal stem cells (MSCs), generated from images obtained by optical phase-contrast microscopy and used to quantify the structural responses of the studied cells to externally applied mechanical loads. Mechanical loading has been shown to affect cell morphology and structure, phenotype, motility and other biological functions. Cells experience mechanical loads naturally, yet under prolonged or sizable loading, damage and cell death may occur, which motivates research regarding the structural behavior of loaded cells. For example, near the weight-bearing boney prominences of the buttocks of immobile persons, tissues may become highly loaded, eventually leading to massive cell death that manifests as pressure ulcers. Cell-specific computational models have previously been developed by our group, allowing simulations of cell deformations under compressive or stretching loads. These models were obtained by reconstructing specific cell structures from series of 2D fluorescence, confocal image-slices, requiring cell-specific fluorescent-staining protocols and costly (confocal) microscopy equipment. Alternative modeling approaches represent cells simply as half-spheres or half-ellipsoids (i.e. idealized geometries), which neglects the curvature details of the cell surfaces associated with changes in concentrations of strains and stresses. Thus, we introduce here for the first time an optical image-based FE modeling, where loads are simulated on reconstructed 3D geometrical cell models from a single 2D, phase-contrast image. Our novel modeling method eliminates the need for confocal imaging and fluorescent staining preparations (both expensive), and makes cell-specific FE modeling affordable and accessible to the biomechanics community. We demonstrate the utility of this cost-effective modeling method by performing simulations of compression of MSCs embedded in a gel.     

Finite element modeling of α-helices and tropocollagen molecules referring to spike of SARS-CoV-2

The newly developed finite element (FE) modeling at the atomic scale was used to predict the static and dynamic response of the α-helix (AH) and tropocollagen (TC) protein fragments, the main building blocks of the spike of the SARS-CoV-2. The geometry and morphology of the spike’s stalk and its connection to the viral envelope were determined from the combination of most recent molecular dynamics (MD) simulation and images of cryoelectron microscopy. The stiffness parameters of the covalent bonds in the main chain of the helix were taken from the literature. The AH and TC were modeled using both beam elements (wire model) and shell elements (ribbon model) in FE analysis to predict their mechanical properties under tension. The asymptotic stiffening features of AH and TC under tensile loading were revealed and compared with a new analytical solution. The mechanical stiffnesses under other loading conditions, including compression, torsion, and bending, were also predicted numerically and correlated with the results of the existing MD simulations and tests. The mode shapes and natural frequencies of the spike were predicted using the built FE model. The frequencies were shown to be within the safe range of 1–20 MHz routinely used for medical imaging and diagnosis by means of ultrasound. These results provide a solid theoretical basis for using ultrasound to study damaging coronavirus through transient and resonant vibration at large deformations.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Load transfer analysis of the distal radius from in-vivo high-resolution CT-imaging.

Prevention of osteoporotic bone fractures requires accurate diagnostic methods to detect the increase in bone fragility at an early stage of osteoporosis. However, today's bone fracture risk prediction, primarily based on bone density measurement, is not sufficiently precise. There is increasing evidence that, in addition to bone density, also the bone microarchitecture and its mechanical loading conditions are important factors determining the fracture risk. Recently, it has been shown that new high-resolution imaging techniques in combination with new computer modeling techniques based on the finite-element (FE) method can account for these additional factors. These techniques might provide information that is more relevant for the prediction of bone fracture risk. So far, however, these new imaged-based FE techniques have not been feasible in-vivo. The objectives of this study were to quantify the load transfer through the trabecular network in a distal radius using a computer model based on in-vivo high-resolution images and to determine if common regions of fractures can be explained as a result of high tissue loading in these regions. The left distal radius and the two adjacent carpal bones of a healthy volunteer were imaged using a high-resolution three-dimensional CT system providing an isotropic resolution of 165 microm. The bone representation was converted into a FE-model that was used to calculate stresses and strains in the trabecular network. The two carpal bones were loaded using different load ratios (for each load case 1000 N in total) representing impact forces on the hand either in near-neutral position or ulnar/radial deviation. The load transfer through the trabecular network of the radius was characterized by the tissue strain energy density (SED) distribution for all load cases. It was found that the distribution of the tissue loading depends on the ratio of the forces acting on the carpal bones. For all load cases the higher SED values (on average: 0.02 +/- 0.08 (S.D.) N mm(-2)) are found in a 10 mm region adjacent to the articular surface which corresponds well with the region where Colles- or Chauffeur-fractures occur. We expect that, eventually, this new approach can lead to a better prediction of the fracture risk than methods based on bone density alone since it accounts for the bone microstructure as well as its loading conditions.     

Towards validation of computational analyses of peri-implant displacements by means of experimentally obtained displacement maps

Micro-finite element (μFE) analysis has recently been introduced for the detailed quantification of the mechanical interaction between bone and implant. The technique has been validated at an apparent level. The aim of this study was to address the accuracy of μFE analysis at the trabecular level. Experimental displacement fields were obtained by deformable image registration, also known as strain mapping (SM), of dynamic hip screws implanted in three human femoral heads. In addition, displacement fields were calculated using μFE analysis. On a voxel-by-voxel basis, the coefficients of determination (R(2)) between experimental and μFE-calculated displacements ranged from 0.67 to 0.92. Linear regression of the mean displacements over nine volumes of interest yielded R(2) between 0.81 and 0.84. The lowest R(2) values were found in regions of very small displacements. In conclusion, we found that peri-implant bone displacements calculated with μFE analysis correlated well with displacements obtained from experimental SM.     

3D real-time micromechanical compressive behaviour of bone-cement interface: experimental and finite element studies

The integrity of bone-cement interface is essential for the long-term stability of cemented total joint arthroplasty. Although several studies have been carried out on bone-cement interface at continuum level, micromechanics of the interface has been studied only recently for tensile and shear loading cases. Fundamental studies of bone-cement interface at microstructural level are critical to the understanding of the failure processes of the interface, where multiple factors may contribute to failure. Here we present a micromechanical study of bone-cement interface under compression, which utilised in situ mechanical testing, time-lapsed microcomputed tomography (CT) and finite element (FE) modelling. Bovine trabecular bone was used to interdigitate with bone cement to obtain bone-cement interface samples, which were tested in step-wise compression using a custom-made loading stage within the μCT chamber. A finite element model was built from the CT images of one of the tested samples and loaded similarly as in the experiment. The simulated stress-displacement response fell within the range of the experimental responses, and the predicted local strain distribution correlated well with the failure pattern in the subject-specific experimental model. Damage evolution with load in the samples was monitored both experimentally and numerically. The results from the FE simulations further revealed the development of damage in the regions of interest during compression, which may be useful towards a micromechanics understanding of the failure processes at bone-cement interface.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

HR-pQCT-based homogenised finite element models provide quantitative predictions of experimental vertebral body stiffness and strength with the same accuracy as μFE models

This study validated two different high-resolution peripheral quantitative computer tomography (HR-pQCT)-based finite element (FE) approaches, enhanced homogenised continuum-level (hFE) and micro-finite element (μFE) models, by comparing them with compression test results of vertebral body sections. Thirty-five vertebral body sections were prepared by removing endplates and posterior elements, scanned with HR-pQCT and tested in compression up to failure. Linear hFE and μFE models were created from segmented and grey-level CT images, and apparent model stiffness values were compared with experimental stiffness as well as strength results. Experimental and numerical apparent elastic properties based on grey-level/segmented CT images (N=35) correlated well for μFE (r2=0.748/0.842) and hFE models (r2=0.741/0.864). Vertebral section stiffness values from the linear μFE/hFE models estimated experimental ultimate apparent strength very well (r2=0.920/0.927). Calibrated hFE models were able to predict quantitatively apparent stiffness with the same accuracy as μFE models. However, hFE models needed no back-calculation of a tissue modulus or any kind of fitting and were computationally much cheaper.     

Topographical brain electrical activity mapping on an IBM-compatible personal computer

A new system for computing brain electrical activity maps on a standard IBM-compatible computer has been developed. The EEG is recorded using a CED1401 intelligent laboratory interface and stored in the computer; colour coded maps are generated using softwave developed in Turbo PASCAL and displayed on the EGA graphics screen. An acceptable computation time of 1.2 s for a 64 × 64 map displayed as a 128 × 128 pixel image has been achieved by incorporating assembly language routines and a maths coprocessor. The system may be readily upgraded as improved hardware becomes available and further software can be added. In addition to triggering auditory, visual and somatosensory stimulators, it provides the potential for the generation of complex stimuli for cognitive experiments by means of mass RAM and digital-to-analogue converters.     

A lower extremity model for muscle-driven simulation of activity using explicit finite element modeling

A key strength of computational modeling is that it can provide estimates of muscle, ligament, and joint loads, stresses, and strains through non-invasive means. However, simulations that can predict the forces in the muscles during activity while maintaining sufficient complexity to realistically represent the muscles and joint structures can be computationally challenging. For this reason, the current state of the art is to apply separate rigid-body dynamic and finite-element (FE) analyses in series. However, the use of two or more disconnected models often fails to capture key interactions between the joint-level and whole-body scales. Single framework MSFE models have the potential to overcome the limitations associated with disconnected models in series. The objectives of the current study were to create a multi-scale FE model of the human lower extremity that combines optimization, dynamic muscle modeling, and structural FE analysis in a single framework and to apply this framework to evaluate the mechanics of healthy knee specimens during two activities. Two subject-specific FE models (Model 1, Model 2) of the lower extremity were developed in ABAQUS/Explicit including detailed representations of the muscles. Muscle forces, knee joint loading, and articular contact were calculated for two activities using an inverse dynamics approach and static optimization. Quadriceps muscle forces peaked at the onset of chair rise (2174 N, 1962 N) and in early stance phase (510 N, 525 N), while gait saw peak forces in the hamstrings (851 N, 868 N) in midstance. Joint forces were similar in magnitude to available telemetric patient data. This study demonstrates the feasibility of detailed quasi-static, muscle-driven simulations in an FE framework.     

A visco-hyperelastic constitutive approach for modeling polyvinyl alcohol sponge

This study proposes the quasi-linear viscoelastic (QLV) model to characterize the time dependent mechanical behavior of poly(vinyl alcohol) (PVA) sponges. The PVA sponges have implications in many viscoelastic soft tissues, including cartilage, liver, and kidney as an implant. However, a critical barrier to the use of the PVA sponge as tissue replacement material is a lack of sufficient study on its viscoelastic mechanical properties. In this study, the nonlinear mechanical behavior of a fabricated PVA sponge is investigated experimentally and computationally using relaxation and stress failure tests as well as finite element (FE) modeling. Hyperelastic strain energy density functions, such as Yeoh and Neo-Hookean, are used to capture the mechanical behavior of PVA sponge at ramp part, and viscoelastic model is used to describe the viscose behavior at hold part. Hyperelastic material constants are obtained and their general prediction ability is verified using FE simulations of PVA tensile experiments. The results of relaxation and stress failure tests revealed that Yeoh material model can define the mechanical behavior of PVA sponge properly compared with Neo-Hookean one. FE modeling results are also affirmed the appropriateness of Yeoh model to characterize the mechanical behavior of PVA sponge. Thus, the Yeoh model can be used in future biomechanical simulations of the spongy biomaterials. These results can be utilized to understand the viscoelastic behavior of PVA sponges and has implications for tissue engineering as scaffold.