JP2011019845A - Fatigue degree measuring device - Google Patents

Fatigue degree measuring device Download PDF

Info

Publication number
JP2011019845A
JP2011019845A JP2009169603A JP2009169603A JP2011019845A JP 2011019845 A JP2011019845 A JP 2011019845A JP 2009169603 A JP2009169603 A JP 2009169603A JP 2009169603 A JP2009169603 A JP 2009169603A JP 2011019845 A JP2011019845 A JP 2011019845A
Authority
JP
Japan
Prior art keywords
head
seat
acceleration signal
fatigue
acceleration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2009169603A
Other languages
Japanese (ja)
Inventor
Hitoshi Uchida
仁 内田
Shogo Kida
正吾 木田
Yutaka Muraki
豊 村木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzuki Motor Corp
Original Assignee
Suzuki Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Suzuki Motor Corp filed Critical Suzuki Motor Corp
Priority to JP2009169603A priority Critical patent/JP2011019845A/en
Publication of JP2011019845A publication Critical patent/JP2011019845A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To measure the degree of fatigue in a continuously seated state on a seat while reducing a load to an occupant. <P>SOLUTION: This fatigue degree measuring device includes a seat acceleration sensor 10 installed at the seat 30 measuring vertical seat acceleration to generate a seat acceleration signal Xn; a head acceleration sensor 12 installed at the occupant's head 32 measuring vertical head acceleration of the head 32 to generate a head acceleration signal Yn; a band pass filter 14 allowing passing of a frequency band BD predetermined on the seat acceleration signal Xn and the head acceleration signal Yn; and a fatigue degree computing part 16 storing the seat acceleration signal Xn and the head acceleration signal Yn for a predetermined fixed time Cn and computing the degree of fatigue TR of the occupant based on a statistical value ST of the seat acceleration signal Xn and head acceleration signal Yn of the fixed time Cn. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、乗員の疲労度を評価する技術分野に関連し、特に、信号処理を用いて疲労度を測定する技術に関する。   The present invention relates to a technical field for evaluating the degree of fatigue of an occupant, and more particularly to a technique for measuring the degree of fatigue using signal processing.

自動車の乗員は、自動車の乗車時間の増加に伴って徐々に疲労を蓄積する。疲労の蓄積の度合いが小さいほど、快適な自動車であり自動車としての商品性は高いといえるため、開発段階で運転者の疲労の度合いを評価することは、シートやサスペンションの開発にとって重要な技術である。
このために、従来より、筋電図や心電図を計測し、また、血液中や唾液中の化学物質を計測して疲労状態を評価する手法が多く行われてきた。
A car occupant gradually accumulates fatigue as the car ride time increases. The smaller the degree of fatigue accumulation, the more comfortable the car and the higher the merchantability of the car, so evaluating the driver's degree of fatigue at the development stage is an important technology for the development of seats and suspensions. is there.
For this reason, conventionally, many methods have been used to measure an electromyogram or an electrocardiogram and to measure a fatigue state by measuring a chemical substance in blood or saliva.

特許文献1には、疲労感に影響する筋疲労部位を特定することを目的として、乗員の腰背部に筋電位測定用の電極を装着し、着座前及び単位時間おきに(段落0026)、錘を持たせて負荷を与えた状態の筋電位信号(段落0044)について、平均パワースペクトルを算出し(段落0046)、疲労の主観的評価との相関を求める手法が開示されている。
特許文献2には、走行状態のシートに着座した人体の疲労度を測定することを目的として、採尿によるアドレナリン分泌量を測定し(段落0018)、シートのクッション特性及び振動特性を説明変数とし、特に振動特性として頭部や腰部の上下加速度や加速度変化率を採用して(段落0019)、重回帰分析により説明変数を選択(段落0022)する手法が開示されている。この特許文献2には、アドレナリン分泌量に影響する説明変数として、背中付近のたわみ量、大腿部付近のたわみ量、頭部の2から4 [Hz] での左右加速度変化率、シートの1から2 [Hz] での上下加速度が開示されている(段落0023)。
特許文献3には、疲労への交感神経機能の影響を反映させた疲労度を測定することを目的として、心拍数や脈波等の生体信号のパワー値を交感神経活動による代償分を補正し、補正したパワー値に基づいて代謝量を疲労度として算出する手法が開示されている。
In Patent Document 1, for the purpose of specifying a muscle fatigue site that affects fatigue feeling, an electrode for measuring myoelectric potential is attached to the occupant's lower back, and the weight is measured before sitting and every unit time (paragraph 0026). A technique is disclosed in which an average power spectrum is calculated (paragraph 0046) for a myoelectric signal (paragraph 0044) in a state in which a load is applied and a correlation with a subjective evaluation of fatigue is obtained.
In Patent Document 2, for the purpose of measuring the fatigue level of the human body seated on the seat in the running state, the amount of adrenaline secretion by urine collection is measured (paragraph 0018), and the cushion characteristics and vibration characteristics of the seat are used as explanatory variables. In particular, a method is disclosed in which vertical accelerations and acceleration change rates of the head and waist are adopted as vibration characteristics (paragraph 0019) and explanatory variables are selected by multiple regression analysis (paragraph 0022). In this patent document, as explanatory variables affecting the amount of adrenaline secretion, the amount of deflection near the back, the amount of deflection near the thigh, the rate of change in lateral acceleration at 2 to 4 [Hz] of the head, and 1 of the seat And the vertical acceleration at 2 [Hz] is disclosed (paragraph 0023).
In Patent Document 3, for the purpose of measuring the degree of fatigue reflecting the effect of sympathetic nerve function on fatigue, the power value of a biological signal such as a heart rate and a pulse wave is corrected for the compensation due to sympathetic nerve activity. A method for calculating the metabolic rate as the degree of fatigue based on the corrected power value is disclosed.

特開2002-224072号公報Japanese Patent Laid-Open No. 2002-224072 特開2003-118458号公報Japanese Patent Laid-Open No. 2003-118458 特開2009-22610号公報JP 2009-22610

上記特許文献1では、筋電図を測定するために人体に電極を貼り付け、電極から信号処理装置までの配線が必要となり、非接触で運転動作を拘束せずに測定をすることができない。
上記特許文献2では、尿を採取してアドレナリン分泌量を測定しなければならないため、シートに着座した状態を継続させた状態での疲労度を測定することができない。その他、血液や唾液を採取する従来例では、測定される乗員の負担が大きく、多数の乗員についてのデータ収集が困難となる。
上記特許文献3では、心拍数や脈波は多様な要因に応じて変動するため、単独で疲労を算出することが難しい。特に、長時間の着座に関する疲労度の測定で誤差が大きくなってしまう。
In Patent Document 1, an electrode is attached to a human body in order to measure an electromyogram, and wiring from the electrode to the signal processing device is required, and measurement cannot be performed without restraining driving operation without contact.
In Patent Document 2, since urine must be collected and the amount of adrenaline secretion must be measured, it is not possible to measure the degree of fatigue in a state where the seated state is continued. In addition, in the conventional example in which blood or saliva is collected, the burden on the occupant to be measured is large, and data collection for a large number of occupants becomes difficult.
In Patent Document 3, since the heart rate and the pulse wave fluctuate according to various factors, it is difficult to calculate fatigue alone. In particular, the error increases when measuring the degree of fatigue related to sitting for a long time.

[課題1]このように、上記従来例では、シートに連続的に着座した状態の疲労度を乗員へ低負担で測定することができない、という不都合があった。
[課題2]さらに、上記従来例では、シートに連続的に着座した状態の疲労度に基づいて当該シートの良否を判定することができない、という不都合があった。
[Problem 1] As described above, the conventional example described above has a disadvantage in that it is impossible to measure the degree of fatigue in a state where the seat is continuously seated on the occupant with a low burden.
[Problem 2] Furthermore, the conventional example described above has a disadvantage that the quality of the seat cannot be determined based on the degree of fatigue in the state of being continuously seated on the seat.

[発明の目的]本発明の目的は、乗員へ低負荷としつつ、シートに連続的に着座した状態での疲労度を測定することにある。   [Object of the invention] An object of the present invention is to measure the degree of fatigue in a state of being continuously seated on a seat while reducing the load on an occupant.

[着眼点]本発明の発明者は、継続的な振動にさらされる人体の伝達関数を想定し、頭部の振動を当該伝達関数の出力と位置づけ、疲労との相関を検討するのが良い、という点に着目した。そして、頭部の振動に関する信号の信号処理を工夫することで、上記課題を解決できるのではないか、との着想に至った。   [Aspects] The inventor of the present invention assumes a transfer function of the human body that is exposed to continuous vibration, positions the head vibration as the output of the transfer function, and examines the correlation with fatigue. We focused on this point. And it came to the idea that the said subject could be solved by devising the signal processing of the signal regarding the vibration of the head.

