stanford-oval/wikipedia huggingface dataset:
キム・ポリング
主な戦績

70kg級の戦績
2006年 - ドイツカデ国際　2位
2007年 - ドイツカデ国際　優勝
2007年 - ヨーロッパカデ柔道選手権大会　2位
2008年 - ドイツジュニア国際　3位
2008年 - 世界ジュニア　7位
2009年 - ヨーロッパジュニア　2位
2010年 - ベルギー国際　ジュニアの部　優勝　シニアの部　2位
2010年 - ロシアジュニア国際　優勝
2010年 - ワールドカップ・タリン　優勝
2010年 - ドイツジュニア国際　優勝
2010年 - ヨーロッパジュニア　優勝
2010年 - ワールドカップ・バーミンガム　優勝
2010年 - 世界ジュニア　優勝
2011年 - ベルギー国際　2位
2011年 - ワールドカップ・オーバーヴァルト　3位
2011年 - ワールドカップ・マドリード　優勝
2011年 - ユニバーシアード　優勝
2011年 - グランプリ・青島　5位
2012年 - ワールドカップ・ソフィア　3位
2012年 - ワールドカップ・ワルシャワ　2位
2013年 - ヨーロッパオープン・ソフィア　優勝
2013年 - グランドスラム・パリ　優勝
2013年 - グランプリ・デュッセルドルフ　優勝
2013年 - ヨーロッパ選手権　優勝
2013年 - ワールドマスターズ 優勝
2013年 - 世界選手権　3位
2013年 - グランプリ・アブダビ　優勝
2013年 - グランドスラム・東京 2位
2014年 - ヨーロッパ選手権　優勝
2014年 - グランドスラム・バクー　優勝
2014年 - 世界選手権　5位
2015年 - グランドスラム・バクー　2位
2015年 - ワールドマスターズ　優勝
2015年 - ヨーロッパ競技大会　優勝
2015年 - 世界選手権　7位
2015年 - グランドスラム・パリ　3位
2015年 - グランドスラム・アブダビ　2位
2016年 - グランプリ・トビリシ　優勝
2016年 - ヨーロッパ選手権　5位
2016年 -　ワールドマスターズ　優勝
2017年 - グランドスラム・アブダビ　2位
2017年 - グランプリ・ハーグ 優勝
2018年 -　グランプリ・チュニス　優勝
2018年 -　グランドスラム・パリ　3位
2018年 - ヨーロッパ選手権　優勝
2019年 - グランプリ・ザグレブ　2位
2019年 -　グランドスラム・アブダビ　優勝
2019年 -　グランドスラム・大阪　2位
2019年 -　ワールドマスターズ　優勝
2021年 -　ワールドマスターズ　3位
2021年 -　グランドスラム・テルアビブ　5位
2021年 -　グランドスラム・アンタルヤ　優勝
2021年 -　ヨーロッパ選手権　5位
2022年 -　ワールドマスターズ　5位
2023年 -　グランドスラム・パリ　3位
2023年 -　世界団体　3位
2023年 -　グランドスラム・アスタナ　3位
2023年 -　グランドスラム・東京　7位
2024年 - グランドスラム・アスタナ　3位
2024年 - 世界選手権　7位
2024年 -　世界団体　3位
(出典、JudoInside.com)
線形多自由度系の振動

振動工学における線形多自由度系の振動（せんけいたじゆうどけいのしんどう）は、線形な特性を持ち、さらに2以上の自由度を持つ系で起きる振動である。運動方程式は一般的に連立2階常微分方程式となり、行列およびベクトルで表現される。
線形多自由度系の振動では、固有モードという多自由度系特有の概念が現れ、自由度の数だけ固有モードと固有振動数の組が存在する。固有モードの直交性によって、減衰の無い系であれば固有モードごとの1自由度系の問題に帰着でき、振動解析を容易化できるのが特徴である。この手法を利用した振動解析手法はモード解析と呼ばれる。減衰のある系でも、比例粘性減衰という仮定を導入することによって、同様なことが可能となる。
モード解析手法は、振動実験結果から振動特性を同定するのにも使われる。有限要素法による連続体の振動計算においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析を容易にする。
背景

振動の問題を分類すると、あるいは振動現象を模した動力学モデルを分類すると、いくつかの視点が存在する。
物体の運動を表すのに必要な座標あるいは変位の数を、自由度という。自由度が1の系を1自由度系という。ある1つの質点がばねとダッシュポットを介して地面に固定され、質点が上下方向のみ動く例を考える。三角関数で表現される励振力がこの質点に加わるとき、この質点の運動方程式は次のような形で与えられる。
()
ここで、m は質量、c は粘性減衰係数、k はばね定数、f0 は外力の振幅、ω は外力の角振動数、t は時間である。x は静的釣り合い位置からの変位で、上部の "˙" は時間微分 (d/dt) を示す。式1.1は、振動問題を考える際の最も基礎的な式となる。1自由度系の振動の問題は、振動現象を理解する上で不可欠な様々な概念を内包しており、1自由度系のモデルが振動問題を扱う基本モデルといえる。
一方、自由度が2以上の系は多自由度系と呼ばれる。振動問題をモデル化する際には1自由度や2自由度にモデル化すれば解析が容易であり、モデル化時の自由度はできるだけ最小限とするのが大切である。しかし、現実の構造物は複雑で、すべてを1自由度系として扱うことができない。実際の機械や建物などでは、1自由度系のモデルではその特性を説明しきれず、多自由度系としての取り扱いが必要になることも多い。
また、系の構成要素が全て線形であるとき、その系を線形系という。系が線形であれば、出力と入力の間には重ね合わせの原理が成り立ち、出力は入力に単純に比例する。振動系が線形であれば、式1.1のように慣性力、減衰力、復元力が、それぞれ加速度、速度、変位に比例する。比例係数が時間変化する場合を除き、線形系の振動問題については理論的手法がほぼ確立できている。
一方で、線形ではない系を非線形系という。摩擦、がたつき、大変位時の材料特性などが非線形の要因となる。厳密に考えると、実際の機械や構造物はなんらかの非線形特性を持ち、実際の振動現象のほとんどは非線形系といえる。非線形系の厳密解を得ることはほとんどできず、一般的には非線形系の問題に対しては近似解法や数値計算に頼らざるを得ない。しかし、線形系の振動として取り扱うことによって、十分に多くの問題を解決することもできる。非線形系であっても、安定な平衡状態周りの微小変位運動であれば線形の理論が当てはまる。特に線形系であれば、多自由度系であっても後述のように1自由度系の問題に帰着できるという顕著な特性がある。
運動方程式 > 一般の多自由度系

