There may be different causes of failures in bone; however, their origin generally lies at the lowest level of structural hierarchy, i.e., at the mineral-collagen composite. Any change in the nanostructure affects the affinity or bonding effectiveness between and within the phases at this level, and hence determines the overall strength and quality of bone. In this study, we propose a novel concept to assess change in the nanostructure and thereby change in the bonding status at this level by revealing change in the orientation distribution characteristics of mineral crystals. Using X-ray diffraction method, a parameter called Degree of Orientation (DO) has been quantified. The DO accounts for the azimuthal distribution of mineral crystals and represents their effective amount along any direction. Changes in the DOs in cortical bone samples from bovine femur with different preferential orientations of mineral crystals were estimated under external loads. Depending on the applied loads, change in the azimuthal distribution of the DOs and the degree of reversibility of the crystals was observed to vary. The characteristics of nanostructural change and thereby possible affect on the strength of bone was then predicted from the reversible or irreversible characteristics of distributed mineral crystals. Significant changes in the organization of mineral crystals were observed; however, variations in the applied stresses and elastic moduli were not evinced at the macroscale level. A novel concept to assess the alteration in nanostructure on the basis of mineral crystals orientation distribution has been proposed. The importance of nanoscale level information obtained noninvasively has been emphasized, which acts as a precise tool to estimate the strength and predict the possible fracture risks in bone.

Among the available tools for the early diagnosis of breast cancer, the elastographic technique based on ultrasounds has many advantages such as the noninvasive measure, the absence of ionizing effects, the high tolerability by patients, and the wide diffusion of the ecographic machines. However this diagnostic procedure is strongly affected by many subjective factors and is considered not reliable enough even to reduce the number of biopsies used to identify the nature of lesions. Therefore in the literature experimental and numerical simulations on physical and virtual phantoms are presented to test and validate procedures and algorithms and to interpret elastosonographic results. In this work, first a description of the elastographic technique and a review of the principal finite element (FE) models are provided and second diagnostic indexes employed to assess the nature of a lump mass are presented. As advances in FE simulations of elastosonography, axisymmetric phantom, and anthropomorphic models are described, which, with respect to the literature, include some features of breast mechanics. In particular deterministic analyses were used to compare the various details of virtual elastograms and also to investigate diagnostic indexes with respect to the regions where strains were considered. In order to improve the reliability of the elastosonographic procedure, univariate and multivariate sensitivity analyses, based on a probabilistic FE approach, were also performed to identify the parameters that mostly influence the deformation contrast between healthy and cancerous tissues. Moreover, synthetic indicators of the strain field, such as the strain contrast coefficient, were evaluated in different regions of interest in order to identify the most suitable for lesion type assessment. The deterministic analyses show that the malignant lesion is characterized by a uniform strain inside the inclusion due to the firmly bonding condition, while in the benign inclusion (loosely bonded) a strain gradient is observed independently from the elastic modulus contrast. The multivariate analyses reveal that the strain contrast depends linearly on the relative stiffness between the lesion and the healthy tissue and not linearly on the interface friction coefficient. The anthropomorphic model shows other interesting features, such as the layer or curvature effects, which introduce difficulties in selecting a reference region for strain assessment. The results show that a simple axisymmetric model with linear elastic material properties can be suitable to simulate the elastosonographic procedure although the breast curvature and layer distinction play a significant role in the strain assessment.

The objective of the present study was to assess the influence of various clinically relevant scenarios on the strain distribution in the biomechanical surrounding of five different dental implant macrogeometries. The biomechanical environment surrounding an implant, i.e., the cortical and trabecular bone, was modeled along with the implant. These models included two different values of the study parameters including loading conditions, trabecular bone elastic modulus, cortical/trabecular bone thickness ratio, and bone loss for five implant designs. Finite element analysis was conducted on the models and strain in the bones surrounding the implant was calculated. Bone volumes having strains in four different windows of 0 – 200 μ ε , 200 – 1000 μ ε , 1000 – 3000 μ ε , and > 3000 μ ε were measured and the effect of each biomechanical variable and their two-way interactions were statistically analyzed using the analysis of variance method. This study showed that all the parameters included in this study had an effect on the volume of bones in all strain windows, except the implant design, which affected only the 0 – 200 μ ε and > 3000 μ ε windows. The two-way interaction results showed that interactions existed between implant design and bone loss, and loading condition, bone loss in the 200 – 1000 μ ε window, and between implant design and loading condition in the 0 – 200 μ ε window. Within the limitations of the present methodology, it can be concluded that although some unfavorable clinical scenarios demonstrated a higher volume of bone in deleterious strain levels, a tendency toward the biomechanical equilibrium was evidenced regardless of the implant design.

