Search

This position paper reviews current research topics in the study of athletic activities, which biomechanists would consider to be important or contentious issues, and focuses on critically evaluating what needs to be done in future research. It concludes that there remain many unresolved issues in the mechanics of athletic activities, many of which overlap with other disciplines. These issues relate to injury mechanisms, the control and co-ordination of movement, and ways of providing biomechanical feedback to enhance performance in athletic activities. Research to address these important issues will increasingly become more question than discipline orientated, must focus more on mechanisms than description, and will involve teams of researchers interacting on interdisciplinary problems.

In the paper by Martens et al. (J. Biomechonics 19.443454. 1986). the loading arrangement for impact bending strength testing of whole human femoral bones is presented diag ramatically in Fig. 1. The authors do not state and we cannot understand how their loading arrangement produced equal forces at the two loading and two support points (as they assumed) since the bones tested were not perfectly uniform along their length in either geometry or mechanical properties. Whilst theabovementioned points might haveaccounted for the 5/33 instances of distal failure reported, it was suggested by the authors that the failures outside the central portion were attributable to greater shear force in this region. We feel that the simple model used for cross-section in the estimate of shear stress did not take into account the considerable cancellous bone content of the distal femur, as shown in Fig. 3b. Furthermore, we wonder about the accuracy of their measurements of the total force applied to the bones. Figure 2 shows the maximum force 6289 N ofone test applied in a time interval of 8.6 ms; the deceleration of the loading fixtures would therefore have been very high. The load cell is shown in Fig. I to be measuring the total force applied to the bones minus the (unknown) inertia force required to decelerate the fixtures; i.e. the load applied to the bone might have been significantly underestimated.

Integrative biomechanics uses biomechanics knowledge and methods at multiple scales and among biological entities to address fundamental and clinical problems at the tissue and organ level. Owing to the large ranges of scale involved, integrative biomechanics is intrinsically multidisciplinary, extending from molecular biophysics to contemporary engineering descriptions of kinematics and bulk constitutive properties. Much of this integration is accomplished through multiscale models of the interactions of interest. Applications can range from the development of new biological knowledge to the creation of new technologies for clinical application. In this white paper, the historical background of, and the rationale behind, integrative biomechanics are reviewed, followed by a sampling of clinical advances that were developed using the integrative approach. Refinements of many of these advances are still needed, and unsolved problems remain, in genomic applications, developing improved interventional procedures and protocols, and personalized medicine. Challenges to achieve these goals include the need for better models and the acquisition and organization of the data needed to parameterize, validate and apply them. These challenges will be overcome, because the advances in characterizing disease risk, personalization of care, and therapeutics that will follow, demand that we continue to move forward in this exciting field.

In simulated or actual microgravity, human and animal postural muscles undergo substantial atrophy: after about 270 days, the muscle mass attains a constant value of about 70% of the initial one. Most animal studies reported preferential atrophy of slow twitch fibres whose mechanical properties change towards the fast type. However, in humans, at the end of a 42-days bed rest study, a similar atrophy of slow and fast fibres was observed. After microgravity, the maximal force of several muscle groups showed a substantial decrease (6–25% of pre-flight values). The maximal power during very short “explosive” efforts of 0.25–0.30 s showed an even greater fall, being reduced to 65% after 1 month and to 45% (of pre-flight values) after 6 months. The maximal power developed during 6–7 s “all-out” bouts on an isokinetic cycloergometer was reduced to a lesser extent, attaining about 75% of pre-flight values, regardless of the flight duration. In these same subjects, the muscle mass of the lower limbs declined by only 9–13%. Thus, a substantial fraction of the observed decreases of maximal power is probably due to a deterioration of the motor co-ordination brought about by the absence of gravity. To prevent this substantial decay of maximal absolute power, we propose that explosive exercise be added to the daily in-flight training schedule. We also describe a system aimed at reducing cardiovascular deconditioning wherein gravity is simulated by the centrifugal acceleration generated by the motion of two counter rotating bicycles ridden by the astronauts on the inner wall of a cylindrical space module. Finally, cycling on circular or elliptical tracks may be useful to reduce cardiovascular deconditioning in permanently manned lunar bases. Indeed, on the curved parts of the path, a cyclist generates an outward acceleration vector (ac). To counterbalance ac, the cyclist must lean inwards, so that the vectorial sum of ac plus the lunar gravity tends to the acceleration of gravity prevailing on Earth.