多自由度系では多数の変数や係数を扱うため、行列とベクトルを使って運動方程式を記す。自由度が n の線形多自由度系の一般的な運動方程式は、次のように書き表される。
()
ここで、ẍ, ẋ, x, f は n 次元縦ベクトル、M, C, K は n 次正方行列で、以下のように表される。
()
()
()
x は変位ベクトルと呼ばれる。その成分の x は、静的な釣り合いの位置を原点とした各自由度の変位あるいは一般化座標を表している。ẋ は速度ベクトル、ẍ は加速度ベクトルと呼ばれる。成分の ẋ は x の時間1階微分すなわち各自由度の速度を意味し、ẍ は x の時間2階微分すなわち各自由度の加速度を意味する。
f は外力ベクトルと呼ばれる。成分の f は各一般化座標に対応して作用する外力で、一般的には時間の関数である。
M は質量行列や慣性行列、質量マトリックスや慣性マトリックスと呼ばれる。K は剛性行列や剛性マトリックスと呼ばれる。C は減衰行列や減衰マトリックスと呼ばれる。各成分の m, c, k は質量、粘性減衰係数、剛性（ばね定数）を表す。ここでの m は慣性係数ともいい、いわゆる質量だけでなく慣性モーメントなども含む。また、k は復元係数ともいい、通常のばね定数の他に回転ばね定数なども含む。
線形多自由度系の一般基礎式2.1は、線形1自由度系の基礎式1.1と形式は同じで、変位と外力がベクトルに、質量とばね定数と粘性減衰係数が行列に置き換わった式となる。ただし、対象の系の規模が大きくなり、複雑なものとなると、ニュートンの運動法則から運動方程式を導出するのは容易ではない。実際に多自由度系の運動方程式を立てる際は、エネルギーのスカラー量から形式的に運動方程式を導けるラグランジュの運動方程式が便利で、間違いを犯しにくく、多用されている。
運動方程式 > 2自由度系の例

建築物の例では、複数の階を持つ多層構造物が多自由度系の問題となる。柱と床から構成されるラーメン構造の建物が振動する場合を考える。柱は床に比べて柔らかいので、構造物が揺れるとき各階の床は水平方向に揺れ、柱は水平方向のばねとして働くと見なせる。これは、水平方向のみに動く2つの質点をばねで連結した2自由度系モデルと等価となる。2層構造物における、1層目の床質量を m1、2層目の床質量を m2、1層目の水平変位を x1、2層目の水平変位を x2、基礎と1層目の間の等価ばね定数を k1、1層目と2層目の間の等価ばね定数を k2 とすると、このモデルの運動方程式は
()
となる。行列表示すれば、
()
となり、この式では、 {\begin{pmatrix}m_{1}&0\\0&m_{2}\end{pmatrix}}が質量行列、 {\begin{pmatrix}k_{1}+k_{2}&-k_{2}\\-k_{2}&k_{2}\end{pmatrix}}が剛性行列である。
並進運動と回転運動が組み合わさった2自由度系の例として、自動車の上下振動とピッチング振動の簡易モデルがある。自動車を側面から見て、上下方向変位とピッチング回転だけの動きを考える。前位側のタイヤとサスペンションを1つのばねと見なして、後位側のタイヤとサスペンションも同様に1つのばねと見なし、剛体を前位と後位を2つのばねが支えているモデルを考える。前位側と後位側のばねのばね定数をそれぞれ k1、k2 とする。剛体の重心位置から前位側のばねまでの距離を l1、前位側のばねまでの距離を l2 とする。剛体の質量を m、重心位置周りのピッチング方向の慣性モーメントを IG とする。この剛体の重心の上下運動 x とその周りのピッチング運動 θ の運動方程式は次のようになる。
()
運動方程式 > 連成項

自動車の上下振動とピッチング振動の簡易モデルである式2.7において、もし k1l1 = k2l2 であれば、x と θ は互いに独立した運動方程式となる。この条件では運動方程式は
()
となり、2つの式はそれぞれ単独で解くことができる。すなわち、x の振動と θ の振動が、互いに干渉すること無しに独立して起こる状態になっている。このような状態を非連成という。運動方程式が単独では解けない連立方程式になっているとき、方程式は連成しているといい、式2.7の (k1l1 − k2l2)x と (k1l1 − k2l2)θ は連成項と呼ばれる。k1l1 = k2l2 のときに、これら連成項が 0 になり、x と θ が独立した運動になる。連成状態にある系を連成系、逆に各自由度の運動が完全に独立している系を非連成系と呼ぶ。
一般に、質量行列 M（式2.2）、減衰行列 C（式2.3）、剛性行列 K（式2.4）における非対角成分の存在は、振動が連成していることを示している。もし系が非連成系であれば、これら全ての行列は、非対角成分が全て 0 の対角行列となる。特に、質量行列が非対角成分によって振動が連成していることを動連成や動的連成といい、質量行列の非対角成分を動連成項という。また、剛性行列が非対角行列であれば静連成や静的連成といい、剛性行列の非対角成分を静連成項という。多自由度系の問題は概して複雑で、個々のパラメータが結果に与える影響の見通しを立てることが難しい。多自由度系では一般的に連成が存在し、これが多自由度系の解析を難しいものにしている。動力学的な設計を行う上では、系を非連成化して影響をわかりやすくすることが有効となる。
不減衰自由振動 > 振動数方程式

継続的な外力が作用せず、外乱だけが与えられて起こる振動を自由振動という。自由振動の振動は、系自体が持つ特性によって決まり、系の動特性を知る上で重要な振動の形態である。自由振動は減衰要素が存在しない場合と存在する場合に分かれる。減衰要素と励振力が存在しない場合の式2.1は
()
となる。ここで、0 は下記のようなゼロベクトルである。
さらに、M と K は正定値であると仮定する。すなわち、M と K は成分が全て実数の対称行列で、それらの二次形式は常に正となる。実際に、質量行列と剛性行列が正定値行列であることは線形多自由度系の一般的な特徴であり、多くの振動系でこの仮定は満たされる。この条件を満たす M と K を前提にすれば、振動系の具体的な構成に依存しない一般性の高い議論を展開できる。
式3.1は定数係数の線形常微分方程式であるため、その解法に従って解を
()
という形式で仮定できる。ここで e はネイピア数、λ は未知定数、 u は n 個の未知定数 u から成る下記のような縦ベクトルである。
式3.1に対して仮定として与えられる解には、他に、最初から単振動を仮定して複素指数関数や三角関数の形式もある。
式3.2を式3.1に代入して整理すると下記のような式になる。
()
この式が恒等的に成り立つには下記の条件が満たされる必要がある。
()
数学的には、この形式の方程式は一般化固有値問題として知られ、λ2 は固有値、u は固有ベクトルと呼ばれる。u = 0 であれば式3.4の条件は満たされるが、これは最初の釣り合いの位置からそのまま静止しているだけ状態を意味する解である。よって、ここでは u = 0 は興味の対象外で、それ以外の式3.4を満たす λ と u について知りたい。u ≠ 0 で、なおかつ式3.4が満たされる条件は、U の係数行列 λ2 M + K の逆行列が存在しないことである。したがって、係数行列の行列式が 0 であればこの条件が満たされる。したがって、λ が、
()
を満たすとき、式3.4が u = 0 以外の解を持つ。ここで、det( ) は行列式を表す。式3.5の行列式を展開すると、 λ2 についての n 次多項式になる。したがって、原理的には λ2 の値を求めることができる。ただし、この多項式を解析的に解くことができるのはせいぜい2自由度あるいは3自由度までで、それ以上の自由度の系になると数値解析で固有値を計算する。ヤコビ法などを使って数値計算するときは、コレスキー分解を使い、一般化固有値問題形式の3.4を標準的な固有値問題の形へ変換する。
上記のとおり、M と K は正定値であると仮定した。このとき、式3.5の解は全て負の実数となる（解を複素指数関数と仮定した場合は全て正の実数）。したがって、n 個の λ2 の解を
()
とおくことができ、平方根を取って 2n 個のλ の値
()
が得られる。ここで j は虚数単位である。n　個の ω は固有角振動数または固有円振動数と呼ばれ、値が小さいものから順に1次、2次、…、n 次の固有角振動数と呼ぶ。特に、最も値が小さい1次の固有角振動数は基本振動数と呼ばれる。式3.5は角振動数を求める式であるため振動数方程式と呼ばれる。あるいは、式3.5は M や K といった系の特性によって構成される式であることから特性方程式とも呼ぶ。
不減衰自由振動 > 固有モード

