脉搏波本构关系实验研究的探索

王礼立; 王晖; 丁圆圆; 陈霞波; 杨黎明; 龚文波; 浣石; 缪馥星

doi:10.11883/bzycj-2022-0434

脉搏波本构关系实验研究的探索

doi: 10.11883/bzycj-2022-0434

1.
宁波大学冲击与安全工程教育部/浙江省重点实验室，浙江宁波 315211
2.
宁波市中医院王晖工作室，浙江宁波 315000
3.
季华实验室，广东佛山 528200

基金项目: 国家自然科学基金（11872218, 11572161）；浙江省国医名师王晖传承工作室建设项目（GZS2020012）；宁波市公益类科技计划（202002N3144）

详细信息

作者简介:
王礼立（1934－　），男，教授，博士生导师，wanglili@nbu.edu.cn

中图分类号: O347.4
计量
- 文章访问数: 297
- HTML全文浏览量: 166
- PDF下载量: 60
- 被引次数: 2
出版历程
- 收稿日期: 2022-10-10
- 修回日期: 2022-11-14
- 网络出版日期: 2022-11-21
- 刊出日期: 2022-12-08

Exploration of experimental study on constitutive relations of pulse waves

1.
Key Laboratory of Impact and Safety Engineering, Ministry of Education, Ningbo University, Ningbo 315211, Zhejiang, China
2.
Wang Hui Workroom, Ningbo Hospital of Traditional Chinese Medicine, Ningbo 315000, Zhejiang, China
3.
Jihua Laboratory, Foshan 528200, Guangdong, China

摘要

摘要: 脉搏波本构关系决定着脉搏波的传播特征。如何通过实验研究来确定脉搏波本构关系，以及如何通过这些方法从现有文献数据来获得脉搏波本构关系，是当前研究的核心之一。本文中探索了3个可行途径：(1)由实测脉搏波波速对压力的关系C(p)进行反分析（无创法）；(2)直接对脉搏波p-V本构关系进行实测（有创法）；(3)由一系列实测脉搏波波形进行Lagrange反分析（无创法）。采用上述方法，根据现有文献数据，发现由C(p)关系的Rogers-Huang简化式可推得指数型p(V)本构关系；由MK-Hughes式可推得对数型p(V)本构关系。脉搏波传播特性随非线性本构参数发生显著变化。按中医体质分类观点，相应的脉搏波本构关系原则上也有不同类型，因人而定。在这个意义上，脉搏波的Lagrange反分析具有广阔发展前景，但它对正确选择测点和提高测量敏感性和精度等方面提出了更高要求。
- 脉搏波 /
- 本构关系 /
- 实验研究 /
- 反分析 /
- 中医
Abstract: The pulse wave constitutive relation determines propagation characteristics of pulse wave. How to determine it through experimental studies, and how to obtain it from the existing literature data through these methods, is one of the cores of the current research. Three feasible approaches were explored: (1) inverse analysis of the relationship C(p) (non-invasive method), (2) direct measurement of pulse wave p-V relationship (invasive method), and (3) Lagrange inverse analysis of a series of measured pulse wave (non-invasive method). By using the above methods and according to the existing literature data, it is found that the exponential p(V) constitutive relation can be deduced from the Rogers-Huang simplified formula of C(p) relation. The logarithmic p(V) constitutive relation can be deduced from the MK-Hughes equation. Pulse wave propagation characteristics vary significantly with nonlinear constitutive parameters. According to the viewpoint of body constitution classification in traditional Chinese medicine (TCM), the corresponding constitutive relation of pulse wave also has different types in principle, depending on people. In this sense, the Lagrange inverse analysis of pulse wave has a broad development prospect, but it puts forward higher requirements on the correct selection of measurement points and improving the sensitivity and accuracy of measurement.
- pulse wave /
- constitutive relation /
- experimental study /
- inverse analysis /
- traditional Chinese medicine (TCM)

HTML全文

由连续介质力学波动理论^[1]，按中医整体观思路，我们为脉搏波系统建立了一个等价的一维波力学模型^[2-3]。问题的控制方程组由如下3个守恒方程（动量守恒、质量守恒和能量守恒）加上系统本构方程（广义的状态方程）组成。

动量守恒方程为：

$\rho_{0} \frac{\partial v}{\partial t}=-\frac{\partial p}{\partial X}$

(1)

质量守恒方程为：

$\frac{{\partial v}}{{\partial X}} = {\rho _0}\frac{{\partial V}}{{\partial t}}$

(2)

能量守恒方程为：

${\rho _0}\frac{{\partial e}}{{\partial t}} = - p\frac{{\partial v}}{{\partial X}}\tag{3a}$

如果把式(2)代入式(3a)可以得到能量守恒方程的另一种形式：