Understanding the mechanical response of the brain to external loadings is of critical importance in investigating the pathological conditions of this tissue during injurious conditions. Such injurious loadings may occur at high rates, for example among others, during road traffic or sport accidents, falls, or due to explosions. Hence, investigating the injury mechanism and design of protective devices for the brain requires constitutive modeling of this tissue at such rates. Accordingly, this paper is aimed at critically investigating the physical background for viscohyperelastic modeling of the brain tissue with scrutinizing the elastic fields pertinent to large, time dependent deformations, and developing a fully nonlinear multimode Maxwell model that can mathematically explain such deformations. The proposed model can be calibrated using the simple monotonic uniaxial deformation of the sample extracted from the tissue, and does not require additional information from relaxation or creep experiments. The performance of the proposed model is examined using the experimental results of two different studies, which reveals a desirable agreement. The usefulness, limitations, and future developments of the proposed model are discussed in this paper.

Ultrasound is a popular and affordable imaging modality, but the nature of freehand ultrasound operation leads to unknown applied loads at non-quantifiable angles. The purpose of this paper was to demonstrate an instrumentation strategy for an ultrasound system to measure probe forces and orientation during freehand imaging to characterize the interaction between the probe and soft-tissue as well as enhance repeatability. The instrumentation included a 6-axis load cell, an inertial measurement unit, and an optional sensor for camera-based motion capture. A known method for compensation of the ultrasound probe weight was implemented, and a novel method for temporal synchronization was developed. While load and optical sensing was previously achieved, this paper presents a strategy for potential instrumentation on a variety of ultrasound machines. A key feature was the temporal synchronization, utilizing the electrocardiogram (EKG) feature built-in to the ultrasound. The system was used to perform anatomical imaging of tissue layers of musculoskeletal extremities and imaging during indentation on an in vivo subject and an in vitro specimen. The outcomes of the instrumentation strategy were demonstrated during minimal force and indentation imaging. In short, the system presented robust instrumentation of an existing ultrasound system to fully characterize the probe force, orientation, and optionally its movement during imaging while efficiently synchronizing all data. Researchers may use the instrumentation strategy on any EKG capable ultrasound systems if mechanical characterization of soft tissue or minimization of forces and deformations of tissue during anatomical imaging are desired.

This paper describes two new methods for computational fracture analysis of human femur using Quantitative Computed Tomography (QCT) voxel-based finite element (FE) simulation. The paper also reports comprehensive mechanical testing for validation of the methods and evaluation of the required material properties. The analyses and tests were carried out on 15 human femurs under 11 different stance-type loading orientations. Several classical forms of subcapital, transcervical, basicervical, and intertrochanteric fractures plus a specific type of subtrochanteric fracture were created and analyzed. A new procedure was developed for prediction of the strengths and the fracture initiation patterns using a FE-based linear scheme. The predicted and observed fracture patterns were in correspondence, and the FE predictions of the fracture loads were in very good agreement with the experimental results. Moreover, the crack initiation and growth behaviors of two subtrochanteric fractures were successfully simulated through a novel implementation of the cohesive zone model (CZM) within a nonlinear FE analysis scheme. The CZM parameters were obtained through a series of experimental tests on different types of specimens and determination of a variety of material properties for different anatomic regions and orientations. The presented results indicated that the locations and patterns of crack initiation, the sequences of crack growth on different paths, and the compatibility of growth increments agreed very well with the observed specifications. Also, very good agreements were achieved between the measured and simulated fracture loads.

Understanding passive skeletal muscle mechanics is critical in defining structure–function relationships in skeletal muscle and ultimately understanding pathologically impaired muscle. In this systematic review, we performed an exhaustive literature search using PRISMA guidelines to quantify passive muscle mechanical properties, summarized the methods used to create these data, and make recommendations to standardize future studies. We screened over 7500 papers and found 80 papers that met the inclusion criteria. These papers reported passive muscle mechanics from single muscle fiber to whole muscle across 16 species and 54 distinct muscles. We found a wide range of methodological differences in sample selection, preparation, testing, and analysis. The systematic review revealed that passive muscle mechanics is species and scale dependent—specifically within mammals, the passive mechanics increases non-linearly with scale. However, a detailed understanding of passive mechanics is still unclear because the varied methodologies impede comparisons across studies, scales, species, and muscles. Therefore, we recommend the following: smaller scales may be maintained within storage solution prior to testing, when samples are tested statically use 2–3 min of relaxation time, stress normalization at the whole muscle level be to physiologic cross-sectional area, strain normalization be to sarcomere length when possible, and an exponential equation be used to fit the data. Additional studies using these recommendations will allow exploration of the multiscale relationship of passive force within and across species to provide the fundamental knowledge needed to improve our understanding of passive muscle mechanics.