得られた固有値あるいは固有角振動数の値を式3.4に代入すれば、0 以外の u の解が得られる。ここで得られる u の各成分は一意な値を持たず、定まるのは各成分の互いの比 u1 : u2 : … : un だけである。一つの固有値あるいは固有角振動数に対応して一つの u が定まり、u は n 個存在する。このような固有値と固有ベクトルの組は固有ペアと呼ばれる。r 次の固有角振動数 ωr に対応する u を ur を表現すれば、 ur は
()
というベクトルである。ただし、上記のとおり、各成分 u1r, u2r, …, unr は互いの比の値を表している。
式3.8で表されるベクトルが式3.4における固有ベクトルであり、振動工学では固有モード、振動モード、固有振動モード、基準振動モード、モードベクトルなどと呼ぶ。n 自由度系には n 個の固有角振動数があり、固有角振動数それぞれに対応する形で n 個の固有モードが存在している。各自由度の振幅比を決める固有モードは、固有角振動数が「振動の速さ」を表しているのに対して、「振動の形」を表していると言える。自由度の数の分だけ固有角振動数が存在し、それら固有角振動数に対応して固有モードが存在していることが、線形多自由度系の特有な性質といえる。
各々の固有モードを、対応する固有角振動数が小さい順に1次固有モード、2次固有モード、…、r 次固有モード、…、n 次固有モードと呼ぶ。特に、基本振動数（1次の固有角振動数）に対応する固有モードは基本モードと呼ばれる。固有モード  ur を次数が低い順に並べて作る、下記のような行列をモード行列やモードマトリックスという。
()
不減衰自由振動 > 固有モードの直交性と正規化

多自由度系の固有モードには直交性という重要な性質がある。r 次の固有モード ur と s 次の固有モード us について考える。ここで、r と s は任意だが、r ≠ s である。M と K は上記のとおり対称行列とする。このとき、  ur, us, M, K には次のような関係がある。
()
()
ここで、ur⊤ は ur の転置行列を表している。これらの式は ur と Mus の内積および ur と Mus の内積が零であることを示しており、M と K に関して ur と us が直交していることを意味している。これが質量行列および剛性行列を介した固有モードの直交性である。別の言い回しでは ur と us が一般直交性を有しているともいう。
固有モードの直交性の物理的な意味は、次数の異なる固有モードの振動の間で力学的エネルギーの移動が全く起きないことを表している。ある固有モードの振動において、運動エネルギーと復元力によるポテンシャルエネルギーは時間的に変化しているが、それらの和の力学的エネルギーは時間に依らず一定に保たれている。そのため、自由振動中にある固有モードの振動が他の固有モードに移り変わったり、他の固有モードが新たに誘起されたりすることはない。
一方で、r = s、すなわち同じ次数同士の固有モードの場合は、上記の式でも左辺は 0 とはならず、ある定数となる。これらの定数を Mr と Kr と表せば、
()
()
と表される。固有モード ur は成分間の比だけを持ち、定まった値を持たないので、Mr と Kr の値もこの段階では定まらない。固有モードの絶対値が決まった後に、Mr と Kr の値も定まる。
定数 Mr と Kr は正の値であり、それぞれを（r 次の）モード質量およびモード剛性という。モード行列を使って式3.12と式3.13を表せば、
()
()
となる。モード質量とモード剛性は、同じ次数の固有角振動数 ωr と下記のような関係がある。
()
固有モードは比の関係を定めるベクトルであった。絶対的な大きさを持たない固有モードに大きさを定める方法として、次のような方法がある。
下記のように、式3.12の右辺 Mr の値が 1 となるように定める。ここで、ūr は、等式を満たすに大きさが決定された固有モードを意味している。
 {\boldsymbol {\overline {u}}}_{r}^{\top }{\boldsymbol {M}}{\boldsymbol {\overline {u}}}_{r}=1(3.17)
下記のように、同じ次数の固有モード同士の内積が 1 となるように定める。
 {\boldsymbol {\overline {u}}}_{r}^{\top }{\boldsymbol {\overline {u}}}_{r}=1(3.18)
固有モードの成分の内、絶対値が最大のものを 1 とおいて定める。
1番目の自由度に対応する固有モード成分 (u1r) を 1 とおいて定める。
以上のように固有モードの大きさを一意に定めることを正規化と呼び、正規化された固有モード ūr を正規固有モードや正規化モードと呼ぶ。1番目の手法による正規固有モードは、特に質量正規固有モードやM-正規固有モードという名で呼ばれる。式3.17を満たすようにしておくと利点が多く、質量正規固有モードはよく使われる。正規固有モード ūr を並べて作る下記のような行列を、正規モード行列という。
()
質量正規固有モードを採用すると、式3.17に対して剛性行列は固有角振動数と
()
という関係を持つ。
不減衰自由振動 > モード座標と一般解

