Numerical simulation of the mechanics of a yeast cell under high hydrostatic pressure

The mechanical effects of the compression of a yeast cell (Saccharomyces cerevisiae) under high hydrostatic pressure used for the processing of food and food ingredients are modelled and simulated with the finite-element method. The cell model consists of a cell wall, cytoplasm a lipid filled vacuole and the nucleus. Material parameters have been taken from literature or have been derived from thermodynamic relationships of water and lipids under high hydrostatic pressure. The model has been validated for a pressure load up to 250 MPa. Comparison of the volume reduction to in situ experimental observations reveals very good agreement. Dimensional analysis of the governing equations shows that transient pressure application in a high-pressure food process does not enhance structural inactivation (mechanical damage), unless pressure oscillation frequencies of 700 MHz are applied. The deformation of the cell under pressure deviates strongly from isotropic volume reduction. Especially, organelle membranes exhibit large effective strain values. Hydrostatic stress conditions are preserved in the interior part of the cell. A pressure load of 400 MPa, which is critical upon disruption of cell organelle membranes, generates an effective strain up to 80%. In the cell wall, the stress state is heterogeneous. Von-Mises stress reaches the critical value upon failure of the cell wall of 70+/-4 MPa at a pressure load between 415 and 460 MPa.     

Design and application of an oscillatory compression device for cell constructs

Mechanical compression has been shown to impact cell activity; however a need for a single device to perform a broader range of parametric studies exists. We have developed an oscillatory displacement controlled device to uniaxially strain cell constructs under both static and dynamic compression and used this device to investigate gene expression in cell constructs. The device has a wide stroke (0.25-4 mm) and frequency range (0.1-3 Hz) and several loading waveforms are possible. Alginate cellular constructs with embedded equine chondrocytes were tested and viability was maintained for the 24 h test period. Off-line mechanical testing is described and a modulus value of 18.2 +/- 1.3 kPa found for alginate disks which indicates the level of stress achieved with this deformation profile. Static (15% strain) and dynamic (15% strain, 1 Hz, triangle waveform) testing of chondrocyte constructs was performed and static compression showed significantly higher collagen II expression than dynamic using quantitative RT-PCR. In contrast, differences in matrix metalloproteinase-3 (MMP-3) expression were statistically insignificant. These studies indicate the utility of our device for studying cell activity in response to compression and suggest further studies regarding how the load and strain spectrum impact chondrocyte activity.     

An experimental method for measuring force on the spinal facet joint: description and application of the method.

A technique is described for measuring load magnitude and resultant load contact location in the facet joint in response to applied loads and moments, and the technique applied to the canine lumbar spine motion segment. Due to the cantilever beam geometry of the cranial articular process, facet joint loads result in surface strains on the lateral aspect of the cranial articular process. Strains were quantified by four strain gages cemented to the bony surface of the process. Strain measured at any one gage depended on the loading site on the articular surface of the caudal facet and on the magnitude of the facet load. Determination of facet loads during in vitro motion segment testing required calibration of the strains to known loads of various magnitudes applied to multiple sites on the caudal facet. The technique is described in detail, including placement of the strain gages. There is good repeatability of strains to applied facet loads and the strains appear independent of load distribution area. Error in the technique depends on the location of the applied facet loads, but is only significant in nonphysiologic locations. The technique was validated by two independent methods in axial torsion. Application of the technique to five in vitro canine L2-3 motion segments testing resulted in facet loads (in newtons, N) of 74+ / -23 N (mean + / -STD) in 2 newton-meter, Nm, extension, to unloaded in flexion. Lateral bending resulted in loads in the right facet of 40+ / -32 N for 1 Nm right lateral bending and 54+ / -29 N for 1 Nm left lateral bending. 4 Nm Torsion with and without 100 N axial compression resulted in facet loads of 92+ / -27 N and 69+ / -19 N, respectively. The technique is applicable to dynamic and in vivo studies.     