${\text{d}}e = - p\;{\text{d}}V\text{，}\;\;或 e = - \int_{{V_0}}^V {p\;{\text{d}}V} \tag{3b}$

本构方程为：

$p = p\left( V \right)$

(4)

式(1)～(4)中：X为物质坐标，t为时间；p为压力，V为比容（与密度互为倒数），v为质点速度，e为比内能，均为(X，t)函数。本构方程 $p = p\left( V \right)$ 一般为非线性函数。

对上述脉搏波系统的控制方程组（式(1)～(4)）进行分析，并展开进一步研究^[2-6]后，可知：

(1)脉搏现象包含着血液的流动和在血液中传播的携带能量（内能e+动能）及整体信息的脉搏波。前者是人眼易见的实体血液本身的物质流；后者是人眼不易见、以波形式传播的能量信息流，可以分别理解为中医的“血”和“气”。中医的“气”是什么？目前并无共识。按上述分析并联系中医气血观可知：血是气传播的物质载体（媒介），气是以波形式传播、携带能量和生命体整体信息、推动血运行的波；即脉搏波不只是反映心脏、血液和血管等循环系统本身局部实体器官的信息。血液物质流与脉搏波能量信息流两者既互有区别，又密切相关耦合。这与中医“气为血之帅，血为气之母”的观点一致。

(2)脉搏波所包含的压力波p(X, t)、质点速度波v(X, t)、比容波V(X, t)和内能波e(X, t)等不同形式的扰动，以波速C传播：

$C = {V_0}\sqrt { - \left( {\frac{{\partial p}}{{\partial V}}} \right)} = \sqrt {\frac{{K\left( p \right)}}{{{\rho _0}}}}$

(5)

式中：K(p)= $- {\text{d}}p/\left( {{\text{d}}V/{V_0}} \right)$ 为非线性体积压缩模量。可见，波速C由非线性本构方程 $p = p\left( V \right)$ 的局部斜率（dp/dV）决定（见图1），随p而变，不是恒值。要注意区分脉搏波传播速度C和血液质点的流动速度v，前者（10⁰～10¹ m/s量级）以比后者（10⁻¹ m/s量级）快得多的速度，把心脏施加给循环血液的扰动（脉动载荷），携带着能量和信息，由近及远地传递到生命体的各部分，不到1 s即已传遍全身。因此，心脏实际上并非像目前教科书所述那样，简单地理解成把血液以流速v泵向全身的泵，而是脉搏波发生器^[6]。

图 1 波速C取决于本构方程p-V曲线的局部斜率（dp/dV）

Figure 1. The wave velocity C depends on the local slope (dp/dV) of the p-V curve

下载: 全尺寸图片幻灯片

(3)式(1)～(4)共同组成脉搏波系统的控制方程组，其中的3个守恒方程是普适的，反映了各种脉搏波的共性方面；而非普适的系统本构方程则反映了不同脉搏波的不同特性。可见，不同脉搏波系统的传播特性主要由脉搏波系统的本构关系决定，而各种病态可以看作偏离常态本构关系的不同反映。联想到近年来引人注目的中医体质学^[7-8]，就体质的重要性和地位而言，相当于前述量化的脉搏波系统中的本构关系。

需要强调，脉搏波本构关系的研究涉及对脉搏波性质的理解。目前很多研究者（例如费兆馥^[9]、Luo等^[10]）把脉搏波理解为脉搏压力波（纵波），也有研究者（例如Xue等^[11]）把脉搏波理解为沿固体血管传播的径向位移波（横波），而实际上脉搏波是流-固耦合和纵波-横波耦合的复杂波^[5]，其波速依赖于2个无量纲参数：血液-血管模量比K_b (p)/E_v(p)和血管径厚比D(p)/h₀，前者是本构关系的反映。对于这样的脉搏波系统的本构关系，迄今未见有关研究报导。

既然脉搏波的传播特性取决于本构关系，问题的研究核心之一就在于如何通过实验研究来确定脉搏波本构关系，以及如何采用这些方法从现有文献数据来获得脉搏波本构关系，这正是本文中研究的目的。

怎么来确定脉搏波系统本构方程呢？目前可从以下3个途径来进行探索：(1)由实测脉搏波波速对压力的关系C(p)进行反分析；(2)对脉搏波p-V本构关系进行实测；(3)由一系列实测脉搏波波形进行Lagrange反分析。下面分别加以探讨。

1. 由实测脉搏波波速与脉压的关系 $C(p)$ 进行反分析

近年来，有关脉搏波传播速度C与血压p的关联性C(p)之研究，在西医界倍受重视。