Osteoporosis-related vertebral body fractures involve large compressive strains of trabecular bone. The small strain mechanical properties of the trabecular bone such as the elastic modulus or ultimate strength can be estimated using the volume fraction and a second order fabric tensor, but it remains unclear if similar estimations may be extended to large strain properties. Accordingly, the aim of this work is to identify the role of volume fraction and especially fabric in the large strain compressive behavior of human trabecular bone from various anatomical locations. Trabecular bone biopsies were extracted from human T12 vertebrae ( n = 31 ) , distal radii ( n = 43 ) , femoral head ( n = 44 ) , and calcanei ( n = 30 ) , scanned using microcomputed tomography to quantify bone volume fraction ( B V / T V ) and the fabric tensor ( M ) , and tested either in unconfined or confined compression up to very large strains ( ∼ 70 % ) . The mechanical parameters of the resulting stress-strain curves were analyzed using regression models to examine the respective influence of B V / T V and fabric eigenvalues. The compressive stress-strain curves demonstrated linear elasticity, yielding with hardening up to an ultimate stress, softening toward a minimum stress, and a steady rehardening followed by a rapid densification. For the pooled experiments, the average minimum stress was 1.89 ± 1.77 MPa , while the corresponding mean strain was 7.15 ± 1.84 % . The minimum stress showed a weaker dependence with fabric as the elastic modulus or ultimate strength. For the confined experiments, the stress at a logarithmic strain of 1.2 was 8.08 ± 7.91 MPa , and the dissipated energy density was 5.67 ± 4.42 MPa . The latter variable was strongly related to the volume fraction ( R 2 = 0.83 ) but the correlation improved only marginally with the inclusion of fabric ( R 2 = 0.84 ) . The influence of fabric on the mechanical properties of human trabecular bone decreases with increasing strain, while the role of volume fraction remains important. In particular, the ratio of the minimum versus the maximum stress, i.e., the relative amount of softening, decreases strongly with fabric, while the dissipated energy density is dominated by the volume fraction. The collected results will prove to be useful for modeling the softening and densification of the trabecular bone using the finite element method.

The glenoid labrum is an integral component of the glenohumeral capsule’s insertion into the glenoid, and changes in labrum geometry and mechanical properties may lead to the development of glenohumeral joint pathology. The objective of this research was to determine the effect that changes in labrum thickness and modulus have on strains in the labrum and glenohumeral capsule during a simulated physical examination for anterior instability. A labrum was incorporated into a validated, subject-specific finite element model of the glenohumeral joint, and experimental kinematics were applied simulating application of an anterior load at 0 deg, 30 deg, and 60 deg of external rotation and 60 deg of glenohumeral abduction. The radial thickness of the labrum was varied to simulate thinning tissue, and the tensile modulus of the labrum was varied to simulate degenerating tissue. At 60 deg of external rotation, a thinning labrum increased the average and peak strains in the labrum, particularly in the labrum regions of the axillary pouch (increased 10.5% average strain) and anterior band (increased 7.5% average strain). These results suggest a cause-and-effect relationship between age-related decreases in labrum thickness and increases in labrum pathology. A degenerating labrum also increased the average and peak strains in the labrum, particularly in the labrum regions of the axillary pouch (increased 15.5% strain) and anterior band (increased 10.4% strain). This supports the concept that age-related labrum pathology may result from tissue degeneration. This work suggests that a shift in capsule reparative techniques may be needed in order to include the labrum, especially as activity levels in the aging population continue to increase. In the future validated, finite element models of the glenohumeral joint can be used to explore the efficacy of new repair techniques for glenoid labrum pathology.