固有モード（固有ベクトル）は互いに直交なので、これらを使った線形結合で変位ベクトル x を表すことができる。つまり、固有モード U を用いて、x を
()
と表すことができる。ここで、q は時間の未知関数で、モード座標、基準座標、正規座標、規準座標、主座標などと呼ばれる。特に正規化した固有モードを使っているモード座標を指して、正規座標と呼ぶこともある。q は qr を 1 次から n 次まで並べた次のような縦ベクトルである。
モード座標 q に対して、元の座標 x を物理座標と呼ぶ。3.21は、正規固有モードを介して物理座標をモード座標に変換していることを意味する。モード座標 qr は、x に対して r 次固有モードūr が寄与する度合いを表しているとも言える。式3.21は展開定理とも呼ばれる。
式3.21を運動方程式3.1へ代入して、左から転置したモード行列 U⊤ を掛けると次のようになる。
()
固有モードの直交性（式3.14と式3.15）を上式に当てはめると、
()
というような q に関する n 個の運動方程式が得られる。それぞれの式の両辺を Mr で割ると、下記のような形になる。
()
式3.24は非連成化されており、各式はそれぞれ独立している。そのため、一つ一つの式は1自由度系の不減衰自由振動と同じであるから、q の各解は以下のようになる。
()
ここで、cr と θr (r = 1, 2, …, n) は初期条件で決まる定数で、cr が振幅、θr が初期位相を意味する。x の一般解は、式3.25を式3.21に代入して
()
と得られる。すなわち、n 自由度系における各自由度の自由振動は、n 個の調和振動の重ね合わせ（和）となっており、それら調和振動のそれぞれの振動数は 1 次から n 次までの固有角振動数となっている。固有モードは、対応する固有角振動数を持つ調和振動成分の各自由度間の振幅比を定めている。位相角を陽に表さずに、余弦関数と正弦関数の和
()
や、指数の正負が異なる複素指数関数の和
()
といった形で一般解を表すこともできる。ar, br, dr, gr も初期条件で決まる定数で、式3.27と式3.27は式3.26へ式変形可能な同値な式である。
以上のような、物理座標をモード座標へ変換し、固有モードの直交性を利用して多自由度系の問題を1自由度系の重ね合わせの問題に帰着させ、解析を行う手法をモード解析（特にモード重畳法）という。後述のように、モード解析は特に有限要素法へ適用することで有効性を発揮する。実物の振動特性を求める上でもモード解析の理論を適用することで各特性を同定でき、この手法は実験モード解析と呼ばれる。
減衰自由振動 > 一般の粘性減衰

減衰行列 C が存在する線形多自由度系の振動について考える。外力が無い場合の n 自由度系の運動方程式は以下のようになる。
()
基礎の上でばねと減衰器が一組となって質点と連結し、直列に連なった3自由度減衰系の典型的な例では、M, C, K は次のような行列となる。
()
()
()
式4.1の解を式3.2と同じように仮定し、式3.2を式4.1に代入して整理すると、
()
となり、u = 0 以外の解を持つという条件から
()
という特性方程式が得られる。式4.6を展開すると、λ に関する 2n 次多項式となる。この多項式を解いて得た解 λ1, λ2, …,  λ2n を式4.5に代入し、u を定め、u に適当な正規化を行う。正規化された各ベクトルを ū1, ū2, …, ū2n と表すと、式4.1の一般解は
()
となる。ここで、a1, a2, …, a2n は初期条件で決まる定数である。2n 個の λ は、解が減衰振動であれば n 組の互いに共役な複素数になる。このとき、固有モードも複素数となり、複素固有モードと呼ばれる。
減衰自由振動 > 比例粘性減衰

式4.1のように一般的な減衰行列 C が運動方程式に存在する場合、正規モード行列によって M と K は対角化できるが、C も同時に対角化することはできない。そのため、不減衰自由振動で可能だったモード座標に変換しての非連成化が不可能となり、モード解析の利点を活かすことができなくなる。そこで、C が下記のような比例粘性減衰として与えられると仮定し、実際の振動解析が行われることも多い。
()
ここで、α と β は定数で、比例減衰定数と呼ばれる。比例粘性減衰を有する系を比例粘性減衰系などと呼ぶ。比例粘性減衰が成り立つと仮定すれば減衰行列を対角化できる。
式4.8の形で与えられる比例粘性減衰は、特にレイリー減衰やレイリー型減衰と呼ばれる。αM のみで仮定されるものは質量比例型減衰、βK のみで仮定されるものは剛性比例型減衰などと呼ぶ。
実際の減衰が比例粘性減衰になっていることはまれであり、比例粘性減衰はあくまでも近似的なものである。しかし、比例粘性減衰の仮定を導入することで、不減衰系と同じ取り扱いが可能となり、モード解析の手法が適用可能になる。もし減衰が全体に分布しているような構造であれば、適当な比例減衰定数を設定すれば、実際の現象を実用問題ないレベルで再現できるという一定の妥当性もある。減衰は摩擦・材料減衰・流体粘性など様々な要因で起こるため、そもそも減衰の適切な定式化自体が難しいといった事情もある。式2.3のように速度の比例定数として与えられる一般の粘性減衰も、多種多様な発生機構によって減衰が起きるという実情に起因して厳密には成立しない。
式4.8を式4.1に代入した場合、
()
という運動方程式になる。解を式3.2のように仮定して上式に代入し、整理すると、
()
となる。さらに、
()
とおけば、式4.10は
()
となる。式4.12は、（複素指数関数や三角関数で解を仮定したときの）不減衰振動の固有値問題における λ を γ に置き換えただけの式になる。したがって、比例粘性減衰を仮定した減衰振動の固有モードは、同一の質量行列と剛性行列を有する不減衰振動の固有モードと同じである。
一方、比例粘性減衰を仮定した減衰振動の固有角振動数は、同一の質量行列と剛性行列の不減衰振動の固有角振動数よりも小さくなる。モード行列 U で減衰行列を対角化すると
()
となる。この行列の成分 αMr + βKr (r = 1, 2, …, n) をモード減衰やモード減衰係数と呼び、 Cr などで表す。ここで、Mr はモード質量、Kr はモード剛性である。さらに、
()
という量を導入して、これでモード減衰を割った量
()
をモード減衰比と呼ぶ。多自由度系の減衰系では固有モードごとに減衰の効果が異なっており、モード減衰比が固有モードごとの減衰の効果の程度を表している。ζr ≥ 1 ならば過減衰の状態であり、その固有モードの振動は起こらない。ζr < 1 ならば減衰振動となり、その固有角振動数は
()
で与えられる。ωd,r を減衰固有角振動数と呼ぶ。以上のように、線形1自由度系の減衰振動の考え方が固有モードごとの振動にも当てはまる。正規座標へ変換を行い、モード減衰比と固有角振動数を用いて運動方程式を表すと、下記のように表現できる。
()
ζr < 1 であれば、qr の一般解は積分定数を ar と br として、
()
となる。物理座標 x へ逆変換すれば、比例粘性減衰系の自由振動の一般解は下記のようになる。
()
強制振動 > 不減衰系

