The gastrointestinal (GI) tract is a continuous channel through the body that consists of the esophagus, the stomach, the small intestine, the large intestine, and the rectum. Its primary functions are to move the intake of food for digestion before storing and ultimately expulsion of feces from the rectum through the anal sphincter. The mechanical behavior of GI tissues thus plays a crucial role for GI function in health and disease. The mechanical properties are typically characterized by a constitutive biomechanical model, which is a mathematical representation of the relation between load and deformation in a tissue. Hence, validated biomechanical constitutive models are essential to characterize and simulate the mechanical behavior of the GI tract under physiological and pathological conditions. Numerous constitutive models have consequently been proposed over the past three decades, mainly inspired by work done in cardiovascular tissues. Here, a comprehensive review of these constitutive models is provided. This review is limited to studies where a model of the strain energy function is proposed to characterize the stress-strain relation of a GI tissue. Several needs are identified for more advanced modeling of GI biomechanics including: 1) Microstructural models that provide actual structure-function relations; 2) Validation of coupled electro-mechanical models accounting for active muscle contractions; 3) Human data under physiological and pathological conditions to develop and validate models. The findings from this review provide guidelines for using existing constitutive models as well as perspective and directions for future studies aimed at establishing new constitutive models for GI tissues.

This work presents a linear smoothing scheme over high-order triangular elements in the framework of a cell-based strain smoothed finite element method for two-dimensional nonlinear problems. The main idea behind the proposed linear smoothing scheme for strain-smoothed finite element method (S-FEM) is no subdivision of finite element cells to sub-cells while the classical S-FEM needs sub-cells. Since the linear smoothing function is employed, S-FEM is able to use quadratic triangular or quadrilateral elements. The modified smoothed matrix obtained node-wise is evaluated. In the same manner with the computation of the strain-displacement matrix, the smoothed stiffness matrix and deformation graident are obtained over smoothing domains. A series of benchmark tests are investigated to demonstrate validity and stability of the proposed scheme. The validity and accuracy are confirmed by comparing the obtained numerical results with the standard FEM using 2nd-order triangular element and the exact solutions.

The design of thin-walled structures is commonly based on the solutions of linear boundary-valued problems, formulated within well-developed theories for elastic plates and shells. However, in modern appliances, especially in MEMS design, it is necessary to take into account non-linear mechanical effects that become decisive for flexible elements. Among substantial non-linear effects, that significantly change the deformation properties of thin plates, are the effects of residual stresses caused by incompatibility of deformations, which inevitably arise during the manufacture of ultra-thin elements. The development of new methods of mathematical modeling of residual stresses and incompatible finite deformations in plates is the subject of the paper. To this end the local unloading hypothesis is used. It makes possible to define smooth fields of local deformations (inverse implant field) for the mathematical formalization of incompatibility. The main outcomes are field equations, natural boundary conditions and conservation laws, derived from least action principle and variational symmetries with account of implant field. The derivations were carried out in the framework of elasticity theory for simple materials and, over and above, within the Cosserat’s theory of a two-dimensional continuum. As illustrative examples the distributions of incompatible deformations in a circular plate are considered.

Polydimethylsiloxane (PDMS) is one of the most popular elastomers and has been used in different fields, especially in biomechanics research. Among the many interesting features of this material, its hyperelastic behaviour stands out, which allows the use of PDMS in various applications, like the ones that mimic soft tissues. However, the hyperelastic behaviour is not linear and needs detailed analysis, especially the characterization of shear strain. In this work, two approaches, numerical and experimental, were proposed to characterize the effect of shear strain on PDMS. The experimental method was implemented as a simple shear testing associated with 3D digital image correlation and was made using two specimens with two thicknesses of PDMS (2 and 4 mm). A finite element software was used to implement the numerical simulations, in which four different simulations using the Mooney-Rivlin, Yeoh, Gent, and Polynomial hyperelastic constitutive models were performed. These approaches showed that the maximum value of shear strain occurred in the central region of the PDMS, and higher values emerged for the 2 mm PDMS thickness. Qualitatively, in the central area of the specimen, the numerical and experimental results have similar behaviours and the values of shear strain are close. For higher values of displacement and thicknesses, the numerical simulation results move further away from experimental values.