Endovascular coil embolization is now widely used to treat cerebral aneurysms (CA) as an alternative to surgical clipping. It involves filling the aneurysmal sac with metallic coils to reduce flow, induce clotting, and promote the formation of a coil/thrombus mass which protects the aneurysm wall from hemodynamic forces and prevents rupture. However, a significant number of aneurysms are incompletely coiled leading to aneurysm regrowth and/or recanalization. Computational models of aneurysm coiling may provide important new insights into the effects of intrasaccular coil and thrombus on aneurysm wall stresses. Porcine blood and platinum coils were used to construct an in vitro coil thrombus mass (CTM) for mechanical testing. A uniaxial compression test was performed with whole blood clots and CTM, with coil packing densities (CPDs) of 10%, 20%, and 30% to obtain compressive stress/strain responses. A fourth-order polynomial mechanical response function was fit to the experimentally obtained stress/strain responses for each CPD in order to represent their mechanical properties for computational simulations. Patient-specific three-dimensional (3D) geometries of three aneurysms with simple geometry and four with complex geometry were reconstructed from digital subtraction angiography (DSA) images. The CPDs were digitally inserted in the aneurysm geometries and finite element modeling was used to determine transmural peak/mean wall stress (MWS) with and without coil packing. Reproducible stress/strain curves were obtained from compression testing of CTM and the polynomial mechanical response function was found to approximate the experimental stress/strain relationship obtained from mechanical testing to a high degree. An exponential increase in the CTM stiffness was observed with increasing CPD. Elevated wall stresses were found throughout the aneurysm dome, neck, and parent artery in simulations of the CAs with no filling. Complete, 100% filling of the aneurysms with whole blood clot and CPDs of 10%, 20%, and 30% significantly reduced MWS in simple and complex geometry aneurysms. Sequential increases in CPD resulted in significantly greater increases in MWS in simple but not complex geometry aneurysms. This study utilizes finite element analysis to demonstrate the reduction of transmural wall stress following coil embolization in patient-specific computational models of CAs. Our results provide a quantitative measure of the degree to which CPD impacts wall stress and suggest that complex aneurysmal geometries may be more resistant to coil embolization treatment. The computational modeling employed in this study serves as a first step in developing a tool to evaluate the patient-specific efficacy of coil embolization in treating CAs.

Currently available knee joint kinematic tracking systems fail to nondestructively capture the subtle variation in joint and soft tissue kinematics that occur in native, injured, and reconstructed joint states. Microcomputed tomography (CT) imaging has the potential as a noninvasive, high-resolution kinematic tracking system, but no dynamic simulators exist to take advantage of this. The purpose of this work was to develop and assess a novel micro-CT compatible knee joint simulator to quantify the knee joint's kinematic and kinetic response to clinically (e.g., pivot shift test) and functionally (e.g., gait) relevant loading. The simulator applies closed-loop, load control over four degrees-of-freedom (DOF) (internal/external rotation, varus/valgus rotation, anterior/posterior translation, and compression/distraction), and static control over a fifth degree-of-freedom (flexion/extension). Simulator accuracy (e.g., load error) and repeatability (e.g., coefficient of variation) were assessed with a cylindrical rubber tubing structure and a human cadaveric knee joint by applying clinically and functionally relevant loads along all active axes. Micro-CT images acquired of the joint at a loaded state were then used to calculate joint kinematics. The simulator loaded both the rubber tubing and the cadaveric specimen to within 0.1% of the load target, with an intertrial coefficient of variation below 0.1% for all clinically relevant loading protocols. The resultant kinematics calculated from the acquired images agreed with previously published values, and produced errors of 1.66 mm, 0.90 mm, 4.41 deg, and 1.60 deg with respect to anterior translation, compression, internal rotation, and valgus rotation, respectively. All images were free of artifacts and showed knee joint displacements in response to clinically and functionally loading with isotropic CT image voxel spacing of 0.15 mm. The results of this study demonstrate that the joint-motion simulator is capable of applying accurate, clinically and functionally relevant loads to cadaveric knee joints, concurrent with micro-CT imaging. Nondestructive tracking of bony landmarks allows for the precise calculation of joint kinematics with less error than traditional optical tracking systems.