Effects of cross-sectional area on human bite studied with raw carrot and surimi gel

The effects of the cross-sectional area of food samples on bite force with molar teeth were investigated using raw carrots and surimi gels. We evaluated human bite force for food samples with different sizes between the upper and lower molars using a multiple-point sheet sensor and electromyography (EMG). The bite force curve and EMG clearly showed textural characteristics of the carrot and gel. In particular, the first peak in the bite curves corresponded to breaking point in the compression test. With increasing cross-sectional area of both foodstuffs, the bite force and contact area increased and the average stress to which the specimen was subjected (mean stress) tended to decrease, while the stress produced between the teeth and the specimen (active stress) did not change. Chewing rhythm and EMG activities were not greatly influenced by sample size. These findings suggest that higher bite force might cause difficulty in biting food with a larger cross-sectional area.     

Tibiofemoral joint contact force in deep knee flexion and its consideration in knee osteoarthritis and joint replacement

The aim of the study was to estimate the tibiofemoral joint force in deep flexion to consider how the mechanical load affects the knee. We hypothesize that the joint force should not become sufficiently large to damage the joint under normal contact area, but should become deleterious to the joint under the limited contact area. Sixteen healthy knees were analyzed using a motion capture system, a force plate, a surface electromyography, and a knee model, and then tibiofemoral joint contact forces were calculated. Also, a contact stress simulation using the contact areas from the literature was performed. The peak joint contact forces (M +/- SD) were 4566 +/- 1932 N at 140 degrees in rising from full squat and 4479 +/- 1478 N at 90 degrees in rising from kneeling. Under normal contact area, the tibiofemoral contact stresses in deep flexion were less than 5 MPa and did not exceed the stress to damage the cartilage. The contact stress simulation suggests that knee prosthesis having the contact area smaller than 200 mm2 may be problematic since the contact stress in deep flexion would become larger than 21 MPa, and it would lead damage or wear of the polyethylene.     

Compressive and tensile mechanical properties of the porcine nasal septum

The expanding nasal septal cartilage is believed to create a force that powers midfacial growth. In addition, the nasal septum is postulated to act as a mechanical strut that prevents the structural collapse of the face under masticatory loads. Both roles imply that the septum is subject to complex biomechanical loads during growth and mastication. The purpose of this study was to measure the mechanical properties of the nasal septum to determine (1) whether the cartilage is mechanically capable of playing an active role in midfacial growth and in maintaining facial structural integrity and (2) if regional variation in mechanical properties is present that could support any of the postulated loading regimens. Porcine septal samples were loaded along the horizontal or vertical axes in compression and tension, using different loading rates that approximate the in vivo situation. Samples were loaded in random order to predefined strain points (2–10%) and strain was held for 30 or 120 seconds while relaxation stress was measured. Subsequently, samples were loaded until failure. Stiffness, relaxation stress and ultimate stress and strain were recorded. Results showed that the septum was stiffer, stronger and displayed a greater drop in relaxation stress in compression compared to tension. Under compression, the septum displayed non-linear behavior with greater stiffness and stress relaxation under faster loading rates and higher strain levels. Under tension, stiffness was not affected by strain level. Although regional variation was present, it did not strongly support any of the suggested loading patterns. Overall, results suggest that the septum might be mechanically capable of playing an active role in midfacial growth as evidenced by increased compressive residual stress with decreased loading rates. However, the low stiffness of the septum compared to surrounding bone does not support a strut role. The relatively low stiffness combined with high stress relaxation under fast loading rates suggests that the nasal septum is a stress dampener, helping to absorb and dissipate loads generated during mastication.     

股骨头负重区松质骨的压缩力学特性

本文对19具新鲜青年(30左右)尸体股骨负重区的松质骨进行了压缩测试,测算并提供了一组松质骨弹性模量,极限应力,应变,屈服载荷及能量吸收值。试验表明松质骨力学性能主要缓解外力冲击加压。     

