A physically motivated constitutive model for 3D numerical simulation of skeletal muscles

A chemo-mechanical constitutive model for muscle activation in bat wing skins

Birds, bats and insects have evolved unique wing structures to achieve a wide range of flight capabilities. Insects have relatively stiff and passive wings, birds have a complex and hierarchical feathered structure and bats have an articulated skeletal system integrated with a highly stretchable skin. The compliant skin of the wing distinguishes bats from all other flying animals and contributes to bats' remarkable, highly manoeuvrable flight performance and high energetic efficiency. The structural and functional complexity of the bat wing skin is one of the least understood although important elements of the bat flight anatomy. The wing skin has two unusual features: a discrete array of very soft elastin fibres and a discrete array of skeletal muscle fibres. The latter is intriguing because skeletal muscle is typically attached to bone, so the arrangement of intramembranous muscle in soft skin raises questions about its role in flight. In this paper, we develop a multi-scale chemo-mechanical constitutive model for bat wing skin. The chemo-mechanical model links cross-bridge cycling to a structure-based continuum model that describes the active viscoelastic behaviour of the soft anisotropic skin tissue. Continuum models at the tissue length-scale are valuable as they are easily implemented in commercial finite element codes to solve problems involving complex geometries, loading and boundary conditions. The constitutive model presented in this paper will be used in detailed finite element simulations to improve our understanding of the mechanics of bat flight in the context of wing kinematics and aerodynamic performance.

Mechanical effects of needle texture on acupoint tissue

OBJECTIVE

This study aims to clarify how the stimulation of acupuncture points is achieved by needles with different surface texture during acupuncture; it also seeks to lessen injury at the insertion site and increase the therapeutic efficacy of acupuncture, by simulating the mechanical effects of various needle surface patterns on Zusanli (ST36) without changing the radius of acupuncture needles.

METHODS

Five acupuncture needle models with different surface patterns, including the smooth needle, the lined needle, the ringed needle, the left-hand threaded needle and the right-hand threaded needle, and a layered model of the Zusanli acupoint were used to investigate how to reduce tissue damage and increase stimulation during acupuncture treatment. Puncturing of the skin as well as lifting-inserting and twisting needle manipulations were simulated using these models, and the degree of damage and force of stimulation caused by the acupuncture needles with different surface patterns during acupuncture were compared.

RESULTS

The smooth needle and the lined needle caused the least tissue damage during insertion, while the left-hand threaded and the right-hand threaded needles caused the most damage. The ringed needle, the left-hand threaded needle and the right-hand threaded needle stimulated the acupoint tissue more during lifting-inserting manipulations, while the lined needle and the smooth needle produced less stimulation. The stimulation of the lined needle on the acupoint tissue was the largest during twisting manipulation, whereas the left-hand threaded needle and the right-hand threaded needle had smaller effects. In lifting-inserting and twisting manipulations, both the left-hand threaded needle and right-hand threaded needle provided more stimulation, but the torsion direction in which they produced better stimulation was the opposite.

CONCLUSION

According to the simulation results, the ringed pattern enhances stimulation best in the lifting-inserting manipulation, whereas the lined pattern enhances stimulation best in the twisting manipulation. Both the right-hand and left-hand thread patterns have certain enhancing effects in these two operations. Taking the geometric properties of the pattern into account, the left-hand thread pattern and the right-hand thread pattern have the geometric characteristics of both the lined pattern and the ringed pattern. To conclude, a pattern perpendicular to the movement direction during the acupuncture manipulation creates more stimulation. These results have significance for future needle design. Please cite this article as: Sun MZ, Wang X, Li YC, Yao W, Gu W. Mechanical effects of needle texture on acupoint tissue. J Integr Med. 2023; Epub ahead of print.

Peak ependymal cell stretch overlaps with the onset locations of periventricular white matter lesions