任意の励振力を受ける減衰の無い系の運動方程式は以下のように表される。
()
例えば、励振力が角振動数 Ω の余弦関数で与えられるとすれば、
()
となる。ここで、fa は、以下のような各自由度に対する励振力の振幅値の縦ベクトルである。
式5.2に対して特解を x = u cosΩt と仮定し、式5.2に代入すれば
()
となるので、u および x は、
()
()
となる。ここで "−1" は逆行列を意味する。したがって、逆行列 (M − Ω2K)−1 を計算すれば式5.5から x の値が分かる。しかし、この逆行列を解析的に解くことは困難で、数値計算を行うにしても自由度の数が増えると膨大な計算量になる。もし励振力の角振動数 Ω を変えると、そのたびに逆行列を計算する必要がある。そのため、実際に x の解を求めるために行われるのは、下記のようなモード解析による手法である。
解を得るために、x を複素数に拡張し、励振力 f を複素指数関数 fae jΩt の形で与えるとする。この場合、計算して解が得られた後に実部あるいは虚部を取ることで、実際の解が得られる。運動方程式5.1の右辺を 0 としたときのモード行列 U が、事前に求められているとする。モード座標への変換式3.21を運動方程式5.1へ適用して、左から U⊤ を掛け、式3.14と式3.15の対角化を適用する。すると、励振力を受ける不減衰系の運動方程式は
()
という独立・非連成の n 個の運動方程式に帰着する。ただし、右辺の Fr は次のような値である。
()
式5.6は線形1自由度系と同じなので、強制振動を表す特解は、
()
となる。式3.21を使ってモード座標を物理座標へ逆変換し、固有角振動数を使って整理すれば強制振動の解は次のようになる。
()
式5.9から分かることは、励振力の角振動数 Ω が系の固有角振動数 ω1, ω2, …, ωn のどれかに近いとき、係数の極限が 0 となってその振動成分が極めて大きくなり、共振が起こるという点である。Ω が固有角振動数のいずれかに一致するとき、x の振幅は無限大に発散する。
強制振動 > 比例粘性減衰系

系が比例粘性減衰系であれば、励振力を受ける場合でもモード座標変換によって独立した1自由度系に分解でき、大自由度の問題も扱いが容易になる。一般の運動方程式2.1に、レイリー減衰の式4.8およびモード座標への変換式3.21を代入する。左から転置した質量正規固有モードの正規モード行列 Ū⊤ を掛けて対角化する。さらに、2ζrωr = α + βωr2 という関係を用いれば、運動方程式は下記のような独立・非連成の n 個の方程式に帰着する。
()
励振力 f が fa cosΩt で与えられるとする。線形1自由度系と同様の解法によって、正規座標 qr の強制振動の解が下記のように求まる。
()
()
解を複素数で表現した場合は、
()
となり、物理座標の解は
()
となる。Ω/ωr は強制振動比などと呼ばれる。通常の固有モードの場合は次式となる。
()
強制振動 > 周波数応答関数

強制振動の特解から、周波数応答関数あるいは伝達関数が求められる。周波数応答関数は、振動系への入力が調和振動の形で与えられるときの、出力と入力の振幅比および位相差を周波数の関数として表す。比例減衰系において、角振動数 Ω の調和外力が q 番目の自由度のみに加わるとする。このときの p 番目の自由度の振動応答は
()
と表される。この場合の周波数応答関数 G(Ω) は[変位]/[力]の単位を持ち、コンプライアンスとも呼ばれる。q 番目加振・p 番目応答の周波数応答関数と p 番目加振・q 番目応答の周波数応答関数は、相反定理によって互いに等しく、Gpq(Ω) = Gqp(Ω) の関係である。不減衰系の場合は、周波数応答関数は次式となる。
()
これらの周波数応答関数によって、系の周波数特性が把握できる。線形多自由度系の場合は、固有角振動数と固有モードを求めれば、調和外力による強制振動に対する周波数特性も同時に理解することができる。横軸に振動数や強制振動比を取り、縦軸に調和外力に対する振幅比や位相差を図示したものを応答曲線や共振曲線と呼ぶ。縦軸をデシベルにして図示したものは特にボード線図と呼ばれる。系の振動特性を振動実験から同定する実験モード解析では、実験値に上記の周波数応答関数でカーブフィッティングして系の各種パラメータを推定する。
上述のように、励振力の振動数が系の固有角振動数に一致すると、振動は極大化する。一方、多自由度系では共振して振幅が発散する現象だけでなく、励振力が特定の角振動数のときに周波数応答関数が極小になる現象もある。このような現象を反共振と呼ぶ。ただし、共振が系全体で起こるのに対して、反共振は一組の加振点・応答点ごとにしか起こらない。周波数応答関数の絶対値を縦軸にした応答曲線上では、共振点は曲線の鋭い山のように現れ、反共振は曲線の鋭い谷のように現れる。反共振は、代表的な制振器である動吸振器の原理として活用される。振動を抑えたい対象にばね・ダンパを介して付加質量を取り付けることによって振動を抑制でき、地震や強風に対する建築構造物の防振や回転機械の防振などに使われる。
強制振動 > 刺激係数

刺激係数は、地震のように構造物の基礎が揺れている場合に、その励振が各モードに対してどのぐらい強く寄与するのかを表す。構造物を多質点系でモデル化して、基礎の変位を y として、各質点の基礎からの相対変位を x とする。各質点は　ẍi + ÿ の加速度を受けるため、運動方程式は
()
となる。ここで 1 は成分が全て 1 の縦ベクトルである。式5.18において、C をレイリー減衰で表し、さらに x をモード座標に置き換え、左から転置したモード行列を掛けて整理すれば、下記のような各モード座標についての式になる。
()
さらに、上式の ÿ の係数を次のように変形する。
()
この βr を r 次の刺激係数と呼び、r 次モードの運動方程式は次のようになる。
()
つまり、βr は励振加速度が各モードに対してどのぐらい寄与しているか、あるいは各モードは基礎励振に対してどのぐらい影響を受けやすいかを表している。さらに、r 次の刺激係数と固有モードの積を r 次の刺激関数と呼び、刺激関数の総和は次のように 1 に等しい。
()
すなわち、対称の振動系に 1 という外力ベクトルが加わったときに各モードが受ける度合いを、刺激関数は表している。
有限要素法による連続体の振動への応用