Using optical full-field measurement based on digital image correlation to measure strain on a tree subjected to mechanical load

Key message

The results indicate the usability of DIC technique in tree biomechanics is advantageous compared to standard procedure since it reveals the induced strain in full-field way, accurately and noninvasively.

Abstract

The goal of this study was to analyze and measure the displacement and strain response behavior of the tree (Juglans regia L.) subjected to a mechanical load. The analysis addresses issues of tree biomechanics that call for a high-quality deformation data measured on trees’ surfaces in a noninvasive way. For these purposes, the 3D digital image correlation (3D-DIC) was proposed and tested against standard extensometers. The measured tree was subjected to a bending by a pulling test to induce strain on the tree stem. The DIC technique successfully provided strain field on a tree bark despite its highly complicated geometry. Vertical averaging of the strain field obtained by the DIC revealed that a tree stem behaves according to beam theory exhibiting compression and tension parts. The absolute values of strain measured by both techniques agree with each other in order of magnitude, DIC returning lower values (approximately 21.1 and 40.8 % for compression and tension part, respectively). The results reveal necessity for future study of stress/strain transfer from xylem to bark. The sensitivity analysis shows that the computed strain highly depends on subset size used in the DIC computation. Based on the projection error and strain, the optimal subset size was found out to be between 21 and 31 pixels. The measurement proved that the DIC method can be successfully used in tree biomechanics. In general, our results and the character of DIC technique that allows testing at different scales depending on imaging methods indicate its big potential in plant biomechanics.     

Stress intensity factors for a vertical surface crack in polyethylene subject to rolling and sliding contact.

Pitting wear is a dominant form of polyethylene surface damage in total knee replacements, and may originate from surface cracks that propagate under repeated tribological contact. In the present study, stress intensity factors, KI and KII, were calculated for a surface crack in a polyethylene-CoCr-bone system in the presence of rolling or sliding contact pressures. Variations in crack length and load location were studied to determine probable crack propagation mechanisms and modes. The crack tip experienced a wide range of mixed-mode conditions that varied as a function of crack length, load location, and sliding friction. Positive KI values were observed for shorter cracks in rolling contact and for all crack lengths when the sliding load moved away from the crack. KII was greatest when the load was directly adjacent to the crack (g/a = +/- 1), where coincidental Mode I stresses were predominantly compressive. Sliding friction substantially increased both KImax and KIImax. The effective Mode I stress intensity factors, Keff, were greatest at g/a = +/- 1, illustrating the significance of high shear stresses generated by loads adjacent to surface cracks. Keff trends suggest mechanisms for surface pitting by which surface cracks propagate along their original plane under repeated reciprocating rolling or sliding, and turn in the direction of sliding under unidirectional sliding contact.     

A myoelectrically based dynamic three-dimensional model to predict loads on lumbar spine tissues during lateral bending.

This work describes a dynamic model of the low back that incorporates extensive anatomical detail of the musculo-ligamentous-skeletal system to predict the load time histories of individual tissues. The dynamic reaction moment about L4/L5 was determined during lateral bending from a linked-segment model. This reaction moment was partitioned into restorative components provided by the disc, ligament strain, and active-muscle contraction using a second model of the spine that incorporated a detailed representation of the anatomy. Muscle contraction forces were estimated using both information from surface electromyograms, collected from 12 sites, and consideration of the modulating effects of muscle length, cross-sectional area and passive elasticity. This modelling technique is sensitive to the different ways in which individuals recruit their musculature to satisfy moment constraints. Time histories of muscle forces are provided. High muscle loads are consistent with the common clinical observation of muscle strain often produced by load handling. Furthermore, the coactivation measured in muscles on both sides of the trunk suggests that muscles are recruited to satisfy the lateral bending reaction torque in addition to performing other mechanical roles such as spine stabilization. If an estimate of the intervertebral joint compression is desired for assessment of lateral bends in industry, then a single equivalent lateral muscle with a moment arm of approximately 3.0-4.0 cm would conservatively capture the effects of muscle co-contraction quantified in this study.     