西医界把脉搏波波速PWV（pulse wave velocity，即本文式(5)所示的C）看作反映动脉弹性和硬化程度的指标。2007年欧洲高血压指南将主动脉PWV高于12 m/s列为高血压的危险因素之一^[12]。《中国高血压防治指南（2018年修订版）》也将PWV列为评估高血压患者心血管风险的重要指标^[13]。进一步，人们把实测的C(p)作为无创伤、无袖带、连续监测血压的重要途径之一，加以更广泛研究。通常从Moens-Korteweg (MK)公式^[14-15]和 Hughes公式^[16-17]出发，建立脉搏波波速C与血压p的关系。

MK公式^[14-15]为：

${C_{{\text{MK}}}} = {\left( {\frac{{E{h_0}}}{{{\rho _{\text{b}}}{D_0}}}} \right)^{1/2}} \tag{6a}$

Hughes公式^[16-17]为：

$E = {E_0}\exp (\zeta p) \tag{6b}$

两者结合，即有：

$C_{{\text{MKH}}}^2 = {A_{{\text{MKH}}}}\exp (\zeta p),\;\;\;\; {A_{{\text{MKH}}}} = \frac{{{E_0}{h_0}}}{{{\rho _{\text{b}}}{D_0}}} \tag{6c}$

式中：E、h₀和D₀(=2R₀)分别为血管弹性模量、初始壁厚和直径，ρ_b为血液密度，E₀为血压为零时的血管弹性模量，ζ为血管的材料系数。

Ma等^[18]指出，MK公式^[14-15]包含2个基本假设：(1)血管是弹性薄壁管（h₀/D₀<10）；(2)随血压的变化，动脉的厚度和半径保持不变。然而，对于人类动脉，这2个假设并不成立；而Hughes公式^[16-17]则完全是经验的，缺乏任何理论基础。Ma等^[18]基于Fung的超弹性模型^[19]，建立了一个由以下无量纲公式表征的、无需上述假设和经验表达式的分析模型（以下简称RH模型）：

$\dfrac{C}{{\sqrt {\dfrac{{{c_F}\exp ({a_2}E_{ZZ}^2)}}{{{\rho _{\text{b}}}}}} }} = f\left( {\dfrac{p}{{{c_F}\exp ({a_2}E_{ZZ}^2)}},{a_1},\dfrac{{{h_0}}}{{{R_0}}}} \right) \tag{7a}$

式中：c_F、a₁和a₂为材料常数，E_ZZ为血管轴向Green应变。在人体血压范围5～20 kPa (38～150 mmHg)内，式(7a)可用如下简化公式近似（以下简称RH简化式）：

$p = \alpha {C^2} + \beta ,\;\;\;\;{\text{or}}\;{C^2} = \frac{{p - \beta }}{\alpha } \tag{7b}$

式中：α和β为依赖于材料和血管几何参数的常数。

其实，按照我们的整体观脉搏波模型给出的波速C表达式（式(5))，可以直接由C(p)反演系统本构关系。事实上，式(5)可改写为：

$\frac{{{\text{d}}p}}{\left({C{{\left( p \right)}}}\right)^2} = - \frac{{{\text{d}}V}}{{V_0^2}} \tag{8}$

一旦C(p)已知，式(8)给出p-V关系的一阶常微分方程，不难由此解出脉搏波系统本构关系。

以RH简化式（式(7b)）为例，把此式代入式(8)后积分，设p=p₀时V=V₀，可解得：

$p = \beta {\text{ + }}\left( {{p_0} - \beta } \right)\exp \left( {\frac{{V_0 - V}}{{\alpha V_0^2}}} \right) \tag{9a}$

或把V表示为p的函数，则有：

$V = {V_0}{\text{ + }}\alpha V_0^2\ln \frac{{{p_0} - \beta }}{{p - \beta }} \tag{9b}$

可见，RH简化式实际上对应于指数形式的p(V)关系或对数形式的V(p)关系。具体的p-V曲线特性取决于本构参数α和β。能不能利用文献数据获得具体的脉搏波本构关系？下面讨论2个实例。

(1)实例1

Ma等^[18]指出，在人类血压范围（5～20 kPa）内，若取：

$\alpha = 0.18\;{\rm{kPa}}\cdot {\rm{s}}^2\cdot {\rm{m}}^{-2}\text{，}\qquad \beta= 2.7\;{\rm{kPa}} \tag{10a}$

RH简化式（式(7b)）能与RH模型（式(7a)）高度一致，如图2(a)所示。

图 2 RH简化式（式(7b)）与RH模型（式(7a)）和文献实测舒张压数据的比较

Figure 2. Comparison of the RH simplified formula (Eq.(7b)) with the RH model (Eq. (7a)) and the measured diastolic blood pressure (DBP) literature data

下载: 全尺寸图片幻灯片

把式(10a)代入式(9a)，取p₀=5 kPa，V₀=0.995×10⁻³ m³/kg，或即ρ₀=1.005×10³ kg/m³，可以得出对应的p-V关系式为：

$p = 2.7 + 2.3\exp \left( { - \left( {V - 0.995} \right)/0.178\;2} \right) \tag{10b}$