多自由度系では、質点または剛体から成る系を想定し、剛性、減衰、慣性などの振動特性も局所的に集中して系内で点在しているというモデルを考えていた。一方で、実際の物体は、物体自体が変形する。実際の物体は連続体としての性質を有しており、質量、剛性などは連続的に系内に分布しているモデルとなる。実際の問題では、多自由度系に近似して取り扱っても十分な場合も多いが、振動時の機械・構造物の各部の変形や応力といったものを知るには連続体として取り扱う必要ある。
連続体の振動の運動方程式は時間と空間に関する偏微分方程式で記述され、厳密解が得られることは限られる。実際の複雑な形状の構造物で連続体の振動を扱うには、実用的には有限要素法という手法が用いられる。有限要素法では、対象の連続体を小さな有限要素に分割し、連続体を多自由度系に置き換えて解を計算する。有限要素法においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析が精度良くかつ容易にできるようになる。
有限要素法による振動解析では、立てられた振動数方程式の数値計算を行い、まず固有振動数と固有モードを得る。次いで、固有モードの直交性を利用して周波数応答関数を得て、応答解析を行う。もし、定式化された運動方程式が1万自由度だとしたら、解くべき方程式は1万元の2次連立微分方程式となり、コンピュータを用いても計算に長時間を要する。しかし、モード解析手法を用いれば、1自由度系の解を1万回解いて、重ね合わせるだけで解が得られる。さらに、対象が大規模自由度になったとしても、自由度の分だけ現れるモードを全て計算する必要性もない。実用的に興味のある外力振動数を含む次数までモードの重ね合わせでも、十分な精度の振動応答解析が可能となる。上記の1万自由度の例えで言えば、1自由度系の解を1万回解く必要もなく、もっと少ない回数の計算で事足りるようになる。これらの長所によって、モード解析手法は有限要素法による振動解析で絶大な威力を発揮し、数十万規模の自由度を扱うような有限要素法計算であっても、モード解析手法の適用によって特段の支障なく計算が可能となる。
係数行列が非対称行列の場合

質量行列 M（式2.2）、減衰行列 C（式2.3）、あるいは剛性行列 K （式2.4）が正定値の条件を満たさない場合、すなわち実対称行列ではなく、非対称行列であるとき、その系では不安定振動が起こることがある。このような条件では、式4.6で表される特性方程式の固有値 λ に、実部が正の固有値が含まれることがありえる。固有値に実部を正とする複素数が含まれるとき、時間とともに振幅が大きくなっていく振動が起こる。このようなメカニズムは、自励振動が起こりえる系で平衡点から振動が成長するか否かを考察する上で基礎となる。自励振動は1自由度系でも起きる現象だが、係数行列が非対称であることによって引きこされる種類の自励振動は、多自由度系特有のものである。
例として、次のような2自由度不減衰系を考える。
()
ただし、k12 ≠ k21 で、剛性行列は非対称行列である。さらに、k12 と k21 のどちらかが正でどちらかが負であるような異符号の関係にあるとき、固有値は
()
という形の複素数となる。λr とλi は
()
()
で与えられ、ここで、ωtr = k11/m1 + k22/m2, ωdiff = k11/m1 − k22/m2, ωsk = (k12/m1)(k21/m2) である。λr は発散率と呼ばれ、自励振動の強さを表す。
このような係数行列の非対称性によって起きる自励振動の事例は機械振動の中で多く見られ、クーロン摩擦による摩擦振動や滑り軸受で起こるオイルホイップなどがある。
参照文献

末岡 淳男・金光 陽一・近藤 孝広、2000、『機械振動学』初版、朝倉書店〈基礎機械工学シリーズ 6〉 ISBN 4-254-23706-5
横山 隆・日野 順市・芳村 敏夫、2015、『基礎振動工学』第2版、共立出版 ISBN 978-4-320-08211-3
長松 昌男・長松 昭男、2018、『実用モード解析入門』初版、コロナ社 ISBN 978-4-339-08227-2
下郷 太郎・田島 清灝、2002、『振動学』初版、コロナ社〈機械系 大学講義シリーズ 11〉 ISBN 4-339-04045-2
入江 敏博・小林 幸徳、2006、『機械振動学通論』第3版、朝倉書店 ISBN 978-4-254-23116-8
吉川 孝雄・松井 剛一・石井 徳章、1987、『機械の力学』初版、コロナ社 ISBN 978-4-339-04273-3
安田 仁彦、2012、『振動工学 ―基礎編』改訂版、コロナ社 ISBN 978-4-339-04624-3
藤田 勝久、2016、『振動工学 ―振動の基礎から実用解析入門まで』新装版、森北出版 ISBN 978-4-627-66542-2
日本機械学会（編）、2004、『機械力学』第1版、丸善〈機械工学便覧 基礎編 α2〉 ISBN 4-88898-116-7
砂子田 勝昭・伊藤 智博・鄭 萬溶・平元 和彦、2012、『わかりやすい振動工学』初版、共立出版 ISBN 978-4-320-08185-7
宮本 裕司・永野 正行・藤谷 秀雄・吉村 智昭、宮本 裕司（編）、2014、『建築振動を学ぶ ―地震から免震・制震まで』初版、理工図書 ISBN 978-4-8446-0824-0
平井 一男・水田 洋司、2018、『耐震工学入門』第3版・補訂版、森北出版 ISBN 978-4-627-46454-4
背戸 一登・丸山 晃市、2002、『振動工学 ―解析から設計まで』第1版、森北出版 ISBN 4-627-66451-6
質点がばねとダッシュポットを介して基礎に固定されたモデル
2層構造物の揺れのモデリング
自動車の上下振動とピッチング振動の簡易モデル
式2.5で表現される2層構造物振動モデルの固有モードの例。k1 = k2, m1 = m2 の場合を示しており、各固有モードの成分の内の一番大きな絶対値が1となるようにして値を定めている。
式2.7で表現される自動車の上下・ピッチング振動モデルの振動の例。図はピッチング θ の振動を示している。パラメータは m = 1600 kg, IG = 2500 kg-m2, k1 = 3500 N/m, k2 = 4100 N/m, l1 = 1.4 m, l2 = 1.6 m で、初期条件は上下変位速度が ẋ = 1 m/s で他は全て 0 という条件の例。
基礎の上でばねと減衰器が一組となって質点と連結し、直列に連なった3自由度減衰系の例
式4.2,4.3,4.4によって与えられた3自由度減衰系の周波数応答関数の応答曲線の例。図は質点 m3 が加振されたときの質量 m2 の応答の例。山になっている点が共振点で、谷になっている点が反共振点。
静脈内区域麻酔

| 静脈内区域麻酔 |
| --- | --- |
 治療法
 旧式のターニケットを使用した静脈内区域麻酔（IVRA）。最新の電子機器では、IVRAに特化した安全機能が充実している。
MeSH:  D000771
 [ウィキデータを編集]
テンプレートを表示
静脈内区域麻酔（じょうみゃくないくいきますい、英: intravenous regional anesthesia: IVRA）またはBier（ビーア）ブロックは、局所麻酔薬を静脈内に注入し、標的部位の循環を遮断する四肢の麻酔法である。この技術は、通常、四肢から体幹に向けて血液を強制的に排出する標的領域の駆血後、続いて血流を安全に停止するために空気圧ターニケットを使用する。麻酔薬は静脈から四肢に注入され、ターニケットが目的の領域内に麻酔薬を保持している間に、周囲の組織に拡散させる。
日本語名称は日本麻酔科学会用語集第五版では、静脈内区域麻酔とされるが、他に静脈内局所麻酔(法)、局所静脈内麻酔(法)などとも呼ばれる。
歴史