Pregnant vehicle occupants experience relatively large acceleration when the vehicle passes a speed-bump. In this paper, the effect of such sudden acceleration on a pregnant uterus is investigated. A biomechanical model representing the fundamental dynamic behaviors of a pregnant uterus has been developed. The model relates to the 32nd week of gestation when the fetus is in head-down, occipito-anterior position. Considering the drag and squeeze effects of the amniotic fluid, we derive a comprehensive differential equation that represents the interaction of the uterus and fetus. Solving the governing equation, we obtain the system response to different speed-bump excitations. Using the fetal head injury criterion (HIC = 390), we evaluate the model response. Three risk zones (Low, Medium, and High) are introduced, and the effects of excitation characteristics on HIC are investigated. HIC enhances, sub-exponentially, as the excitation amplitude (width) increases (decreases). Three risk-bounds, corresponding to 25%, 75%, and 100% risk of injury, are developed in the “width-amplitude” and the “frequency-amplitude” planes. Considering a typical speed-bump of width and excitation amplitude of 0.5 m and 0.12 m, respectively, the driver should not hit the speed-bump at 42 km/h or more. We advise hitting such speed-bumps under 25 km/h, based on this paper's findings. According to the risk-bounds, the injury risk of an arbitrary speed-bump excitation, at any desired vehicle speed, can be determined. The findings can help to understand how a pregnant uterus and fetus are subjected to risk caused by a vehicle passing a speed-bump and to expand our knowledge to improve safety during pregnancy.

Spinal loads are recognized to play a causative role in back disorders and pain. Knowledge of lumbar spinal loads is required in proper management of various spinal disorders, effective risk prevention and assessment in the workplace, sports and rehabilitation, realistic testing of spinal implants as well as adequate loading in in vitro studies. During the last few decades, researchers have used a number of techniques to estimate spinal loads by measuring in vivo changes in the intradiscal pressure, body height, or forces and moments transmitted via instrumented vertebral implants. In parallel, computational models have been employed to estimate muscle forces and spinal loads under various static and dynamic conditions. Noteworthy is the increasing growth in latter computational investigations. This paper aims to review, compare and critically evaluate the existing literature on in vivo measurements and computational model studies of lumbar spinal loads to lay the foundation for future biomechanical studies. Towards this goal, the paper reviews in separate sections models dealing with static postures (standing, sitting, lying) as well as slow and fast dynamic activities (lifting, sudden perturbations and vibrations). The findings are helpful in many areas such as work place safety design and ergonomics, injury prevention, performance enhancement, implant design and rehabilitation management.

The paper focuses on the connection between elastic and electrical properties of cortical bone. Both these properties are governed by microstructure that consists of several pore systems filled with mechanically soft and electrically conductive tissue. Microstructural changes induced by aging, various diseases, microgravity conditions etc. lead to variation in both properties. The paper address the problem of evaluation of the changes in mechanical performance (decrease in Young’s moduli) via monitoring electrical conductivity. The theoretical results are verified experimentally.

The pathological changes associated with the development of atherosclerotic plaques within arterial vessels result in significant alterations to the mechanical properties of the diseased arterial wall. There are several methods available to characterise the mechanical behaviour of atherosclerotic plaque tissue, and it is the aim of this paper to review the use of uniaxial mechanical testing. In the case of atherosclerotic plaques, there are nine studies that employ uniaxial testing to characterise mechanical behaviour. A primary concern regarding this limited cohort of published studies is the wide range of testing techniques that are employed. These differing techniques have resulted in a large variance in the reported data making comparison of the mechanical behaviour of plaques from different vasculatures, and even the same vasculature, difficult and sometimes impossible. In order to address this issue, this paper proposes a more standardised protocol for uniaxial testing of diseased arterial tissue that allows for better comparisons and firmer conclusions to be drawn between studies. To develop such a protocol, this paper reviews the acquisition and storage of the tissue, the testing approaches, the post-processing techniques and the stress-strain measures employed by each of the nine studies. Future trends are also outlined to establish the role that uniaxial testing can play in the future of arterial plaque mechanical characterisation.