式中：p的单位为kPa，V的单位为10⁻³ m³/kg。相应的曲线如图3(a)所示。

图 3 不同情况下RH简化式（式(7b)对应的p-V曲线

Figure 3. p-V curves for the RH simplified formula (Eq. (7b)) under different conditions

下载: 全尺寸图片幻灯片

(2)实例2

Ma等^[18]还指出，若取：

$\alpha = 0.046\;{\rm{kPa}}\cdot {\rm{s}}^2\cdot {\rm{m}}^{-2}\text{，}\qquad \beta= 5.1\;{\rm{kPa}} \tag{11a}$

RH简化式(7b)能与文献报道的人体舒张压实测数据较好吻合，如图2(b)所示。

把式(11a)代入式(9a)，取p₀=6 kPa，V₀=0.995×10⁻³ m³/kg (ρ₀=1.005×10³ kg/m³)，可以得出对应的p-V关系式为：

$p = 5.1 + 0.9\exp \left( { - \left( {V - 0.995} \right)/0.045\;54} \right) \tag{11b}$

式中：p的单位为kPa，V的单位为10⁻³ m³/kg。相应的曲线如图3(b)所示。

对比图3中两曲线后可见，图3(b)曲线陡峭得多，意味着波速C随脉压p升高而迅速增大。由此可见，不同的本构参数α和β对于本构关系p-V曲线有明显影响，必然会影响脉搏波波速C以及不同形式脉搏波的波形传播特征，从而影响脉搏波能量和信息的传递。这为我们下一步定量化研究中医脉诊的脉象特征与系统本构关系p-V曲线特征之间的联系，包括研究常态及病态脉象与系统本构关系p-V曲线之间的定量联系等等，提供了一个良好开端和基础。

2. 对脉搏波 $p{\text{-}}V$ 关系进行实测

对脉搏波p-V本构关系最直接的实验研究方法无疑是对其进行实测。为说明这一点，先回顾一下式(6)所示的MK公式和Hughes公式。这一在西医界获得广泛认可和采用的公式，从我们的整体观脉搏波模型的观点来重新审视的话，MK公式只是本模型的式(5)在“血管为弹性薄壁管”特定假设下的简化特例，而Hughes公式实质是对狗进行体内（in vivo）+体外（in vitro）的脉搏p-V关系实测后所得出的。

事实上，由弹性力学^[20]知内压p下薄壁圆管的环应力σ_θ和环应变ε_θ分别为：

${\sigma _\theta } = \frac{{p{R_0}}}{{{h_0}}}\text{，}\;\;\; {\varepsilon _\theta } = \frac{{R - {R_0}}}{{{R_0}}}$

微分后，计及Hooke定律 ${\text{d}}{\sigma _\theta } = E{\text{d}}{\varepsilon _\theta }$ ，有：

${\text{d}}R = \frac{{R_0^2}}{{E{h_0}}}{\text{d}}p \tag{12a}$

另外，从几何关系可得薄壁血管内腔的体积压缩变形的变化（−dV_v/V₀）为：

$- \frac{{{\text{d}}{V_{\text{v}}}}}{{{V_0}}} = - {\text{d}}\left( {\frac{{{\text{π }}R_0^2 - {\text{π }}{R^2}}}{{{\text{π }}R_0^2}}} \right)\left/{V_0} = \frac{{2R}}{{{V_0}R_0^2}}{\text{d}}R\right. \tag{12b}$

由式(12a)和(12b)消去dR，可得：

$- {V_0}\frac{{{\text{d}}p}}{{{\text{d}}{V_{\text{v}}}}} = \frac{{E{h_0}}}{{2{\rho _0}R}} = {K_{\text{v}}} \tag{12c}$

式中：K_v(p)为血管涨缩所对应的等效体积压缩模量。对比式(6c)可知，K_MKH=K_v= $C_{{\rm{MK}}}^2$ 。将式(12c)代入式(5)的K，与MK公式（式(6a)）对比后可知，式(6a)只是式(5)的简化特例；并且MK公式只考虑了总容积变化dV中薄壁血管壁涨缩所作的贡献dV_v，而忽略了血液本构容积变化的贡献dV_b（相当于血液不可压缩假定）^[5]。

至于Hughes公式，回顾其原文^[16-17]可知，Hughes等在狗的活体体内实验中，用三轴超声导管测量了主动脉内外径随血压的变化，如图4(a)所示。由式(12)的分析可知，由此实际上已得出比容随压力的变化V(p)，即相应的活体p-V本构关系。而Hughes等在死狗体外实验中，则直接对p-V本构曲线进行了测量，如图4(b)所示。

4(a) Hughes等^[16-17]在狗的活体体内实验中，用三轴超声导管测得的主动脉内外径随血压的变化