[課題解決手段1]実施例1に対応する第1群の本発明は、シートに設置され当該シートの垂直方向のシート加速度を測定してシート加速度信号Xnを生成するシート加速度センサと、乗員の頭部に設置され当該頭部32の垂直方向の頭部加速度を測定して頭部加速度信号Ynを生成する頭部加速度センサと、前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させるバンドパスフィルタと、前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間蓄積すると共に、当該一定時間の前記シート加速度信号及び頭部加速度信号の統計値に基づいて、前記乗員の疲労度を算出する疲労度算出部とを備えた、という構成を採っている。
これにより、上記課題1を解決した。
[Problem Solving Means 1] A first group of the present invention corresponding to the first embodiment includes a seat acceleration sensor that is installed on a seat and generates a seat acceleration signal Xn by measuring a seat acceleration in the vertical direction of the seat; A head acceleration sensor installed on the head and measuring the head acceleration in the vertical direction of the head 32 to generate a head acceleration signal Yn, and the seat acceleration signal Xn and the head acceleration signal Yn are predetermined. A band-pass filter that passes the frequency band, and the seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined time, and based on the statistical values of the seat acceleration signal and the head acceleration signal for the predetermined time. And a fatigue level calculating unit for calculating the fatigue level of the occupant.
Thereby, the said subject 1 was solved.

[課題解決手段2] 実施例2に対応する第2群の本発明は、第1群と同様に、シートに設置され当該シートの垂直方向のシート加速度を測定してシート加速度信号Xnを生成するシート加速度センサと、乗員の頭部に設置され当該頭部の垂直方向の頭部加速度を測定して頭部加速度信号Ynを生成する頭部加速度センサと、前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させるバンドパスフィルタと、前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間蓄積すると共に、当該一定時間の前記シート加速度信号及び頭部加速度信号の統計値に基づいて、前記乗員の疲労度TRを算出する疲労度算出部とを備えている。
そして、課題解決手段2は、取り替え可能にシートを保持するシート保持部と、前記シートを加振する加振部とを備えた、という構成を採っている。
これにより、上記課題2を解決した。
[Problem Solving Means 2] As in the first group, the second group of the present invention corresponding to Example 2 measures the vertical seat acceleration of the seat and generates the seat acceleration signal Xn. A seat acceleration sensor, a head acceleration sensor installed on the head of an occupant to measure a head acceleration in the vertical direction of the head to generate a head acceleration signal Yn, the seat acceleration signal Xn, and the head acceleration signal A band-pass filter that passes a predetermined frequency band for Yn, the seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined time, and the seat acceleration signal and head acceleration for the predetermined time are stored. And a fatigue level calculation unit that calculates the fatigue level TR of the occupant based on the statistical value of the signal.
The problem solving means 2 has a configuration in which a sheet holding unit that holds a replaceable sheet and a vibration unit that vibrates the sheet are provided.
Thereby, the said subject 2 was solved.

本発明は、本明細書の記載及び図面を考慮して各請求項記載の用語の意義を解釈し、各請求項に係る発明を認定すると、各請求項に係る発明は、上記背景技術等との関連において次の有利な効果を奏する。   The present invention interprets the meaning of the terms described in each claim in consideration of the description of the present specification and the drawings, and certifies the invention according to each claim. There are the following advantageous effects in relation to

[発明の作用効果1] 課題解決手段1の疲労度測定装置は、バンドパスフィルタが、シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させ、疲労度算出部が、前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間蓄積すると共に、当該一定時間の前記シート加速度信号Xn及び頭部加速度信号Ynの統計値に基づいて、前記乗員の疲労度を算出することで、シートに着座後、時間経過に応じた疲労の蓄積に伴って大きくなる頭部加速度信号の特徴を捉えて、着座状態のまま時間的制約もない状態で疲労度を算出することができる。   [Functional Effect 1 of Invention] In the fatigue level measuring apparatus of Problem Solving Means 1, the bandpass filter passes a predetermined frequency band for the seat acceleration signal Xn and the head acceleration signal Yn, and the fatigue level calculation unit includes: The seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined time, and the fatigue level of the occupant is determined based on the statistical values of the seat acceleration signal Xn and the head acceleration signal Yn for the predetermined time. By calculating the degree of fatigue after sitting on the seat, capturing the characteristics of the head acceleration signal that increases as fatigue accumulates over time, and calculating the degree of fatigue in the seated state without any time constraints Can do.

[発明の作用効果2] 課題解決手段2のシート評価装置は、 加振部が、評価対象のシートを加振し、疲労度算出部が、この加振に応じたシート加速度信号Xn及び頭部加速度信号Ynに応じて疲労度TRを算出するため、評価対象のシートでの疲労度を定量的に測定することができる。このように、シートに連続的に着座した状態の疲労度に基づいて当該シートの良否を判定することができる。   [Functional Effect 2 of Invention] In the sheet evaluation apparatus of Problem Solving Means 2, the vibration unit vibrates the evaluation target sheet, and the fatigue calculation unit calculates the seat acceleration signal Xn and the head corresponding to the vibration. Since the fatigue level TR is calculated according to the acceleration signal Yn, the fatigue level of the evaluation target sheet can be quantitatively measured. Thus, the quality of the seat can be determined based on the degree of fatigue in the state where the seat is continuously seated.

本発明の実施形態の構成例を示すブロック図である。(実施例1から2)It is a block diagram which shows the structural example of embodiment of this invention. (Examples 1 and 2) 本実施形態のハードウエア資源の構成例を示すブロック図である。(実施例1から2)It is a block diagram which shows the structural example of the hardware resource of this embodiment. (Examples 1 and 2) 乗員の頭部の振動の三次元での各成分の比率を示すグラフ図である。(実施例1から2)It is a graph which shows the ratio of each component in three dimensions of a passenger | crew's head vibration. (Examples 1 and 2) 頭部の振動のZ軸成分である頭部加速度信号の一例を示すグラフ図である。(実施例1から2)It is a graph which shows an example of the head acceleration signal which is a Z-axis component of the vibration of a head. (Examples 1 and 2) 実施例1の処理工程の一例を示すフローチャートである。(実施例1から2)3 is a flowchart illustrating an example of a processing process according to the first embodiment. (Examples 1 and 2) 加速度信号の一例を示す説明図である。(実施例1から2)It is explanatory drawing which shows an example of an acceleration signal. (Examples 1 and 2) 伝達比変化量と主観的疲労度との相関関係の一例を示すグラフ図である。(実施例1から2)It is a graph which shows an example of the correlation between the transmission ratio change amount and the subjective fatigue level. (Examples 1 and 2) 各時点tnと、伝達比Tnと、伝達比変化量との関係の一例を示す説明図である。(実施例1から2)It is explanatory drawing which shows an example of the relationship between each time tn, transmission ratio Tn, and transmission ratio change amount. (Examples 1 and 2) 実施例2の構成例を示すブロック図である。(実施例2)10 is a block diagram illustrating a configuration example of Example 2. FIG. (Example 2) 伝達比変化量の時間変化の一例を示すグラフ図である。(実施例1から2)It is a graph which shows an example of the time change of the transmission ratio change amount. (Examples 1 and 2)

発明を実施するための形態として、2つの実施例を開示する。実施例1は疲労度測定装置であり、実施例2はシート評価装置である。実施例1から2までを含めて実施形態という。
本実施形態は、自動車車室内等、エンジンや路面からの入力により、継続的に振動にさらされる状況にある人体の疲労の程度(疲労度TR)を、人体の頭部32の振動に基づいて測定しようとするものである。
Two embodiments are disclosed as modes for carrying out the invention. Example 1 is a fatigue measurement apparatus, and Example 2 is a sheet evaluation apparatus. Embodiments including Examples 1 and 2 are referred to as embodiments.
In the present embodiment, the degree of fatigue (fatigue degree TR) of a human body that is continuously exposed to vibration by an input from an engine or road surface such as in a car interior is based on the vibration of the head 32 of the human body. It is to be measured.

人体の頭部32は、骨盤から脊柱を通して人体最上部に支持されている。しかし脊柱そのものは、多くの関節が積層された柔軟な構造で、そのままでは柔らかいバネを垂直に立て、その上に重りを置いたような状況になり、頭部32を安定して支持することはできない。このため人体では、柔らかな脊柱の周囲を囲んだ多くの筋肉が発生する張力によって、姿勢を保持して頭部32を安定に保つような構造・性質を持っている。
しかし筋肉の状態は、ストレスや疲労の蓄積などによって変化するため、人体における頭部32が支持される状態も、筋肉の状態変化の影響を受けると考えられる。
実際に自動車運転時における人体へシート30から入力される振動と、頭部32の振動の状態を測定し、長時間の自動車走行にともなう変化を観測すると、一定の周波数範囲における振動のシート30から頭部32への伝達比率が変化する傾向が見られた。
本実施形態は、このことを利用して自動車運転のような人体が振動暴露される環境下において、人体の筋肉状態の変化、具体的には疲労の状態変化を、振動計測の結果化から推測・定量化しようとするものである。
The human head 32 is supported from the pelvis through the spine to the top of the human body. However, the spine itself is a flexible structure in which many joints are stacked, and as it is, a soft spring stands upright and a weight is placed on it, and the head 32 is supported stably. Can not. For this reason, the human body has a structure and properties that maintain the posture and keep the head 32 stable by the tension generated by many muscles surrounding the soft spine.
However, since the muscle state changes due to accumulation of stress, fatigue, etc., it is considered that the state in which the head 32 is supported in the human body is also affected by the change in the muscle state.
When the vibration input from the seat 30 to the human body during actual driving of the vehicle and the vibration state of the head 32 are measured and the changes associated with the long-time driving of the vehicle are observed, There was a tendency for the transmission ratio to the head 32 to change.
In the present embodiment, this is used to infer from the result of vibration measurement the change in the muscle state of the human body, specifically the change in the state of fatigue, in an environment where the human body is exposed to vibration such as driving a car.・ I want to quantify.