The liver harvested from a donor must be preserved and transported to a suitable recipient immediately for a successful liver transplantation. In this process, the preservation period is the most critical, since it is the longest and most tissue damage occurs during this period due to the reduced blood supply to the harvested liver and the change in its temperature. We investigate the effect of preservation period on the dynamic material properties of bovine liver using a viscoelastic model derived from both impact and ramp and hold experiments. First, we measure the storage and loss moduli of bovine liver as a function of excitation frequency using an impact hammer. Second, its time-dependent relaxation modulus is measured separately through ramp and hold experiments performed by a compression device. Third, a Maxwell solid model that successfully imitates the frequency- and time-dependent dynamic responses of bovine liver is developed to estimate the optimum viscoelastic material coefficients by minimizing the error between the experimental data and the corresponding values generated by the model. Finally, the variation in the viscoelastic material coefficients of bovine liver are investigated as a function of preservation period for the liver samples tested 1 h, 2 h, 4 h, 8 h, 12 h, 24 h, 36 h, and 48 h after harvesting. The results of our experiments performed with three animals show that the liver tissue becomes stiffer and more viscous as it spends more time in the preservation cycle.

The impact and penetration of high speed particles with the human skin is of interest for targeted drug delivery by transdermal powder injection. However, it is often difficult to perform penetration experiments on dermal tissue using micron scale particles. To address this, a finite element model of the impact and penetration of a 2 μ m gold particle into the human dermis was developed and calibrated using experiments found in the literature. Using dimensional analysis, the model was linked to a larger scale steel ball-gelatin system in order to extract key material parameters for both systems and perform impact studies. In this manner, an elastic modulus of 2.25 MPa was found for skin, in good agreement with reported values from the literature. Further gelatin experiments were performed with steel, polymethyl methacrylate, titanium, and tungsten carbide balls in order to determine the effects of particle size and density on penetration depth. Both the finite element model and the steel-gelatin experiments were able to predict the penetration behavior that was found by other investigators in the study of the impact of typical particles used for vaccine delivery into the human dermis. It can therefore be concluded that scaled up systems utilizing ballistic gelatins can be used to investigate the performance of transdermal powder injection technology.

In the past decade, novel tools for surgical planning and disease diagnosis have been developed to detect the liver’s mechanical properties. Some tools utilize superficial indentation type measurements to determine the elasticity of the liver parenchyma and to assume material homogeneity. In fact, the liver is a soft tissue covered with a connective sheathing that is called Glisson’s capsule. This article quantifies the effect of this capsule on the measured or “effective” elastic modulus obtained by indentation with a spherical geometry. Two sets of parametric computational studies in which the Glisson capsule thickness and elasticity were varied, demonstrated the relationship between the measured elastic response and the underlying parenchymal elastic response. Previously reported in vivo indentation data on the human liver were utilized to determine the elasticity of its parenchyma. The results indicated a linear relationship between the effective (measured) elastic response and the underlying parenchyma for the Glisson capsule thicknesses considered. When previously published human liver indentation data were analyzed, the measured elastic modulus was nearly 6.9% greater than the underlying parenchyma elastic modulus. Although the analyzed data were obtained from a single liver and yet displayed a significant variation, the Glisson capsule elasticity induced a significant but systematic error as well. The Glisson capsule thickness error was negligible for capsule parameters associated with a normal liver. Based on this work, an emphasis on the Glisson capsule’s contribution to the mechanical response of the liver would enhance the clinical potential of indentation-based novel tools for liver care.

As an alternative to purely mechanical methods, optical tracking of passive osmotic swelling was used to assess mechanical properties of the porcine lens capsule. A simple model was developed accounting for the permeability of the lens fiber cells and capsule to water, the concentration of fixed charges in the fiber cells, and the capsule’s resistance to the swelling of fiber cells. Fitting the model solution to experimental data provided an estimate of the elastic modulus of the lens capsule under the assumption of linear isotropic elasticity. The calculated elastic modulus at a fixed charge density of 20 mol m − 3 was 2.0 ± 0.5 MPa ( mean ± 95 % confidence interval; n = 15 ) for 0.1% saline solution, 0.64 ± 0.3 MPa ( n = 10 ) for 0.2% saline solution, and 0.28 ± 0.5 MPa ( n = 6 ) for 0.5% saline solution. These values are comparable to previously reported moduli of elasticity for the porcine lens capsule at small strains ( < 10 % ) , and the slight increase with hypotonicity is consistent with the nonlinear mechanical behavior of the lens capsule. Although limited by being a single measurement on a heterogeneous tissue, osmotic swelling provides a quantitative assessment of the stiffness of the lens capsule without requiring dissection or manipulation of the lens. Thus, the new method could be useful for small animal models.

