非侵入式人体软组织损伤风险监测系统

本发明涉及人体软组织监测领域，具体涉及一种非侵入式人体软组织损伤风险监测系统。

背景技术：

由糖尿病及相关综合征引起的脊髓神经损伤(sci)、神经源性疾病、其他行动受限疾病会导致患有这些病症的人群在站立或者坐着时，软组织的骨突出部位长时间地受压迫作用，并且软组织无法感知所受内部应力，导致血流障碍而产生由深及浅的坏死。当体表皮肤组织发现病变时，患者软组织其实已经完全坏死，不可治愈或者患者虽然皮肤仍完整但软组织损伤已出现。现有的医疗条件，更多的是有创的医疗手段或者是对体表组织的监测，而缺乏对软组织的监测。从医学角度分析，软组织不同于体表组织的一大重要特征就是肌肉组织更加容易受到损伤，肌肉线条更加敏感，肌肉强度比较脆弱，这也决定了软组织不可能像体表组织那样接受相同的监测方法。通过非侵入式医疗手段对坐骨结节或足部跟骨处表面压力进行实时测量，再通过算法实现对软组织内部应力状态评估，达到监测的目的。

非侵入式医疗手段可以让患者获得无痛治疗和更好的舒适度，并降低住院率。非侵入式医疗手段可以准确、及时且按需实现患者监测。因此，对于软组织损伤的监测成为一种必要。

技术实现要素：

本发明是为了解决上述问题而进行的，目的在于提供一种非侵入式人体软组织损伤风险监测系统。

本发明提供了一种非侵入式人体软组织损伤风险监测系统，用于通过赫兹模型对人体软组织的损伤风险进行监测，具有这样的特征，包括：信息采集装置，用于通过医学影像和软组织测试装置，获取人体预定位置的软组织在不同压痕深度z和不同杨氏模量下的形态参数，软组织测试装置用于对杨氏模量进行控制变量；信息处理装置，包括第一信息处理部、第二信息处理部、第三信息处理部以及第四信息处理部，第一信息处理部用于对形态参数进行有限元分析，得到赫兹模型在不同的大变形条件下的校正因子α，第二信息处理部用于对赫兹模型进行修正，得到修正后赫兹模型，第三信息处理部用于通过修正后赫兹模型计算得到软组织内部应力σzz(t)，第四信息处理部用于将软组织内部应力σzz(t)与软组织的损伤阈值进行比较，对软组织目前的生存状态进行评价以及预测软组织未来出现损伤的风险；控制装置，用于控制信息采集装置进行信息采集，以及控制信息处理装置进行信息处理；其中，第一信息处理部包括第一信息处理部第一单元、第一信息处理部第二单元、第一信息处理部第三单元以及第一信息处理部第四单元，第一信息处理部第一单元用于通过设定压痕深度z和杨氏模量的条件构建一个可近似代表预定位置的软组织结构的三维有限元模型，从三维有限元模型的应力分布图中读出软组织的模拟内部应力σfe，第一信息处理部第二单元用于通过赫兹模型及其公式计算原始内部应力σh，第一信息处理部第三单元用于计算得到校正因子α，校正因子α等于模拟内部应力σfe与原始内部应力σh之间的商，第一信息处理部第四单元用于改变压痕深度z和杨氏模量的条件设定，通过第一信息处理部第一单元、第一信息处理部第二单元以及第一信息处理部第三单元得到多组对应不同压痕深度z和不同杨氏模量的校正因子α，并拟合出校正因子α在不同压痕深度z和不同杨氏模量下的对应关系。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，预定位置为坐骨结节或足部跟骨。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，医学影像为超声、ct或mri扫描。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，形态参数包括预定位置的曲率半径rit、弹性模量eit、泊松比vit、软组织原始厚度h，预定位置对应的肌肉组织的曲率半径rst、弹性模量est、泊松比vst，以及预定位置的末端半球骨的半径a。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，其中，第一信息处理部第一单元包括第一模块、第二模块以及第三模块，第一模块用于针对预定位置的末端半球骨搭建一个半球模型，半球模型垂直于刚性界面，半球模型和刚性界面的材料属性采用ogden模型，第二模块用于对半球模型逐步施加非线性几何形式的载荷，并对半球模型的表面的节点进行位移计算，并构建应力分布图，第三模块用于从应力分布图中读出软组织的模拟内部应力σfe。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，原始内部应力σh的计算公式为其中f为预定位置的表面压力，数值为人体所受重力的一半。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，修正后的赫兹模型及其计算公式如下，其中git和gst分别是预定位置和预定位置对应的肌肉组织的剪切模量。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，第三信息处理部计算软组织内部应力σzz(t)的具体过程为，对辅助形态参数a进行计算，辅助形态参数a的计算公式为对曲率进行计算，曲率的计算公式为对与时间有关的印迹半径a(t)进行计算，计算与时间有关的印迹半径a(t)的公式为其中f(t)是与时间相关的体表压力，与时间相关的体表压力f(t)通过设置在预定位置的体表接触区的压力传感器获得，对当量圆半径r进行计算，计算当量圆半径r的公式为对辅助形态参数进行计算，辅助形态参数的计算公式为，对与时间有关的肌骨接触区域上的应力峰值p0(t)进行计算，应力峰值的计算公式为对软组织内部应力σzz(t)进行计算，软组织内部应力σzz(t)的计算公式为：

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，软组织内部应力σzz(t)的计算方式在以下情况时，有例外的计算方法：z＝0,r＝0,σzz(t)＝p0(t)，z＝0,r＜a,z＝0,r≥a,σzz(t)＝0。