Understanding the pressure distributions at the limb-socket interface is essential to the design and evaluation of prosthetic components for lower limb prosthesis users. Force sensing resistors (FSRs) are employed in prosthetics research to measure pressure at this interface due to their low cost, thin profile, and ease of use. While FSRs are known to be sensitive to many sources of error, few studies have systematically quantified these errors using test conditions relevant to lower limb prosthetics. The purpose of this study was to evaluate FSR accuracy for the measurement of lower limb prosthetics interface pressures. Two FSR models (Flexiforce A201 and Interlink 402) were subjected to a series of prosthetic-relevant tests. These tests included: (1) static compression, (2) cyclic compression, and (3) a combined static and cyclic compression protocol mimicking a variable activity (walk–sit–stand) procedure. Flexiforce sensors outperformed Interlink sensors and were then subjected to two additional tests: (4) static curvature and (5) static shear stress. Results demonstrated that FSRs experienced significant errors in all five tests. We concluded that: (1) if used carefully, FSRs can provide an estimate of prosthetic interface pressure, but these measurements should be interpreted within the expected range of possible measurement error given the setup; (2) FSRs should be calibrated in a setup that closely matches how they will be used for taking measurements; and (3) both Flexiforce and Interlink sensors can be used to estimate interface pressures; however, in most cases Flexiforce sensors are likely to provide more accurate measurements.

This paper presents in vivo mechanical characterization of the muscularis, submucosa, and mucosa of the porcine stomach wall under large deformation loading. This is particularly important for the development of gastrointestinal pathology-specific surgical intervention techniques. The study is based on testing the cardiac and fundic glandular stomach regions using a custom-developed compression ultrasound elastography system. Particular attention has been paid to elucidate the heterogeneity and anisotropy of tissue response. A Fung hyperelastic material model has been used to model the mechanical response of each tissue layer. A univariate analysis comparing the initial shear moduli of the three layers indicates that the muscularis (5.69 ± 4.06 kPa) is the stiffest followed by the submucosa (3.04 ± 3.32 kPa) and the mucosa (0.56 ± 0.28 kPa). The muscularis is found to be strongly distinguishable from the mucosa tissue in the cardiac and fundic regions based on a multivariate discriminant analysis. The cardiac muscularis is observed to be stiffer than the fundic muscularis tissue (shear moduli of 7.96 ± 3.82 kPa versus 3.42 ± 2.96 kPa), more anisotropic (anisotropic parameter of 2.21 ± 0.77 versus 1.41 ± 0.38), and strongly distinguishable from its fundic counterpart. The results are consistent with the tissue morphology and are in accordance with our previous ex vivo tissue study. Finally, a univariate comparison of the in vivo and ex vivo initial shear moduli for each layer shows that the muscularis and submucosa tissues are softer while in vivo, but the mucosa tissue is stiffer while in vivo. The results concerning the mechanical properties highlight the inhomogeneity and anisotropy of multilayer stomach tissue.

Osteoarthritis (OA) is a degenerative joint disease that affects millions of people worldwide, yet its disease mechanism is not clearly understood. Animal models have been established to study disease progression by initiating OA through modified joint mechanics or altered biological activity within the joint. However, animal models often do not have the capability to directly relate the mechanical environment to joint damage. This review focuses on a novel in vivo approach based on controlled, cyclic tibial compression to induce OA in the mouse knee. First, we discuss the development of the load-induced OA model, its different loading configurations, and other techniques used by research laboratories around the world. Next, we review the lessons learned regarding the mechanobiological mechanisms of load-induced OA and relate these findings to the current understanding of the disease. Then, we discuss the role of specific genetic and cellular pathways involved in load-induced OA progression and the contribution of altered tissue properties to the joint response to mechanical loading. Finally, we propose using this approach to test the therapeutic efficacy of novel treatment strategies for OA. Ultimately, elucidating the mechanobiological mechanisms of load-induced OA will aid in developing targeted treatments for this disabling disease.