<1 疲労度測定装置>
<1.1 頭部z軸の統計値>
まず、本実施形態の実施例1を開示する。実施例1は、乗員に低負荷な測定で、かつ、長時間シート30に着座した状態そのものの疲労度TRを測定するために、シート加速度と頭部32加速度とを測定し、信号処理により疲労度TRを算出しようとするものである。
<1 Fatigue measuring device>
<1.1 Head z-axis statistics>
First, Example 1 of this embodiment is disclosed. In the first embodiment, the seat acceleration and the head 32 acceleration are measured in order to measure the fatigue level TR in a state where the occupant is seated on the seat 30 for a long time, and the fatigue is obtained by signal processing. The degree TR is to be calculated.

実施例1の疲労度測定装置は、その主要な要素として、シート加速度センサ10と、頭部加速度センサ12と、バンドパスフィルタ14と、疲労度算出部16とを備えている。
図1に示す例では、シート加速度センサ10は、シート30に設置され当該シート30の垂直方向のシート加速度を測定してシート加速度信号Xnを生成する。シート加速度は、路面、ブレーキ及びエンジンの影響による振動の加速度である。頭部加速度センサ12は、乗員の頭部32に設置され当該頭部32の垂直方向の頭部加速度を測定して頭部加速度信号Ynを生成する。頭部加速度は、シート加速度による人体の加速と、このシートの振動に対抗して姿勢を維持した際のシート加速度より低減される頭部32の加速度である。バンドパスフィルタ14は、前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯BDを通過させることで、ノイズ成分を除去する。
The fatigue level measuring apparatus according to the first embodiment includes a seat acceleration sensor 10, a head acceleration sensor 12, a bandpass filter 14, and a fatigue level calculation unit 16 as main elements.
In the example shown in FIG. 1, the seat acceleration sensor 10 is installed on the seat 30 and measures the seat acceleration in the vertical direction of the seat 30 to generate the seat acceleration signal Xn. The seat acceleration is the acceleration of vibration due to the influence of the road surface, the brake and the engine. The head acceleration sensor 12 is installed on the occupant's head 32 and measures the head acceleration in the vertical direction of the head 32 to generate a head acceleration signal Yn. The head acceleration is the acceleration of the head 32 that is reduced by the acceleration of the human body due to the seat acceleration and the seat acceleration when the posture is maintained against the vibration of the seat. The band pass filter 14 removes a noise component by passing a predetermined frequency band BD for the seat acceleration signal Xn and the head acceleration signal Yn.

疲労度算出部16は、前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間Cn蓄積すると共に、当該一定時間Cnの前記シート加速度信号Xn及び頭部加速度信号Ynの統計値STに基づいて、前記乗員の疲労度TRを算出する。   The fatigue degree calculation unit 16 accumulates the seat acceleration signal Xn and the head acceleration signal Yn for a predetermined time Cn, and the statistical value ST of the seat acceleration signal Xn and the head acceleration signal Yn for the predetermined time Cn. Based on the above, the fatigue level TR of the occupant is calculated.

図2に示す例では、疲労度測定装置は、ハードウエア資源として、頭部32及びシート30の加速度をそれぞれ測定する加速度センサ40と、この加速度センサ40からの信号(シート加速度信号Xn及び頭部加速度信号Yn)をそれぞれ増幅するアンプ42と、増幅された信号をデジタル信号に変換するA/Dコンバーター44と、このデジタル信号を処理する記録分析装置36とを備えている。   In the example illustrated in FIG. 2, the fatigue measurement apparatus includes, as hardware resources, an acceleration sensor 40 that measures accelerations of the head 32 and the seat 30, and signals from the acceleration sensor 40 (the seat acceleration signal Xn and the head). The amplifier 42 amplifies the acceleration signal Yn), the A / D converter 44 that converts the amplified signal into a digital signal, and the recording analyzer 36 that processes the digital signal.

記録分析装置36は、例えばコンピュータであり、デジタル信号を一時的に記憶する一次メモリー46と、このデジタル信号を処理するCPU等の演算装置48と、演算装置48の主記憶となる主記憶装置50と、不揮発にデータを記憶するハードディスク等の補助記憶装置52とを備えている。A/Dコンバーター44でデジタル信号に変換されたシート加速度信号Xn及び頭部加速度信号Ynは、補助記憶装置52に格納される。演算装置48は、疲労度測定用プログラムを実行することで、図1に示す疲労度算出部16として機能し、一定時間Cn分の統計値STを算出し、この統計値STに基づいて疲労度TRを算出する。   The recording / analyzing device 36 is, for example, a computer, and includes a primary memory 46 that temporarily stores digital signals, an arithmetic device 48 such as a CPU that processes the digital signals, and a main storage device 50 that is a main memory of the arithmetic device 48. And an auxiliary storage device 52 such as a hard disk for storing data in a nonvolatile manner. The seat acceleration signal Xn and the head acceleration signal Yn converted into digital signals by the A / D converter 44 are stored in the auxiliary storage device 52. The computing device 48 functions as the fatigue level calculation unit 16 shown in FIG. 1 by executing the fatigue level measurement program, calculates a statistical value ST for a predetermined time Cn, and based on this statistical value ST, the fatigue level Calculate TR.

図3に、運転中の人体の頭部32の振動の割合を示す。自動車運転時にあっては、頭部32に伝わる走行中の振動は主に垂直方向のZ軸成分が大きく観測され、例えば、振動成分を三次元に分割するとZ軸成分は75%となる。そして、このZ軸成分の振動の変化が疲労度と良い相関を示す。
このため、実施例1では、第1に、Z軸方向(垂直方向)の振動のみを用いて評価を行う。この例では、図1に示すシート加速度センサ10及び頭部加速度センサ12は、垂直成分であるZ軸成分の加速度を検出する。
FIG. 3 shows the vibration ratio of the head 32 of the human body during driving. When driving an automobile, the Z-axis component in the vertical direction is mainly observed in the vibration during traveling transmitted to the head 32. For example, when the vibration component is divided into three dimensions, the Z-axis component becomes 75%. The change in vibration of the Z-axis component shows a good correlation with the fatigue level.
For this reason, in Example 1, first, evaluation is performed using only vibration in the Z-axis direction (vertical direction). In this example, the seat acceleration sensor 10 and the head acceleration sensor 12 shown in FIG. 1 detect acceleration of a Z-axis component that is a vertical component.

なお、頭部32のピッチ方向回転成分は、頭部32のできるだけ上の部分について、その前後方向の動きを計測することで近似できるため、ヘッドレストなど頭部32の近傍の基準点と頭部32の距離を光学的手段などで計測することにより、非接触で頭部32の動きを測定するようにしても良い。   Note that the pitch direction rotation component of the head 32 can be approximated by measuring the movement in the front-rear direction of the portion as high as possible of the head 32, so the reference point near the head 32 such as a headrest and the head 32. The movement of the head 32 may be measured in a non-contact manner by measuring the distance by optical means or the like.

実施例1では、第2に、シート加速度信号Xn及び頭部加速度信号Ynについて、予め定められた周波数帯BDのみ抽出すると良い。
図4を参照すると、バンドパスフィルタ14が通過させる周波数帯BDは、0.5 [Hz] から20 [Hz] とすることが望ましい。0.5 [Hz] 以下の領域と、20 [Hz] 以上の領域とは、計測誤差などが大きいため、0.5 [Hz] から 20 [Hz]の帯域での比較が有効である。
In the first embodiment, secondly, it is preferable to extract only the predetermined frequency band BD for the seat acceleration signal Xn and the head acceleration signal Yn.
Referring to FIG. 4, the frequency band BD that the band pass filter 14 passes is desirably 0.5 [Hz] to 20 [Hz]. Comparison between 0.5 [Hz] and 20 [Hz] is effective because the measurement error is large between the region below 0.5 [Hz] and the region above 20 [Hz].

そして、シート30の座面から一定の振動を人体に継続して入力した際、頭部32の振動についてZ軸方向の加速度を観察すると、初期の時点t1での頭部32の頭部加速度信号Y1のパワーは図4の点線で示す状態であり、所定時間経過して疲労が蓄積された時点t2では、頭部加速度信号Y2は図4の実線で示す状態に変化する。
図4に示す例では、周波数帯BDにて示す範囲にて時点t2での実線で示す頭部加速度信号Y2のパワーがY1よりも大きい。
When a constant vibration is continuously input from the seating surface of the seat 30 to the human body, when the acceleration in the Z-axis direction is observed for the vibration of the head 32, the head acceleration signal of the head 32 at the initial time t1. The power of Y1 is in the state indicated by the dotted line in FIG. 4, and at time t2 when fatigue is accumulated after a predetermined time has elapsed, the head acceleration signal Y2 changes to the state indicated by the solid line in FIG.
In the example shown in FIG. 4, the power of the head acceleration signal Y2 indicated by the solid line at the time point t2 is larger than Y1 in the range indicated by the frequency band BD.