Alterations to joint tissues, including subchondral bone, occur with osteoarthritis. A microindentation technique was developed to determine the local compressive modulus of subchondral bone. This test, in conjunction with a cartilage indentation test at the same location, could evaluate changes of these material properties in both tissues. The accuracy of the technique was determined by applying it to materials of known moduli. The technique was then applied to rat tibial plateaus to characterize the local moduli of the subchondral bone. An established nanoindentation method was adopted to determine the modulus of subchondral bone following penetration of the overlying articular cartilage. Three cycles of repeated loadings were applied (2.452 N, 30 s hold). The slope of the load-displacement response during the unloading portion of the third cycle was used to measure the stiffness. Indentation tests were performed on two polyurethane foams and polymethyl-methacrylate for validation ( n = 15 ) . Regression analysis was used to compare the moduli with reference values. Subchondral bone moduli of tibial plateaus from Sprague-Dawley rats ( n = 5 ) were measured for central and posterior locations of medial and lateral compartments. An analysis of variance was used to analyze the effects of compartment and test location. The measured moduli of the validation materials correlated with the reference values ( R 2 = 0.993 , p = 0.05 ). In rat tibial plateaus, the modulus of the posterior location was significantly greater than the center location ( 4.03 ± 1.00 GPa and 3.35 ± 1.16 GPa respectively, p = 0.03 ). The medial compartment was not different from the lateral compartment. This method for measuring the subchondral bone in the same location as articular cartilage allows studies of the changes in these material properties with the onset and progression of osteoarthritis.

A continuum model based on the contact mechanics theory was developed and used for evaluating virus indentation forces at the early stage of membrane invagination, as well as the work of the virus indentation forces and virus-cell contact pressures in a receptor-mediated endocytosis, depending on the virus size and virus/cell stiffnesses. The model indicated that early virus indentation forces are in the order of 1–10 pN and for a given extent of virus engulfment, they increase linearly with the elastic modulus of the host cell and also with the square of the virus radius. The work of invagination at the initial phase of virus endocytosis is in the order of tens of zeptojoules and peak virus-cell contact pressures at this stage are in the order of hundreds of Pascals to several kPa. For a given extent of virus engulfment, peak and average contact pressures increase linearly with the elastic modulus of the host cell but interestingly, they are negligibly affected by the virus size. The present model may be useful in the fields of cellular biomechanics, virology and nanodrug delivery to evaluate mechanical factors during the early phase of membrane invagination.

Elastography, which is the imaging of soft tissues on the basis of elastic modulus (or, more generally, stiffness) has become increasingly popular in the last decades and holds promise for application in many medical areas. Most of the attention has focused on inhomogeneous materials that are locally isotropic, the intent being to detect a (stiff) tumor within a (compliant) tissue. Many tissues of mechanical interest, however, are anisotropic, so a method capable of determining material anisotropy would be attractive. We present here an approach to determine the mechanical anisotropy of inhomogeneous, anisotropic tissues, by directly solving the finite element representation of the Cauchy stress balance in the tissue. The method divides the sample domain into subdomains assumed to have uniform properties and solves for the material constants in each subdomain. Two-dimensional simulated experiments on linear anisotropic inhomogeneous systems demonstrate the ability of the method, and simulated experiments on a nonlinear model demonstrate the ability of the method to capture anisotropy qualitatively even though only a linear model is used in the inverse problem. As with any inverse problem, ill-posedness is a serious concern, and multiple tests may need to be done on the same sample to determine the properties with confidence.