Implantable MEMS compressive stress sensors: Design,fabrication and calibration with application to the disc annulus

Physiological stresses are fundamental to biomechanical testing, mechanobiological analyses, implant design, and tissue engineering. The purpose of this study was to design, fabricate, and evaluate compressive stress sensors packaged for extended, in vivo implantation in the annulus of the intervertebral disc. A commercial microelectromechanical systems (MEMS) pressure sensor die was selected as the active element for a custom stress sensor. The sensor die was modified and packaged to protect the electrical system from the biochemical and biomechanical environment. Completed sensors were calibrated under hydrostatic pressure and solid contact compression. Calibrations were performed before and after 8 weeks of in vivo implantation in a porcine disc. For the two reported sensors, stress and voltage were linearly correlated over a range of 0–1.8 MPa with less than 5% change in sensitivity. Sensitivity to solid contact stress was within 10% of that from hydrostatic pressure. In contrast to most previous studies, in which disc pressure was measured in the fluidic nucleus pulposus, these sensors may be used to measure in vivo dynamic compressive stresses in the annulus at magnitudes typical of the musculoskeletal system in a large animal over a relatively long post-operative time.     

The mechanics of the first bite

An analysis of the action of the incisor teeth in humans is presented in terms of the fracture of food particles. It is predicted that the resistance of foods with an essentially linear elastic response to an initial bite by the incisors will depend on the square root of the product of two food properties, Young's modulus and toughness. This quantity should be approximately equal to the product of the stress at cracking during a bite, and the square root of the length of a notch or indentation from which that crack initiates. As a test of the theory, the relationship between in vivo stresses and the depth of incisal penetration, measured during bites on seven 'snack' foods by 10 subjects, and food properties established from mechanical testing, was investigated. Theory and experiment were found to be in excellent agreement. A dimensionless index of the efficiency of incision is suggested, relating fracture performance by subjects to values from a testing machine. This appears to have a high level of inter-subject discrimination with efficiencies varying about threefold. The method appears to have potential applications in dentistry, food science and studies of human and primate evolution.     

Adhesive strength of single muscle cells to basement membrane at myotendinous junctions

Whole muscles loaded to failure frequently fail at or near myotendinous junctions. The present investigation was directed toward determining the breaking stress and failure site of intact and injured myotendinous junction preparations consisting of muscle cells dissected free from surrounding parallel structures but still attached to tendon collagen fibers. These tests show that the breaking stress for intact myotendinous units is 2.7 x 10(5) N/m2, expressed relative to cell cross-sectional area. Failure occurs immediately external to the junction membrane between the cell membrane and lamina densa of the basement membrane. Site and stress at failure are independent of strain and strain rate over a biologically relevant range. Breaking stress in the plane of the membrane, corrected for membrane folding, is 1.2 X 10(4) N/m2. This value is not significantly greater than stress at maximum isometric tension for these cells at these sarcomere lengths. After compression injury, cells fail within the compression site at significantly lower stress (1.9 X 10(5) N/m2). These findings suggest that, in muscle strain injuries that occur under conditions simulated here, failure occurs at myotendinous junctions unless the muscle has suffered previous compression injury leading to failure within the muscle.     

Investigation of impact loading rate effects on the ligamentous cervical spinal load-partitioning using finite element model of functional spinal unit C2–C3