4(a). Aortic inner and outer diameters along with aortic pressure recorded simultaneously and continuously with a triaxial ultrasonic catheter in the in-vivo measurements for dogs by Hughes et al^[16-17]

下载: 全尺寸图片幻灯片

4(b) Hughes等^[16-17]在死狗体外直接测得的p-V曲线

4(b). The p-V curve directly measured in the in-vitro measurements for dogs by Hughes et al^[16-17]

下载: 全尺寸图片幻灯片

Hughes等^[16-17]的实验研究结果当时以E(p) 的形式（式6(b)）表述，而没有以更具实质意义的p-V本构关系的形式表述，显然是受到当时已被广泛认可的MK公式的影响。其实，由Hughes公式（式(6b)）和式(12c)可导出：

${\text{d}}V = - \frac{D}{{{E_0}{h_0}}}\exp \left( { - \zeta p} \right){\text{d}}p$

(13)

设p=p₀时V=V₀，并设D=D₀（即设压力对管径的影响暂可忽略），式(13)积分后有：

$V = {\eta _0} + {\eta _1}\exp \left( { - \zeta p} \right),\;\;\;\;{\text{or}}\;p = - \frac{1}{\zeta }\ln \left( {\frac{{V - {\eta _0}}}{{{\eta _1}}}} \right)$

(14)

${\eta _0} = {V_0} - \frac{{{D_0}}}{{{E_0}{h_0}\zeta }}\exp \left( { - \zeta {p_0}} \right)，\;\;\;{\eta _1} = \frac{{{D_0}}}{{{E_0}{h_0}\zeta }}$

由此可见，MK-Hughes公式（式(6)）实际上对应于式(14)所示的指数函数形式的V(p)，或对数函数形式的p(V)本构关系，式中的η₀和η₁依赖于p₀、V₀、D₀/h₀、E₀和本构参数ζ。

Hughes等的贡献常常被归结于提出了Hughes公式（式(6b)）。其实，更有价值的贡献在于对脉搏波p-V本构关系进行体内+体外的直接实测。Hughes等是开创这方面研究的先驱，虽然尚限于狗的实验。

把上述脉搏波p-V关系实测法与C(p)关系反分析法两者相比较，前者属于直接实测法，其误差一般可小于后者的间接反演法；但前者属于有创法，而后者属于无创法，因而更便于推广应用。不过，不论是基于RH简化式反演的p-V本构关系（式(9)和图3），还是基于Hughes等实测结果建立的p-V本构关系（式(14)），都只是提供了一个基于大量实测数据的统计平均结果。从中医脉诊学的观点来看，脉象因人、因时、因地而异，换句话说，脉搏波的本构关系必定因人、因时、因地而异；甚至于同一个人因测量部位、测量时辰、病况变化而异。原则上，如同中医体质分类那样，本构关系也有不同类型，需要因人而定。

如何能够针对个人、无创地由实测数据来确定其脉搏波p-V关系？本文中建议采用脉搏波反分析法，即通过实测一系列波形进行Lagrange反分析^[2-3,21]，下面作进一步具体讨论。

3. 由系列实测波形进行Lagrange反分析

由连续介质力学波动理论^[1-3]，Lagrange反分析基于动量守恒方程（式(1)）和质量守恒方程（式(2)）。前者建立了压力p(X, t)的偏导数和质点速度v(X, t)的偏导数之间的关系，而后者建立了比容V(X, t)的偏导数和质点速度v(X, t)的偏导数之间的关系。通过求解v(X, t)作为过渡，可以求取脉搏波的p-V本构方程（式(4)）。具体步骤如下：

(1)在设定的测点设置由n(n≥3)个相邻的测压元件组成的阵列式压力传感器。一旦实测得到p (X_i, t)，i=1, 2 , ··· , n，就可由式(1)求得 $\partial v/\partial t$ ；应用初始条件 $\left.v\right|_{t=0}=v_{0}$ ，即可积分求得v (X, t)。

(2)求得v (X, t)后，就可由式(2)求得 $\partial V/\partial t$ ；应用初始条件 $\left.V\right|_{t=0}=V_{0}$ ，即可积分求得V (X, t) 。

(3)由 p(X, t)和V(X, t)消去t，即得p-V本构方程。

根据实测的一系列脉搏波信息去反求系统本构方程，力学上称为解“反问题”（反演）。其实，中医的切脉相当于根据脉诊信息对病况做出诊断，以便对症下药，予以诊治，以数理语言来表述，同样是在解“反问题”。两者在思维方式上是相通的。

为正确而成功地施行脉搏波的反分析，有必要强调以下2点：(1)正确选取测点；(2)满足必需的测量分辨率和精度。