止血帯と注射による局所麻酔は、1908年にアウグスト・ビーア（英語版）によって初めて紹介された。彼はエスマルヒ駆血帯（英語版）を使用して腕を駆血し、2つの止血帯の間にプロカインを注射して、その部位に麻酔効果と鎮痛効果を急速にもたらした。効果的であることが証明されたものの、C.McK. ホームズが1963年に再導入するまでは一般的では無かった。今日、この手技は、経済性、迅速な回復、信頼性、簡便性から一般的になっている。
方法

プロトコルは、地域の標準的な手順と手術対象の四肢によって異なる。一般的な手技は、ビーアが行ったように弾性包帯(エスマルヒ駆血帯（英語版）)を使用して四肢を駆血することから始め、血液を心臓に向かって近位に絞り出す。次に、空気圧ターニケットを四肢に装着し、動脈圧より30mmHg高く膨張させてすべての血管を閉塞し、次に弾性包帯を除去する。高用量の局所麻酔薬、通常はアドレナリンを含まないリドカインまたはプリロカインを、駆血した四肢のできるだけ遠位にゆっくりと注入する。静脈は麻酔薬で満たされ、約6～8分後に麻酔薬が局所組織に浸透する。その後、その領域の手術、または整復（英語版）を開始する。この時点で、患部が活発な血流から隔離されていることが重要である。鎮痛効果は、使用する麻酔薬の投与量と種類にもよるが、最大2時間持続する。低血圧、痙攣、不整脈、死に至る可能性のある血液中の麻酔薬の過剰摂取を避けるためには、待機時間と患部の血流遮断が重要である。ブピバカインやエチドカインなどの心毒性のある局所麻酔薬は、厳禁である。
安全性

IVRAの安全性と有効性は、臨床文献で十分に確立されている。ただし、ブピバカインやエチドカインなどの心毒性のある局所麻酔薬は禁忌である。IVRAを遠位肢、特に前腕に適用する場合は、患者に止血帯使用の禁忌がある場合（例えば、低酸素分圧による大量溶血や血流制限による溶血クリーゼのリスクがある鎌状赤血球症など）を除き、処置時間を短くする（最長2時間）ことが望ましい。IVRAに関連する合併症のシステマティックレビューでは、1964年から2005年の間に64例が報告されており、他の手技と比較して良好な結果が得られている。IVRAに関連する合併症の多くは、麻酔薬の種類、機器の不適切な使用や選択、技術的なミスが大きな要因となっている。現在では、患者の安全性を向上させるために、様々な安全対策がとられている。
安全性 > 装置

麻酔科医や外科医からの報告では、機器の適切な選択、検査、メンテナンスが重要な安全対策として挙げられている。最も安全な止血帯装置には、各チャネルの独立した四肢閉塞圧力測定などのIVRA固有の機能や、ヒューマンエラーを減らす専用の安全ロックアウトと組み合わせた二連式止血帯カフが必要である。さらに、IVRAプロトコルには、手術前に、手動または自動にかかわらず、機器の定期的な予防保守および性能テストの手順を含める必要がある。
安全性 > 補助薬

補助薬は、麻酔作用を促進し、副作用を最小限に抑えることにより、IVRAの安全性を向上させる。たとえば、ベンゾジアゼピンとフェンタニルは、痙攣発作を予防するために、神経ブロックを改善するために、それぞれ追加されることがよくある。
安全性 > 使用時の安全対策

標準化された実践への順守を含む改善されたプロトコルも、合併症の可能性と影響を改善するのに役立つ。例えば、四肢保護パッドと密着可能な止血帯を適用することで、組織の損傷を防ぎ、過度ではなく十分なターニケット圧力により、麻酔薬を四肢内に確実に留まらせることができる。速すぎるカフの開放やカフの膨張不足を避けるように注意する必要がある。万一、合併症が発生した場合には、常に生体情報モニタを装着しておき、蘇生薬や蘇生器具をすぐに使用できるようにしておくと、迅速な対応が可能になる。
アウグスト・ビーア
東京都道224号神戸山多幸線

東京都道224号神戸山多幸線（とうきょうとどう224ごう こうべやまたこうせん）は、東京都神津島村字神戸山と同字多幸を連絡する都道（一般都道）である。
神津本道の標識
本線から太平洋を望む
路線データ

起点 東京都神津島村字神戸山
終点 東京都神津島村字多幸三浦港
延長 12,267m
面積 148,002m2
通称 神津本道（本線の全区間など）東京都通称道路名設定公告整理番号119＝1984年5月1日公告
支線
神津島空港アクセス道路
神津島漁協前 − 神津島港
都道認定要件 地方開発のため特に必要な道路
沿革

2017年4月1日 - 神戸山多幸線が路線認定、長浜多幸線が路線廃止
2020年4月1日 - 村道森田線を都道に編入
風早郡

令制国一覧 > 南海道 > 伊予国 > 風早郡
日本 > 四国地方 > 愛媛県 > 風早郡
愛媛県風早郡の位置
風早郡（かざはやぐん）は、愛媛県（伊予国）にあった郡。
古代

風速国造の領域を中心として建てられた。
古代 > 式内社

『延喜式』神名帳に記される郡内の式内社。
| 神名帳> 社名 | 神名帳> 読み | 神名帳> 格 | 神名帳> 付記 | 比定社> 社名 | 比定社> 所在地 | 比定社> 備考 | 集成 |
| --- | --- | --- | --- | --- | --- | --- | --- |
風早郡 2座（並小）表 話  歴
| 國津比古命神社 | クニツヒコノ | 小 |  | 國津比古命神社 | 愛媛県松山市八反地 |  |  |
| 櫛玉比賣命神社 | クシタマヒメノ- | 小 |  | 櫛玉比賣命神社 | 愛媛県松山市高田 |  |  |
凡例を表示
明治維新後

1878年（明治11年）に行政区画として発足した当時の郡域は、松山市の一部（小川、光洋台、西谷、大西谷、菅沢町、神次郎町、上総町、梅木町、九川以北および興居島・釣島を除く忽那諸島）にあたる。
明治維新後 > 沿革

