Search

Skeletal muscles are living tissues that can undergo large deformations in short periods of time and that can be activated to produce force. In this paper we use the principles of continuum mechanics to propose a dynamic, fully non-linear, and three-dimensional model to describe the deformation of these tissues. We model muscles as a fibre-reinforced composite and transversely isotropic material. We introduce a flexible computational framework to approximate the deformations of skeletal muscle to provide new insights into the underlying mechanics of these tissues. The model parameters and mechanical properties are obtained through experimental data and can be specified locally. A semi-implicit in time, conforming finite element in space scheme is used to approximate the solutions to the governing nonlinear dynamic model. We provide a series of numerical experiments demonstrating the application of this framework to relevant problems in biomechanics, and also discuss questions around model validation.

The propagation of internal gravity waves in stratified media, such as those found in ocean basins and lakes, leads to the development of geometrical patterns called "attractors". These structures accumulate much of the wave energy and make the fluid flow highly singular. In more analytical terms, the cause of this phenomenon has been attributed to the presence of a continuous spectrum in some nonlocal zeroth-order pseudo-differential operators. In this work, we analyze the generation of these attractors from a numerical analysis perspective. First, we propose a high-order pseudo-spectral method to solve the evolution problem (whose long-term behaviour is known to be not square-integrable). Then, we use similar tools to discretize the corresponding eigenvalue problem. Since the eigenvalues are embedded in a continuous spectrum, we compute them using viscous approximations. Finally, we explore the effect that the embedded eigenmodes have on the long-term evolution of the system., Comment: 23 pages, 16 figures, MATLAB codes available at http://www.github.com/javieralmonacid/zeroth-order-operators

The purpose of this study was to report numerical validation of a 3D finite element model of contracting muscle. The model was based on continuum theory for fibre-reinforced composite materials. Here we simulated contractions for an idealized medial gastrocnemius muscle in man, using the model. Simulations were performed to test the force-length relation of the whole muscle, to evaluate the changes in internal fascicle geometry during contractions, and to assess the importance of material formulations for the aponeurosis and tendon. The simulation results were compared to previously published experimental values. The force-length profile for the whole muscle showed a realistic profile. As the muscle contracted the fascicles curved into S-shaped trajectories and curled around 3D paths, both of which matched previous experimental findings. As the fascicles shortened they increased in their cross-sectional area, but this increase was asymmetric with the smaller increase occurring within the fascicle-plane: the Poisson's ratio in this plane matched that previously shown from ultrasound imaging. The distribution of strains in the aponeurosis and tendon was shown to be a function of their material properties. This study demonstrated that the model could replicate realistic patterns of whole muscle-force, and changes to the internal muscle geometry, and so will be useful for testing mechanisms that affect the structural changes within contracting muscle.

We present degrees of freedom to accompany the approximation spaces already presented in a companion paper and thus complete the definition of families of high-order conforming finite elements on pyramids for the spaces of the de Rham complex. We prove that the elements are unisolvent; are compatible with conventional tetrahedral and hexahedral elements; satisfy a commuting diagram property and contain high-degree polynomials. We also tabulate shape functions for each element., Comment: The two parts of this paper have been re-combined into one manuscript so part II has been withdrawn. The full version is now available at http://arxiv.org/abs/math/0610206