全部

力学与实践 ›› 2018, Vol. 40 ›› Issue (3): 300-307.doi: 10.6052/1000-0879-17-286

• 应用研究 • 上一篇    下一篇

离散梯度法在基于图像的计算生物力学中的应用1)

钱璟*, 黄冠江*, 张莹雪,†,2   

  1. *深圳浥清科技有限公司,广东深圳518112;
    † 华侨大学信息科学与工程学院,福建厦门361021;
  • 收稿日期:2017-08-18 修回日期:2018-01-10 出版日期:2018-06-15 发布日期:2018-07-17
  • 作者简介:2) 张莹雪,讲师,主要研究方向为图像处理。E-mail: zyx@hqu.edu.cn

DISCRETE GRADIENT METHOD FOR IMAGE-BASED BIOMECHANICAL STRESS ANALYSIS 1)

QIAN Jing*, HUANG Guanjiang*, ZHANG Yingxue†,†,2   

  1. *Shenzhen Yiqing Technology Ltd, Shenzhen 518112, Guangdong, China;
    † Information Science and Engineering School, Huaqiao University, Xiamen 361021, Fujian, China;
  • Received:2017-08-18 Revised:2018-01-10 Online:2018-06-15 Published:2018-07-17
  • Contact: 1) 福建省高校产学合作资助项目(2017H6011)

摘要:

介绍一种可用于计算生物力学的离散梯度方法,此方法可利用离散的点云模型直接进行数值模拟分析而不需要传统的几何模型。将离散梯度法应用于点云模型需要首先确定模型中点之间的相邻关系和每个点所分配的材料体积,然后通过用广义的有限差分的形式定义了梯度插值向量,并以此向量来近似函数在每个离散点上的梯度。从弱形式出发,推导建立了适用于弹性固体大变形问题的求解器,并具有和有限元法中双线性四边形单元一致的准确性和收敛性。着重描述了一种可以从医学图像中快速提取材料点并建立点云模型的方法,以及利用三角划分和重心划分确定材料点之间的相邻关系和每个材料点体积的具体过程。通过腹主动脉瘤膨胀的静力学模拟分析,展示了离散梯度法的实用性和准确性。该算法实现了基于医学图像进行生物力学分析的过程自动化,为病体特异性的研究和治疗提供便利和实用的工具。

关键词:

生物力学|数值模拟|离散梯度|图像处理|点云模型

Abstract:

This paper proposes a discrete gradient method for the stress analysis for biological systems. A point-cloud is taken as the geometric input instead of the conventional CAD model. Before applying the discrete gradient method on the point-cloud model, the neighboring relationship is defined among points as well as the volume occupied by each point. The gradient interpolation vectors, which can approximate the gradient of a function, are defined for each point in the form of the generalized finite difference. The pointwise strain is calculated by using the discrete differentials involving the nodal displacements of a set of neighboring points. A mechanical solver using the discrete gradient method for finite strain elasticity is developed in weak form. It can be shown that this solver retains the similar accuracy and convergence rate of bilinear quadrilateral finite elements, with a locking-free behavior, and is more tolerant to the mesh distortion. An e±cient method is developed to extract the point-cloud model from medical images. Since a material constituent comprises pixels within a certain range of gray-scale values, the pixels within given thresholds are isolated and an initial point-cloud is formed. The physical coordinates are inferred from the image resolution. The Delaunay tessellation and the barycentric subdivision are utilized to provide the neighboring relation and the point volume. A static analysis of the abdominal aortic aneurysm inflation is carried out to demonstrate the usefulness of the method. Despite the use of the tessellation, the method is not element-based because the element-wise assumed solution is never constructed. A distinct feature of this method is that the entire process is completed with a minimal user interference or even fully automatically. This is significant for applications where a timely analysis is desired.

Key words:

biomechanics|stress analysis|discrete gradient method|medical images processing|point-cloud

中图分类号: 

  • Q66