In vitro electromechanical and biomechanical testing of articular cartilage provide critical information about the structure and function of this tissue. Difficulties obtaining fresh tissue and lengthy experimental testing procedures often necessitate a storage protocol, which may adversely affect the functional properties of cartilage. The effects of storage at either 4 ° C for periods of 6 days and 12 days, or during a single freeze-thaw cycle at − 20 ° C were examined in young bovine cartilage. Non-destructive electromechanical measurements and unconfined compression testing on 3 mm diameter disks were used to assess cartilage properties, including the streaming potential integral (SPI), fibril modulus (Ef), matrix modulus (Em), and permeability ( k ) . Cartilage disks were also examined histologically. Compared with controls, significant decreases in SPI (to 32.3 ± 5.5 % of control values, p < 0.001 ), Ef (to 3.1 ± 41.3 % of control values, p = 0.046 ), Em (to 6.4 ± 8.5 % of control values, p < 0.0001 ), and an increase in k (to 2676.7 ± 2562.0 % of control values, p = 0.004 ) were observed at day 12 of refrigeration at 4 ° C , but no significant changes were detected at day 6. A trend toward detecting a decrease in SPI (to 94.2 ± 6.2 % of control values, p = 0.083 ) was identified following a single freeze-thaw cycle, but no detectable changes were observed for any biomechanical parameters. All numbers are mean ± 95 % confidence interval. These results indicate that fresh cartilage can be stored in a humid chamber at 4 ° C for a maximum of 6 days with no detrimental effects to cartilage electromechanical and biomechanical properties, while one freeze-thaw cycle produces minimal deterioration of biomechanical and electromechanical properties. A comparison to literature suggested that particular attention should be paid to the manner in which specimens are thawed after freezing, specifically by minimizing thawing time at higher temperatures.

Creep is an active form of time-dependent viscoelastic deformation that occurs in bone tissue during daily life. Recent findings indicate bone mineralization, which is involved in determining the elastic and plastic properties of bone matrix, can also contribute in controlling its viscoelastic property. Nanoindentation viscosity was used as a direct measure for the capacity of a material to resist viscous-like flow under loading. The objectives of this study were to examine (1) whether the nanoindentation viscosity obtained using the traditional viscoelastic Voigt model can describe creep response of bone matrix and (2) how the nanoindentation viscosity is related to contact hardness and elastic modulus. The Voigt model accurately described the creep behavior of bone matrix ( r 2 > 0.96 , p < 0.001 ) . The nanoindentation viscosity had strong relationships with nanoindentation contact hardness ( r 2 = 0.94 , p < 0.001 ) and modulus ( r 2 = 0.83 , p < 0.001 ) independent of tissue ages of osteonal bone matrix. The strong positive relationships of nanoindentation viscosity with contact hardness and modulus can be interpreted as increases in the mineral portion of bone matrix may limit the interfibril motion of collagen while enhancing the mechanical stability of bone. We suggest that previous nanoindentation results can be reanalyzed to characterize the viscoelastic creep using the Voigt model.

Age-dependent outcomes following traumatic brain injury motivate the study of brain injury biomechanics in experimental animal models at different stages of development. Finite element models of the rat brain are used to better understand the mechanical mechanisms behind these age-dependent outcomes; however, age- and region-specific rat brain tissue mechanical properties are required for biofidelity in modeling. Here, we have used the atomic force microscope (AFM) to measure region-dependent mechanical properties for subregions of the cortex and hippocampus in P10, P17, and adult rats. Apparent elastic modulus increased nonlinearly with indentation strain, and a nonlinear Ogden hyperelastic model was used to fit the force-deflection data. Subregional heterogeneous distributions of mechanical properties changed significantly with age. Apparent elastic modulus was also found to increase overall with age, increasing by > 100 % between P10 and adult rats. Unconfined compression tests ( ε = − 0.3 ) were performed on whole slices of the hippocampus and cortex of P10, P17, and adult rats to verify the mechanical properties measured with the AFM. Mean apparent elastic modulus at an indentation strain of 30% from AFM measurements for each region and age correlated well with the long-term elastic modulus measured from 30% unconfined compression tests (slope not significantly different from 1, p > 0.05 ). Protein, lipid, and sulfated glycosaminoglycan content of the brain increased with age and were positively correlated with tissue stiffness, whereas water content decreased with age and was negatively correlated with tissue stiffness. These correlations can be used to hypothesize mechanistic models for describing the mechanical behavior of brain tissue as well as to predict relative differences between brain tissue mechanical properties of other species, at different ages, and for different regions based on differences in tissue composition.