「旧高旧領取調帳」に記載されている明治初年時点での支配は以下の通り。（85村）
| 知行 | 知行 | 村数 | 村名 |
| --- | --- | --- | --- |
| 幕府領 | 松山藩預地 | 2村 | 粟井村、小浜村 |
| 藩領 | 伊予松山藩 | 78村 | 土手内村、辻村、中西内村、中西外村、八反地村、寺谷村、神田村、波田村、猪木村、滝本村、猿川村、中村、米之野村、庄府村、儀式村、小山田村、猿川原村、尾儀原村、才之原村、庄村、萩原村、上難波村、中通村、浅海原村、浅海本谷村、下難波村、北条村、中須賀村、片山村、夏目村、苞木村、鹿峰村、河原村、久保村、常竹村、安岡村、和田村、鴨之池村、磯河内村、小川村、大西谷村、西谷村、本谷村、小川谷村、平林村、麓村、客村、菅沢村、神次郎村、城山村、恩地村、柳谷村、牛谷村、大河内村、横谷村、善応寺村、高山村、佐古村、常保免村、宮内村、別府村、閏谷村、九川村、梅木村、上総村、小屋村、野忽那村、畑里村、饒村、吉木村、熊田村、長師村、宮野村、神浦村、二神村、津和地村、院内村、安居島村 |
| 藩領 | 伊予大洲藩 | 4村 | 睦月村、宇和間村、上怒和村、元怒和村 |
| 幕府領・藩領 | 松山藩預地・大洲藩 | 1村 | 大浦村 |
慶応4年1月27日（1868年2月20日） - 戊辰戦争により松山藩預地が大洲藩預地となる(「旧高旧領取調帳」では高知藩預処となっているが、国立公文書館所蔵「府藩県高帳」では風間郡内高754石9斗8升5合が大洲藩御預所となっており、「旧高旧領取調帳」の誤記と考えられる)。
明治4年
1月19日（1871年3月9日） - 大洲藩預地(風間郡内大浦村外二箇村)が倉敷県の管轄となる。
7月14日（1871年8月29日） - 廃藩置県により、藩領が松山県、大洲県の管轄となる。
9月14日（1871年10月27日） - 倉敷県の管轄地域が丸亀県の管轄となる。
11月15日（1871年12月26日） - 第1次府県統合により、丸亀県の管轄地域および旧松山藩領が松山県、旧大洲藩領が宇和島県の管轄となる。
明治5年
2月9日（1872年3月17日） - 松山県の管轄地域が石鉄県の管轄となる。
6月23日（1872年7月28日） - 宇和島県の管轄地域が神山県の管轄となる。
明治6年（1873年）2月20日 - 全域が愛媛県の管轄となる。
明治11年（1878年）12月16日 - 郡区町村編制法の愛媛県での施行により、行政区画としての風早郡が発足。「風早野間郡役所」が北条村に設置され、野間郡とともに管轄。
明治17年（1884年）11月 - 「風早和気温泉久米郡役所」が温泉郡松山榎町に設置され、同郡および和気郡・久米郡とともに管轄。
明治22年（1889年）12月15日 - 町村制の施行により、下記の各村が発足。全域が現・松山市。（12村）
浅海村 ← 浅海本谷村、浅海原村
難波村 ← 下難波村、中通村、上難波村、庄村［大部分］
立岩村 ← 才之原村、滝本村、猪木村、猿川村、中村、米之野村、庄府村、儀式村、小山田村、猿川原村、尾儀原村、萩原村、庄村［一部］
正岡村 ← 八反地村、中西内村、中西外村、寺谷村、院内村、波田村、神田村
北条村 ← 北条村、辻村、土手内村、安居島村
河野村 ← 別府村、宮内村、善応寺村、横谷村、閏谷村、九川村、常保免村、佐古村、高山村、牛谷村、大河内村、中須賀村、片山村、夏目村
粟井村 ← 小川村、磯河内村、鴨之池村、和田村、河原村、安岡村、常竹村、本谷村、西谷村、大西谷村、客村、麓村、平林村、小川谷村、久保村、鹿峰村、苞木村
五明村 ← 菅沢村、神次郎村、城山村、柳谷村、恩地村、小屋村、梅木村、上総村
睦野村 ← 睦月村、野忽那村
東中島村 ← 大浦村、小浜村、長師村、宮野村、神浦村
西中島村 ← 粟井村、畑里村、饒村、吉木村、熊田村、宇和間村
神和村 ← 元怒和村、津和地村、二神村、上怒和村
明治30年（1897年）4月1日 - 郡制の施行により、温泉郡・久米郡・風早郡・和気郡および伊予郡・下浮穴郡の各一部の区域をもって、改めて温泉郡が発足。同日風早郡廃止。
21.浅海村 22.難波村 23.立岩村 24.正岡村 25.北条村 26.河野村 27.粟井村 28.五明村 29.睦野村 30.東中島村 31.西中島村 32.神和村（紫：松山市 1 - 9は温泉郡 11 - 15は久米郡 41 - 50は和気郡 51 - 56・61 - 68は下浮穴郡 71 - 80は伊予郡）
イローナ・スルピアネク

イローナ・スルピアネク (Ilona Slupianek、1956年9月24日 - )は、旧東ドイツの陸上競技選手。女子砲丸投の元世界記録保持者であり、1980年モスクワオリンピック女子砲丸投の金メダリストである。
主な実績

| 年 | 大会 | 場所 | 種目 | 結果 | 記録 |
| --- | --- | --- | --- | --- | --- |
| 1977 | IAAFワールドカップ | デュッセルドルフ(西ドイツ) | 砲丸投 | 1位 | 20.93m |
| 1978 | ヨーロッパ陸上選手権 | プラハ(チェコスロバキア) | 砲丸投 | 金 | 21.41m |
| 1979 | ユニバーシアード | メキシコシティ(メキシコ) | 砲丸投 | 金 | 20.48m |
| 1979 | IAAFワールドカップ | モントリオール(カナダ) | 砲丸投 | 1位 | 20.98m |
| 1980 | オリンピック | モスクワ(ソビエト連邦) | 砲丸投 | 金 | 22.41m |
| 1981 | IAAFワールドカップ | ローマ(イタリア) | 砲丸投 | 1位 | 20.60m |
| 1982 | ヨーロッパ陸上選手権 | アテネ(ギリシャ) | 砲丸投 | 金 | 21.59m |
| 1983 | 世界陸上選手権 | ヘルシンキ(フィンランド) | 砲丸投 | 銅 | 20.56m |
| 1987 | 世界室内陸上選手権 | インディアナポリス(アメリカ合衆国) | 砲丸投 | 銀 | 20.28m |
世界記録

22.36m – 1980年5月2日
22.45m – 1980年5月10日
井狩春男

井狩 春男（いかり はるお、1945年3月13日 - ） は、日本のエッセイスト・評論家。埼玉県出身。
経歴・人物

中央大学第二文学部仏文科中退。その後、書籍取次店の鈴木書店に入社。店売で商品知識を身につけ、やがて対出版社の窓口である仕入部に配属される。出版界初の手書きによる情報誌「日刊まるすニュース」を一人で担当する。
1985年、『返品のない月曜日』によりデビュー。
新風舎出版賞の審査委員長も務め、さらに自身の著書『本は読むより書く方が10倍楽しい』でも同賞に応募することを勧めているなど、新風舎との関わりが深かった。