Processing math: 70%

EI / SCOPUS / CSCD 收录

中文核心期刊

侧倾−垂向双轴激励下的坐姿人体生物动力学特性

林泽峰, 林杰威, 郑旭, 刘迟, 邱毅

林泽峰, 林杰威, 郑旭, 刘迟, 邱毅. 侧倾−垂向双轴激励下的坐姿人体生物动力学特性[J]. 声学学报, 2024, 49(2): 183-193. DOI: 10.12395/0371-0025.2023197
引用本文: 林泽峰, 林杰威, 郑旭, 刘迟, 邱毅. 侧倾−垂向双轴激励下的坐姿人体生物动力学特性[J]. 声学学报, 2024, 49(2): 183-193. DOI: 10.12395/0371-0025.2023197
LIN Zefeng, LIN Jiewei, ZHENG Xu, LIU Chi, QIU Yi. Biodynamic characteristics of seated human body under roll-vertical vibration[J]. ACTA ACUSTICA, 2024, 49(2): 183-193. DOI: 10.12395/0371-0025.2023197
Citation: LIN Zefeng, LIN Jiewei, ZHENG Xu, LIU Chi, QIU Yi. Biodynamic characteristics of seated human body under roll-vertical vibration[J]. ACTA ACUSTICA, 2024, 49(2): 183-193. DOI: 10.12395/0371-0025.2023197
林泽峰, 林杰威, 郑旭, 刘迟, 邱毅. 侧倾−垂向双轴激励下的坐姿人体生物动力学特性[J]. 声学学报, 2024, 49(2): 183-193. CSTR: 32049.14.11-2065.2023197
引用本文: 林泽峰, 林杰威, 郑旭, 刘迟, 邱毅. 侧倾−垂向双轴激励下的坐姿人体生物动力学特性[J]. 声学学报, 2024, 49(2): 183-193. CSTR: 32049.14.11-2065.2023197
LIN Zefeng, LIN Jiewei, ZHENG Xu, LIU Chi, QIU Yi. Biodynamic characteristics of seated human body under roll-vertical vibration[J]. ACTA ACUSTICA, 2024, 49(2): 183-193. CSTR: 32049.14.11-2065.2023197
Citation: LIN Zefeng, LIN Jiewei, ZHENG Xu, LIU Chi, QIU Yi. Biodynamic characteristics of seated human body under roll-vertical vibration[J]. ACTA ACUSTICA, 2024, 49(2): 183-193. CSTR: 32049.14.11-2065.2023197

侧倾−垂向双轴激励下的坐姿人体生物动力学特性

详细信息
    通讯作者:

    邱毅, yiqiu@zju.edu.cn

  • 中图分类号: 43.40

  • PACS: 
    • 43.40  (结构声学与振动)