We study apparent elastic moduli of trabecular bone, which is represented, for simplicity, by a two- or three-dimensional periodic cellular network. The term “apparent” refers to the case when the region used in calculations (or specimen size) is smaller than a representative volume element and the moduli depend on the size of that region and boundary conditions. Both the bone tissue forming the network and the pores (represented by a very soft material) are assumed, for simplicity, as homogeneous, linear elastic, and isotropic. In order to investigate the effects of scale and boundary conditions on the moduli of these networks we vary the specimen size and apply four different boundary conditions: displacement, traction, mixed, and periodic. The analysis using periodic boundary conditions gives the effective moduli, while the displacement, traction, and mixed boundary conditions give apparent moduli. The apparent moduli calculated using displacement and traction boundary conditions bound the effective moduli from above and below, respectively. The larger is the size of the region used in our calculations, the closer are the bounds. Our choice of mixed boundary conditions gives results that are very close to those obtained using periodic boundary conditions. We conduct this analysis computationally using a finite element method. We also investigate the effect of mismatch in elastic moduli of bone tissue and soft fill, trabecular bone structure geometry, and bone tissue volume fraction on the apparent elastic moduli of idealized periodic models of trabecular bone. This study gives guidance on how the size of the specimen and boundary conditions (used in experiments or simulations) influence elastic moduli of cellular materials. This approach is applicable to heterogeneous materials in general.

Advanced solid freeform fabrication (SFF) techniques have been an interest for constructing tissue engineered polymeric scaffolds because of its repeatability and capability of high accuracy in fabrication resolution at the scaffold macro- and microscales. Among many important scaffold applications, hydrogel scaffolds have been utilized in tissue engineering as a technique to confide the desired proliferation of seeded cells in vitro and in vivo into its architecturally porous three-dimensional structures. Such fabrication techniques not only enable the reconstruction of scaffolds with accurate anatomical architectures but also enable the ability to incorporate bioactive species such as growth factors, proteins, and living cells. This paper presents a bioprinting system designed for the freeform fabrication of porous alginate scaffolds with encapsulated endothelial cells. The bioprinting fabrication system includes a multinozzle deposition system that utilizes SFF techniques and a computer-aided modeling system capable of creating heterogeneous tissue scaffolds. The manufacturing process is biologically compatible and is capable of functioning at room temperature and relatively low pressures to reduce the fluidic shear forces that could deteriorate biologically active species. The deposition system resolution is 10 μ m in the three orthogonal directions X Y Z and has minimum velocity of 100 μ m / s . The ideal concentrations of sodium alginate and calcium chloride were investigated to determine a viable bioprinting process. The results indicated that the suitable fabrication parameters were 1.5% ( w / v ) sodium alginate and 0.5% ( w / v ) calcium chloride. Degradation studies via mechanical testing showed a decrease in the elastic modulus by 35% after 3 weeks. Cell viability studies were conducted on the cell encapsulated scaffolds for validating the bioprinting process and determining cell viability of 83%. This work exhibits the potential use of accurate cell placement for engineering complex tissue regeneration using computer-aided design systems.

During temporomandibular joint (TMJ) function, the mandibular condylar cartilage plays a prime role in the distribution and absorption of stresses generated over the condyle. Biomechanical characterization of the tissue under compression, however, is still incomplete. The present study investigates the regional variations in the elastic and equilibrium moduli of the condylar cartilage under high strains using unconfined compression and stress relaxation, with aims to facilitate future tissue engineering studies. Porcine condylar cartilages from five regions (anterior, central, lateral, medial, and posterior) were tested under unconfined compression. Elastic moduli were obtained from the linear regions of the stress-strain curves corresponding to the continuous deformation. Equilibrium moduli were obtained from the stress relaxation curves using the Kelvin model. The posterior region was the stiffest, followed by the middle (medial, central, and lateral) regions and the anterior region, respectively. Specifically, in terms of the equilibrium modulus, the posterior region was 1.4 times stiffer than the middle regions, which were in turn 1.7 times stiffer than the anterior region, although only the difference between anterior and posterior regions was statistically significant. No significant differences in stiffness were observed among the medial, central, lateral, and posterior regions. A positive correlation between the thickness and stiffness of the cartilage was observed, reflecting that their regional variations may be related phenomena caused in response to cartilage loading patterns. Condylar cartilage was less stiff under compression than in tension. In addition, condylar cartilage under compression appears to behave in a manner similar to the TMJ disc in terms of the magnitude of moduli and drastic initial drop in stress after a ramp strain.