さらに、実施例1では、第3に、図4に点線及び実線で示す周波数成分等を得るため、頭部加速度信号Ynを一定時間Cn蓄積し、その蓄積した信号の周波数を算出すると良い。この一定時間Cnとしては、数分から十数分程度が望ましい。   Furthermore, in the first embodiment, thirdly, in order to obtain frequency components and the like indicated by the dotted line and the solid line in FIG. 4, it is preferable to accumulate the head acceleration signal Yn for a certain period of time Cn and calculate the frequency of the accumulated signal. The fixed time Cn is preferably about several minutes to several tens of minutes.

図4に点線で示すように、着席直後の時点t1ではシート30の振動を人体で減衰させ頭部32の位置を安定させている状態では、頭部加速度信号Y1の振動成分は小さい。その後、図4に実線で示すように、所定の時間が経過して乗員に疲労が蓄積されると、頭部32の位置を安定させる人体の機能が弱まり、頭部32の振動が大きくなる。すると、時点t2では、図4に示すように、頭部加速度信号Y2のパワーは時点t1の頭部加速度信号Y1のパワーよりも大きくなる。   As indicated by a dotted line in FIG. 4, at the time t1 immediately after sitting, in the state where the vibration of the seat 30 is attenuated by the human body and the position of the head 32 is stabilized, the vibration component of the head acceleration signal Y1 is small. Thereafter, as shown by a solid line in FIG. 4, when fatigue is accumulated in the occupant after a predetermined time has elapsed, the function of the human body that stabilizes the position of the head 32 is weakened, and the vibration of the head 32 is increased. Then, at time t2, as shown in FIG. 4, the power of the head acceleration signal Y2 becomes larger than the power of the head acceleration signal Y1 at time t1.

この時間経過に伴う頭部加速度信号Ynの特徴を信号処理により効果的で安定して計算するために、実施例1では、第4に、頭部加速度信号Yn単体ではなく、シート加速度信号Xnと頭部加速度信号Ynとに基づいた処理をする。すなわち、まず、頭部加速度信号Y1の振動成分には、シート30の振動による成分と、シート30の振動に抗して頭部32の位置を安定させようとする人体の筋力による減衰を受けた振動成分とが含まれている。そして、この状態は、シート30の振動を入力とし、頭部32の振動を出力とした伝達系と考えることができる。さらに、頭部加速度信号Y1を単体で取り扱うのではなく、シート加速度信号Xnと関連させて取り扱うことで、頭部32を安定させる筋力の疲労の程度をより直接的に定量化することができる。シート加速度信号Xnと頭部加速度信号Ynとに基づいた処理の例としては、頭部32を維持する筋力の機能を伝達関数として捉えると、両信号を除する計算が適用できる。また、頭部加速度信号Y1に残存したシート加速度信号Ynの量を検出するには、両信号の減算としても良い。   In order to effectively and stably calculate the characteristics of the head acceleration signal Yn with the passage of time by signal processing, in the first embodiment, fourthly, not the head acceleration signal Yn alone but the seat acceleration signal Xn Processing based on the head acceleration signal Yn is performed. That is, first, the vibration component of the head acceleration signal Y1 was attenuated by a component caused by the vibration of the seat 30 and by a muscle strength of the human body that tries to stabilize the position of the head 32 against the vibration of the seat 30. Contains vibration components. This state can be considered as a transmission system in which the vibration of the seat 30 is an input and the vibration of the head 32 is an output. Furthermore, the degree of muscle strength fatigue that stabilizes the head 32 can be more directly quantified by handling the head acceleration signal Y1 in association with the seat acceleration signal Xn instead of handling it alone. As an example of processing based on the seat acceleration signal Xn and the head acceleration signal Yn, if the function of the muscular strength that maintains the head 32 is regarded as a transfer function, a calculation that removes both signals can be applied. In order to detect the amount of the seat acceleration signal Yn remaining in the head acceleration signal Y1, both signals may be subtracted.

実施例1では、第5に、シート加速度信号Xnと頭部加速度信号Ynとに基づいた処理に際して、両信号の統計値STを用いる。両信号には、一定時間Cn分のデータが含まれており、全てのデータを両信号間で除算や減算しても値が不安定となってしまう。このため、実施例1では、両信号についてそれぞれ統計値STを算出しておき、この両者の統計値STを除算又は減算する。統計値STの例としては、図4に示す各時点tnでのパワーの差が良好に現れる統計値STとすることが望ましい。例えば、一定時間Cnの全体について、唯一の平均値や中央値とする。平均値にしても、加算平均や、二乗平均平方根を採用することができる。また、一定時間Cnを予め定められた周波数の区間で区切り、各区間について統計値STを算出し、周波数の区間毎に両信号の除算や減算をしても良い。   In the first embodiment, fifthly, the statistical value ST of both signals is used in the processing based on the seat acceleration signal Xn and the head acceleration signal Yn. Both signals contain data for a certain time Cn, and even if all data are divided or subtracted between the two signals, the values become unstable. For this reason, in the first embodiment, the statistical value ST is calculated for both signals, and the statistical value ST of both is divided or subtracted. As an example of the statistical value ST, it is desirable to use the statistical value ST in which the difference in power at each time point tn shown in FIG. For example, a single average value or median value is set for the entire fixed time Cn. As the average value, an addition average or a root mean square can be adopted. Alternatively, the predetermined time Cn may be divided into predetermined frequency sections, the statistical value ST may be calculated for each section, and the two signals may be divided or subtracted for each frequency section.

このように、第1に、シート30及び頭部32の垂直方向(Z軸方向)の振動を測定し、第2に、予め定められた周波数帯BD内の信号を取り出し、第3に、一定時間Cn分の蓄積をし、第4に、シート加速度信号Xnと頭部加速度信号Ynとに基づいた処理をし、そして、第5に、シート加速度信号Xn及び頭部加速度信号Ynのそれぞれの統計値STを算出して統計値STの除算や減算をすることで、乗員の疲労度TRを良好に算出することができる。そして、統計値STの除算や減算の値を、時点t1, t2, tnで比較することで、疲労の大きさに応じた疲労度TRを算出することができる。   Thus, firstly, the vibration in the vertical direction (Z-axis direction) of the seat 30 and the head 32 is measured, and secondly, a signal in a predetermined frequency band BD is taken out, and thirdly, constant Accumulate time Cn, and fourth, process based on the seat acceleration signal Xn and the head acceleration signal Yn, and fifth, the statistics of the seat acceleration signal Xn and the head acceleration signal Yn, respectively. By calculating the value ST and dividing or subtracting the statistical value ST, the occupant's fatigue level TR can be calculated satisfactorily. Then, by comparing the division and subtraction values of the statistical value ST at the time points t1, t2, and tn, the fatigue level TR corresponding to the magnitude of fatigue can be calculated.

実施例1では、疲労度の指標算出に使用する頭部振動として、頭部Z軸方向並進加速度の0.5-20Hzの領域を利用しているが、頭部振動のピッチ方向回転成分(頭が頷く運動をする方向)も頭部Z軸並進の代替信号として使用可能で、頭部Z軸方向並進振動の測定と類似した結果を得ることができる。   In the first embodiment, a 0.5-20 Hz region of translational acceleration in the head Z-axis direction is used as the head vibration used for calculating the fatigue index. The direction of movement) can also be used as an alternative signal for the head Z-axis translation, and a result similar to the measurement of translational vibration in the head Z-axis direction can be obtained.

・1.1頭部z軸の統計値の効果
上述のように、バンドパスフィルタ14が、シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯BDを通過させ、疲労度算出部16が、前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間Cn蓄積すると共に、当該一定時間Cnの前記シート加速度信号Xn及び頭部加速度信号Ynの統計値STに基づいて、前記乗員の疲労度TRを算出することで、シート30に着座後、時間経過に応じた疲労の蓄積に伴って大きくなる頭部加速度信号Ynの特徴を捉えて、着座状態のまま時間的制約もない状態で疲労度TRを算出することができる。さらに、両信号の一定時間Cnの統計値STがあれば疲労度TRを算出できるため、着座状態のまま連続して同一の測定条件下で複数回の測定をすることができる。
1.1 Effect of Head Z-Axis Statistics As described above, the bandpass filter 14 passes the predetermined frequency band BD for the seat acceleration signal Xn and the head acceleration signal Yn, and the fatigue level calculation unit 16 The seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined time Cn, and the occupant is based on the statistical value ST of the seat acceleration signal Xn and the head acceleration signal Yn for the predetermined time Cn. By calculating the degree of fatigue TR, the head acceleration signal Yn, which increases as the fatigue accumulates over time after the seat 30 is seated, is captured, and there is no time restriction in the seated state. It is possible to calculate the fatigue level TR. Furthermore, since the fatigue level TR can be calculated if there is a statistical value ST of the constant time Cn of both signals, a plurality of measurements can be performed continuously under the same measurement conditions while in the sitting state.

<1.2伝達比と伝達比変化量>
実施例1では、疲労及び疲労の蓄積の程度と極めて相関の高い測定をすることができる。この例では、主要な要素として、疲労度算出部16が、パワー算出処理18と、伝達比算出処理20とを備える。また、好ましい例では、伝達比変化量算出処理22を備えると良い。
<1.2 Transmission ratio and transmission ratio change>
In Example 1, it is possible to perform measurement that is highly correlated with fatigue and the degree of fatigue accumulation. In this example, the fatigue level calculation unit 16 includes a power calculation process 18 and a transmission ratio calculation process 20 as main elements. In a preferred example, a transmission ratio change amount calculation process 22 may be provided.