The cervical spine functions as a complex mechanism that responds to sudden loading in a unique manner, due to intricate structural features and kinematics. The spinal load-sharing under pure compression and sagittal flexion/extension at two different impact rates were compared using a bio-fidelic finite element (FE) model of the ligamentous cervical functional spinal unit (FSU) C2–C3. This model was developed using a comprehensive and realistic geometry of spinal components and material laws that include strain rate dependency, bone fracture, and ligament failure. The range of motion, contact pressure in facet joints, failure forces in ligaments were compared to experimental findings. The model demonstrated that resistance of spinal components to impact load is dependent on loading rate and direction. For the loads applied, stress increased with loading rate in all spinal components, and was concentrated in the outer intervertebral disc (IVD), regions of ligaments to bone attachment, and in the cancellous bone of the facet joints. The highest stress in ligaments was found in capsular ligament (CL) in all cases. Intradiscal pressure (IDP) in the nucleus was affected by loading rate change. It increased under compression/flexion but decreased under extension. Contact pressure in the facet joints showed less variation under compression, but increased significantly under flexion/extension particularly under extension. Cancellous bone of the facet joints region was the only component fractured and fracture occurred under extension at both rates. The cervical ligaments were the primary load-bearing component followed by the IVD, endplates and cancellous bone; however, the latter was the most vulnerable to extension as it fractured at low energy impact.     

Shear wave elastography can assess the in-vivo nonlinear mechanical behavior of heel-pad

This study combines non-invasive mechanical testing with finite element (FE) modelling to assess for the first time the reliability of shear wave (SW) elastography for the quantitative assessment of the in-vivo nonlinear mechanical behavior of heel-pad. The heel-pads of five volunteers were compressed using a custom-made ultrasound indentation device. Tissue deformation was assessed from B-mode ultrasound and force was measured using a load cell to calculate the force – deformation graph of the indentation test. These results were used to design subject specific FE models and to inverse engineer the tissue’s hyperelastic material coefficients and its stress – strain behavior. SW speed was measured for different levels of compression (from 0% to 50% compression). SW speed for 0% compression was used to assess the initial stiffness of heel-pad (i.e. initial shear modulus, initial Young’s modulus). Changes in SW speed with increasing compressive loading were used to quantify the tissue’s nonlinear mechanical behavior based on the theory of acoustoelasticity. Statistical analysis of results showed significant correlation between SW-based and FE-based estimations of initial stiffness, but SW underestimated initial shear modulus by 64%(±16). A linear relationship was found between the SW-based and FE-based estimations of nonlinear behavior. The results of this study indicate that SW elastography is capable of reliably assessing differences in stiffness, but the absolute values of stiffness should be used with caution. Measuring changes in SW speed for different magnitudes of compression enables the quantification of the tissue’s nonlinear behavior which can significantly enhance the diagnostic value of SW elastography.     

Detecting coral bleaching using high-resolution satellite data analysis and 2-dimensional thermal model simulation in the Ishigaki fringing reef,Japan

In 2007, high-temperature-induced mass coral mortality was observed in a well-developed fringing reef area on the southeastern coast of Ishigaki Island, Japan. To analyze the response of the corals to thermal stress, the coral cover was examined using Quickbird data, taken across the reef flat just before and after the bleaching event and performing a reef scale horizontal 2-dimensional thermal model simulation. The Quickbird data consisted of multispectral (MSS) imagery, which had a spatial resolution of 2.4 m, and panchromatic (PAN)-fused multispectral imagery, which had a 0.6-m spatial resolution. The observed changes in coral cover implied that the delineation of partially bleached coral was more precise with PAN + MSS. The classification accuracy achieved using PAN + MSS (93%) was superior to that obtained using MSS (88%). The in situ water temperature observations and 2-dimensional thermal model simulation results indicated that the water temperature fluctuated greatly in the inner reef area in late July 2007. Different thermal stress indices, including daily average temperature, daily maximum excess temperature, and daily accumulated temperature, were examined to define a suitable index that represented the severity of the thermal stress on coral cover. The results suggested that the daily accumulated temperature that occurred during the maximum sea surface temperature period of the bleaching season provided the best predictor of bleaching. The changes in water temperature, bathymetry, and coral patch size affected the severity of bleaching; therefore, the spatial dependence of these variables was examined using Moran’s I and Lagrange multiplier tests. An investigation of the effect of coral patch sizes on coral bleaching indicated that large coral patches were less affected than the small patches, which were more likely to suffer bleaching and coral mortality.     