Contains fulltext : 137140.pdf (Publisher’s version ) (Open Access) Deformations of the atherosclerotic vascular wall induced by the pulsating blood can be estimated using ultrasound strain imaging. Because these deformations indirectly provide information on mechanical plaque composition, strain imaging is a promising technique for differentiating between stable and vulnerable atherosclerotic plaques. This paper first explains 1-D radial strain estimation as applied intravascularly in coronary arteries. Next, recent methods for noninvasive vascular strain estimation in a transverse imaging plane are discussed. Finally, a compounding technique that our group recently developed is explained. This technique combines motion estimates of subsequently acquired focused ultrasound images obtained at various insonification angles. However, because the artery moves and deforms during the multi-angle acquisition, errors are introduced when compounding. Recent advances in computational power have enabled plane wave ultrasound acquisition, which allows 100 times faster image acquisition and thus might resolve the motion artifacts. In this paper the performance of strain imaging using plane wave compounding is investigated using simulations of an artery with a vulnerable plaque and experimental data of a two-layered vessel phantom. The results show that plane wave compounding outperforms 0 degrees focused strain imaging. For the simulations, the root mean squared error reduced by 66% and 50% for radial and circumferential strain, respectively. For the experiments, the elastographic signal-to-noise and contrast-to-noise ratio (SNRe and CNRe) increased with 2.1dB and 3.7dB radially, and 5.6dB and 16.2dB circumferentially. Because of the high frame rate, the plane wave compounding technique can even be further optimized and extended to 3D in future.

For nearly a century, articular cartilage has been known for its exceptional tribological properties. For nearly as long, there have been research efforts to elucidate the responsible mechanisms for application toward biomimetic bearing applications. It is now widely accepted that interstitial fluid pressurization is the primary mechanism responsible for the unusual lubrication and load bearing properties of cartilage. Although the biomechanics community has developed elegant mathematical theories describing the coupling of solid and fluid (biphasic) mechanics and its role in interstitial lubrication, quantitative gaps in our understanding of cartilage tribology have inhibited our ability to predict how tribological conditions and material properties impact tissue function. This paper presents an analytical model of the interstitial lubrication of biphasic materials under migrating contact conditions. Although finite element and other numerical models of cartilage mechanics exist, they typically neglect the important role of the collagen network and are limited to a specific set of input conditions, which limits general applicability. The simplified approach taken in this work aims to capture the broader underlying physics as a starting point for further model development. In agreement with existing literature, the model indicates that a large Peclet number, Pe, is necessary for effective interstitial lubrication. It also predicts that the tensile modulus must be large relative to the compressive modulus. This explains why hydrogels and other biphasic materials do not provide significant interstitial pressure under high Pe conditions. The model quantitatively agrees with in-situ measurements of interstitial load support and the results have interesting implications for tissue engineering and osteoarthritis problems. This paper suggests that a low tensile modulus (from chondromalacia or local collagen rupture after impact, for example) may disrupt interstitial pressurization, increase shear stresses, and activate a condition of progressive surface damage as a potential precursor of osteoarthritis.

The relative rigid body motions between the femur and the tibia (termed tibiofemoral kinematics) during flexion activities provide an objective measure of knee function. However, to obtain physiologic tibiofemoral kinematics, coordinate systems must meet certain requirements which neither have been explicitly stated nor in fact satisfied in any previous publication known to the author. Thus the primary purpose of this paper was to explain these requirements. The only valid method for obtaining physiologic tibiofemoral kinematics is to insure that coordinate systems meet four necessary requirements which are 1) the joint coordinate system must follow the convention of Grood and Suntay ( Grood and Suntay, 1983 ), 2) the body-fixed axes must correspond to the functional FE axis of the femur and the functional IE axis of the tibia, 3) the origins of the Cartesian coordinate systems, which serve as a stepping stone for computing translations, must lie on the functional body-fixed axes, and 4) the reference position between the tibia and femur must correspond to neutral with the knee in extension and the change in the components of the vector locating the origins of the Cartesian coordinate systems from the neutral position must be used to determine translations by projecting changes in vector components along the joint coordinate system axes. Neither the paper by Grood and Suntay nor the ISB recommendation ( Wu et al., 2002 ) which adopted the joint coordinate system of Grood and Suntay explains these requirements. Indeed meeting these requirements may not be possible following the ISB recommendation thus indicating the need for revision to this recommendation.