Articular cartilage focal defects are common soft tissue injuries potentially linked to osteoarthritis development. Although several defect characteristics likely contribute to osteoarthritis, their relationship to local tissue deformation remains unclear. Using finite element models with various femoral cartilage geometries, we explore how defects change cartilage deformation and joint kinematics assuming loading representative of the maximum joint compression during the stance phase of gait. We show how defects, in combination with location-dependent cartilage mechanics, alter deformation in affected and opposing cartilages, as well as joint kinematics. Small and average sized defects increased maximum compressive strains by approximately 50% and 100%, respectively, compared to healthy cartilage. Shifts in the spatial locations of maximum compressive strains of defect containing models were also observed, resulting in loading of cartilage regions with reduced initial stiffnesses supporting the new, elevated loading environments. Simulated osteoarthritis (modeled as a global reduction in mean cartilage stiffness) did not significantly alter joint kinematics, but exacerbated tissue deformation. Femoral defects were also found to affect healthy tibial cartilage deformations. Lateral femoral defects increased tibial cartilage maximum compressive strains by 25%, while small and average sized medial defects exhibited decreases of 6% and 15%, respectively, compared to healthy cartilage. Femoral defects also affected the spatial distributions of deformation across the articular surfaces. These deviations are especially meaningful in the context of cartilage with location-dependent mechanics, leading to increases in peak contact stresses supported by the cartilage of between 11-34% over healthy cartilage.

An analytical theory for the unconfined creep behavior of a cylindrical inclusion (simulating a soft tissue tumor) embedded in a cylindrical background sample (simulating normal tissue) is presented and analyzed in this paper. Both the inclusion and the background are considered as fluid-filled, porous materials, each of them being characterized by a set of mechanical properties. Specifically, in this paper, the inclusion is considered to be less permeable than the background. The cylindrical sample is compressed using a constant pressure within two frictionless plates and is allowed to expand in an unconfined way along the radial direction. Analytical expressions for the effective Poisson's ratio (EPR) and fluid pressure inside and outside the inclusion are derived and analyzed. The theoretical results are validated using finite element models (FEMs). Statistical analysis shows excellent agreement between the results obtained from the developed model and the results from FEM. Thus, the developed theoretical model can be used in medical imaging modalities such as ultrasound poroelastography to extract the mechanical parameters of tissues and/or to better understand the impact of different mechanical parameters on the estimated displacements, strains, stresses, and fluid pressure inside a tumor and in the surrounding tissue.

An analytical theory for the unconfined creep behavior of a cylindrical inclusion (simulating a soft tissue tumor) embedded in a cylindrical background sample (simulating normal tissue) is presented and analyzed in this paper. Both the inclusion and the background are considered as fluid-filled, porous materials, each of them being characterized by a set of mechanical parameters. Specifically, in this derivation, the inclusion is assumed to have significantly higher interstitial permeability than the background. The formulations of the effective Poisson's ratio (EPR) and fluid pressure in the inclusion and in the background are derived for the case of a sample subjected to a creep compression. The developed analytical expressions are validated using finite element models (FEM). Statistical comparison between the results obtained from the developed model and the results from FEM demonstrates accuracy of the proposed theoretical model higher than 99.4%. The model presented in this paper complements the one reported in the companion paper (Part I), which refers to the case of an inclusion having less interstitial permeability than the background.

Aging and degeneration of the intervertebral disk are noted by changes in tissue composition and geometry, including a decrease in nucleus pulposus (NP) area. The NP centroid is positioned slightly posterior of the disk's centroid, but the effect of NP size and location on disk joint mechanics is not well understood. We evaluated the effect of NP size and centroid location on disk joint mechanics under dual-loading modalities (i.e., compression in combination with axial rotation or bending). A finite element model (FEM) was developed to vary the relative NP area (NP:Disk area ratio range = 0.21–0.60). We also evaluated the effect of NP position by shifting the NP centroid anteriorly and posteriorly. Our results showed that compressive stiffness and average first principal strains increased with NP size. Under axial compression, stresses are distributed from the NP to the annulus, and stresses were redistributed toward the NP with axial rotation. Moreover, peak stresses were greater for disks with a smaller NP area. NP centroid location had a greater impact on intradiscal pressure during flexion and extension, where peak pressures in the posterior annulus under extension was greater for disks with a more posteriorly situated NP. In conclusion, the findings from this study highlight the importance of closely mimicking NP size and location in computational models that aim to understand stress/strain distribution during complex loading and for developing repair strategies that aim to recapitulate the mechanical behavior of healthy disks.