首先，测点取何处？原则上凡能测量到脉搏波的地方都可以，当然以脉搏波信号越强越好。联系到至今仍指导着中医临床实践的“寸口三部九候诊法”，研究者们显然会更多关注桡动脉切脉部位“寸、关、尺”3处。初步容易设想，能不能直接对寸、关、尺3部的脉搏波进行Lagrange反分析^[2-3]。但实践发现并不可行，这一方面由于寸、关、尺3部的间距约为10 mm，达不到反分析差分数值分析的精度要求；另一方面，也是更重要的，由于按照中医寸口脉分候脏腑学说，左手寸、关、尺分候心、肝、肾，右手寸、关、尺分候肺、脾、命门（右肾），则显然在寸、关、尺3处分别确定的本构关系既有基本共性，又分别包含不同的分候特性，不能等同视之。事实上，脉诊实践中常常出现关部脉搏波最强，寸部次之，尺部最弱的情况。由于Lagrange反分析适用于波强单调变化的简单波，这样的强弱起伏的脉搏波分布显然不符合这一要求。所以就寸口脉而言，应该对寸、关、尺分别测量，在每个测点采用阵列式传感器测取系列波形后进行反分析。这样才能更好结合中医脉诊特点，以促进传统中医脉诊的定量化发展。

其次，在进行Lagrange反分析时都要经历数值微分和数值积分运算，从而存在一个如何满足对于测量分辨率和计算精度等要求的问题。下面通过一个实例来加以讨论。

Hu等^[22]采用图5(a)所示的电容型阵列传感器测量了脉搏压力波的时空分布。阵列传感器长10 mm，宽7.5 mm，含4个长度方向和3个宽度方向共12个微传感器。典型的“关”部的测量结果如图5(b)所示，图中纵坐标是无量纲压力波幅值，横坐标是采样数。把图5(b)的中间一列前3行（即图中Ch2、Ch5、Ch8）的实测无量纲压力波形换算为p-t图，则如图5(c)所示。图中波形沿脉搏波传播方向依次递减，可满足Lagrange反分析的要求。然而，实验中采用的扫描速率为100 Hz，相当于采样时间即时间分辨率为10 ms，其精度对于Lagrange反分析是远远不够的。由于微传感器中心间距为2.5 mm（参看图5(a)），如果波速按照10 m/s计，对应的时间差为0.25 ms。因此，传感器的时间测量精度至少应该是25 μs，才能保证两位有效数的测量精度。

5(a) Hu等^[22]的实验中采用的电容型阵列传感器

5(a). The capacitor array sensor used in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

5(b) Hu等^[22]的实验中“关”部各通道的实测波形

5(b). The waveform of each channel at Guan pulse taking position in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

5(c) 对应于Hu等^[22]的实验中电容型阵列传感器通道2、5、8的p-t图的电信号

5(c). Electrical signals corresponding to the p-t waveforms measured at Ch2, Ch5 and Ch8 of the capacitor array sensor in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

Hu等^[22]的实验虽然提供了有益的信息，但因其扫描速率不够高，尚难以据此进行Lagrange反分析，由此足以说明测量精度的重要性。今后需要研制具有微秒级别时间分辨率的毫米级阵列传感器，以推动Lagrange反分析法在脉搏波本构关系研究中的应用。

这是目前尚无由脉搏波Lagrange反分析来确定本构关系之实例的主要原因，也是我们正在进行中的研究工作，将另文讨论。

4. 结　论

脉搏波控制方程组的本构关系决定着脉搏波的传播特征，但迄今缺乏这方面的有关研究报导。如何通过实验研究来确定脉搏波本构关系，以及如何通过这些方法从现有文献数据来获得脉搏波本构关系，是当前研究的核心之一。

本文中探索了3个可行途径：(1)由实测脉搏波波速C对压力p的关系C(p)进行反分析（无创法)；(2)直接对脉搏波p-V关系进行实测（有创法）；(3)由一系列实测脉搏波波形进行Lagrange反分析（无创法）。

采用上述方法，根据现有文献数据，成功地获得了不同类型的脉搏波本构关系。由C(p)关系的RH简化式实际上可推得指数型p(V)本构关系。由MK-Hughes式实际上可推得对数型p(V)本构关系。脉搏波传播特性随非线性本构参数发生显著变化。这一研究还为今后从大量已有文献数据探索其背后更具实质意义的本构关系提供了一个新思路、新途径。

按照中医体质分类观点，脉搏波本构关系必定也有不同类型，具体因人、因时、因地而异；并且对个人还因测量部位、测量时辰、病况变化而异，原则上不存在单一普适的关系。在这个意义上，脉搏波的Lagrange反分析具有广阔发展前景，但它对正确选择测点和提高测量敏感度和精度等方面提出了更高要求。

图 1 波速C取决于本构方程p-V曲线的局部斜率（dp/dV）

Figure 1. The wave velocity C depends on the local slope (dp/dV) of the p-V curve