Deep and periventricular white matter hyperintensities (dWMH/pvWMH) are bright appearing white matter tissue lesions in T2-weighted fluid attenuated inversion recovery magnetic resonance images and are frequent observations in the aging human brain. While early stages of these white matter lesions are only weakly associated with cognitive impairment, their progressive growth is a strong indicator for long-term functional decline. DWMHs are typically associated with vascular degeneration in diffuse white matter locations; for pvWMHs, however, no unifying theory exists to explain their consistent onset around the horns of the lateral ventricles. We use patient imaging data to create anatomically accurate finite element models of the lateral ventricles, white and gray matter, and cerebrospinal fluid, as well as to reconstruct their WMH volumes. We simulated the mechanical loading of the ependymal cells forming the primary brain-fluid interface, the ventricular wall, and its surrounding tissues at peak ventricular pressure during the hemodynamic cycle. We observe that both the maximum principal tissue strain and the largest ependymal cell stretch consistently localize in the anterior and posterior horns. Our simulations show that ependymal cells experience a loading state that causes the ventricular wall to be stretched thin. Moreover, we show that maximum wall loading coincides with the pvWMH locations observed in our patient scans. These results warrant further analysis of white matter pathology in the periventricular zone that includes a mechanics-driven deterioration model for the ventricular wall.

Skeletal muscle: Modeling the mechanical behavior by taking the hierarchical microstructure into account

Skeletal muscles ensure the mobility of mammals and are complex natural fiber-matrix-composites with a hierarchical microstructure. In this work, we analyze the muscle's mechanical behavior on the level of fascicles and muscle fibers. We introduce continuum mechanics hyperelastic material models for the connective tissue endomysium and the embedded muscle fibers. The coupled electrical, chemical and mechanical processes taking place in activated contracting muscle fibers are captured including the temporal change of the activation level and the spatial propagation of the activation potential in fibers. In our model, we investigate the material behavior of fascicle, fiber and endomysium in the fiber direction and examine interactions between muscle fiber and endomysium by considering the temporal and spatial change of muscle fiber activation. In addition, a loading case of normal and shear forces is applied to analyze the fiber lifting force and the lifting height of unipennate muscles with different pennation angles. Moreover, the development of local stresses and strains in fibers and endomysium for different strains are studied. The simulation results allow to identify regions in high risk of damage. Optimal arrangements of unipennate muscle microstructure are found for either very small or very large pennation angles.

A comprehensive and volumetric musculoskeletal model for the dynamic simulation of the shoulder function

We present a volumetric and extensive finite element model of the shoulder usable in the context of inverse control, in which the scapula is left unconstrained on the ribcage. Such a model allows for exploring various shoulder movements, which are essential for making patient-specific decisions. The proposed model consists of 23 volumetric muscles parts modelled using the finite element method. The glenohumeral, acromioclavicular and sternoclavicular joints are modelled with soft ball-socket constraints. The musculoskeletal model can be controlled by a tracking-based algorithm, finding the excitations values in the muscles needed to follow some target points. The moment arms obtained during abduction and rotation are compared with the literature, which includes results from cadaveric data and a fine FE model of the rotator cuff and the deltoid. We simulated the paralysis of serratus anterior, a main reason of scapular winging, and compared it with its physiological counterpart. A deficiency in the range of motion as well as a reduction in upward rotation were observed, which both corroborate clinical observations. This is one of the most comprehensive model of the shoulder, which can be used to study complex pathologies of the shoulder and their impact on functional outcome such as range-of-motion.

Experimental and Numerical Characterization of the Mechanical Masseter Muscle Response During Biting