パワー算出処理18は、前記シート加速度信号Xnの前記一定時間Cnの二乗平均平方根XRnと、前記頭部加速度信号Ynの前記一定時間Cnの二乗平均平方根YRnとを算出する。
パワー算出処理18は、計測開始からある時点t1で一定時間Cn計測したシート加速度信号Xn、頭部加速度信号Ynに、周波数帯BDを0.1 [Hz] から20 [Hz] としたバンドパスフィルタ14を適用した後の信号について、信号の個々の値をxn, ynとした際、そのシグナルパワーに相当するXR1, YR1を次式(1)及び(2)で算出する。
$
図2に示すA/Dコンバーター44は、一定の間隔で加速度信号をA/D変換し、サンプリングする。すると、シート加速度信号Xn及び頭部加速度信号Ynは、それぞれ、N個の数値データとして収集される。
The power calculation process 18 calculates the root mean square XRn of the fixed time Cn of the seat acceleration signal Xn and the root mean square YRn of the fixed time Cn of the head acceleration signal Yn.
The power calculation process 18 uses a bandpass filter 14 that changes the frequency band BD from 0.1 [Hz] to 20 [Hz] on the seat acceleration signal Xn and head acceleration signal Yn measured for a certain time Cn at a certain time t1 from the start of measurement. With respect to the signal after application, when the individual values of the signal are xn and yn, XR1 and YR1 corresponding to the signal power are calculated by the following equations (1) and (2).
$
The A / D converter 44 shown in FIG. 2 performs A / D conversion on the acceleration signal at regular intervals and samples it. Then, the seat acceleration signal Xn and the head acceleration signal Yn are collected as N pieces of numerical data, respectively.

X1 = {x1, x2, x3, ...xN} ... 0.1 [Hz] から20 [Hz] のシート加速度信号Xn
Y1 = {y1, y2, y3, ...yN} ... 0.1 [Hz] から20 [Hz] の頭部加速度信号Yn
X1 = {x1, x2, x3, ... xN} ... Seat acceleration signal Xn from 0.1 [Hz] to 20 [Hz]
Y1 = {y1, y2, y3, ... yN} ... Head acceleration signal Yn from 0.1 [Hz] to 20 [Hz]

Figure 2011019845
Figure 2011019845

パワー算出処理18は、同様に、次のある時点t2での二乗平均平方根XR2, YR2を算出する。そして、測定時点tnが多数ある際には、各測定時点tnの一定時間Cnの個々の信号の値から二乗平均平方根を算出する。   Similarly, the power calculation process 18 calculates the root mean square XR2 and YR2 at the next certain time point t2. Then, when there are a large number of measurement time points tn, the root mean square is calculated from the values of the individual signals for a fixed time Cn at each measurement time point tn.

伝達比算出処理20は、次式(3)に従い、前記シート加速度信号Xnの二乗平均平方根XRnで前記頭部加速度信号YRnの二乗平均平方根YRnを除することで、前記シート30の振動の前記頭部32の振動への伝達比Tnを算出する。すなわち、伝達比算出処理20は、入力に相当するXRnで、出力に相当するYRnを除して比を算出し、これをある時点での振動の伝達比T1とする。この伝達比は、時点tnにおける乗員の疲労度TRを定量化した値である。   The transmission ratio calculation process 20 divides the root mean square YRn of the head acceleration signal YRn by the root mean square XRn of the seat acceleration signal Xn in accordance with the following equation (3), whereby the head of the vibration of the seat 30 is calculated. A transmission ratio Tn to vibration of the unit 32 is calculated. That is, the transmission ratio calculation process 20 calculates the ratio by dividing XRn corresponding to the input and YRn corresponding to the output, and this is set as the vibration transmission ratio T1 at a certain time. This transmission ratio is a value obtained by quantifying the occupant fatigue level TR at the time point tn.

Figure 2011019845
Figure 2011019845

好ましい例では、さらに、伝達比変化量算出処理22が、前記一定時間Cnを超える時間間隔の2つの時点tnについて、前記伝達比Tnを算出させると共に、次式(4)に従って、当該2つの時点tnの各伝達比Tnの差の絶対値を前記疲労度TRと相関する伝達比変化量DTnとして算出する。すなわち、時点t1から所定の時間経過後の時点t2において、同様にシート加速度信号X2及び頭部加速度信号Y2を用いて伝達比T2を計算し、このT2と時点t1で算出したT1との差の絶対値を、実施例1での伝達比変化量とする。   In a preferred example, the transmission ratio change amount calculation process 22 further calculates the transmission ratio Tn for two time points tn in the time interval exceeding the certain time Cn, and the two time points according to the following equation (4): The absolute value of the difference between the transmission ratios Tn of tn is calculated as the transmission ratio change amount DTn correlated with the fatigue level TR. That is, at time t2 after a predetermined time has elapsed from time t1, similarly, the transmission ratio T2 is calculated using the seat acceleration signal X2 and the head acceleration signal Y2, and the difference between T2 and T1 calculated at time t1 is calculated. The absolute value is the transmission ratio change amount in the first embodiment.

Figure 2011019845
Figure 2011019845

図5を参照すると、疲労度測定装置は、まず、シート30に設置され当該シート30の垂直方向(Z軸方向)のシート加速度を測定して図6に示すシート加速度信号Xnを補助記憶装置50に記録する(ステップS01,シート加速度記録工程)。このステップS01と並行して、乗員の頭部32に設置され当該頭部32の垂直方向(Z軸方向)の頭部加速度を測定して図6に示す頭部加速度信号Ynを補助記憶装置50に記録する(ステップS02,頭部加速度記録工程)。そして、前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯BDを通過させることで、ノイズ成分を除去する(ステップS03,S04,バンドパスフィルタ工程)。周波数帯BDとしては、上述のように、0.5 [Hz] から 20 [Hz]とすると良い。   Referring to FIG. 5, the fatigue level measuring apparatus is first installed on the seat 30 to measure the seat acceleration in the vertical direction (Z-axis direction) of the seat 30, and the seat acceleration signal Xn shown in FIG. (Step S01, sheet acceleration recording step). In parallel with this step S01, the head acceleration in the vertical direction (Z-axis direction) installed on the head 32 of the passenger is measured, and the head acceleration signal Yn shown in FIG. (Step S02, head acceleration recording step). Then, the noise component is removed by passing a predetermined frequency band BD with respect to the seat acceleration signal Xn and the head acceleration signal Yn (steps S03, S04, band pass filter step). As described above, the frequency band BD is preferably 0.5 [Hz] to 20 [Hz].

そして、疲労度算出部16は、図6に示す予め定められた一定時間Cnの前記シート加速度信号Xn及び頭部加速度信号Ynの統計値STに基づいて、前記乗員の疲労度TRを算出する(疲労度算出工程)。この疲労度算出工程では、まず、前記一定時間Cnを超える時間間隔の2つの時点tnについて、前記シート加速度信号Xnの統計値STで前記頭部加速度信号Ynの統計値STを除することで、前記シート30の振動の前記頭部32の振動への伝達比Tnを算出する(ステップS05,伝達比算出工程)。統計値STとしては、例えば、二乗平均平方根XRn, YRnとすると良い。   Then, the fatigue level calculation unit 16 calculates the occupant fatigue level TR based on the statistical value ST of the seat acceleration signal Xn and the head acceleration signal Yn for a predetermined time Cn shown in FIG. Fatigue degree calculation process). In this fatigue level calculation step, first, for two time points tn of the time interval exceeding the predetermined time Cn, the statistical value ST of the head acceleration signal Yn is divided by the statistical value ST of the seat acceleration signal Xn, A transmission ratio Tn of the vibration of the seat 30 to the vibration of the head 32 is calculated (step S05, transmission ratio calculation step). The statistical value ST may be, for example, root mean square XRn, YRn.

この伝達比Tnについて、運転前の時点t1と所定の時間経過後の時点t2の値を比較すると、図7に示すように、筋肉状態の変化や姿勢の崩れを原因としてその比率が減少または増大する。すなわち、前記2つの時点tnの各伝達比Tnの差の絶対値を前記疲労度TRと相関する伝達比変化量DTnとして算出する(ステップS06,伝達比変化量算出工程)。
なお、2つの条件、例えば、2種類のシート30の比較やサスペンションの比較をするためには、伝達比Tnの差の絶対値ではなく、単純に差として、この差の値が2種類の測定結果で異なるか否かという判定をしても良い。
As for the transmission ratio Tn, when the value of the time point t1 before driving and the value of the time point t2 after elapse of a predetermined time are compared, as shown in FIG. 7, the ratio decreases or increases due to a change in muscle state or a collapse of posture. To do. That is, the absolute value of the difference between the transmission ratios Tn at the two time points tn is calculated as the transmission ratio change amount DTn correlated with the fatigue level TR (step S06, transmission ratio change amount calculation step).
It should be noted that in order to compare two conditions, for example, two types of seats 30 and suspensions, it is not an absolute value of a difference in transmission ratio Tn, but simply a difference, and this difference value is measured in two types. It may be determined whether the results are different.