A fluid–structure interaction (FSI) model of a left anterior descending (LAD) coronary artery was developed, incorporating transient blood flow, cyclic bending motion of the artery, and myocardial contraction. The three-dimensional (3D) geometry was constructed based on a patient's computed tomography angiography (CTA) data. To simulate disease conditions, a plaque was placed within the LAD to create a 70% stenosis. The bending motion of the blood vessel was prescribed based on the LAD spatial information. The pressure induced by myocardial contraction was applied to the outside of the blood vessel wall. The fluid domain was solved using the Navier–Stokes equations. The arterial wall was defined as a nonlinear elastic, anisotropic, and incompressible material, and the mechanical behavior was described using the modified hyper-elastic Mooney–Rivlin model. The fluid (blood) and solid (vascular wall) domains were fully coupled. The simulation results demonstrated that besides vessel bending/stretching motion, myocardial contraction had a significant effect on local hemodynamics and vascular wall stress/strain distribution. It not only transiently increased blood flow velocity and fluid wall shear stress, but also changed shear stress patterns. The presence of the plaque significantly reduced vascular wall tensile strain. Compared to the coronary artery models developed previously, the current model had improved physiological relevance.

Application of tibiofemoral compression force (TCF) has been shown to produce anterior cruciate ligament (ACL) injuries in a laboratory setting. A new robotic testing methodology was utilized to predict ACL forces generated by TCF without directly loading the ligament. We hypothesized that ACL force, directly recorded by a miniature load cell during an unconstrained test, could be predicted by measurements of anterior tibial restraining force (ARF) recorded during a constrained test. The knee was first flexed under load control with 25 N TCF (tibial displacements and rotations unconstrained) to record a baseline kinematic pathway. Tests were repeated with increasing levels of TCF, while recording ACL force and knee kinematics. Then tests with increasing TCF were performed under displacement control to reproduce the baseline kinematic pathway (tibia constrained), while recording ARF. This allowed testing to 1500 N TCF since the ACL was not loaded. TCF generated ACL force for all knees (n = 10) at 50 deg flexion, and for eight knees at 30 deg flexion (unconstrained test). ACL force (unconstrained test) and ARF (constrained test) had strong linear correlations with TCF at both flexion angles (R 2 from 0.85 to 0.99), and ACL force was strongly correlated with ARF at both flexion angles (R 2 from 0.76 to 0.99). Under 500 N TCF, the mean error between ACL force prediction from ARF regression and measured ACL force was 4.8±7.3 N at 30 deg and 8.8±27.5 N at 50 deg flexion. Our hypothesis was confirmed for TCF levels up to 500 N, and ARF had a strong linear correlation with TCF up to 1500 N TCF.

This study formulates a finite element algorithm for frictional contact of solid materials, accommodating finite deformation and sliding. The algorithm uses a penalty method regularized with an augmented Lagrangian scheme to enforce contact constraints in a nonmortar surface-to-surface approach. Use of a novel kinematical approach to contact detection and enforcement of frictional constraints allows solution of complex problems previously requiring mortar methods or contact smoothing algorithms. Patch tests are satisfied to a high degree of accuracy with a single-pass penalty method, ensuring formulation errors do not affect the solution. The accuracy of the implementation is verified with Hertzian contact, and illustrations demonstrating the ability to handle large deformations and sliding are presented and validated against prior literature. A biomechanically relevant example addressing finger friction during grasping demonstrates the utility of the proposed algorithm. The algorithm is implemented in the open source software febio , and the source code is made available to the general public.

A wide range of materials have been used for the development of intervertebral cages. Poly(propylene fumarate) (PPF) has been shown to be an excellent biomaterial with characteristics similar to trabecular bone. Hydroxyapatite (HA) has been shown to enhance biocompatibility and mechanical properties of PPF. The purpose of this study was to characterize the effect of PPF augmented with HA (PPF:HA) and evaluate the feasibility of this material for the development of cervical cages. PPF was synthesized and combined with HA at PPF:HA wt:wt ratios of 100:0, 80:20, 70:30, and 60:40. Molds were fabricated for testing PPF:HA bulk materials in compression, bending, tension, and hardness according to ASTM standards, and also for cage preparation. The cages were fabricated with and without holes and with porosity created by salt leaching. The samples as well as the cages were mechanically tested using a materials testing frame. All elastic moduli as well as the hardness increased significantly by adding HA to PPF (p < 0.0001). The 20 wt % HA increased the moduli significantly compared to pure PPF (p < 0.0001). Compressive stiffness of all cages also increased with the addition of HA. HA increased the failure load of the porous cages significantly (p = 0.0018) compared with nonporous cages. PPF:HA wt:wt ratio of 80:20 proved to be significantly stiffer and stronger than pure PPF. The current results suggest that this polymeric composite can be a suitable candidate material for intervertebral body cages.

