Department of Orthopaedic Surgery, University of Pennsylvania, Philadelphia, 19104-6081, USA.
Biomech Model Mechanobiol. 2010 Oct;9(5):651-8. doi: 10.1007/s10237-010-0194-x. Epub 2010 Feb 11.
Modeling of connective tissues often includes collagen fibers explicitly as one of the components. These fibers can be oriented in many directions; therefore, several studies have considered statistical distributions to describe the fiber arrangement. One approach to formulate a constitutive framework for distributed fibers is to express the mechanical parameters, such as strain energy and stresses, in terms of angular integrals. These integrals represent the addition of the contribution of infinitesimal fractions of fibers oriented in a given direction. This approach leads to accurate results; however, it requires lengthy calculations. Recently, the use of generalized structure tensors has been proposed to represent the angular distribution in the constitutive equations of the fibers. Although this formulation is much simpler and fewer calculations are required, such structure tensors can only be used when all the fibers are in tension and the angular distribution is small. However, the amount of error introduced in these cases of non-tensile fiber loading and large angular distributions have not been quantified. Therefore, the objective of this study is to determine the range of values of angular distribution for which acceptable differences (less than 10%) between these two formulations are obtained. It was found, analytically and numerically, that both formulations are equivalent for planar distributions under equal-biaxial stretch. The comparison also showed, for other loading conditions, that the differences decrease when the fiber distribution is very small. Differences of less than 10% were usually obtained when the fiber distribution was very low (κ ≈ 0.03; κ ranges between 0 and 1/3, for aligned and isotropic distributed fibers, respectively). This range of angular distribution greatly limits the types of tissue that can be accurately analyzed using generalized structure tensors. It is expected that the results from this study guide the selection of a proper approach to analyze a particular tissue under a particular loading condition.
结缔组织的建模通常明确地将胶原蛋白纤维作为其中一个组成部分。这些纤维可以有许多方向的取向;因此,有几项研究考虑了统计分布来描述纤维排列。一种为分布纤维制定本构框架的方法是根据角积分来表示机械参数,如应变能和应力。这些积分表示在给定方向上取向的纤维的无穷小分数的贡献的总和。这种方法可以得到准确的结果;然而,它需要冗长的计算。最近,已经提出了广义结构张量的使用,以在纤维的本构方程中表示角分布。虽然这种表述形式简单得多,并且需要的计算也少,但只有当所有纤维都处于拉伸状态且角分布较小时,才能使用这种结构张量。然而,在这些非拉伸纤维加载和大角度分布的情况下引入的误差量尚未量化。因此,本研究的目的是确定在这两种表述形式之间可以获得可接受差异(小于 10%)的角分布范围。分析和数值研究发现,在等双轴拉伸下,两种表述形式对于平面分布都是等效的。比较还表明,对于其他加载条件,当纤维分布非常小时,差异会减小。当纤维分布非常低时(κ≈0.03;κ在 0 到 1/3 之间变化,分别对应于取向和各向同性分布的纤维),通常可以获得小于 10%的差异。这种角度分布范围极大地限制了可以使用广义结构张量准确分析的组织类型。预计本研究的结果将指导在特定加载条件下选择分析特定组织的适当方法。