この伝達比Tnは、乗員によって増大する場合と、減少する場合とがある。しかし、総じて、人体の骨格がとる姿勢やそれをささえる筋肉の状態が、初期(疲労していない状態)から変化した大きさに応じて変化すると考えられる。そして、伝達比変化量DTnは、図7に示すように、人間が感じる主観的な疲労度と良い相関を示す。   This transmission ratio Tn may be increased or decreased by the occupant. However, as a whole, it is considered that the posture of the human skeleton and the state of the muscles that support it change according to the magnitude that has changed from the initial stage (the state that is not fatigued). As shown in FIG. 7, the transmission ratio change amount DTn shows a good correlation with the subjective fatigue level felt by humans.

実施例1では、このことを利用して、 0.5 [Hz] から 20 [Hz]の周波数帯BDにおける垂直方向のシート30の座面の振動エネルギーと乗員の頭部32の振動エネルギー比率の変化を観測することで、人体の疲労度TRを推定する。   In the first embodiment, by utilizing this, the change in the vibration energy ratio of the seat surface of the seat 30 in the vertical direction and the vibration energy ratio of the occupant's head 32 in the frequency band BD of 0.5 [Hz] to 20 [Hz] is obtained. By observing it, the human fatigue level TR is estimated.

図8に示すように、実施例1記載の手法では、長時間の着座について同一条件下で複数回の測定を安定して行うことができ、そして、時点t1, t2, t3についてシート加速度信号X1, X2, X3及び頭部加速度信号Y1, Y2, Y3を得て、伝達比T1, T2, T3を算出し、伝達比変化量DT1_2, DT2_3のみならず、DT1_3を算出することもできる。   As shown in FIG. 8, in the method described in the first embodiment, a plurality of measurements can be stably performed under the same conditions for a long time seating, and the seat acceleration signal X1 at time points t1, t2, and t3. , X2, X3 and head acceleration signals Y1, Y2, Y3 are obtained and transmission ratios T1, T2, T3 are calculated, and not only transmission ratio change amounts DT1_2, DT2_3 but also DT1_3 can be calculated.

図5に示す処理工程は、図2に示す演算装置48が疲労度測定用プログラムを実行することで、実現することができる。この疲労度測定用プログラムは、図1のパワー算出処理18、伝達比算出処理20、伝達比変化量算出処理22の機能を発揮するためのプログラムコードを備えると良い。   The processing steps shown in FIG. 5 can be realized by the arithmetic unit 48 shown in FIG. 2 executing the fatigue degree measurement program. This fatigue degree measurement program may include program codes for exhibiting the functions of the power calculation process 18, the transmission ratio calculation process 20, and the transmission ratio change amount calculation process 22 of FIG.

・1.2 伝達比と伝達比変化量の効果
上述のように、伝達比算出処理20が、統計値STとして伝達比Tnを算出し、伝達比変化量算出処理22が、時点t1の伝達比T1と時点t2の伝達比T2との差(の絶対値)DT1_2を算出し、これを疲労度TRとするため、図7に示すように主観的疲労度との相関の高い疲労度TRを着座状態のままで、乗員への負荷が小さい状態で、かつ、長時間連続的な測定環境を構築しつつ、複数の時点Tn間で蓄積された疲労度TRを定量化することができる。
1.2 Effect of Transmission Ratio and Transmission Ratio Change As described above, the transmission ratio calculation process 20 calculates the transmission ratio Tn as the statistical value ST, and the transmission ratio change calculation process 22 calculates the transmission ratio T1 at the time point t1. To calculate the difference (absolute value) DT1_2 from the transmission ratio T2 at time t2 and set this as the fatigue level TR, the fatigue level TR having a high correlation with the subjective fatigue level is set as shown in FIG. The fatigue level TR accumulated between a plurality of time points Tn can be quantified while the measurement environment is continuously constructed for a long time with a small load on the occupant.

このように、実施例1では、人体の頭部32の垂直方向並進加速度について、 0.5 [Hz] から 20 [Hz] の周波数帯を抽出して、数秒から数十秒間の二乗平均平方根(統計値ST)を算出する。そして、人体が着座するシート30の座面についても、同様に0.5-10Hzの周波数帯で垂直方向並進加速度の二乗平均の平方根を算出する。さらに、算出した二乗平均値について、人体の頭部32の値をシート30の値で除して、信号の伝達比Tnを算出する。そして、一定時間経過後の信号伝達比T2と、初期の伝達比T1との差の絶対値DT1_2を算出して、その大きさから人体の疲労度TRを推測する。   As described above, in the first embodiment, the frequency band from 0.5 [Hz] to 20 [Hz] is extracted for the vertical translational acceleration of the human head 32, and the root mean square (statistical value) for several seconds to several tens of seconds is extracted. ST) is calculated. Similarly, for the seating surface of the seat 30 on which the human body is seated, the root mean square of the vertical translational acceleration is calculated in the frequency band of 0.5-10 Hz. Further, with respect to the calculated mean square value, the value of the human head 32 is divided by the value of the seat 30 to calculate a signal transmission ratio Tn. Then, the absolute value DT1_2 of the difference between the signal transmission ratio T2 after a predetermined time and the initial transmission ratio T1 is calculated, and the fatigue level TR of the human body is estimated from the magnitude.

このため、人体の疲労の状態について、手間がかかり比較的誤差が大きくなる主観的官能評価を用いず、外部から電気信号として測定できる加速度データを利用して推測できる。そして、0.5 [Hz] から 20 [Hz] の周波数帯に限定して利用することで、ノイズなどの影響を受けにくくして疲労による人体振動の変化を計測できる。
さらに、車室内での人体の頭部32の振動は垂直並進方向の成分が多く、この方向のみを計測・評価することで、3軸全てを計測・評価するのと同等な効果を簡便に得られるまた、人体の頭部32の振動は、車体側の振動の大きさに影響を受けるため、シート30の座面振動との比率を用いることで、車両走行状態の観測値への影響を最小限とする。
そして、シート30の性質や個人差などによる疲労状態の違いにより、単純な数値の増減では疲労状態評価がやり難い。初期の伝達比と、一定時間経過後の伝達比の変化量絶対値をとることで、主観による疲労感の大小と相関のよい物理値を得ることができる。
For this reason, the state of fatigue of the human body can be inferred using acceleration data that can be measured as an electrical signal from the outside without using subjective sensory evaluation, which is laborious and relatively large in error. By limiting the frequency range from 0.5 [Hz] to 20 [Hz], it is possible to measure changes in human vibration due to fatigue without being affected by noise.
Furthermore, the vibration of the human head 32 in the passenger compartment has many components in the vertical translation direction. By measuring and evaluating only this direction, the same effects as measuring and evaluating all three axes can be obtained easily. In addition, since the vibration of the head 32 of the human body is affected by the magnitude of the vibration on the vehicle body side, the influence on the observed value of the vehicle running state is minimized by using the ratio with the seat surface vibration of the seat 30. Limit.
Further, due to the difference in the fatigue state due to the nature of the sheet 30 and individual differences, it is difficult to evaluate the fatigue state by simply increasing or decreasing the numerical value. By taking the initial transmission ratio and the absolute value of the amount of change in the transmission ratio after a lapse of a certain time, a physical value having a good correlation with the subjective fatigue level can be obtained.

<2 シート評価装置>
次に、実施例2を開示する。
実施例2のシート評価装置は、その主要な要素として、実施例1の主要な要素に加え、シート保持部24と、加振部26とを備えている。
図8に示す例では、さらに、スクリーン28を備えると共に、図2に示すアンプ42及び記録分析装置36を備えている。
<2 Sheet evaluation device>
Next, Example 2 is disclosed.
The sheet evaluation apparatus according to the second embodiment includes a sheet holding unit 24 and a vibration unit 26 as main elements in addition to the main elements according to the first embodiment.
In the example shown in FIG. 8, the screen 28 is further provided, and the amplifier 42 and the recording analyzer 36 shown in FIG.

シート保持部24は、取り替え可能にシート30を保持する。四輪車用のシート30を取付可能な機構を有すると良い。加振部26は、前記シート30を加振する。この加振は、実際の四輪車の走行状態をシミュレートするための加振であり、予め設定された走行環境と車速とに応じた振動をシート30に与えると良い。スクリーン28に、走行環境の画像を記憶し乗員の操作に応じて走行環境を表示制御するコンピュータを併設し、スクリーン28に走行環境を表示することで、擬似的な運転環境を提供すると良い。   The sheet holding unit 24 holds the sheet 30 in a replaceable manner. It is preferable to have a mechanism capable of attaching a seat 30 for a four-wheel vehicle. The vibration unit 26 vibrates the sheet 30. This excitation is an excitation for simulating the actual traveling state of the four-wheeled vehicle, and it is preferable to give the seat 30 vibration according to a preset traveling environment and vehicle speed. It is preferable to provide a pseudo driving environment by storing an image of the driving environment on the screen 28 and providing a computer that displays and controls the driving environment according to the operation of the occupant and displaying the driving environment on the screen 28.