Biodynamic characteristics of seated human body under roll-vertical vibration

  • 摘要:

    研究了侧倾−垂向双轴激励下的坐姿人体动力学特性。测取了8名被试者在单轴随机侧倾振动下的动力学响应, 通过人椅交界面处的角加速度和力矩定义了视在惯量, 识别了人体在侧倾方向上位于1 Hz和2.5 Hz方向的两阶共振。随后引入垂向振动, 对侧倾−垂向双轴激励下的人体响应开展频响函数分析与偏相干分析, 分析了人体垂向响应与侧倾响应的耦合关系。研究发现, 侧倾振动在0~10 Hz测试频率范围内对人体动力学响应特性几乎没有影响, 垂向振动在两阶侧倾共振频率附近对包括模态频率、视在惯量幅值等侧倾动力学特性几乎没有影响。最后, 面向侧倾与垂向振动激励, 建立了人体动力学模型并识别了其中的动力学参数, 实现了单一模型对双轴激励下人体响应的预测。

    Abstract:

    This paper presents a study on the biodynamic characteristics of seated human body under roll-vertical vibration. Initially, the dynamic responses of eight subjects under single-axis random roll vibration was measured. The apparent inertia was introduced, defined as the frequency response function between angular acceleration and torque measured at the interface between human body and seat. Two-order resonances of human body in the roll direction at 1 Hz and 2.5 Hz were identified. Subsequently, vertical vibrations were combined with the roll ones to investigate the human response under dual-axis roll-vertical vibration. Frequency response function analysis and partial coherence analysis were conducted, and the coupling between the vertical and roll responses of human body was discussed. The study revealed that roll vibration had minimal influences on the dynamic response of seated human body in the vertical direction within 0−10 Hz. Similarly, vertical vibration had negligible impact on human response in the roll direction around the first two resonance frequencies, including modal frequencies and associated modulus of roll apparent inertia. Finally, a dynamic model was established to represent the biodynamic characteristics of human body under combined roll and vertical vibration, and dynamic parameters of the model were identified.

  • 振动广泛存在于人们的日常生活中, 以交通活动为代表, 振动暴露是影响交通运输装备乘坐舒适性的重要因素。由于舒适性感受具有主观性, 获取具有代表性的主观感受指标, 有赖于一定数量的人群参与主观评价实验。相关测试的时间和经济成本十分可观, 甚至根本难以在产品设计与开发的过程中开展。因此, 如果能够建立主观感受与振动加速度等客观参数之间的联系, 以此开展振动舒适性设计, 就可以在充分考虑人体实际感受的同时, 兼顾工程应用中的可行性。而这一目标的实现是以了解、掌握人体在振动暴露下的动力学响应特性为前提的。

    振动激励下人体客观动力学响应特性的诸多研究成果有力推动了振动舒适性设计的进步, 以单轴平动, 尤其是单轴垂向振动激励下的研究成果最为丰富。垂向振动激励下测取的视在质量表明, 人体在4~8 Hz附近存在一阶共振, 部分研究也在8~10 Hz与10~16 Hz附近观察到较弱的二阶、三阶共振[1-3]。在水平方向的类似研究通过前后、侧向视在质量识别了人体在1 Hz附近与2.5 Hz附近的两阶共振[4]。部分研究还讨论了关于激励强度[5-6]、被试者条件(如坐姿[7-8]、性别[9]等)对人体动力学特性的影响。基于上述研究成果, 通过集中质量法[10-12]、多体动力学[13]和有限元方法[14]形成了多种复杂程度与适应性不同的人体动力学模型, 实现了振动激励下人体响应的预测。然而, 在实际振动环境中, 振动绝不会是单一方向的。单轴激励下, 人体的交叉轴视在质量有时也会产生明显的幅值[15-16]。这意味着某一方向的振动激励也会引起人体在其他方向上的响应, 多个研究都观察到明显的垂向−前后方向交叉轴视在质量, 且在5 Hz全身垂向模态频率附近出现峰值[17-18]。当多轴激励共同作用时, 人体的同轴响应与交叉轴响应可能发生耦合, 基于此, 有学者开展了多轴振动激励下的人体生物动力学特性研究以及人体生物动力学模型研究。垂向−前后双轴激励下, 同轴和交叉轴视在质量与单轴激励下表现出基本一致的频率分布特征, 部分实验中会观察到因新激励轴引入而出现的峰值频率前移, 如Mansfield[19]观察到相比单轴激励, 双轴平动激励下的同轴、交叉轴视在质量峰值出现在更低的频率范围; 类似地, Qiu[18,20]观察到了随前后方向激励强度增加, 垂向同轴视在质量峰值频率前移的现象。进一步引入侧向激励后就可以对三轴平动激励下的人体动力学特性进行分析。从交叉轴视在质量上看, 垂向与侧向激励下的人体响应耦合并不强, 因此三轴平动耦合下的人体动力学特性与垂向−前后双轴激励下的表现具有一致性[21-22]。Mansfield等[19,23]比较了平动三个方向: 前后、侧向、垂向单轴激励下和两两组合形成的双轴激励及三轴激励下的视在质量, 也得到了类似的结论, 且认为人体共振频率随激励强度的增加而降低的非线性特征是对于总输入振动量而言的, 而非单一轴。基于人体动力学试验测试结果, 有研究进一步建立了面向多轴激励的人体动力学模型[24-25], 模型通过使用不同方向上的弹簧、阻尼单元反映了人体在对应方向上的动力学特征, 并使用多方向上测取的视在质量进行模型参数识别。

    上述研究中, 单轴人体动力学特性的研究以平动方向为主, 尤其是垂向激励下的研究较为充分。相比之下, 对旋转激励下的人体客观动力学特性仍缺乏了解, 而转动激励引起的人体主观不适感已被证明是十分明显的, 尤其是侧倾激励。Pason等[26]对于振动主观舒适度的研究表明, 同幅值的侧倾振动将比俯仰振动导致坐姿人体更强烈的不舒适感, Beard等[27-28]则确认了侧倾、侧向激励引起的主观不舒适感受。对单轴转动激励下的人体动力学响应研究不够充分, 也限制了对转动与平动激励叠加作用的分析和考虑, 使得转动−平动多轴激励对人体动力学特征的影响效果尚不明确。

    本文考虑单轴侧倾振动、侧倾−垂向双轴振动激励, 通过对垂向视在质量与侧倾视在惯量的测量和分析, 识别转动激励下的人体生物动力学特性, 并探究人体垂向响应与侧倾响应间是否存在耦合作用, 从转动激励与多轴激励的角度对现有人体动力学理论进行扩充。在此基础上, 本文将建立面向侧倾、垂向双轴激励的人体动力学模型, 在转动(侧倾)方向上扩展现有人体动力学模型, 同时提供多轴激励下人体动力学建模方法。本文所开展的研究有助于从机理上解释人体主观感受的成因, 为振动舒适性评价、预测和设计提供指导。

    人体动力学测试实验系统如图1所示, 包括六自由度运动平台, 坐垫式加速度传感器、测力台、数采系统及计算机、刚性座椅与其他附件等。

    图  1  人体动力学测试系统

    运动平台可以提供六轴振动激励, 其动平台直径1 m; 各轴运动范围为垂向±16 cm、前后和侧向±30 cm、侧倾和俯仰±9°、横摆±12°; 由伺服电机驱动, 能够在垂向峰−峰值8 cm的运动条件下提供±10 m/s2的振动加速度。

    人椅交界面所受到的转矩即为侧倾激励下的人体响应, 通过安装在刚性座椅上表面的测力台进行测量。测力台四角分别安装有四枚精确校准的力传感器, 所测得的电荷信号经电荷放大器输入采集系统。测力台基于L3010型力传感器(熙源电子科技有限公司, 扬州, 中国)定制, 并配有专用电荷放大器。角加速度通过两个布置于人椅交界面的两个中心相距18 cm的356B41型坐垫式加速度传感器(PCB Piezotronics, Inc., Depew, USA)测得。在振动台表面还布置有两个356A26型加速度传感器, 与座椅表面加速度传感器配合使用, 用于检查座椅刚性。

    测试数据通过Smicenter Scadas Mobil (Siemens, Munich, Germany)采集。全部实验中采样频率均为512 Hz, 抗混叠滤波器设置为40 Hz。

    试验使用一高0.52 m, 宽0.47 m, 进深0.41 m的刚性座椅, 从而将运动平台产生的振动原本地传递至人体。刚性座椅在10 Hz以内不存在模态。座椅配有安全带、急停按钮以及靠背, 靠背仅作安全保护之用, 防止被试者后仰跌落, 测试期间不允许被试者倚靠。此外, 还为被试者提供高度可选的脚托, 以保证大腿与座椅表面基本处于平行状态, 即大腿与座椅表面的接触程度为平均[17]

    八名年龄22~33岁之间的男性被试者参与了实验, 平均体重72.1 kg (55~90 kg), 平均身高176.1 cm (168~188 cm)。被试者在实验过程中被要求全程系紧安全带, 并手握急停按钮, 同时保持正坐姿态、不倚靠靠背, 身体自然坐直, 即上身既不紧绷挺直, 也不过于放松至曲背。同时要求被试者在全部测试中尽量保持身体坐姿不变, 并避免主动的肌肉动作来抑制身体摆动。参与本测试的全部被试者均签署了知情同意书,该试验由浙江大学科技伦理委员会审批批准。

    测试使用6组定义于0.5~10 Hz的随机激励信号, 方向分别为单轴垂向、单轴侧倾、双轴垂向−侧倾, 强度如表1所示。其中, 激励信号1, 2, 3用于分析单轴侧倾振动下的人体动力学响应; 激励信号4, 5, 6用于分析侧倾−垂向振动下的人体动力学响应, 激励信号6是激励信号4和5的组合。考虑到平直路面行驶的车辆, 其驾驶室地板发生的侧倾振动一般并不十分强烈, 较平动方向振动引起的不舒适感受也较弱, 因此测试面向轨道车辆、工程机械及重型卡车等应用场景, 选取强度对应于工程机械或重型卡车在颠簸路面如比利时路、破损水泥路等工况下车厢内侧倾振动强度。在信号频率范围内, 其具有大致平坦的功率谱分布, 适用于响应特性的测量。每组信号长12秒, 包含最前一秒与最后一秒使用余弦信号进行了调制, 以保护运动平台启停时不产生冲击。

    表  1  侧倾−侧向激励下人体响应特性测试用激励信号
    激励信号垂向
    (m/s2 r.m.s.)
    侧倾
    (rad/s2 r.m.s.)
    101.0
    201.5
    302.0
    402.0
    52.00
    62.02.0
    下载: 导出CSV 
    | 显示表格

    人体全身客观响应由人椅交界面处测得的频响函数表示。对于垂向振动, 即视在质量, 定义为

    M(f)=GaVFV(f)GaV(f), (1)

    其中, GaV(f)是人椅交界面处垂向振动加速度的自谱; GaVFV(f)是人椅交界面处垂向力与垂向加速度的互谱。类似地, 本文定义用于侧倾方向人体动力学响应特性分析的视在惯量, 即人椅交界面处的转矩和角加速度之间的频响函数, 表示为

    I(f)=GaRMR(f)GaR(f), (2)

    其中, GaR(f)是人椅交界面处侧倾方向角加速度的自谱, GaRMR(f)是人椅交界面处侧倾方向转矩和角加速的互谱。与视在质量类似, 视在惯量具有动态转动惯量的物理意义, 可以用来描述转动激励下人体的动力学特性。考虑到本文使用的多轴激励使用信号时长12秒, 所计算的频响函数包括垂向视在质量和视在惯量统一使用 0.25 Hz频域分辨率。

    对于多轴振动激励下的人体动力学响应测试, 使用多输入单输出方法(MISO)计算视在质量(或视在惯量)并开展偏相干分析。以垂向响应作为输出和侧倾响应作为输出的偏相干分析模型如图2所示。

    图  2  侧倾−垂向振动激励下的人体动力学响应特性分析MISO模型 (下标v, r分别表示垂向、侧倾振动输入, 下标V, R分别表示垂向、侧倾响应输出) (a) 垂向响应输出; (b) 侧倾响应输出

    应用MISO方法的视在质量(视在惯量)计算方法为

    MFAA(f)=GaAFA(1GaAaBGaBFAGaBaBGaAFA)GaAaA(1γ2), γ2 = |GaAaB|2GaAaAGaBaB, (3)

    其中, GaAaBA方向和B方向加速度(角加速度)的互谱, GaAaAGaBaB 分别是A方向和B方向加速度(角加速度)的自谱, GaAFAGaBFA 分别是A方向、B方向加速度(角加速度)与A方向力(力矩)的互谱, γA方向和B方向振动加速度的相干函数。A方向振动激励与A方向响应的常相干函数表示为

    γ2AA(f)=|GaAFA(f)|2GaAaA(f)GFAFA(f), (4)

    B方向振动激励与A方向响应的偏相干函数表示为

    γ2BA.A(f)=|GaBFA.aA(f)|2GaBaB.aA(f)GFBFB(f), (5)

    其中, GaAFAA方向振动加速度(角加速度)与A方向力(转矩)的互谱, GaBFA.aAA方向条件输入振动激励与B方向振动激励的互谱, GaBaB.FAB方向条件输入振动激励与A方向响应的互谱。

    侧倾激励下八名被试者于刚性座椅上的动力学响应特性测试结果即视在惯量如图3所示。对于大多数被试者, 都可以在1 Hz附近观察到视在惯量的第一个峰值, 也可视作坐姿人体侧倾激励下的一阶共振。在3 Hz附近, 可以观察到人体第二阶峰值, 二阶峰值的幅值要远低于一阶峰值。推测一阶模态可能来自于侧倾激励下由脊椎控制的全身摆动, 二阶模态则来自于上身肌肉−骨骼组织与腰椎及坐骨的异相运动。

    图  3  单轴侧倾激励下八名被试者的视在惯量 (a) 1.0 rad/s2 r.m.s.; (b) 1.5 rad/s2 r.m.s.; (c) 2.0 rad/s2 r.m.s.

    图3中可以观察到明显的被试者间差异性(Inter-Subject Variability)。对于不同的被试者, 一阶共振频率存在于0.75 Hz 至1.25 Hz范围内。二阶共振频率的被试者间差异性更为明显, 尤其是对于视在惯量的幅值。二阶共振频率处视在惯量的幅值可能与被试者体重具有相关性。被试者1、被试者2、被试者8具有全部被试者中最重的坐重, 他们三人也在二阶共振频率处产生了较为显著的响应, 如图4所示。这一体重与视在惯量幅值的关系暗示了侧倾激励下坐姿人体的二阶模态很可能不是由脊椎摆动等形式的全身振动所控制, 而是与身体局部的肌肉和脂肪组织有关。此外, 受脊椎控制的一阶模态在其模态频率处, 视在惯量幅值也未呈现相似的规律。相比一阶模态, 视在惯量的二阶峰值较低、共振区更宽, 这对应于骨盆较上身脊椎−肌肉结构更大的阻尼[29]。利用平均差异系数(Average Coefficient of Variation, ACV)来描述被试者间差异性。平均差异系数定义为各频率下八名被试者视在惯量幅值标准差与均值之比的平均值[30], 表示为

    图  4  强度为1.5 rad/s2 r.m.s.单轴侧倾激励下八名被试者的视在惯量
    ACV=1nnj=11NNk=1(Mkj¯Mj)2¯Mj (6)

    其中, n为0.5~10 Hz内的谱线数, 本文采用0.25 Hz频率分辨率, 故n = 39; N为被试者数目, 本文中N = 8; Mkj为第k名被试者在第j个频率处的视在惯量幅值; ¯Mj为八名被试者在第j个频率处的视在质量幅值均值。三种振动强度下的平均差异系数分别为0.27, 0.23, 0.22。这表明被试者间差异性在低振动强度下更为明显。

    对八名被试者的视在惯量幅值取中位数, 结果如图5所示。激励强度对一阶、二阶响应峰值频率和峰值频率处视在惯量幅值影响的统计分析如表2所示。在侧倾激励下, 激励强度对模态频率的影响并不明显, 三种测试强度下坐姿人体视在惯量的一阶峰值均出现在1 Hz左右。这与垂向激励下人体表现出的非线性特征(即人体共振频率随振动强度增加而前移[1])有所不同。这一现象意味着在实验所涉及的强度范围内, 人体侧倾方向的非线性特征并不强烈, 同时也为使用线性弹簧、阻尼所构成的动力学系统对人体侧倾激励下的全身响应特性进行建模提供了理论基础。同时, 激励强度对一阶峰值频率处的视在惯量幅值的影响也较为轻微。相比一阶模态, 二阶模态的非线性特征相对明显: 随振动强度的增加, 视在惯量二阶峰值频率前移, 峰值频率处幅值增加。相对明显的非线性特征也是肌肉和脂肪组织参与支配二阶模态的佐证之一。

    图  5  单轴侧倾激励下八名被试者的视在惯量中位数
    表  2  激励强度对共振频率和共振频率处视在惯量幅值影响的Wilcoxon 符号秩检验结果
    激励强度 1 rad/s2 r.m.s. 1.5 rad/s2 r.m.s.
    1.5 rad/s2 r.m.s. ns/ns/ns/*
    2.0 rad/s2 r.m.s. ns/ns/ns/ns ns/ns/ns/*
    a) ns表示影响不显著, 即p>0.05; * 表示p<0.05; **表示 p<0.01; ***表示 p<0.005, Wilcoxon;
    b) 由“/”分隔的显著性水平结果, 依次为激励强度对一阶共振频率/一阶共振频率处幅值/二阶共振频率/二阶共振频率处幅值的影响。
    下载: 导出CSV 
    | 显示表格

    八名被试者单、双轴激励下的垂向同轴视在质量如图6所示, 视在质量中位数如图7所示。八名被试者中, 尽管存在被试者间差异性, 但是除一名被试者一阶共振频率较低、共振区明显较窄外, 其余被试者表现出了一致的振动响应规律。单轴垂向视在质量在5 Hz附近观察到了一阶共振峰, 部分被试者在7 Hz附近观察到二阶共振峰。这一结果与文献报告的垂向人体动力学特性一致[1,14]。侧倾振动引入后, 人体垂向动力学特性未发生明显改变, 一阶共振出现的频率、视在质量的幅值均无明显变化。对八名被试者的视在惯量一阶峰值频率、峰值频率处幅值开展统计分析, Wilcoxon符号秩检验的结果均表现出影响不显著(p>0.05)。

    图  6  八名被试者的垂向同轴视在质量 (a) 单轴垂向振动激励; (b) 侧倾−垂向双轴激励
    图  7  八名被试者单、多轴垂向同轴视在质量中位数

    八名被试者单、双轴激励下的侧倾同轴视在惯量如图8所示, 视在惯量中位数结果如图9所示。与前文测试结果一致, 在1 Hz附近观察到明显的一阶共振峰值, 部分被试者在3~4 Hz附近可以观察到二阶峰值。相比垂向视在质量, 侧倾视在惯量的被试者间差异性较为明显。被试者间差异在单轴侧倾激励下的人体动力学响应测试中也有所表现(见2.1节)。垂向振动引入后, 即垂向−侧倾振动激励相比单轴侧倾振动激励下, 视在惯量峰值频率、峰值频率处的幅值均未出现明显变化。Wilcoxon符号秩检验证实了垂向激励的引入对于一阶共振频率及共振频率处的视在惯量幅值都没有显著影响(p>0.05)。宽频范围内垂向激励的引入对视在惯量的影响, 将通过3.2节偏相干分析予以进一步讨论。

    图  8  八名被试者侧倾同轴视在惯量 (a) 单轴侧倾激励; (b)侧倾−垂向双轴激励
    图  9  八名被试者单轴侧倾、侧倾−垂向激励下的侧倾同轴视在惯量中位数

    通过比较垂向−侧倾与单轴垂向、单轴侧倾条件下垂向视在质量与侧倾视在惯量的变化(见图6图9), 可以发现, 侧倾振动激励对于人体垂向一阶模态频率及视在惯量幅值的影响并不显著; 与之类似, 垂向激励对于侧倾响应也影响甚微。由于3.1节所应用的统计检验方法仅能考察一阶模态频率及模态频率处的视在质量(惯量)幅值所受影响, 本节通过偏相干分析在整个实验频率范围考察激励轴是否对人体振动响应产生了贡献。

    图10所示为双轴激励下的偏相干分析结果。图10(a)对应图2(a)垂向响应输出MISO模型。在该模型中, 垂向响应的常相干函数几乎与多重相干函数重合。也就是说, 尽管振动输入为垂向、侧倾两激励轴组合而成, 但垂向响应几乎全部由垂向振动导致, 这意味着侧倾与侧向激励的出现在整个实验测试的频率范围内基本不会对全身垂向响应产生影响。

    图  10  垂向−侧倾激励下人体动力学响应的偏相干分析结果 (a) 以垂向响应作为输出; (b)以侧倾响应作为输出

    图10(b)对应图2(b)以侧倾响应为输出的MISO模型。从图中可以观察到, 以侧倾响应作为输出的模型在3 Hz以上频率范围多重相干函数亦处于较低水平。这可能与转动响应的测试方法有关。垂向−侧倾这一测试工况中本身存在垂向激励, 即人体也必然产生明显的垂向响应; 3~10 Hz频率范围远离侧倾共振频率, 而正好覆盖垂向的多个模态频率尤其是4~5 Hz的全身模态频率, 这就造成了这一频率范围内人体侧倾响应比垂向响应微弱。而人体的侧倾响应又是通过测力台上两力传感器在垂直方向采集的数据计算得到, 从而导致了在该频率范围所测得的侧倾响应信号多重相干函数较低。理论上讲, 将人体侧倾响应和垂向响应分别测量有助于提升精度, 然而在人椅交界面上通过两个设备严格区分侧倾响应与垂向响应实际上是难以做到的。使用运动分析系统等非测力设备监测人体响应可能有助于区分两种响应, 相关测试手段在三平动轴激励下的人体运动特性分析测试中有所报告[31]。在5 Hz附近, 图10(b)中可以看到垂向偏相干函数值与侧倾常相干函数值接近, 但考虑到较低的多重相干函数值, 并不能轻易认定在这一频率范围, 垂向激励对于侧倾响应是有明显贡献的。在0.5~3 Hz的低频范围, 由于处于人体侧倾共振频率附近, 侧倾响应十分明显, 多重相干函数依然处在很高的水平, 数据具有较高的分析价值。在这一频率范围, 侧倾常相干函数基本与多重相干函数重合, 垂向偏相干函数处于极低水平。也就是说, 在人体两阶侧倾模态频率附近, 垂向振动对侧倾响应几乎没有贡献。

    通过集中质量法建立侧倾−垂向方向的坐姿人体动力学模型。由于垂向激励下的人体响应不与侧倾方向耦合, 因此, 反映垂向−侧倾激励下的人体动力学模型只需满足单轴激励下的人体动力学特性并合理设计拓扑即可, 模型垂向自由度无需与侧倾自由度耦合。首先建立模型的侧倾方向自由度, 随后附加垂向自由度, 建立面向侧倾、垂向方向振动响应预测的坐姿人体动力学模型, 模型结构如图11所示。

    图  11  侧倾−垂向人体动力学模型

    首先建立侧倾方向自由度。假设脊椎的摆动和弯曲可以通过有限个刚体结构的转动来实现, 而不需要任何的平动自由度。考虑到脊椎的层状结构, 脊椎在受侧倾激励下的横向弯曲是通过两个代表躯干结构的串联刚体反映的[32-33], 两个串联刚体称为上躯干m2和下躯干m1。需要说明的是, 由于惯性参数并不是完全给定而是通过参数识别过程进行调整, 因此上、下躯干这两个模型结构与人体结构并非完全对应。在上、下躯干的刚体之间以及躯干和座椅之间, 通过两组线性扭簧−阻尼单元连接。考虑到视在惯量测试中并没有观察到非线性现象, 这两组弹簧−阻尼单元可以在一定程度上反映椎间盘和大腿−臀部组织在侧倾方向的动力学特征。每一个躯干结构的回转中心即弹簧阻尼单元的连接位置为其几何下边缘中点。考虑到人体是关于中矢面对称的, 躯干结构的重心定义于其几何中心。考虑到二阶侧倾模态与上身肌肉−骨骼结构的水平运动有关, 在代表上躯干的刚体中心, 布置一具有水平自由度的刚体m4表示上身肌肉组织, 其通过水平方向平动弹簧、阻尼单元与上躯干刚体连接。侧倾方向共包含三个自由度, 分别为下躯干侧倾方向转动, 广义坐标定义为θ1; 上躯干侧倾方向转动, 广义坐标定义为θ2; 上身肌肉组织侧向平动, 广义坐标定义为x4; 输入侧倾振动激励的位移定义为θb

    随后, 建立模型的垂向自由度。与侧倾自由度一致, 模型包含上、下两个串联的集中质量(m2+m4m1), 同时还有一集中质量与其并联(m3), 形成在垂直方向上混联的三个自由度。这一模型引入的垂向自由度与ISO 5982-2019[34]所推荐的反映人体垂向动力学特性的三自由度混联模型接近, 同时与侧倾方向三自由度匹配。在此基础上, 还需对模型垂向自由度进行说明: (1) ISO 5982-2019[34]所推荐模型的各部分并不与人体某部分解剖学结构对应, 因此尽管在侧倾方向上, 本文所给出模型的三个自由度有大致的物理意义, 对应下躯干、上躯干、胸部肌骨结构, 但在垂向方向模型结构与身体结构的对应关系则不明确。(2) 模型垂向第二自由度的质量由对应单轴侧倾人体模型中m2m4共同组成。模型的垂向第一自由度(对应刚体m1)广义坐标定义为z1, 第二自由度(对应刚体m2+m4)定义为z2, 第三自由度定义为z3, 输入垂向振动激励的位移定义为zb

    由此, 可以建立人体模型的动力学方程。模型侧倾方向动力学方程为

    (γ1m1h2r1+m2(h1+hr1)2)¨θ1+m2(h1+hr1)hr2¨θ2+[k1+k2+k4(h1+hr1)2]θ1+[k2+k4(h1+hr1)hr2]θ2k4(h1+hr1)x4 + [c1+c2+c4(h1+hr1)2]˙θ1+[c2+c4(h1+hr1)hr2]˙θ2c4(h1+hr1)˙x4=k1θb+c1˙θb,
    (γ2m2hr22+m2hr22)¨θ2+m2hr2(h1+hr1)¨θ1+(k2+k4hr22)θ2+[k2+k4hr2(h1+hr1)]θ1k4hr2x4+(c2+c4hr22)˙θ2+[c2+c4hr2(h1+hr1)]˙θ1c4hr2˙x4=0,
    m4¨x4+k4[x4(h1+hr1)θ1hr2θ2]+c4[x4(h1+hr1)˙θ1hr2˙θ2]=0, (7)

    其中, hr1是下躯干质心与座椅连接点到下躯干质心的距离, h1是上下躯干连接点到下躯干质心的距离, hr2 是上躯干质心到上下躯干连接点的距离; ki,i = 1,2ci,i = 1,2分别为转动方向刚度和阻尼; k4c4分别为上身肌肉侧向平动刚度和阻尼; γi,i = 1,2为转动惯量计算修正系数, mi,i = 1,2为上、下躯干质量。

    模型垂直方向动力学方程为

    m1¨z1 + (z1zb)kv1+(˙z1˙zb)cv1(z2z1)kv2(˙z2˙z1)cv2=0,(m2+m4+ml+mr)¨z2+(z2z1)kv2+(˙z2˙z1)cv2 = 0,m3¨z3+(z3zb)kv3 + (˙z3˙zb)cv3 = 0, (8)

    其中, kvi, i=1,2,3、cvi, i=1,2,3为垂向自由度刚度及阻尼; m3为模型中第三垂向自由度质量。

    通过比较试验与模型获得的人椅交界面处侧倾视在惯量, 对前文所建立的人体模型中动力学参数进行识别。

    模型参数识别基于一个最优化问题, 目标函数为

    min (9)

    其中, {M_{{\text{p\_mod}}}}\left( {{f_i}} \right) {M_{{\text{p\_ang}}}}\left( {{f_i}} \right) 分别为模型预测视在惯量幅值和相位; {M_{{\text{m\_mod}}}}\left( {{f_i}} \right) {M_{{\text{m\_ang}}}}\left( {{f_i}} \right) 分别为由实验预测得到的视在惯量幅值和相位; fi为离散频率点, 本文中为0.5~10 Hz范围内以0.25 Hz为间隔的所有频率。目标函数是待识别参数 \lambda 的函数, \lambda 表示为 {\mathbf{\lambda }} = {( {\begin{array}{*{20}{c}} {{m_1}}&{{k_1}}&{{c_1}}&{{m_2}}&{{k_2}}&{{c_2}}&{{\gamma _1}}&{{\gamma _2}}&{{m_4}}&{{k_4}}&{{c_4}} \end{array}} )^{\text{T}}}

    类似地, 模型垂直方向动力学参数由以下最优化问题识别:

    {\text{min }}E({\mathbf{\lambda }}) = \sqrt {\frac{1}{N}\sum\limits_{i = 1}^N {{{\left( {\left| {{M^v}_{{\text{p\_mod}}}\left( {{f_i}} \right)} \right| - \left| {{M^v}_{{\text{m\_mod}}}\left( {{f_i}} \right)} \right|} \right)}^2}} + \frac{1}{N}\sum\limits_{i = 1}^N {{{\left( {\left| {{M^v}_{{\text{p\_ang}}}\left( {{f_i}} \right)} \right| - \left| {{M^v}_{{\text{m\_ang}}}\left( {{f_i}} \right)} \right|} \right)}^2}} } , (10)

    其中, {\lambda _v} = {( {\begin{array}{*{20}{c}} {{m_3}}&{{k_{v1}}}&{{c_{v1}}}&{{k_{v2}}}&{{c_{v2}}}&{{k_{v4}}}&{{c_{v4}}} \end{array}} )^{\text{T}}}

    侧倾−垂向人体动力学模型参数识别结果如表3所示。模型预测的垂向视在质量、侧倾视在惯量与试验实测结果的比较如图12所示, 模型与试验测试所得视在质量(视在惯量)在峰值频率与峰值频率处幅值、相位上对应良好, 垂向、侧倾方向视在质量(视在惯量)的一阶峰值频率和峰值频率处幅值都较为接近。这表明所识别的动力学参数具有一定代表性。由所识别的模型参数可知,在侧倾旋转方向上,人体下躯干即腰部附近的刚度较上躯干高,但阻尼相对更低;胸部水平运动的刚度和阻尼均很低,这符合人体胸部结构在侧倾激励下产生明显响应的现象,对应于侧倾视在惯量上观察到的二阶峰值。

    表  3  垂向−侧倾人体动力学模型动力学参数
    惯量 刚度 阻尼
    m1 27.1 kg k1 4.39 × 102 Nm/rad c1 52.9 Nms/rad
    m2 12.0 kg k2 1.68 × 102 Nm/rad c2 106 Nms/rad
    mf 1.50 kg k4 3.31 × 102 N/m c4 6.28 Ns/m
    m4 2.60 kg kv1 7.21 × 104 N/m cv1 1.02 × 103 Ns/m
    γ1 0.165 kv2 1.60 × 104 N/m cv2 683 Ns/m
    γ2 0.0892 kv3 6.96 × 104 N/m cv3 129 Ns/m
    m3 15.1 kg
    下载: 导出CSV 
    | 显示表格
    图  12  试验测量和模型预测结果 (a) 垂向视在质量; (b) 侧倾视在惯量

    本文从侧倾振动下的人体动力学特性出发, 首先测取了单轴侧倾激励下的视在惯量, 识别了人体转动方向上的模态频率; 随后研究扩展至双轴激励, 进一步测取了侧倾−垂向振动激励下的垂向同轴视在质量和侧倾同轴视在惯量, 并通过偏相干分析探究了人体垂向、侧倾振动响应间的耦合关系; 最后建立了适应于侧倾−垂向激励下响应预测的人体动力学模型。根据本文所开展的侧倾、侧倾−垂向人体动力学响应特性测试结果, 获得了以下结论:

    (1)单轴侧倾激励下, 人体在1 Hz和2.5 Hz附近出现两阶共振。人体动力学特性在侧倾方向上的非线性特征不明显。

    (2)对于垂向−侧倾组合激励下, 侧倾振动激励对于人体垂向动力学特性在0~10Hz频率范围内几乎不存在影响。与之类似, 垂向激励对于两阶全身侧倾模态频率及模态频率下视在惯量幅值的影响也十分微弱。在垂向、侧倾方向的两阶共振频率范围内, 人体侧倾、垂向动力学响应的耦合不明显;

    (3)本文建立了六自由度侧倾、垂向人体动力学模型并识别了模型中的动力学参数, 模型包含三个侧倾方向自由度以及三个垂向自由度, 具有预测侧倾、垂向双轴振动激励下人体动力学响应的能力。

  • 图  12   试验测量和模型预测结果 (a) 垂向视在质量; (b) 侧倾视在惯量

    图  1   人体动力学测试系统

    图  2   侧倾−垂向振动激励下的人体动力学响应特性分析MISO模型 (下标v, r分别表示垂向、侧倾振动输入, 下标V, R分别表示垂向、侧倾响应输出) (a) 垂向响应输出; (b) 侧倾响应输出

    图  3   单轴侧倾激励下八名被试者的视在惯量 (a) 1.0 rad/s2 r.m.s.; (b) 1.5 rad/s2 r.m.s.; (c) 2.0 rad/s2 r.m.s.

    图  4   强度为1.5 rad/s2 r.m.s.单轴侧倾激励下八名被试者的视在惯量

    图  5   单轴侧倾激励下八名被试者的视在惯量中位数

    图  6   八名被试者的垂向同轴视在质量 (a) 单轴垂向振动激励; (b) 侧倾−垂向双轴激励

    图  7   八名被试者单、多轴垂向同轴视在质量中位数

    图  8   八名被试者侧倾同轴视在惯量 (a) 单轴侧倾激励; (b)侧倾−垂向双轴激励

    图  9   八名被试者单轴侧倾、侧倾−垂向激励下的侧倾同轴视在惯量中位数

    图  10   垂向−侧倾激励下人体动力学响应的偏相干分析结果 (a) 以垂向响应作为输出; (b)以侧倾响应作为输出

    图  11   侧倾−垂向人体动力学模型

    表  1   侧倾−侧向激励下人体响应特性测试用激励信号

    激励信号垂向
    (m/s2 r.m.s.)
    侧倾
    (rad/s2 r.m.s.)
    101.0
    201.5
    302.0
    402.0
    52.00
    62.02.0
    下载: 导出CSV

    表  2   激励强度对共振频率和共振频率处视在惯量幅值影响的Wilcoxon 符号秩检验结果

    激励强度 1 rad/s2 r.m.s. 1.5 rad/s2 r.m.s.
    1.5 rad/s2 r.m.s. ns/ns/ns/*
    2.0 rad/s2 r.m.s. ns/ns/ns/ns ns/ns/ns/*
    a) ns表示影响不显著, 即p>0.05; * 表示p<0.05; **表示 p<0.01; ***表示 p<0.005, Wilcoxon;
    b) 由“/”分隔的显著性水平结果, 依次为激励强度对一阶共振频率/一阶共振频率处幅值/二阶共振频率/二阶共振频率处幅值的影响。
    下载: 导出CSV

    表  3   垂向−侧倾人体动力学模型动力学参数

    惯量 刚度 阻尼
    m1 27.1 kg k1 4.39 × 102 Nm/rad c1 52.9 Nms/rad
    m2 12.0 kg k2 1.68 × 102 Nm/rad c2 106 Nms/rad
    mf 1.50 kg k4 3.31 × 102 N/m c4 6.28 Ns/m
    m4 2.60 kg kv1 7.21 × 104 N/m cv1 1.02 × 103 Ns/m
    γ1 0.165 kv2 1.60 × 104 N/m cv2 683 Ns/m
    γ2 0.0892 kv3 6.96 × 104 N/m cv3 129 Ns/m
    m3 15.1 kg
    下载: 导出CSV
  • [1]

    Fairley T E, Griffin M J. The apparent mass of the seated human body: Vertical vibration. J. Biomech., 1989; 22(2): 81−94 DOI: 10.1016/0021-9290(89)90031-6

    [2] 高江华, 侯之超, 何乐, 等. 坐姿人体垂向振动特性及其两自由度模型. 中国科学: 技术科学, 2011; 41(12): 1640−1648 DOI: 10.1007/sl1431-011-4461-6
    [3]

    Kumar V, Saran V H. Biodynamic model of the seated human body under the vertical whole body vibration exposure. Int. J. Acoust. Vib., 2019; 24(4): 657−664 DOI: 10.20855/ijav.2019.24.41335

    [4]

    Fairley T E, Griffin M J. The apparent mass of the seated human body in the fore-and-aft and lateral directions. J. Sound Vib., 1990; 139(2): 299−306 DOI: 10.1016/0022-460X(90)90890-C

    [5]

    Huang Y, Griffin M J. Nonlinear dual-axis biodynamic response of the semi-supine human body during longitudinal horizontal whole-body vibration. J. Sound Vib., 2008; 312(1-2): 273−295 DOI: 10.1016/j.jsv.2007.10.047

    [6]

    Kuang H, Qiu Y, Liu C, et al. A nonlinear biodynamic model of seated human body exposed to vertical vibration with sensitivity analysis. Mech. Syst. Signal Process., 2023; 204: 110758 DOI: 10.1016/j.ymssp.2023.110758

    [7]

    Dong R C, He L, Du W, et al. Effect of sitting posture and seat on biodynamic responses of internal human body simulated by finite element modeling of body-seat system. J. Sound Vib., 2019; 438: 543−554 DOI: 10.1016/j.jsv.2018.09.012

    [8]

    Mansfield N J, Griffin M J. Effects of posture and vibration magnitude on apparent mass and pelvis rotation during exposure to whole-body vertical vibration. J. Sound Vib., 2002; 253(1): 93−107 DOI: 10.1006/jsvi.2001.4251

    [9]

    Nawayseh N, Sinan H A, Alteneiji S, et al. Effect of gender on the biodynamic responses to vibration induced by a whole-body vibration training machine. Proc. Inst. Mech. Eng., Part H: J. Eng. Med., 2019; 233(3): 383−392 DOI: 10.1177/0954411919830122

    [10] 张志飞, 胡正权, 徐中明, 等. 面向振动响应特性的坐姿人体动力学模型. 振动与冲击, 2016; 35(4): 104−109 DOI: 10.13465/j.cnki.jvs.2016.04.017
    [11] 白先旭, 程伟, 徐时旭, 等. 坐姿人体四自由度动力学模型研究——集中参数模型及其在汽车乘坐舒适性研究中的应用. 工程设计学报, 2017; 24(6): 638−647 DOI: 10.3785/j.issn.1006-754X.2017.06.005
    [12]

    Wei L, Griffin M J. Mathematical models for the apparent mass of the seated human body exposed to vertical vibration. J. Sound Vib., 1998; 212(5): 855−874 DOI: 10.1006/jsvi.1997.1473

    [13]

    Zheng G, Qiu Y, Griffin M J. An analytic model of the in-line and cross-axis apparent mass of the seated human body exposed to vertical vibration with and without a backrest. J. Sound Vib., 2011; 330(26): 6509−6525 DOI: 10.1016/j.jsv.2011.06.026

    [14]

    Kitazaki S, Griffin M J. A modal analysis of whole-body vertical vibration, using a finite element model of the human body. J. Sound Vib., 1997; 200(1): 83−103 DOI: 10.1006/jsvi.1996.0674

    [15]

    Nawayseh N, Griffin M J. Non-linear dual-axis biodynamic response to vertical whole-body vibration. J. Sound Vib., 2003; 268(3): 503−523 DOI: 10.1016/S0022-460X(03)00254-2

    [16]

    Huang Y, Griffin M J. Nonlinear dual-axis biodynamic response of the semi-supine human body during vertical whole-body vibration. J. Sound Vib., 2008; 312(1-2): 296−315 DOI: 10.1016/j.jsv.2007.10.046

    [17]

    Nawayseh N, Griffin M J. Tri-axial forces at the seat and backrest during whole-body vertical vibration. J. Sound Vib., 2004; 277(1-2): 309−326 DOI: 10.1016/j.jsv.2003.09.048

    [18]

    Qiu Y, Griffin M J. Biodynamic responses of the seated human body to single-axis and dual-axis vibration. Ind. Health, 2010; 48(5): 615−627 DOI: 10.2486/indhealth.MSWBVI-26

    [19]

    Mansfield N J, Maeda S. Comparison of the apparent masses and cross-axis apparent masses of seated humans exposed to single- and dual-axis whole-body vibration. J. Sound Vib., 2006; 298(3): 841−853 DOI: 10.1016/j.jsv.2006.06.017

    [20]

    Qiu Y, Griffin M J. Biodynamic response of the seated human body to single-axis and dual-axis vibration: Effect of backrest and non-linearity. Ind. Health, 2012; 50(1): 37−51 DOI: 10.2486/indhealth.MS1298

    [21]

    Hinz B, Blüthner R, Menzel G, et al. Apparent mass of seated men — Determination with single-and multi-axis excitations at different magnitudes. J. Sound Vib., 2006; 298(3): 788−809 DOI: 10.1016/j.jsv.2006.06.020

    [22]

    Mandapuram S, Rakheja S, Boileau P E, et al. Apparent mass and head vibration transmission responses of seated body to three translational axis vibration. Int. J. Ind. Ergon., 2012; 42(3): 268−277 DOI: 10.1016/j.ergon.2012.02.002

    [23]

    Mansfield N J, Maeda S. The apparent mass of the seated human exposed to single-axis and multi-axis whole-body vibration. J. Biomech., 2007; 40(11): 2543−2551 DOI: 10.1016/j.jbiomech.2006.10.035

    [24]

    Wu J, Qiu Y. Modelling of seated human body exposed to combined vertical, lateral and roll vibrations. J. Sound Vib., 2020; 485: 115509 DOI: 10.1016/j.jsv.2020.115509

    [25]

    Sun C, Liu C, Zheng X, et al. An analytical model of seated human body exposed to combined fore-aft, lateral, and vertical vibration verified with experimental modal analysis. Mech. Syst. Signal Process., 2023; 200: 110527 DOI: 10.1016/j.ymssp.2023.110527

    [26]

    Parsons K, Griffin M. Vibration and comfort, II. Rotational seat vibration. Ergonomics, 1982; 25(7): 631−644 DOI: 10.1080/00140138208925024

    [27]

    Beard G F, Griffin M J. Discomfort caused by low-frequency lateral oscillation, roll oscillation and roll-compensated lateral oscillation. Ergonomics, 2013; 56(1): 103−114 DOI: 10.1080/00140139.2012.729613

    [28]

    Beard G F, Griffin M J. Discomfort of seated persons exposed to low frequency lateral and roll oscillation: Effect of seat cushion. Appl. Ergon., 2014; 45(6): 1547−1557 DOI: 10.1016/j.apergo.2014.05.002

    [29]

    Millard M, Uchida T, Seth A, et al. Flexing computational muscle: modeling and simulation of musculotendon dynamics. J. Biomech. Eng., 2013; 135(2): 021005 DOI: 10.1115/1.4023390

    [30]

    Nawayseh N, Griffin M J. Tri-axial forces at the seat and backrest during whole-body fore-and-aft vibration. J. Sound Vib., 2005; 281(3-5): 921−942 DOI: 10.1016/j.jsv.2004.02.047

    [31]

    Hinz B, Menzel G, Bluthner R, et al. Seat-to-head transfer function of seated men ― Determination with single and three axis excitations at different magnitudes. Ind. Health, 2010; 48(5): 565−583 DOI: 10.2486/indhealth.MSWBVI-03

    [32]

    Matsumoto Y, Griffin M J. Modelling the dynamic mechanisms associated with the principal resonance of the seated human body. Clinical Biomech., 2001; 16: S31−S44 DOI: 10.1016/S0268-0033(00)00099-1

    [33]

    Pankoke S, Buck B, Woelfel H. Dynamic FE model of sitting man adjustable to body height, body mass and posture used for calculating internal forces in the lumbar vertebral disks. J. Sound Vib., 1998; 215(4): 827−839 DOI: 10.1006/jsvi.1998.1699

    [34]

    ISO 5982: 2019. Mechanical vibration and shock — Range of idealized values to characterize human biodynamic response under whole-body vibration. International Organization for Standardization, 2019

图(12)  /  表(3)
计量
  • 文章访问数:  311
  • HTML全文浏览量:  55
  • PDF下载量:  125
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-08-30
  • 修回日期:  2023-10-27
  • 刊出日期:  2024-03-04

目录

/

返回文章
返回