Data analysis can be the most challenging aspect of a research study. Having been taught statistical techniques that tend to be based on finding significant differences or significant relationships, difficulties arise when trying to determine if a newly developed method is equally as good as the established method (the gold standard). Testing for significant differences is rigorous and it would be rare for researchers to report significant differences without using an appropriate statistical test. Testing for agreement is assessed with far less rigour. Analysis of papers in this journal suggests that testing for agreement is an area that could be improved by a better understanding of statistical methods by biomechanics researchers. This perspectives paper focusses on informing the reader about the assessment of agreement between two methods.

Characterising soft biological tissues outside the linear viscoelastic regime is challenging due to their complex behaviour. In addition, the viscoelastic properties of tissues have been shown to be sensitive to sample preparation and loading regime resulting in inconsistent data varying by orders magnitude in the literature. This paper presents a novel technique to characterise the non-linear behaviour of tissues which uses Fourier Transformation to decompose the stress output waveform under large amplitude oscillatory shear (LAOS) into harmonic contributions. The effect of varying preload, the compressive strain exerted on a liver tissue specimen prior to shear testing to minimise slip, was also investigated. Results showed that in the linear regime, preload affects the viscoelastic response of liver. Histological analysis indicated that there were structural changes as a result of the preload that may be linked to the differences in observed behaviour. Fourier analysis was used to extract the first and third harmonic components of the shear moduli at large strain. At 50% shear strain, a change in the third harmonic component of the shear moduli was accompanied by a marked change in the micro-structural arrangement of the sinusoids. This paper demonstrates a method of efficiently characterising soft biological tissues under large amplitude oscillatory shear under combined loading.

This study examined the effect of friction between the hand and grip surface on a person's grip strategy and force generation capacity. Twelve young healthy adults performed power grip exertions on an instrumented vertical cylinder with the maximum and 50% of maximum efforts (far above the grip force required to hold the cylinder), while normal and shear forces at each phalanx of all five fingers in the direction orthogonal to the gravity were recorded. The cylinder surface was varied for high-friction rubber and low-friction paper coverings. An increase in surface friction by replacing the paper covering with the rubber covering resulted in 4% greater mean phalanx normal force (perpendicular to the cylinder surface) and 22% greater mean phalanx shear force in either the proximal or distal direction of the digits (p

The interaction between the flowing blood and deforming arterial wall is critical in understanding the role of hemodynamic forces such as wall shear stress (WSS) in atherosclerosis. Numerical methods have been extensively used to understand the nature of flow around atherosclerosis susceptible regions of the vascular tree in order to establish the exact role of WSS in atherosclerosis. Unfortunately, most of the numerical studies have been performed on rigid arterial geometries, which do not take into account the effect of the interaction between the flowing blood and the dynamics of the flexible arterial wall. In vivo , blood vessels are continuously deforming with every contraction and relaxation of the heart during the cardiac cycle. This paper forms the first of the two-part paper series discussing the need for fluid-structure interaction (FSI) in hemodynamic WSS analysis. The paper presents a well validated FSI based numerical model, capable of accurately predicting flow transients in arteries. The numerical model is validated using analytical solutions and experiments conducted on polyurethane mock artery, with the numerical predictions, analytical solutions and experimental data comparing very well. Numerical studies are performed using OpenFOAM, a 3D Finite Volume Method(FVM) based C++ library. Deposited by bulk import

Accurate spatial location of joint center (JC) is a key issue in motion analysis since JC locations are used to define standardized anatomical frames, in which results are represented. Accurate and reproducible JC location is important for data comparison and data exchange. This paper presents a method for JC locations based on the multiple regression algorithms without preliminary assumption on the behavior of the joint-of-interest. Regression equations were obtained from manually palpable ALs on each bone-of-interest. Results are presented for all joint surfaces found on the clavicle, scapula and humeral bone. Mean accuracy errors on the JC locations obtained on dry bones were 5.2+/-2.5 mm for the humeral head, 2.5+/-1.1 mm for the humeral trochlea, 2.3+/-0.9 mm for the humeral capitulum, 8.2+/-3.9 mm for the scapula glenoid cavity, 7.2+/-3.2 mm for the scapular aspect of the acromio-clavicular joint, 3.5+/-1.8mm for the clavicular aspect of the sternoclavicular joint and 3.2+/-1.4 mm for the clavicular aspect of the acromio-clavicular joint. In-vitro and in-vivo validation accuracy was 5.3 and 8.5mm, respectively, for the humeral head center location. Regression coefficients for joint radius dimension and joint surface orientation were also processed and reported in this paper.