Characterizing the biomechanical and biotribological properties for articular surfaces in healthy, damaged, and repaired states will both elucidate the understanding of mechanical degradation and lubricating phenomena and enhance the development of functional tissue engineered cartilage and surgical repair techniques. In recent work, a new methodology involving concomitant linear translational and oscillating rotational motion was developed to determine the frictional and wear characteristics of articular cartilage. The impetus of this work was to further characterize the biomechanical characteristics from stress relaxation and dynamic cyclical indentation testing of normal and damaged articular cartilage and to correlate the biotribological characteristic findings with the biomechanical data. Quasilinear viscoelastic (QLV) theory was used to curve fit the stress-relaxation data, while the dynamic data were used both to determine the dynamic properties through fast Fourier transform analysis and to validate the dynamic behavior based on the properties obtained from the QLV theory. Comparisons of the curve-fit parameters showed a significant decrease in pre- versus postwear elastic response, A ( p < 0.04 ) , and viscous response, c ( p < 0.01 ) . In addition, the short term relaxation time, τ 1 ( p < 0.0062 ) , showed a significant decrease between surfaces with and without a defect. The magnitude of the complex modulus from dynamic tests revealed a decrease due to wear, l G l postwear ∕ l G l prewear < 1 ( p < 0.05 ) . The loss factor, tan δ , was generally greater while l G l was less for those specimens experiencing rotation. A linear regression analysis was performed to correlate μ static and μ initial with the curve-fit QLV parameters, A , B , c , τ 1 , and τ 2 . Increasing coefficients of friction correlated with decreases in the elastic response, A , viscous response, c , and the short term relaxation time constant, τ 1 , while B became increasingly nonlinear and τ 2 became shorter postwear. Qualitatively, scanning electron microscopy photographs revealed the mechanical degradation of the tissue surface due to wear. Surfaces with a defect had an increased amount of wear debris, which ultimately contributed to third body wear. Surfaces without a defect had preferentially aligned abrasions, while those surfaces not within the wear path showed no signs of wear. The efficacy of various repair techniques and innovative repair tissue models in comparison to normal and worn articular surface tissue can be determined through experimental designs involving both biomechanical and biotribological parameter characterizations. The development of this comprehensive testing scenario involving both biotribological and biomechanical characteristics is essential to the continued development of potential articular repair tissue.

Freezing of biological tissue is emerging in various biomedical applications. The success of these applications requires precise control of the tissue functionality, which is closely associated with the microstructure of the extracellular matrix (ECM). In the present study, the spatiotemporal effects of freezing on the ECM were experimentally and theoretically investigated by approximating biological tissue as a poroelastic material saturated with interstitial fluid. The experiments with type I collagen gel showed that its matrix underwent two distinct levels of structural changes due to freezing: enlarged pore structure of the matrix and increased collagen fibril diameters. The extent of these changes was augmented as the freezing temperature was lowered. The theoretical model suggested that the interstitial fluid might be transported toward the unfrozen region from the phase change interface due to the volumetric expansion associated with the water-ice phase change, and the transported fluid could interact with the matrix and enlarge its pore structure. The model also illustrated the effects of matrix structural properties on this interaction including initial porosity, hydraulic conductivity, and elastic modulus. These results imply that an identical macroscopic freezing protocol may result in different microstructural alterations of poroelastic materials depending on the structural properties of the matrix. This may be relevant to understanding the tissue-type dependent outcomes of cryomedicine applications and be useful in designing cryomedicine applications for a wide variety of tissues.