Actin stress fibres and cell-cell adhesion molecules in tendons: organisation in vivo and response to mechanical loading of tendon cells in vitro.

Tendons consist of parallel longitudinal rows of cells separated by collagen fibres. The cells are in intimate contact longitudinally within rows, and laterally via sheet-like lateral cell processes between rows. At points of contact, they are linked by gap junctions. Since tendons stretch under load, such cell contacts require protection. Here we describe the organisation of the actin cytoskeleton and actin-based cell-cell interactions in vivo and examine the effect of cyclic tensile loading on tendon cells in vitro. Cells within longitudinal rows contained short longitudinally running actin stress fibres. Each fibre was aligned with similar fibres in the cells longitudinally on either side, and fibres appeared to be linked via adherens junctions. Overall, these formed long oriented rows of stress fibres running along the rows of tendon cells. In culture, junctional components n-cadherin and vinculin and the stress fibre component tropomyosin increased in strained cultures, whereas actin levels remained constant. These results suggest that: (1) cells are linked via actin-associated adherens junctions along the line of principal strain; and (2) under load, cells appear to attach themselves more strongly together, and assemble more of their cytoplasmic actin into stress fibres with tropomyosin. Taken together, this suggests that cell-cell contacts are protected during stretch, and also that the stress fibres, which are contractile, may provide an active mechanism for recovery from stretch. In addition, stress fibres are ideally oriented to monitor tensile load and thus may be important in mechanotransduction and the generation of signals passed via the gap junction network.     

Load sharing between solid and fluid phases in articular cartilage: II--Comparison of experimental results and u-p finite element predictions.

Experimental measurements in conjunction with theoretical predictions were used to determine the extent of load supported by the fluid phase of cartilage at the articular surface. The u-p finite element model was used to simulate the loading of six separate porcine knee joints and to predict surface deformations of the cartilage layer on the lateral femoral condyle. Representative geometry for the condyle, contact pressures, and intrinsic material properties of the cartilage layer were supplied from experimental measures (see Part I). The u-p finite element predictions for surface deformations of the cartilage layer were obtained for several load partitioning states between the solid and fluid phases of cartilage at the articular surface. These were then compared to actual surface deformations obtained experimentally. It appeared from the comparison that approximately 75 percent of the applied load was borne by the fluid phase at the articular surface under this loading regime. This was qualitatively in agreement with the hypothesis that an applied load to articular joints is partitioned at the surface to the two phases according to the surface area ratios of the solid and fluid phases. It appeared that the solid phase was shielded from the total applied stress on the articular surface by the fluid and could be a reason for the excellent durability of the tissue under the demanding conditions in a diarthrodial joint.     

High strain rate response of rabbit femur bones

Strain rate dependence of the mechanical response of hard tissues has led to a keen interest in their dynamic properties. The current study attempts to understand the high strain rate characteristics of rabbit femur bones. The testing was conducted using a split-Hopkinson pressure bar equipped with a high speed imaging system to capture the fracture patterns. The bones were also characterized under quasi-static compression to enable comparison with the high strain rate results. The quasi-static compressive moduli of the epiphyseal and diaphyseal regions were measured to be in the range of 2–3 and 5–7 GPa, respectively. Under high strain rate loading conditions the modulus is observed to increase with strain rate and attains values as high as 15 GPa for epiphyseal and 30 GPa for diaphyseal regions of the femur. The strength at high strain rate was measured to be about twice the quasi-static strength value. A large number of small cracks initiated on the specimen surface close to the incident bar. Coalescence of crack branches leading to fewer large cracks resulted in specimen fragmentation. In comparison, the quasi-static failure was due to shear cracking.