Witzenburg Colleen M, Dhume Rohit Y, Lake Spencer P, Barocas Victor H
IEEE Trans Med Imaging. 2016 Jan;35(1):29-41. doi: 10.1109/TMI.2015.2453316. Epub 2015 Jul 7.
Variations in properties, active behavior, injury, scarring, and/or disease can all cause a tissue's mechanical behavior to be heterogeneous. Advances in imaging technology allow for accurate full-field displacement tracking of both in vitro and in vivo deformation from an applied load. While detailed strain fields provide some insight into tissue behavior, material properties are usually determined by fitting stress-strain behavior with a constitutive equation. However, the determination of the mechanical behavior of heterogeneous soft tissue requires a spatially varying constitutive equation (i.e., one in which the material parameters vary with position). We present an approach that computationally dissects the sample domain into many homogeneous subdomains, wherein subdomain boundaries are formed by applying a betweenness based graphical analysis to the deformation gradient field to identify locations with large discontinuities. This novel partitioning technique successfully determined the shape, size and location of regions with locally similar material properties for: (1) a series of simulated soft tissue samples prescribed with both abrupt and gradual changes in anisotropy strength, prescribed fiber alignment, stiffness, and nonlinearity, (2) tissue analogs (PDMS and collagen gels) which were tested biaxially and speckle tracked (3) and soft tissues which exhibited a natural variation in properties (cadaveric supraspinatus tendon), a pathologic variation in properties (thoracic aorta containing transmural plaque), and active behavior (contracting cardiac sheet). The routine enables the dissection of samples computationally rather than physically, allowing for the study of small tissues specimens with unknown and irregular inhomogeneity.
组织的属性、活性行为、损伤、瘢痕形成和/或疾病的变化都可能导致其力学行为的异质性。成像技术的进步使得能够精确地全场跟踪体外和体内在施加负荷下的变形。虽然详细的应变场能提供一些关于组织行为的见解,但材料属性通常是通过将应力-应变行为与本构方程拟合来确定的。然而,确定异质软组织的力学行为需要一个空间变化的本构方程(即材料参数随位置变化的方程)。我们提出了一种方法,通过计算将样本域分解为许多均匀的子域,其中子域边界是通过对变形梯度场应用基于中介中心性的图形分析来识别具有大不连续性的位置而形成的。这种新颖的划分技术成功地确定了具有局部相似材料属性的区域的形状、大小和位置,这些区域包括:(1) 一系列模拟软组织样本,其各向异性强度、规定的纤维排列、刚度和非线性既有突然变化也有逐渐变化;(2) 双轴测试并进行散斑跟踪的组织类似物(聚二甲基硅氧烷和胶原凝胶);(3) 具有属性自然变化的软组织(尸体冈上肌腱)、具有属性病理变化的软组织(含有透壁斑块的胸主动脉)以及活性行为(收缩心肌片)。该方法能够通过计算而非物理方式对样本进行剖析,从而允许对具有未知和不规则不均匀性的小组织标本进行研究。