下载: 全尺寸图片幻灯片

图 2 RH简化式（式(7b)）与RH模型（式(7a)）和文献实测舒张压数据的比较

Figure 2. Comparison of the RH simplified formula (Eq.(7b)) with the RH model (Eq. (7a)) and the measured diastolic blood pressure (DBP) literature data

下载: 全尺寸图片幻灯片

图 3 不同情况下RH简化式（式(7b)对应的p-V曲线

Figure 3. p-V curves for the RH simplified formula (Eq. (7b)) under different conditions

下载: 全尺寸图片幻灯片

4(a) Hughes等^[16-17]在狗的活体体内实验中，用三轴超声导管测得的主动脉内外径随血压的变化

4(a). Aortic inner and outer diameters along with aortic pressure recorded simultaneously and continuously with a triaxial ultrasonic catheter in the in-vivo measurements for dogs by Hughes et al^[16-17]

下载: 全尺寸图片幻灯片

4(b) Hughes等^[16-17]在死狗体外直接测得的p-V曲线

4(b). The p-V curve directly measured in the in-vitro measurements for dogs by Hughes et al^[16-17]

下载: 全尺寸图片幻灯片

5(a) Hu等^[22]的实验中采用的电容型阵列传感器

5(a). The capacitor array sensor used in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

5(b) Hu等^[22]的实验中“关”部各通道的实测波形

5(b). The waveform of each channel at Guan pulse taking position in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

5(c) 对应于Hu等^[22]的实验中电容型阵列传感器通道2、5、8的p-t图的电信号

5(c). Electrical signals corresponding to the p-t waveforms measured at Ch2, Ch5 and Ch8 of the capacitor array sensor in the experiment by Hu et al^[22]

下载: 全尺寸图片幻灯片

参考文献(22)