Standard registration techniques of bone morphology to motion analysis data often lead to unsatisfactory motion simulation because of discrepancies during the location of anatomical landmarks in the datasets. This paper describes an iterative registration method of a three-dimensional (3D) skeletal model with both 6 degrees-of-freedom joint kinematics and standard motion analysis data. The method is demonstrated in this paper on the lower limb. The method includes two steps. A primary registration allowed synchronization of in vitro kinematics of the knee and ankle joints using flexion/extension angles from in vivo gait analysis. Results from primary registration were then improved by a so-called advanced registration, which integrated external constraints obtained from experimental gait pre-knowledge. One cadaver specimen was analyzed to obtain both joint kinematics of knee and ankle joints using 3D electrogoniometry, and 3D bone morphology from medical imaging data. These data were registered with motion analysis data from a volunteer during the execution of locomotor tasks. Computer graphics output was implemented to visualize the results for a motion of sitting on a chair. Final registration results allowed the observation of both in vivo motion data and joint kinematics from the synchronized specimen data. The method improved interpretation of gait analysis data, thanks to the combination of realistic 3D bone models and joint mechanism. This method should be of interest both for research in gait analysis and medical education. Validation of the overall method was performed using RMS of the differences between bone poses estimated after registration and original data from motion analysis.

Codman's paradox refers to a specific pattern of motion at the shoulder joint. It asked how a mysterious axial rotation about the longitudinal axis of the arm occurred during two or three sequential arm rotations that did not involve rotation about the long-axis. The objective of this paper was to find how the mysterious axial rotation occurred in the Codman's paradox. First, Codman's paradox and Codman's rotation were defined in general situations that involved arm rotations about orthogonal axes starting from the neutral attitude as well as rotations about non-orthogonal axes and starting from non-neutral attitudes. Then a general law of motion was proposed to answer the question of the Codman's paradox, which is stated as when the long-axis of the arm performs a closed-loop motion by three sequential rotations defined as Codman's rotation, it produces an equivalent axial rotation angle about the long-axis. The equivalent axial rotation angle equals the angle of swing-the second rotation in the three sequential long-axis rotations. Validity of the proposed law of motion is demonstrated by computer simulation of various Codman's rotations. Clinical relevance of the proposed law of motion is also discussed in the paper.

The aim of this paper was to introduce the idea of topology optimization in engineering to the simulation of bone morphology. The external shape of bone structure was predicted with the quantitative bone functional adaptation theory. The high-order nonlinear equation of bone remodeling proposed by Zhu et al. (J. Biomech. 35(7)(2002)951) combining with the finite element method was adopted to a rectangular design domain, which occupies a larger space than the external shape of bone structure. It was at this point we imported the idea of topology optimization in engineering. The proximal femur was used here as an example, whose external shape and internal density distribution were simultaneously simulated quantitatively to validate that the external shape of bone structure could be successfully predicted in this way. Then the growth of vertebral body from young to old was simulated numerically in its coronal section to discuss the significance of the prediction of external shape. The study in this paper provides computational basis for further studies on osteophyte formation, osteoporosis, osteoarthritis, bone growth and even bone evolution, etc.

When a mechanical load is applied very rapidly to the thoracic wall, part of the internal damage is suspected to be due to a “high-frequency” injury mechanism, that is, a phenomenon in which waves are involved. This paper addresses a specific high-frequency mechanism for lung injury in which a stress wave is generated through rapid acceleration of the body wall. Displacement-related injuries, which are rather “low-frequency” phenomena, are not considered. The present work was done in the context of assessing behind armor blunt trauma (injury to thoracic organs occurring when a bullet is stopped by a body armor) through mathematical modeling. One aspect of the thorax response to high-speed blunt impact and an associated injury mechanism are investigated based on an idealized model of thorax and a set of computations presented in previous papers. The injury mechanism considered elucidates a possible mathematical relationship between the acceleration at the surface of the thoracic wall and the occurrence of lung injury.

Motion capture systems currently used in biomechanical analysis introduce systematic measurement errors that appear in the form of noise in recorded displacement signals. The noise is unacceptably amplified when differentiating displacements to obtain velocities and accelerations. To avoid this phenomenon, it is necessary to smooth the displacement signal prior to differentiation in order to eliminate the noise introduced by the experimental system. The use of singular spectrum analysis (SSA) is presented in this paper as an alternative to traditional digital filtering methods. SSA is a novel non-parametric technique based on principles of multivariate statistics. The original time series is decomposed into a number of additive time series, each of which can be easily identified as being part of the modulated signal, or as being part of the random noise. Several examples that demonstrate the superiority of this technique over other methods used in biomechanical analysis are presented in this paper.

