Burd H J, Regueiro R A
Department of Engineering Science, Oxford University, Oxford, UK.
Department of Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, Boulder, CO, USA.
Biomech Model Mechanobiol. 2015 Nov;14(6):1363-78. doi: 10.1007/s10237-015-0680-2. Epub 2015 May 9.
An axisymmetric finite element implementation of a previously described structural constitutive model for the human lens capsule (Burd in Biomech Model Mechanobiol 8(3):217-231, 2009) is presented. This constitutive model is based on a hyperelastic approach in which the network of collagen IV within the capsule is represented by an irregular hexagonal planar network of hyperelastic bars, embedded in a hyperelastic matrix. The paper gives a detailed specification of the model and the periodic boundary conditions adopted for the network component. Momentum balance equations for the network are derived in variational form. These balance equations are used to develop a nonlinear solution scheme to enable the equilibrium configuration of the network to be computed. The constitutive model is implemented within a macroscopic finite element framework to give a multiscale model of the lens capsule. The possibility of capsule wrinkling is included in the formulation. To achieve this implementation, values of the first and second derivatives of the strain energy density with respect to the in-plane stretch ratios need to be computed at the local, constitutive model, level. Procedures to determine these strain energy derivatives at equilibrium configurations of the network are described. The multiscale model is calibrated against previously published experimental data on isolated inflation and uniaxial stretching of ex vivo human capsule samples. Two independent example lens capsule inflation analyses are presented.
本文提出了一种轴对称有限元方法,用于实现先前描述的人晶状体囊膜结构本构模型(Burd,《生物力学模型与分子生物力学》8(3):217 - 231,2009年)。该本构模型基于超弹性方法,其中囊膜内的IV型胶原网络由嵌入超弹性基质中的不规则六边形平面超弹性杆网络表示。本文详细说明了该模型以及网络组件所采用的周期性边界条件。以变分形式推导了网络的动量平衡方程。这些平衡方程用于开发一种非线性求解方案,以计算网络的平衡构型。本构模型在宏观有限元框架内实现,以给出晶状体囊膜的多尺度模型。公式中考虑了囊膜起皱的可能性。为实现这一实现,需要在局部本构模型层面计算应变能密度相对于面内拉伸比的一阶和二阶导数。描述了在网络平衡构型下确定这些应变能导数的过程。该多尺度模型根据先前发表的关于离体人囊膜样本的孤立膨胀和单轴拉伸的实验数据进行校准。给出了两个独立的晶状体囊膜膨胀分析示例。