在本发明提供的一种非侵入式人体软组织损伤风险监测系统中，还可以具有这样的特征：其中，第四信息处理部将软组织内部应力σzz(t)与软组织的损伤阈值进行比较的具体过程为：若软组织内部应力σzz(t)达到了软组织的损伤阈值，则判定软组织受到挫伤；若软组织内部应力σzz(t)接近但未达到软组织的损伤阈值，则判定软组织有受到挫伤的风险；若软组织内部应力σzz(t)远未达到软组织的损伤阈值，则判定软组织没有受到挫伤的风险。

发明的作用与效果

根据本发明所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过赫兹模型和医疗影像所得的数据对人体预定位置的末端半球骨进行模拟，并结合设置在预定位置的体表接触区的压力传感器，对人体软组织进行损伤情况、损伤风险的检测，相比传统技术，实现了本发明所涉及的一种非侵入式人体软组织损伤风险监测系统在实现检测功能的情况下，没有对人体造成创伤。

根据本发明所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过赫兹模型对人体预定位置的末端半球骨进行模拟，并构建有限元模型对赫兹模型进行修正，保证了发明所涉及的一种非侵入式人体软组织损伤风险监测系统的准确性。

根据本发明所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过设置在预定位置的体表接触区的压力传感器实时对人体预定位置进行检测，从而实现对预定位置软组织的损伤情况以及损伤风险的实时判断，保证了发明所涉及的一种非侵入式人体软组织损伤风险监测系统的实时性。

附图说明

图1是本实施例中非侵入式人体软组织损伤风险监测系统的结构示意图；

图2是本发明的实施例中非侵入式人体软组织损伤风险监测系统的方法流程图；

图3是本发明的实施例中足部跟骨的医学扫描示意图；

图4是本发明的实施例中第一信息处理部的结构示意图；

图5是本发明的实施例中足部跟骨有限元模型的示意图；以及

图6是本发明的实施例中足部跟骨有限元模型的内部应力分布图。

具体实施方式

为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，以下实施例结合附图对本发明非侵入式人体软组织损伤风险监测系统作具体阐述。

图1是本实施例中非侵入式人体软组织损伤风险监测系统的结构示意图。如图1所示，非侵入式人体软组织损伤风险监测系统100包括信息采集装置10、信息处理装置20以及控制装置30。控制装置30用于控制信息采集装置10进行信息采集，以及控制信息处理装置20进行信息处理。

图2是本实施例中非侵入式人体软组织损伤风险监测系统的方法流程图，非侵入式人体软组织损伤风险监测系统用于通过赫兹模型对人体软组织的损伤风险进行监测。如图2所示，非侵入式人体软组织损伤风险监测系统的运行流程包括以下步骤：

步骤s1，信息采集装置10通过医学影像和软组织测试装置，获取人体预定位置的软组织在不同压痕深度z和不同杨氏模量下的形态参数。

医学影像为超声、ct或mri扫描，本实施例中采用ct。软组织测试装置为一种常见通用的医疗装置，用于对杨氏模量进行控制变量。预定位置为坐骨结节或足部跟骨，本实施例中的预定位置为足部跟骨。

形态参数包括预定位置的曲率半径rit、弹性模量eit、泊松比vit、软组织原始厚度h，预定位置对应的肌肉组织的曲率半径rst、弹性模量est、泊松比vst，以及预定位置的末端半球骨的半径a。

图3是本实施例中足部跟骨的医学扫描示意图。如图3所示，足部跟骨的骨突出部位近似一个半球形，可以作为经典赫兹半球接触模型的简化再现。从图3中可以直接得到的形态参数包括足部跟骨的半球骨末端的半径a，软组织原始厚度h。

信息处理装置20包括第一信息处理部21、第二信息处理部22、第三信息处理部23以及第四信息处理部24。

步骤s2，第一信息处理部21对形态参数进行有限元分析，得到赫兹模型在不同的大变形条件下的校正因子α。

图4是本实施例中的第一信息处理部的结构示意图。如图4所示，第一信息处理部21包括第一信息处理部第一单元211、第一信息处理部第二单元212、第一信息处理部第三单元213以及第一信息处理部第四单元214。