Dentin constitutes the major part of human tooth. It is composed of a large number of tubules with both variational radii and radially parallel pattern. In addition, peritubular dentin surrounds each tubule lumen and has a higher elastic modulus than the matrix of dentin, i.e. intertubular dentin. Considering the above microstructural characteristics, a micromechanics model is used in this paper to evaluate the overall elastic properties of dentin. Five independent effective elastic parameters in transverse isotropic elasticity matrix can be expressed analytically by the material parameters of peri- and intertubular dentin and the volume fraction of tubules. To determine the effectivity of this theoretical model, a finite element (FE) model simulating a longitudinal tooth slice in moire fringe testing of Wang and Weiner (J. Biomech. 31 (1998) 135) was performed. Furthermore, the FE model was developed incorporate modeling of variation of tubule's diameter and softer characteristic of intertubular dentin near the dentin-enamel junction and around the pulp chamber. It turned out that the isoline figure of longitudinal displacement by FE calculation has very similar patterns to the moire fringe results. However, the FE results of displacement by traditional stress-strain models which regard dentin as a homogeneous and isotropic material show an obviously different strain distributions as compared to published moire fringes results. Thus the inhomogeneous and anisotropic model developed in this paper more accurately reflects the true physical nature of human dentin.

Underwater shock can produce very high accelerations, resulting in severe human injuries. In this paper, a shock-structure-human interaction model is proposed to study the biodynamic response of a shipboard sitting subject to ship motion induced by underwater shock (ship shock motion) wherein, the human body is modeled using a lumped parameter system with the parameters obtained from dynamic tensile tests. The results obtained from the human model used in this paper and living human drop test are also compared. Numerical results have revealed the characteristics of human response to ship shock motion. The part in direct contact with the structure (like the pelvis) is much more vulnerable than other parts (like the head). The influences of structural damping and stiffness on the peak loads acting on the human body are investigated. Both damping and stiffness have important influences on the pelvis, but have much less influences on other parts. Injury criteria in the literature are also summarized to facilitate injury assessment.

The purpose of this paper is to introduce some concepts of low back injury for use towards developing better injury risk reduction strategies and advancing rehabilitation of the injured spine. Selected issues in low back injury are briefly reviewed and discussed, specifically, the types of tissue loads that cause low back injury, methods to investigate tissue loading, and issues which are important considerations when formulating injury avoidance strategies such as spine posture, and prolonged loading of tissues over time. Finally, some thoughts on current practice are expressed to stimulate discussion on directions for injury reduction efforts in the future, particularly, the way in which injuries are reported, the use of simple indices of risk such as load magnitude, assessment of the injury and development of injury avoidance strategies. This paper was written for a general biomechanics audience and not specifically for those who are spine specialists.

This paper discusses the manner in which random Gaussian errors affect the determination of body segment kinematics. For the process of modelling rigid body (RB) motion, three types of kinematic errors, input, measured and theoretical, are identified. These correspond to errors in: the determination of three-dimensional observed points, the RB fit of those points, and the estimation of true RB positions, respectively. Of these, the theoretical error is most critical and most pivotal. Accuracy is provided when the theoretical error is minimised, yet only the measured error can be minimised by RB modelling algorithms. In computer simulations one may determine the effect that such manipulations have on theoretical error, yet in most experimental conditions this value may not even be calculated. Fortunately, computer simulations can be performed to determine the inter-relationships between types of RB modelling errors. Such simulations can also be used to investigate the effects of RB shape. In this paper, Monte Carlo simulations were performed on three unit radius RBs; a triangle, a square and a tetrahedron. Although the use of the triangle provided the lowest measured error, this also coincided with the greatest theoretical error. The use of redundant points was found to yield superior theoretical accuracies. A slight advantage was gained with use of the non-planar point arrangement on the tetrahedron, both the measured and theoretical errors were reduced. Finally, the superiority of RB modelling over individual point tracking was reflected in all of the results; between 33 and 50% of the input error was eliminated with the use of RB modelling.