シート加速度センサ10と、頭部加速度センサ12と、バンドパスフィルタ14と、疲労度算出部16(記録分析装置36)の構成は実施例1と同様である。   The configurations of the seat acceleration sensor 10, the head acceleration sensor 12, the bandpass filter 14, and the fatigue level calculation unit 16 (recording analysis device 36) are the same as those in the first embodiment.

図9に示すシート評価装置は、実施例1による測定を行うための加振機能をもったシミュレータである。運転席床面に接続された加振部26によって、シート30を経由して人体に自動車走行を模擬した振動を与えながら、予め定められた時間の運転シミュレーションを行って、シート30の座面と乗員の頭部32に固定した加速度センサ40によって、頭部加速度信号Xnとシート加速度信号Ynとのデータを収集する。頭部加速度信号Xnとシート加速度信号Ynを収集すると、図5に示す手順に従って伝達比Tnを算出する。   The sheet evaluation apparatus shown in FIG. 9 is a simulator having a vibration function for performing measurement according to the first embodiment. A vibration simulation unit 26 connected to the driver's seat floor surface performs a driving simulation for a predetermined period of time while applying vibrations simulating automobile travel to the human body via the seat 30. Data of the head acceleration signal Xn and the seat acceleration signal Yn is collected by the acceleration sensor 40 fixed to the occupant's head 32. When the head acceleration signal Xn and the seat acceleration signal Yn are collected, the transmission ratio Tn is calculated according to the procedure shown in FIG.

具体的な使用例としては、例えば、自動車の乗員用のシート30について、何種類かを比較した場合に、一定時間の走行実験、あるいは自動車の走行振動を模擬したシミュレータ実験を行ったあとに、シート30からの入力振動と、人体の頭部32の振動を用いて伝達比変化量DTnを算出して比較することで、シート30の違いによる乗員の疲労への影響の大小を比較することができ、より疲れない自動車のシート30の開発を支援することが可能となる。   As a specific use example, for example, when comparing several types of seats 30 for an automobile occupant, after performing a driving experiment for a certain period of time or a simulator experiment simulating the driving vibration of an automobile, By calculating and comparing the transmission ratio change amount DTn using the input vibration from the seat 30 and the vibration of the head 32 of the human body, it is possible to compare the magnitude of the influence on the occupant's fatigue due to the difference in the seat 30. This makes it possible to support the development of an automobile seat 30 that is less tired.

図10に、シート評価装置を用いて、2種類のシート30について、実験開始直後の伝達比をT1として、同じシミュレータ実験を行いながら、時間経過に伴って3回の伝達比計測を行い、それぞれについてT1との間で伝達比変化量を算出した例を示す。図10に示すように、2種類のシート30について、双方とも時間経過に比例して疲労度TRが増加すること、点線で示すシート30の方が大きい疲労度をもたらしている。このように2つのシート30の間で伝達比の変化に差があった場合、伝達比変化量の時間経過に伴う増加が少ないシート30は疲労が少なく、変化の多いシート30が疲労の多いシート30であると推測することができる。
また、シート30の軽量化やコストダウンなどの改良を行う際、構造変更によって人体の疲労感を増加させるような、性能の劣化が出ているか否かを評価することが可能となる。
In FIG. 10, using the sheet evaluation apparatus, with respect to two types of sheets 30, the transmission ratio immediately after the start of the experiment is set to T1, and the same simulator experiment is performed, and the transmission ratio is measured three times over time. Shows an example of calculating the transmission ratio change amount with respect to T1. As shown in FIG. 10, for both types of sheets 30, the fatigue level TR increases in proportion to the passage of time, and the sheet 30 indicated by the dotted line brings about a higher level of fatigue. In this way, when there is a difference in the transmission ratio change between the two sheets 30, the sheet 30 in which the amount of change in the transmission ratio change is small with time has little fatigue, and the sheet 30 with much change has a lot of fatigue. It can be estimated that it is 30.
In addition, when performing improvements such as weight reduction and cost reduction of the seat 30, it is possible to evaluate whether or not performance deterioration has occurred that increases the fatigue of the human body due to structural changes.

・2 シート評価装置の効果
上述のように、実施例2では、加振部26が、評価対象のシート30を加振し、疲労度算出部16が、この加振に応じたシート加速度信号Xn及び頭部加速度信号Ynに応じて疲労度TRを算出するため、評価対象のシート30での疲労度TRを定量的に測定することができる。この疲労度TRの評価については、実施例1と同様の効果を奏する他、シート30に連続的に着座した状態の疲労度TRに基づいて、複数のシート30の良否を比較し、判定することができる。
-2 Effect of Sheet Evaluation Device As described above, in Example 2, the vibration unit 26 vibrates the evaluation target sheet 30, and the fatigue level calculation unit 16 performs the sheet acceleration signal Xn according to the vibration. Since the fatigue level TR is calculated in accordance with the head acceleration signal Yn, the fatigue level TR on the evaluation target seat 30 can be quantitatively measured. Regarding the evaluation of the fatigue level TR, in addition to the same effect as in the first embodiment, the quality of the plurality of seats 30 is compared and determined based on the fatigue level TR in the state of being continuously seated on the seat 30. Can do.

10 シート加速度センサ
12 頭部加速度センサ
14 バンドパスフィルタ
16 疲労度算出部
18 パワー算出処理
20 伝達比算出処理
22 伝達比変化量算出処理
24 シート保持部
26 加振部
28 スクリーン
30 シート
32 頭部
36 記録分析装置
40 加速度センサ
42 アンプ
44 A/Dコンバーター
46 一次メモリー
48 演算装置
50 主記憶装置
52 補助記憶装置
Xn, X1, X2, X3 シート加速度信号
Yn, Y1, Y2, Y3 頭部加速度信号
BD 周波数帯
Cn 一定時間
ST 統計値
TR 疲労度
XRn, YRn 二乗平均平方根
Tn, T1, T2, T3 伝達比
tn, t1, t2, t3 時点
DTn, DT1_2, DT1_3, DT2_3 伝達比変化量
DESCRIPTION OF SYMBOLS 10 Sheet acceleration sensor 12 Head acceleration sensor 14 Band pass filter 16 Fatigue degree calculation part 18 Power calculation process 20 Transmission ratio calculation process 22 Transmission ratio change amount calculation process 24 Sheet holding part 26 Excitation part 28 Screen 30 Sheet 32 Head part 36 Recording analysis device 40 Acceleration sensor 42 Amplifier 44 A / D converter 46 Primary memory 48 Arithmetic device 50 Main storage device 52 Auxiliary storage device
Xn, X1, X2, X3 Seat acceleration signal
Yn, Y1, Y2, Y3 Head acceleration signal
BD frequency band
Cn fixed time
ST statistics
TR fatigue
XRn, YRn root mean square
Tn, T1, T2, T3 transmission ratio
tn, t1, t2, t3
DTn, DT1_2, DT1_3, DT2_3 Transmission ratio change

Claims (6)