While the anisotropic behavior of the complex composite myocardial tissue has been well characterized in recent years, the compressibility of the tissue has not been rigorously investigated to date. In the first part of this study, we present experimental evidence that passive-excised porcine myocardium exhibits volume change. Under tensile loading of a cylindrical specimen, a volume change of 4.1 ± 1.95 % is observed at a peak stretch of 1.3. Confined compression experiments also demonstrate significant volume change in the tissue (loading applied up to a volumetric strain of 10%). In order to simulate the multiaxial passive behavior of the myocardium, a nonlinear volumetric hyperelastic component is combined with the well-established Holzapfel–Ogden anisotropic hyperelastic component for myocardium fibers. This framework is shown to describe the experimentally observed behavior of porcine and human tissues under shear and biaxial loading conditions. In the second part of the study, a representative volumetric element (RVE) of myocardium tissue is constructed to parse the contribution of the tissue vasculature to observed volume change under confined compression loading. Simulations of the myocardium microstructure suggest that the vasculature cannot fully account for the experimentally measured volume change. Additionally, the RVE is subjected to six modes of shear loading to investigate the influence of microscale fiber alignment and dispersion on tissue-scale mechanical behavior.

Computational fluid dynamics (CFD) modeling of myocardial bridging (MB) remains challenging due to its dynamic and phasic nature. This study aims to develop a patient-specific CFD model of MB. There were two parts to this study. The first part consisted of developing an in silico model of the left anterior descending (LAD) coronary artery of a patient with MB. In this regard, a moving-boundary CFD algorithm was developed to simulate the patient-specific muscle compression caused by MB. A second simulation was also performed with the bridge artificially removed to determine the hemodynamics in the same vessel in the absence of MB. The second part of the study consisted of hemodynamic analysis of three patients with mild and moderate and severe MB in their LAD by means of the developed in silico model in the first part. The average shear stress in the proximal and bridge segments for model with MB were significantly different from those for model without MB (proximal segment: 0.32 ± 0.14 Pa (with MB) versus 0.97 ± 0.39 Pa (without MB), P < 0.0001 — bridge segment: 2.60 ± 0.94 Pa (with MB) versus 1.50 ± 0.64 Pa (without MB), P < 0.0001). When all three patients were evaluated, increasing the degree of vessel compression shear stress in the proximal segment decreased, whereas the shear stress in the bridge segment increased. The presence of MB resulted in hemodynamic abnormalities in the proximal segment, whereas segments within the bridge exhibited hemodynamic patterns which tend to discourage atheroma development.

In this study, the damage evolution of liver tissue was quantified at the microstructural level under tensile, compression, and shear loading conditions using an interrupted mechanical testing method. To capture the internal microstructural changes in response to global deformation, the tissue samples were loaded to different strain levels and chemically fixed to permanently preserve the deformed tissue geometry. Tissue microstructural alterations were analyzed to quantify the accumulated damages, with damage-related parameters such as number density, area fraction, mean area, and mean nearest neighbor distance (NND). All three loading states showed a unique pattern of damage evolution, in which the damages were found to increase in number and size, but decrease in NND as strain level increased. To validate the observed damage features as true tissue microstructural damages, more samples were loaded to the above-mentioned strain levels and then unloaded back to their reference state, followed by fixation. The most major damage-relevant features at higher strain levels remained after the release of the external loading, indicating the occurrence of permanent inelastic deformation. This study provides a foundation for future structure-based constitutive material modeling that can capture and predict the stress-state dependent damage evolution in liver tissue.

A subset of temporomandibular joint (TMJ) disorders is attributed to joint degeneration. The pig has been considered the preferred in vivo model for the evaluation of potential therapies for TMJ disorders, and practical considerations such as cost and husbandry issues have favored the use of young, skeletally immature animals. However, the effect of growth on the biochemical and biomechanical properties of the TMJ disk and articulating cartilage has not been examined. The present study investigates the effect of age on the biochemical and biomechanical properties of healthy porcine TMJs at 3, 6, and 9 months of age. DNA, hydroxyproline, and glycosaminoglycan (GAG) content were determined and the disks and condyles were tested in uniaxial unconfined stress relaxation compression from 10% to 30% strain. TMJ disks were further assessed with a tensile test to failure technique, which included the ability to test multiple samples from the same region of an individual disk to minimize the intraspecimen variation. No differences in biochemical properties for the disk or compressive properties at 30% stress relaxation in the disk and condylar cartilage were found. In tension, no differences were observed for peak stress and tensile modulus. The collagen content of the condyle was higher at 9 months than 3 months (p < 0.05), and the GAG content was higher at 9 months than 6 months (p < 0.05). There was a trend of increased compressive instantaneous modulus with age. As such, age-matched controls for growing pigs are probably appropriate for most parameters measured.