A mathematical model of lubrication of human synovial joints under squeeze-film conditions is presented in this several-part paper. Squeeze-film action leads to a concentration of hyaluronic-acid-protein macromolecular complex in the synovial fluid between the approaching cartilage surfaces as a result of the diffusion of water and low molecular weight substances through the cartilage surfaces or along the gap. Increasing viscosity of synovial fluid delays the approach of these surfaces and the formation of stable gels then protects cartilage, if sliding motion ensues, before fluid film lubrication is restored. In Part I of the present paper synovial fluid is considered as a mixture of two incompressible fluids. The material parameters of this mixture of fluids are found using previously published experimental results. Squeeze-film analysis is carried out for the axially symmetric synovial film.

There are many diseases which cause detrimental changes in the trabecular structure of cancellous bone, leading to mechanical failure of the tissue. One approach to understanding the mechanisms of these diseases is to create idealized models that recreate the morphology of the tissue. This paper presents a partial development of such a model. Further histological methods must be developed before a complete definition of morphologically valid models is possible. In a histological section of cancellous bone, the orientation and length of the trabecular surfaces determine how a line drawn across the bone section will intersect the bone-marrow interface. The distribution of the average length between intersections for a set of parallel lines is defined as the mean intercept length distribution. In this paper, the average surface morphology and volume of the average structure of cancellous bone is determined from an examination of the mean intercept length. The average structure of cancellous bone contains a repeated structural element (SE). As a result, the basic bone structure is analogous to a brick wall made from many similar bricks. For a group of 107 specimens, a strong relationship between structural element volume (SE.V) and bone volume fraction (BV/TV) is demonstrated, SE.V = 0.017 kappa (BV/TV)-2.05 mm3, R2 = 0.93, with kappa a model-dependent constant. For the same specimens, the structural element surface (SE.S) showed the relationship, SE.S = 0.144 kappa (BV/TV)-1.35, R2 = 0.92. As a result of the inverse square dependence of structural element volume on bone volume fraction, it is predicted that cancellous bone strength is inversely proportional to structural element volume.

Growth across the lifespan of males and females produces change in the form and structure of the human body which must be recognized in biomechanical models. The different morphologies for age span, gender and racial groupings have been identified and quantified by anthropometrists working primarily with measuring sticks and tapes. These parameters are not the parameters used to solve the dynamics problems of human movement. Models of the body composed of joints and rigid links have been formed and the inertial properties of the segments are required. The research which has been conducted to estimate these segment inertias is incomplete, dispersed, often flawed and inconclusive. However, it is essential that we have a reasonable understanding of the differences in inertia between individuals and changes of inertia within the individual regardless of age, sex and race. Segment inertia parameters can provide valuable insight into growth and the differences between individuals and populations. Segment morphology ranges from embryonic to excessive and each stage of growth or development has potential for the investigation of relationships between shape, size, inertia, weight and movement. The purpose of this paper is to examine research into segment inertias which provides useful information about the differences and changes to be expected. Although some guidance is provided, particularly with respect to bench mark studies, the paper is not intended as a discussion of methodologies. The studies which are pertinent to changes and differences in inertia make it clear that there are many interesting differences between age spans and, within these, differences between males and females and between races. If, in place of the expediency of proportions, segment parameters were measured and integrated into all studies, we would have a more complete representation of morphological changes and differences.

The effects and interaction of the anatomical displacements in the human knee are a prerequisite to an accurate assessment and communication of the kinematic data. For the kinematic information to be used to improve diagnosis and treatment, and for better prosthetic design and installation, there must be clear, concise, and universal definitions of the displacements. In general, the displacements are defined as three translations and three rotations. In this paper, anatomic landmarks on the femur and on the tibia are used to define the locations and orientations of the six displacement axes; i.e. three translational and three rotational displacement axes. The most commonly accepted kinematic representation of the knee joint, in the literature, is a special geometry three-cylindric open chain in which the axes of the cylindric joints are defined according to the rotational displacement axes. The sequentially adjacent joint axes are assumed to not only intersect but to intersect at right-angles. The open chain permits a total of six degrees of freedom between a Cartesian reference frame attached to the femur and a Cartesian reference frame attached to the tibia. In this paper, the three rotational axes are shown to be skewed and off-set from each other, therefore, a three-cylindric open chain with skewed joint axes is proposed to measure the six displacements between the two reference frames. The authors believe that the proposed open chain is the most general to date and provides a more realistic representation of the displacements in the knee. To illustrate the significance of the reference frames on the interpretation of measured data, the anterior/posterior drawer is plotted against per cent gait cycle for three existing open chains and the proposed open chain.