シートに設置され当該シートの垂直方向のシート加速度を測定してシート加速度信号Xnを生成するシート加速度センサと、
乗員の頭部に設置され当該頭部の垂直方向の頭部加速度を測定して頭部加速度信号Ynを生成する頭部加速度センサと、
前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させるバンドパスフィルタと、
前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間蓄積すると共に、当該一定時間の前記シート加速度信号及び頭部加速度信号の統計値に基づいて、前記乗員の疲労度を算出する疲労度算出部と、を備えた、
ことを特徴とする疲労度測定装置。
A seat acceleration sensor that is installed on the seat and measures a seat acceleration in the vertical direction of the seat to generate a seat acceleration signal Xn;
A head acceleration sensor that is installed on the head of the occupant and measures a head acceleration in the vertical direction of the head to generate a head acceleration signal Yn;
A bandpass filter that passes a predetermined frequency band for the seat acceleration signal Xn and the head acceleration signal Yn;
The seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined period of time, and the degree of fatigue of the occupant is calculated based on statistical values of the seat acceleration signal and the head acceleration signal for the predetermined period of time. A fatigue degree calculation unit,
A fatigue measurement apparatus characterized by the above.
前記疲労度算出部が、
前記シート加速度信号Xnの前記一定時間の二乗平均平方根と、前記頭部加速度信号Ynの前記一定時間の二乗平均平方根とを算出するパワー算出処理と、
前記シート加速度信号Xnの二乗平均平方根で前記頭部加速度信号の二乗平均平方根を除することで、前記シートの振動の前記頭部の振動への伝達比を算出する伝達比算出処理と、を備えた、
ことを特徴とする請求項1記載の疲労度測定装置。
The fatigue level calculation unit
A power calculation process for calculating the root mean square of the fixed time of the seat acceleration signal Xn and the root mean square of the constant time of the head acceleration signal Yn;
A transmission ratio calculation process for calculating a transmission ratio of the vibration of the seat to the vibration of the head by dividing the root mean square of the head acceleration signal by the root mean square of the seat acceleration signal Xn. The
The fatigue degree measuring apparatus according to claim 1, wherein:
前記疲労度算出部が、
前記一定時間を超える時間間隔の2つの時点について、前記伝達比を算出させると共に、当該2つの時点の各伝達比の差の絶対値を前記疲労度と相関する伝達比変化量として算出する伝達比変化量算出処理を備えた、
ことを特徴とする請求項2記載の疲労度測定装置。
The fatigue level calculation unit
A transmission ratio for calculating the transmission ratio for two time points exceeding the predetermined time and calculating an absolute value of a difference between the transmission ratios at the two time points as a transmission ratio change amount correlated with the fatigue level. With a change amount calculation process,
The fatigue degree measuring apparatus according to claim 2.
シートに設置され当該シートの垂直方向のシート加速度を測定してシート加速度信号Xnを記録するシート加速度記録工程と、
このシート加速度記録工程と並行して、乗員の頭部に設置され当該頭部の垂直方向の頭部加速度を測定して頭部加速度信号Ynを記録する頭部加速度記録工程と、
前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させるバンドパスフィルタ工程と、
予め定められた一定時間の前記シート加速度信号Xn及び頭部加速度信号Ynの統計値に基づいて、前記乗員の疲労度を算出する疲労度算出工程とを備え、
この疲労度算出工程が、
前記一定時間を超える時間間隔の2つの時点について、前記シート加速度信号Xnの統計値で前記頭部加速度信号Ynの統計値を除することで、前記シートの振動の前記頭部の振動への伝達比を算出する伝達比算出工程と、
前記2つの時点の各伝達比の差を前記疲労度と相関する伝達比変化量として算出する伝達比変化量算出工程と、を備えた、
ことを特徴とする疲労度測定方法。
A sheet acceleration recording step that is installed in the sheet and measures the sheet acceleration in the vertical direction of the sheet and records the sheet acceleration signal Xn;
In parallel with this seat acceleration recording step, a head acceleration recording step that is installed on the occupant's head and measures the head acceleration in the vertical direction of the head and records the head acceleration signal Yn,
A bandpass filter step of passing a predetermined frequency band for the seat acceleration signal Xn and the head acceleration signal Yn;
A fatigue degree calculating step of calculating the degree of fatigue of the occupant based on statistical values of the seat acceleration signal Xn and the head acceleration signal Yn for a predetermined period of time;
This fatigue level calculation process
By dividing the statistical value of the head acceleration signal Yn by the statistical value of the seat acceleration signal Xn at two time points that exceed the predetermined time, the vibration of the seat is transmitted to the vibration of the head. A transmission ratio calculation step for calculating the ratio;
A transmission ratio change amount calculating step of calculating a difference between the transmission ratios at the two time points as a transmission ratio change amount correlated with the degree of fatigue.
A method for measuring fatigue level.
請求項4記載の方法を演算装置を用いて実行するための疲労度測定用プログラム。   A fatigue degree measuring program for executing the method according to claim 4 using an arithmetic unit. 取り替え可能にシートを保持するシート保持部と、
前記シートを加振する加振部と、
シートに設置され当該シートの垂直方向のシート加速度を測定してシート加速度信号Xnを生成するシート加速度センサと、
乗員の頭部に設置され当該頭部の垂直方向の頭部加速度を測定して頭部加速度信号Ynを生成する頭部加速度センサと、
前記シート加速度信号Xn及び頭部加速度信号Ynについて予め定められた周波数帯を通過させるバンドパスフィルタと、
前記シート加速度信号Xn及び頭部加速度信号Ynを予め定められた一定時間蓄積すると共に、当該一定時間の前記シート加速度信号及び頭部加速度信号の統計値に基づいて、前記乗員の疲労度TRを算出する疲労度算出部と、を備えた、
ことを特徴とするシート評価装置。
A sheet holding section for holding the sheet in a replaceable manner;
A vibration unit for vibrating the sheet;
A seat acceleration sensor that is installed on the seat and measures a seat acceleration in the vertical direction of the seat to generate a seat acceleration signal Xn;
A head acceleration sensor that is installed on the head of the occupant and measures a head acceleration in the vertical direction of the head to generate a head acceleration signal Yn;
A bandpass filter that passes a predetermined frequency band for the seat acceleration signal Xn and the head acceleration signal Yn;
The seat acceleration signal Xn and the head acceleration signal Yn are accumulated for a predetermined period of time, and the occupant's fatigue level TR is calculated based on statistical values of the seat acceleration signal and the head acceleration signal for the predetermined period of time. A fatigue degree calculation unit,
A sheet evaluation apparatus.
JP2009169603A 2009-07-18 2009-07-18 Fatigue degree measuring device Pending JP2011019845A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009169603A JP2011019845A (en) 2009-07-18 2009-07-18 Fatigue degree measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009169603A JP2011019845A (en) 2009-07-18 2009-07-18 Fatigue degree measuring device

Publications (1)

Publication Number Publication Date
JP2011019845A true JP2011019845A (en) 2011-02-03

Family

ID=43630434

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009169603A Pending JP2011019845A (en) 2009-07-18 2009-07-18 Fatigue degree measuring device

Country Status (1)

Country Link
JP (1) JP2011019845A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101421122B1 (en) 2013-01-30 2014-07-22 국방과학연구소 Prediction Method of Muscle Fatigue and Prediction System of Muscle Fatigue
WO2017068880A1 (en) * 2015-10-19 2017-04-27 株式会社デンソー Motion estimation system, motion estimation method, and wearable device
KR20180041458A (en) * 2016-10-14 2018-04-24 광운대학교 산학협력단 Earable apparatus and method for measuring stress of human
JP2019517847A (en) * 2016-04-29 2019-06-27 フリーア、ロジック、インコーポレイテッドFreer Logic, Inc. Non-Contact Body- and Head-Based Monitoring of Brain Electrical Activity
CN111381229A (en) * 2018-12-28 2020-07-07 松下知识产权经营株式会社 Estimation method, estimation device, and recording medium

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101421122B1 (en) 2013-01-30 2014-07-22 국방과학연구소 Prediction Method of Muscle Fatigue and Prediction System of Muscle Fatigue
WO2017068880A1 (en) * 2015-10-19 2017-04-27 株式会社デンソー Motion estimation system, motion estimation method, and wearable device
JP2017077296A (en) * 2015-10-19 2017-04-27 株式会社デンソー Behavior estimation system, behavior estimation method, and wearable device
JP2019517847A (en) * 2016-04-29 2019-06-27 フリーア、ロジック、インコーポレイテッドFreer Logic, Inc. Non-Contact Body- and Head-Based Monitoring of Brain Electrical Activity
US11517240B2 (en) 2016-04-29 2022-12-06 Freer Logic, Inc. Non-contact body and head based monitoring of brain electrical activity
JP7208012B2 (en) 2016-04-29 2023-01-18 フリーア、ロジック、インコーポレイテッド Devices and methods for monitoring electrical activity produced by the brain and methods for monitoring physiological conditions in humans
KR20180041458A (en) * 2016-10-14 2018-04-24 광운대학교 산학협력단 Earable apparatus and method for measuring stress of human
CN111381229A (en) * 2018-12-28 2020-07-07 松下知识产权经营株式会社 Estimation method, estimation device, and recording medium
US11561279B2 (en) 2018-12-28 2023-01-24 Panasonic Intellectual Property Management Co., Ltd. Radar estimating method, device and medium to extract living body vector information
CN111381229B (en) * 2018-12-28 2024-05-14 松下知识产权经营株式会社 Estimation method, estimation device, and recording medium

Similar Documents

Publication Publication Date Title
EP2116185B1 (en) Biosignal analyzer and method of analyzing biosignal
Demić et al. Some aspects of the investigation of random vibration influence on ride comfort
Van Niekerk et al. The use of seat effective amplitude transmissibility (SEAT) values to predict dynamic seat comfort
EP2113198B1 (en) Load body state judgment device, vehicle seat, and computer program
JP2011019845A (en) Fatigue degree measuring device
Faris et al. Ride quality of passenger cars: an overview on the research trends
JP2009195384A (en) Fatigue analyzer and computer program
JP2009022610A (en) Degree of fatigue computing device and computer program
JP2009097973A (en) Vehicle vibration analysis method and device
JP3029103B1 (en) Ride comfort evaluation method and apparatus for vehicles
Lin et al. Biodynamic response of seated human body to roll vibration: Effect of armrest support
JP2002224072A (en) Muscular fatigue evaluating device, muscular fatigue evaluating method, computer program for executing this method and storage medium for storing this computer program
JP2007209453A (en) Muscle fatigue evaluation system
Patelli et al. Effects of seating on the discomfort caused by mechanical shocks: Measurement and prediction of SEAT values
Zheng et al. Evaluation of sternocleidomastoid muscle activity of a passenger in response to a car's lateral acceleration while slalom driving
Ibicek et al. Quantification of human discomfort in a vehicle using a four-post rig excitation
JP4216546B2 (en) Vehicle occupant fatigue evaluation method
Polanco et al. Effect of rider posture on bicycle comfort
Arun et al. Normalized frontal impact biofidelity kinematic corridors using post mortem human surrogates
Mashino et al. A study on ride comfort control method using heart rate variability in consideration of the physiological state
Taluja et al. Methodology to Quantitatively Evaluate the Secondary Ride Characteristics of a Vehicle
Atsumi et al. Evaluation of vehicle motion sickness due to vehicle vibration
JP2004286502A (en) Method and apparatus for measuring evaluation index of motion sickness in train
KR100784280B1 (en) Method for tire ride performance analysis
Morisaki et al. On the vibration transfer characteristics from the seat of the vehicle to the occupant