步骤s2包括以下子步骤：

步骤s2-1，第一信息处理部第一单元211通过设定压痕深度z和杨氏模量的条件，构建一个可近似代表预定位置的软组织结构的三维有限元模型，从三维有限元模型的应力分布图中读出软组织的模拟内部应力σfe。

第一信息处理部第一单元211包括第一模块2111、第二模块2112以及第三模块2113。

步骤s2-1包括以下子步骤：

步骤s2-1-1，第一模块2111针对预定位置的末端半球骨搭建一个半球模型。图5是本实施例中足部跟骨有限元模型的示意图。如图5所示，半球模型垂直于刚性界面，半球模型和刚性界面的材料属性采用ogden模型。

步骤s2-1-2，第二模块2112对半球模型逐步施加非线性几何形式的载荷，并对半球模型的表面的节点进行位移计算，并构建应力分布图。

图6是本实施例中足部跟骨有限元模型的内部应力分布图，步骤s2-1-2构建应力分布图如图6所示。

步骤s2-1-3，第三模块2113从应力分布图中读出软组织的模拟内部应力σfe。

步骤s2-2，第一信息处理部第二单元212经赫兹模型及其公式计算原始内部应力σh。原始内部应力σh的计算公式为

其中f为预定位置的表面压力，数值为人体所受重力的一半。

步骤s2-3，第一信息处理部第三单元213计算得到校正因子α，校正因子α等于模拟内部应力σfe与原始内部应力σh之间的商，

步骤s2-4，第一信息处理部第四单元214改变压痕深度z和杨氏模量的条件设定，重复步骤s2-1至步骤s2-3，得到多组对应不同压痕深度z和不同杨氏模量的校正因子α，并拟合出校正因子α在不同压痕深度z和不同杨氏模量下的对应关系。

步骤s3，第二信息处理部22对赫兹模型进行修正，得到修正后赫兹模型。修正后的赫兹模型及其计算公式如下：

其中git和gst分别是预定位置和预定位置对应的肌肉组织的剪切模量。

步骤s4，第三信息处理部23通过修正后赫兹模型计算得到软组织内部应力σzz(t)。计算软组织内部应力的具体过程为：

对辅助形态参数a进行计算，辅助形态参数a的计算公式为

对曲率进行计算，曲率的计算公式为

对与时间有关的印迹半径a(t)进行计算，计算与时间有关的印迹半径a(t)的公式为

其中f(t)是与时间相关的体表压力，与时间相关的体表压力f(t)通过设置在预定位置的体表接触区的压力传感器获得。本实施例中的压力传感器采用薄片式柔性压力传感器。

对当量圆半径r进行计算，计算当量圆半径r的公式为

对辅助形态参数进行计算，辅助形态参数的计算公式为，

对与时间有关的肌骨接触区域上的应力峰值p0(t)进行计算，应力峰值的计算公式为

对肌肉内部应力σzz(t)进行计算，软组织内部应力σzz(t)的计算公式为：

以下是计算软组织内部应力σzz(t)的一些特殊情况：

z＝0,r＝0,σzz(t)＝p0(t)，

z＝0,r≥a,σzz(t)＝0。

步骤s5，第四信息处理部24将软组织内部应力σzz(t)与软组织的损伤阈值进行比较，对软组织目前的生存状态进行评价以及预测软组织未来出现损伤的风险。

将软组织内部应力σzz(t)与软组织的损伤阈值进行比较的具体过程为：若软组织内部应力σzz(t)达到了软组织的损伤阈值，则判定软组织受到挫伤；若软组织内部应力σzz(t)接近但未达到软组织的损伤阈值，则判定软组织有受到挫伤的风险；若软组织内部应力σzz(t)远未达到软组织的损伤阈值，则判定软组织没有受到挫伤的风险。

实施例的作用与效果

根据本实施例所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过赫兹模型和医疗影像所得的数据对人体预定位置的末端半球骨进行模拟，并结合设置在预定位置的体表接触区的压力传感器，对人体软组织进行损伤情况、损伤风险的检测，相比传统技术，实现了本发明所涉及的一种非侵入式人体软组织损伤风险监测系统在实现检测功能的情况下，没有对人体造成创伤。

根据本实施例所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过赫兹模型对人体预定位置的末端半球骨进行模拟，并构建有限元模型对赫兹模型进行修正，保证了发明所涉及的一种非侵入式人体软组织损伤风险监测系统的准确性。

根据本实施例所涉及的一种非侵入式人体软组织损伤风险监测系统，本发明通过设置在预定位置的体表接触区的压力传感器实时对人体预定位置进行检测，从而实现对预定位置软组织的损伤情况以及损伤风险的实时判断，保证了发明所涉及的一种非侵入式人体软组织损伤风险监测系统的实时性。

上述实施方式为本发明的优选案例，并不用来限制本发明的保护范围。