[1]	王礼立. 应力波基础 [M]. 2版. 北京: 国防工业出版社, 2005. WANG L L. Foundation of stress waves [M]. 2nd ed. Beijing, China: National Defense Industry Press, 2005.
[2]	WANG L L, WANG H. Mechanics modeling and inverse analyses of pulse wave system from the view-point of traditional Chinese medicine [C]//Proceedings of the 35th International Conference on Ocean, Offshore and Arctic Engineering (OMAE2016). Busan, South Korea: BEXCO Exhibition Convention Center, 2016.
[3]	王礼立, 王晖. 脉搏波系统的力学模型及反演兼对若干中医学问题的讨论 [J]. 力学学报, 2016, 48(6): 1416–1424. DOI: 10.6052/0459-1879-15-322. WANG L L, WANG H. Mechanics modeling and inverse analyses of pulse waves system with discussions on some concepts in the traditional Chinese medicine [J]. Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(6): 1416–1424. DOI: 10.6052/0459-1879-15-322.
[4]	王礼立, 王晖, 杨黎明, 等. 论脉搏波客观化和定量化研究的症结所在 [J]. 中华中医药杂志, 2017, 32(11): 4855–4863. WANG L L, WANG H, YANG L M, et al. Crux of objectification and quantification of pulse waves [J]. China Journal of Traditional Chinese Medicine and Pharmacy, 2017, 32(11): 4855–4863.
[5]	缪馥星, 王晖, 王礼立, 等. 血液-血管耦合特性与脉搏波传播特性的关系 [J]. 爆炸与冲击, 2020, 40(4): 041101. DOI: 10.11883/bzycj-2020-0082. MIAO F X, WANG H, WANG L L, et al. Relationship between the blood-vessel coupling characteristics and the propagation of pulse waves [J]. Explosion and Shock Waves, 2020, 40(4): 041101. DOI: 10.11883/bzycj-2020-0082.
[6]	王晖, 王礼立, 缪馥星, 等. 论心脏功能的“泵说”与“波说” [J]. 爆炸与冲击, 2020, 40(11): 111101. DOI: 10.11883/bzycj-2020-0386. WANG H, WANG L L, MIAO F X, et al. On “pump theory” and “wave theory” of cardiac function [J]. Explosion and Shock Waves, 2020, 40(11): 111101. DOI: 10.11883/bzycj-2020-0386.
[7]	王琦. 中医体质学 [M]. 北京: 人民卫生出版社, 2009. WANG Q. Constitutionology of Chinese Medicine [M]. Beijing, China: People’s Medical Publishing House, 2009.
[8]	王晖. 体质的中医保健 [M]. 浙江宁波: 宁波出版社, 2009. WANG H. Traditional Chinese Medicine health care of body constitutions[M]. Ningbo, Zhejiang, China: Ningbo Press, 2009.
[9]	费兆馥. 现代中医脉诊学 [M]. 2版. 北京: 人民卫生出版社, 2006. FEI Z F. Contemporary sphygmology in traditional Chinese medicine [M]. 2nd ed. Beijing, China: People’s Medical Publishing House, 2006.
[10]	LUO C H, CHUNG Y F, HU C S, et al. Possibility of quantifying TCM finger-reading sensation: I. bi-sensing pulse diagnosis instrument [J]. European Journal of Integrative Medicine, 2012, 4(3): 255–262. DOI: 10.1016/j.eujim.2012.03.003.
[11]	XUE Y, SU Y, ZHANG C, et al. Full-field wrist pulse signal acquisition and analysis by 3D Digital Image Correlation [J]. Optics and Lasers in Engineering, 2017, 98: 76–82. DOI: 10.1016/j.optlaseng.2017.05.018.
[12]	2007 Guidelines for the management of arterial hypertension [J]. Hypertension, 2007, 25: 1105–1187.
[13]	中国高血压防治指南(2018年修订版) [J]. 中国心血管杂志, 2019, 24(1): 24–56. Guidelines for the prevention and treatment of hypertension in China (2018 revision) [J]. Chinese Journal of Cardiovascular Sciences, 2019, 24(1): 24–56.
[14]	MOENS A I. Die pulskurve[M]. Leiden, the Netherlands: Brill, 1878.
[15]	KORTEWEG D J. Ueber die fortpflanzungsgeschwindigkeit des schalles in elastischen röhren [J]. Annalen der Physik Chemie, 1878, 241(12): 525–542. DOI: 10.1002/andp.18782411206.
[16]	HUGHES D J, BABBS C F, GEDDES L A, et al. Measurements of Young’s modulus of elasticity of the canine aorta with ultrasound [J]. Ultrason Imaging, 1979, 1(4): 356–367. DOI: 10.1177/016173467900100406.
[17]	HUGHES D J, FEARNOT N E, BABBS C F, et al. Continuous measurement of aortic radius change in vivo with an intra-aortic ultrasonic catheter [J]. Medical and Biological Engineering and Computing, 1985, 23(3): 197–202. DOI: 10.1007/BF02446857.
[18]	MA Y, CHOI J, HOURLIER-FARGETTE A, et al. Relation between blood pressure and pulse wave velocity for human arteries [J]. Proceedings of the National Academy of Sciences, 2018, 115(44): 11144–11149. DOI: 10.1073/pnas.1814392115.
[19]	FUNG Y C. Biomechanics: circulation [M]. New York, USA: Springer, 1997.
[20]	TIMOSHENKO S P, GOODIER J N. Theory of elasticity [M]. 3rd ed. London, UK: McGraw-Hill Book Company, 1970.
[21]	王礼立, 胡时胜, 杨黎明, 等. 材料动力学 [M]. 合肥: 中国科学技术大学出版社, 2017. WANG L L, HU S S, YANG L M, et al. Dynamics of materials [M]. Hefei, Anhui, China: Press of University of Science and Technology of China, 2017.
[22]	HU C S, CHUNG Y F, YEH C C, et al. Temporal and spatial properties of arterial pulsation measurement using pressure sensor array [J]. Evidence-Based Complementary and Alternative Medicine: eCAM, 2012, 2012(1): 745127–745135. DOI: 10.1155/2012/745127.

施引文献

期刊类型引用(2)

1.	陈明东，黄晓东，李潮. 多学科交叉融合背景下的人体血压和血管弹性实验. 物理实验. 2024(01): 22-27+33 . 百度学术
2.	陈琥颖，龚文波，贾梦娜，丁溢群，陈一夫，廖晟极. 2型糖尿病气虚证和阴虚证常见脉象的脉搏波特征研究. 中医临床研究. 2024(32): 7-11 . 百度学术

其他类型引用(0)

资源附件(0)

访问统计

图(8)

计量

文章访问数: 297
HTML全文浏览量: 166
PDF下载量: 60
被引次数: 2

1. 由实测脉搏波波速与脉压的关系 $C(p)$ 进行反分析
2. 对脉搏波 $p{\text{-}}V$ 关系进行实测
3. 由系列实测波形进行Lagrange反分析
4. 结　论

脉搏波本构关系实验研究的探索

doi: 10.11883/bzycj-2022-0434

作者简介: 王礼立（1934－ ），男，教授，博士生导师，wanglili@nbu.edu.cn

计量

出版历程

Exploration of experimental study on constitutive relations of pulse waves

1. 由实测脉搏波波速与脉压的关系C(p)C(p) 进行反分析

2. 对脉搏波p-Vp{\text{-}}V 关系进行实测

3. 由系列实测波形进行Lagrange反分析

4. 结 论

期刊类型引用(2)

其他类型引用(0)

计量

出版历程

目录

1. 由实测脉搏波波速与脉压的关系C(p)C(p) 进行反分析

2. 对脉搏波p-Vp{\text{-}}V 关系进行实测

3. 由系列实测波形进行Lagrange反分析