Many clinical case series have reported the predisposing factors for C5 palsy and have presented comparisons of the two types of laminoplasty. However, there have been no biomechanical studies focusing on cervical spinal cord and nerve root following laminoplasty. The purpose of this study is to investigate biomechanical changes in the spinal cord and nerve roots following the two most common types of laminoplasty, open-door and double-door laminoplasty, for cervical ossification of the posterior longitudinal ligament (OPLL). A finite element (FE) model of the cervical spine and spinal cord with nerve root complex structures was developed. Stress changes in the spinal cord and nerve roots, posterior shift of the spinal cord, and displacement of the cervical nerve roots were analyzed with two types of cervical laminoplasty models for variations in the degree of canal occupying ratio and shape of the OPLL. The shape and degree of spinal cord compression caused by the OPLL had more influence on the changes in stress, posterior shift of the spinal cord, and displacement of the nerve root than the type of laminoplasty. The lateral-type OPLL resulted in imbalanced stress on the nerve roots and the highest nerve root displacement. Type of laminoplasty and shape and degree of spinal cord compression caused by OPLL were found to influence the changes in stress and posterior displacement of the cervical spinal cord and nerve roots. Lateral-type OPLL might contribute to the development of C5 palsy due to the imbalanced stress and tension on the nerve roots after laminoplasty.

Percutaneous pedicle screw fixation (PPSF) is a well-known minimally invasive surgery (MIS) employed in the treatment of thoracolumbar burst fractures (TBF). However, hardware failure and loss of angular correction are common limitations caused by the poor support of the anterior column of the spine. Balloon kyphoplasty (KP) is another MIS that was successfully used in the treatment of compression fractures by augmenting the injured vertebral body with cement. To overcome the limitations of stand-alone PPSF, it was suggested to augment PPSF with KP as a surgical treatment of TBF. Yet, little is known about the biomechanical alteration occurred to the spine after performing such procedure. The objective of this study was to evaluate and compare the immediate post-operative biomechanical performance of stand-alone PPSF, stand-alone-KP, and KP-augmented PPSF procedures. Novel three-dimensional (3D) finite element (FE) models of the thoracolumbar junction that describes the fractured spine and the three investigated procedures were developed and tested under mechanical loading conditions. The spinal stiffness, stresses at the implanted hardware, and the intradiscal pressure at the upper and lower segments were measured and compared. The results showed no major differences in the measured parameters between stand-alone PPSF and KP-augmented PPSF procedures, and demonstrated that the stand-alone KP may restore the stiffness of the intact spine. Accordingly, there was no immediate post-operative biomechanical advantage in augmenting PPSF with KP when compared to stand-alone PPSF, and fatigue testing may be required to evaluate the long-term biomechanical performance of such procedures.

There is lack of investigation capturing the complex mechanical interaction of tissue-engineered intervertebral disk (IVD) constructs in physiologically relevant environmental conditions. In this study, mechanical characterization of anisotropic electrospinning (ES) substrates made of polycaprolactone (PCL) was carried out in wet and dry conditions and viability of human bone marrow derived mesenchymal stem cells (hMSCs) seeded within double layers of ES PCL were also studied. Cyclic compression of IVD-like constructs composed of an agarose core confined by ES PCL double layers was implemented using a bioreactor and the cellular response to the mechanical stimulation was evaluated. Tensile tests showed decrease of elastic modulus of ES PCL as the angle of stretching increased, and at 60 deg stretching angle in wet, the maximum ultimate tensile strength (UTS) was observed. Based on the configuration of IVD-like constructs, the calculated circumferential stress experienced by the ES PCL double layers was 40 times of the vertical compressive stress. Confined compression of IVD-like constructs at 5% and 10% displacement dramatically reduced cell viability, particularly at 10%, although cell presence in small and isolated area can still be observed after mechanical conditioning. Hence, material mechanical properties of tissue-engineered scaffolds, composed of fibril structure of polymer with low melting point, are affected by the testing condition. Circumferential stress induced by axial compressive stimulation, conveyed to the ES PCL double layer wrapped around an agarose core, can affect the viability of cells seeded at the interface, depending on the mechanical configuration and magnitude of the load.