Predictive simulations of the mastication system would significantly improve our understanding of temporomandibular joint (TMJ) disorders and the planning of cranio-maxillofacial surgery procedures. Respective computational models must be validated by experimental data from in vivo characterization of the mastication system's mechanical response. The present pilot-study demonstrates the feasibility of a combined experimental and numerical procedure to validate a computer model of the masseter muscle. An experimental setup is proposed that provides a simultaneous bite force measurement and ultrasound-based visualization of muscle deformation. The direct comparison of the experimentally observed and numerically predicted muscle response demonstrates the predictive capabilities of such anatomically accurate biting models. Differences between molar and incisor biting are investigated; muscle deformation is recorded for three different bite forces in order to capture the effect of increasing muscle fiber recruitment. The three-dimensional (3D) muscle deformation at each bite position and force-level is approximatively reconstructed from ultrasound measurements in five distinct cross-sectional areas (four horizontal and one vertical cross section). The experimental work is accompanied by numerical simulations to validate the predictive capabilities of a constitutive muscle model previously formulated. An anatomy-based, fully 3D model of the masseter muscle is created from magnetic resonance images (MRI) of the same subject. The direct comparison of experimental and numerical results revealed good agreement for maximum bite forces and masseter deformations in both biting positions. The present work therefore presents a feasible in vivo measurement system to validate numerically predicted masseter muscle contractions during mastication.

Topology optimization of fusiform muscles with a maximum contraction

Understanding the optimal designs in nature is critical in bionics. This paper presents a method for designing the configuration of fusiform muscle with a maximum contractile displacement based on topology optimization methods. A nearly incompressible continuum constitutive model of skeletal muscle is utilized. The contractile displacement from the relaxed state to the contracted state is regarded as the objective function. To handle the numerical difficulties that result from the existence of element density, an energy interpolation equation is employed, and a modification of the constitutive model of skeletal muscle is proposed. Several numerical examples are given to demonstrate the reasonability of the proposed method.

Human skeletal muscle behavior in vivo: Finite element implementation, experiment, and passive mechanical characterization

In this study, useful methods for active human skeletal muscle material parameter determination are provided. First, a straightforward approach to the implementation of a transversely isotropic hyperelastic continuum mechanical material model in an invariant formulation is presented. This procedure is found to be feasible even if the strain energy is formulated in terms of invariants other than those predetermined by the software's requirements. Next, an appropriate experimental setup for the observation of activation-dependent material behavior, corresponding data acquisition, and evaluation is given. Geometry reconstruction based on magnetic resonance imaging of different deformation states is used to generate realistic, subject-specific finite element models of the upper arm. Using the deterministic SIMPLEX optimization strategy, a convenient quasi-static passive-elastic material characterization is pursued; the results of this approach used to characterize the behavior of human biceps in vivo indicate the feasibility of the illustrated methods to identify active material parameters comprising multiple loading modes. A comparison of a contact simulation incorporating the optimized parameters to a reconstructed deformed geometry of an indented upper arm shows the validity of the obtained results regarding deformation scenarios perpendicular to the effective direction of the nonactivated biceps. However, for a valid, activatable, general-purpose material characterization, the material model needs some modifications as well as a multicriteria optimization of the force-displacement data for different loading modes.

Elastic-viscoplastic modeling of soft biological tissues using a mixed finite element formulation based on the relative deformation gradient

The characteristic highly nonlinear, time-dependent, and often inelastic material response of soft biological tissues can be expressed in a set of elastic-viscoplastic constitutive equations. The specific elastic-viscoplastic model for soft tissues proposed by Rubin and Bodner (2002) is generalized with respect to the constitutive equations for the scalar quantity of the rate of inelasticity and the hardening parameter in order to represent a general framework for elastic-viscoplastic models. A strongly objective integration scheme and a new mixed finite element formulation were developed based on the introduction of the relative deformation gradient-the deformation mapping between the last converged and current configurations. The numerical implementation of both the generalized framework and the specific Rubin and Bodner model is presented. As an example of a challenging application of the new model equations, the mechanical response of facial skin tissue is characterized through an experimental campaign based on the suction method. The measurement data are used for the identification of a suitable set of model parameters that well represents the experimentally observed tissue behavior. Two different measurement protocols were defined to address specific tissue properties with respect to the instantaneous tissue response, inelasticity, and tissue recovery.