OKAZAKI MANUFACTURING COMPANY

TD-1116R4

機械工学の基礎解説
保護管の強度計算を理解するために

目次
1 . はじめに2 . 強度計算のための基礎知識
3 . 材料力学
4 . 流体力学
5 . 機械力学
6 . 熱力学【参考文献】【参考資料】【計算例】-------------------------------------

2 page2 page2 page3 page6 page8 page9 page9 page11 page

株式会社岡崎製作所
東京技術部
作成日：2014.10.28
改

-1-

定：

4

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

１． はじめに温度計保護管の強度計算を理解するためには、機械工学系の様々な基礎知識を必要とします。強度計算の解説を読む前に、以下に記す基礎知識を知ることで一層の理解が深まります。既に機械工学系の知識を持っている方は復習と思ってお読み下さい。
２． 強度計算のための基礎知識必要とする分野は以下のようなものです。
（１）材料力学：材料に生ずる応力に関する知識（２）流体力学：流体によって受ける力や、発生する渦流に関する知識（３）機械力学：固有振動数と共振にかかわる知識（４）熱力学：ガスの場合の密度を求めるための知識
３． 材料力学
基本的には材料の破壊にかかわる力学で、静的破壊の他、疲労破壊、クリープ破壊および脆性破壊がある。これらは材料に加わる荷重と発生する応力を評価する方法とも言える。裏返して言えば、材料の持
つ強度の評価である。
（１） 荷重機械および構造物は変形する固体から成り立っている。固体の表面に外部から作用して、それを変形させる力学作用を荷重（load）といい、荷重を作用させることを負荷するという。荷重は大きさと向きを持った力であり、SI 単位ではＮ（ニュートン）。質量に重力の加速度を掛け合わせたものが重量になり、同様の単位になる。γ（重量）＝ρ（密度）･ｇ（重力の加速度）
荷重の分類法は色々あるが、基本的に表-1 に分類できる。
表-1

速度による荷重の分類

静荷重（static load）

動荷重（dynamic load）

静止する荷重繰り返し荷重（repeated load）

大きさが繰り返して周期的に変動

交番荷重（alternate load）

大きさと向きも周期的に変動

ランダム荷重（random load）

大きさがランダムに変動

衝撃荷重（impact load）

急激に作用

その他の分類法として、荷重の分布様式による分類（等分布荷重、不等分布荷重、集中荷重）
および作用による分類（垂直荷重、せん断荷重、引っ張り荷重、圧縮荷重等）がある。
（２） 応力
物体に力が作用すると変形し、この変形に抵抗して物体内部に力を生じる、この力を内力（internal force）という。物体を仮想断面で切り取れば、仮想面上の内力は外力とつり合っており、単位面積あたりの内力を応力（stress）という。応力を仮想面上に垂直な方向と平行な方向に分解して垂直な成分を圧縮応力（normal stress）σ、平行な成分をせん断応力（shearing stress）τといい、仮想面を引っ張る［圧縮する］時の垂直応力を引っ張り応力（tensile stress）［圧縮応力（compressive stress）］と呼ぶ。応力は圧力と同様の単位になるが、

図-１

応力の定義

圧力は物体の表面に均等にかかる荷重の単位面積あたりの値としたもので、圧力に受圧面積を掛けたものが荷重となる。この荷重により材料内部に発生するものが応力である。SI 単位系ではN/mm2 またはPa（パスカル）で表される。

-2-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4（３） はり棒が外力により曲げ作用を受ける場合、これをはり（beam）という。このはりに曲げ作用を及ぼす外力としてははりに垂直や傾斜を持った荷重W、長さ方向に分布した荷重w、および曲げモーメントM などがある。図-1 のような場合、荷重W1 や分布荷重w1 は下向きを正とし、反力R などは上向きを正とする。この場合、はりのある断面mn にはその断面の左側または右側に作用する荷重や反力の総和に等しい力Fとそれらによる偶力Ｍが作用するが、この力Fをせん断力（shearing force）、偶力Mを曲げモーメント（bending moment）という。
W1

m

m

w1

M
x

R1

n

n

R2

F

y

図−２

はりに加わる荷重

はりに生ずる応力は、荷重に比例するが、はりの断面形状により大きく異なる。基本的にはりに生ずる応力 σ＝M/Z で表される。ここで、Zを断面係数という。断面係数は各種便覧等に求め方が掲載されている。また、はりの断面形状によって決まる、断面二次モーメントＩも重要な要素であり、各種計算に使用される。保護管の計算で最も多く使用されるパイプ状断面の断面
二次モーメントと断面係数を表-2 に示す。
表−2
断面形状

パイプ断面

断面二次モーメントI𝜋
𝑑2! − 𝑑1!
64

断面係数

Z

𝜋 𝑑2! − 𝑑1!
32𝑑2（４） 曲げモーメント
一般にはりは支持される方法により、荷重を受け発生する応力は異なる。保護管等の支持は片端固定で他端は自由である。この様な場合荷重を受けることによって発生する曲げモーメントの最大値は根本（固定部分）に発生する。曲げモーメントは荷重が一点集中荷重と仮定すると、M=F×Lで表され、ここでL はモーメント
アームといい、荷重の集中した点から固定点までの距

図−３

片持ちはり

離である。このはりを片持ちはりと呼ぶ。
４． 流体力学
元々は水力学から発展した学問分野で、近世になって航空技術の発展とともに空気力学系が更なる発展を遂げた。物体は固体、液体および気体の三つの形態を取る。その中で液体と気体のように、形を容易に変えるものを通常は流体（Fluid）と呼ぶ。但し、気体の密度は圧力や温度の変化とともに0 近くまで変化するようにモデル化されるが、液体の密度はほとんど変わらない。
通常工業的に扱う流体の基本的な性質は粘性と圧縮性であるが、音速を超える領域以外では圧縮性の影
響は少ない。従って、粘性の影響が最も大きくなる。

-3-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4（１） 流体の性質
流体の性質を表すときにレイノルズ（Re）数がよく取り上げられる。ある物体を対象としたときに、その周囲を流れるときの状態を表すための無次元量で、動粘度（ν）、代表寸法（D）および流体の速度（q）で表現される。
𝑅𝑒 =

𝑞∙𝐷
− − − − − − − (1)𝜈

また、流体の物性を表すものに、密度（ρ）、比体積（v）および比重がある。物質の単位体積あたりの質量を密度（density）、単位質量あたりの体積を比体積（specific volume）といい、それぞれ互いに逆関数の関係にある。水の密度に対する他の物質の密度の比を比重（specific gravity）という。
その他に、前述のRe 数を求める際に必要となった、流体の粘性を表す粘度（µ）と動粘度（ν）
がある。流体に速度勾配があって流れているとき、流れに平行な面に沿ってせん断力 τ が作用する。流体の持つこの性質を粘性という。通常流体は配管中を流れることが多い。その場合、管壁近くでの流速は大きな速度勾配を持つ。従って、そこに生じるせん断力が配管の抵抗となるわけである。また、温度計の保護管のように流体中に固定されて置かれた円柱は、その近傍に大きな速度勾配を生じさせる。その結果流体の粘性により保護管に荷重がかかる。これを抗力とよんでいる。従って、前述のRe 数は同一の流体であっても、配管の抵抗を考慮する場合はD=配管内径であり、保護管の抗力を求める場合は、D=保護管外径として算出するため異なった値になる。
（２） 流体中の物体に働く力
流体中の三次元物体には一般に揚力、抗力、横力、ローリングモーメント、ヨーイングモーメ
ントおよびピッチングモーメントが作用する。ここで、保護管の場合には一端支持の円柱に限定されるため、前者の3 つの力のうち二つを考えれば良く、モーメントは無視できる。それぞれの力は流体の密度と速度に比例する計算式で与えられる。ここで、

𝜌𝐹 = 𝐶! 𝑉 ! 𝐴 − − − − − − − (2)2F ：揚力または抗力（N）
V ：流体の速度（m/s）
A ：保護管の投影面積（m2）Cx ：揚力（CL）または抗力（CD）係数
ρ ：流体密度（kg/m3）
この各係数は、世界中の様々な研究者により研究されている。さらに保護管（流体中の丸棒）の場合、上述の力の他、カルマン渦による変動揚力や変動抗力が生ずることが研究から明らかになっている。以下に、考慮するべき抗力と変動揚力の考え方を述べる。
１） 抗力
流れに垂直におかれた円柱の抗力係数CD はRe 数に対し図-3 のように変化する。高粘度流体では粘度が大きいため、Re数が極めて小さくなり、抗力係数CDが大きくなる。そのため、
攪拌機等において高粘度流体の温度を測定する際は、保護管の曲がりに注意が必要である。
低いRe数範囲では各種近似式が成り立つ。Re = 2×104 〜2×105 ではCD = 1.2 でほぼ一定である。強度計算において、正確な粘度を知る事は極め重要である。

-4-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

Re数と抗力係数CD

図−４２） カルマン渦（変動揚力）

流れにおかれた円柱の後流側には交互に渦が発生する。これをカルマンの渦列といい、その渦数は流体速度 V とストローハル（St）数と円柱の外径 D によって決まる。Re 数とSt数の関係は図−５で表される。ここで、カルマン渦数 fs は式（３）表される。一般的にこのカルマン渦が保護管に強制振動を与えるとされているため、強制振動数ともいわれる。
𝑓! =

図−５

𝑆𝑡 ∙ 𝑉− − − − − − − (3)𝐷

Re数とSt数３） 対称渦（変動抗力）
1995年に発生した「高速増殖炉もんじゅ」の事故以来、話題になった対称渦について簡単に触れたい。詳細は別の書類にまとめるが、その発生原因は二通りあるようである。
まず一つは、円柱の交流にできる渦は交互に発生している間、渦の放出振動数fsと同一の変動揚力とその2倍の振動数に等しい変動抗力が働く。変動揚力は左右の渦一つずつに対し１回発生するが、変動抗力は一つの渦に対し１回発生するからである。従って、円柱が剛体でない場合はその固有振動数fnとfsが一致すると共振現象を生じ変動揚力は極大値となる。
2fs と一致すると、変動抗力方向の振幅が極大となる。
もう一つは、無次元振動数（換算流速） Vr = V/(fn･D) が1.25 〜2.5 近傍で自励振動的（流体−構造連成振動的）な対称渦が発生する現象である。この現象は流れの中にある円柱の先
端部の動きがある程度以上大きくなったときに生じ、機械的な円柱の動きと流れによる影響で相互渦が対称渦に転ずる。これらの現象は、当社が異なった2箇所（機関）において、液体を用いた流体実験を行った結果、Vr = 2 近傍での抗力方向の共振が確認されている。
これらの計算手法はJSME S012に詳細に記載されている。しかしながら、このVr = 2 近傍での抗力方向における共振現象を、ASME PTC 19.3 TW-2010 においては計算時に考慮していないため、注意が必要である。

-5-

OKAZAKI MANUFACTURING COMPANY

図−６

TD-1116R4円柱の振動形態（ASME TW-2010及びJSME S012）

５． 機械力学（機械振動学）
機械工学で振動学の知識が応用されるのは、問題の機械系をいくつかのばねと集中質量、あるいはいく
つかのばねと支持された剛体と見なして、その系の固有振動数を計算するという場合が大部分である。
一般的な１自由度振動には自由振動と強制振動があり、それぞれに減衰の作用しない場合と、減衰が作用する場合の振動がある。
（１）自由振動１） 無減衰自由振動（単振動）
外部から力が作用しない時の振動を自由振動といい、x(t) = A cos（ωt + φ）で表される。変位 x は時間 t とともに正弦的に変化する。この様な振動を単振動という。
A：振幅、ω：固有角振動数、φ：初期位相で、振動数 f、
周期 T および角振動数 ω の間には f = 1/T = ω / (2 π)の関係がある。この様に減衰のない時の振動数を系の固

図−７

単振動

図−８

減衰振動

有振動数（Natural frequency）と呼ぶ。
２） 粘性減衰が作用する自由振動
実際にはすべての振動系には何らかの粘性減衰が作用し、減衰係数比 ζ の大きさにより振動の様子が変化する。ζ ≧1 の場合は無周期運動となり、ζ ＜ 1 の場合は振幅が時間とともに減少する運動、すなわち減衰振動となる。この場
合の振動が下式で表される。

𝑥 = 𝑒 !!"# 𝑥! cos 𝜔𝑡 +𝜁𝜔𝑥! + 𝑣!
sin 𝜔𝑡𝜔− − − − − − − (4)（２）強制振動１） 無減衰強制振動
振動体が外部から作用する周期的外力によって行う定常振動を強制振動という。減衰の作用しない強制振動方程式の解は自由振動解と強制振動解の重ね合わせとなり、外力の振動数と固有振動数が一致すると振幅は時間とともに大きくなる。
２） 粘性減衰が作用する強制振動
粘性減衰が作用する強制振動は自由振動解と強制振動解の重ね合わせとなるが、自由振動
解は時間とともに減衰し、強制振動解だけが残る。従って、振動体の運動は周期的外力と-6-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

同じ振動数を持ち位相がずれただけの単振動である。その振幅は静的変位の倍数として次式で表される。
𝐹! =

𝑋
=
𝑋!

1
1 − 𝑓!

! !

− − − − − − − (5)+ 2𝜁𝑓!

!

ここで、fkは振動数の比でfk=fs/fnで表される。この式中で両者の振動数が一致した値でFM は最大値を示す。この状態を共振状態と呼ぶ。

図−９

振幅倍率曲線（３） 固有振動数
温度計の保護管を対象とする場合、その保持方法と形態から、「片持ちはり」という分類になる。材料力学の項の図-2 に示した形状である。
このときの固有振動数は下式で表される。（ここで使用する単位は、強度計算の解説で使用するものと同一とする）この固有振動数は一様断面はりの横振動数で、振動のモードにより、1 次、2次および3次までの振動数係数 λ が機械工学便覧等に掲載されている。式（６）中の1012は縦弾性係数（N/mm2）と材料密度（kg/m3）の単位を合わせるための乗数である。
𝑓! =

𝜆!
𝐸∙𝐼
∙ 10!" − − − − − − − (6)!
2𝜋𝐿 𝑤 ∙ 𝐴!

-7-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

６． 熱力学
ここで扱う熱力学は、本来の分野のかなり狭い部分のみである。どちらかというと物理に近い分野といえる。一般にプロセス流体を扱う場合、その対象となるものは液体、蒸気あるいは気体（ガス）である。
これらの物質の状態はその種類と温度、圧力によって変化するのが常である。流体の性質を知るためには、温度、圧力、流速、粘度、密度等が重要な要素となる。液体の場合は圧力の変化で密度は大きく変化しないが、気体や蒸気は圧力、温度によってその性質は大きく変化する。
（１） ガスの性質
以降に述べるボイル- シャール（シャルル）の法則に厳密に従うガスを完全ガス（Perfect Gas）
あるいは理想気体（Ideal Gas）と称する。この2 法則はガスの分子が容積を有せず、かつまた分子間に引力が働かないと仮定した場合に成り立つ。ところが実在ガスでは分子が容積を有し、また分子間に引力が存在している。従って完全ガスは正確には理想的ガスにすぎない。しかし、空気、水素、酸素、窒素およびヘリウムなどは普通の状態において完全ガスと見なしても実用上差し支えない。これに対して水、アンモニア、炭酸ガスおよび亜硫酸ガスなどの蒸気は、普通の状態では完全ガスとして取り扱えないが、これらにおいても過熱の度が高まるにつれて次第に完全ガスの性質に近づくのである
１） ボイル- シャールの法則（Boyle and Charles' law : Combined gas law）
気体の状態方程式は絶対圧力をP (Pa)、比体積をv (m3/kg)、密度をρ (kg/m3)、絶対温度T (K)、
気体定数（gas constant）をRとすれば、Pv = RTあるいは

ρ =P/(RT)で表される

２） アボガドロの法則
一般気体定数(universal gas constant) R はアボガドロの法則に導かれて、すべての気体に対して同一となる。すなわち分子量をM とすればR =MR=8314.8 J/(kmol･K)で表される。
３） 完全ガスの密度計算
１）と２）の関係から、完全ガスのモル数（kg/kmol）が判るとガスの密度が計算できる。
R=8314.8/M

であるから

ρ =P/(8314.8/M × T )kg/m3

ここで、P（Pa）、M（kg/kmol）、T（K）の単位系である。
（２） 蒸気の性質
蒸気（水蒸気）は前項で述べたように、一般の完全ガスと同等には扱えず、日本機械学会（JSME）
発行の蒸気表を用いて、各温度・圧力における密度と粘度を求める。蒸気の密度をモル数（H2Oとすれば18）で与えると、完全ガスと見なして計算する事になるため、注意が必要である。蒸気表を読む際の注意は、表に記載の圧力は絶対圧力で表記され、Pa（=1.01972 × 10-5 kg/cm2）表記と言うことである。温度、圧力ともに高い場合には過熱蒸気となり、完全ガスに近づく。
JSMEの蒸気表の付録CDにあるプログラムで、温度・圧力を入力すれば、密度・粘度は自動計算できる。
（３） 粘度と動粘度
粘性(viscosity)はせん断、あるいは角変形に対する流体の抵抗であり、流れの速さの違いをならして一様にしようとする流体の性質をいう。粘性の程度（大きさ）を表すのが粘度で、粘性係数、粘性率と呼ぶこともある。動粘度(kinematic viscosity)と区別するために絶対粘度とも呼ばれる。粘度は記号として µ 又は η が用いられ、SI単位はPa･s（パスカル秒）である。CGS単位系
では P（ポアズ）が用いられた。動粘度は記号として ν が用いられ、粘度 µ との関係は流体密度を ρ とすると、ν＝μ/ρで表される。SI単位はm2/sで表され、ストークス（St）という単位を用いることもある。

-8-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4【参考文献】
１． 機械工学便覧：A4

材料力学

日本機械学会編

２． 機械工学便覧：A5

流体工学

日本機械学会編

３． 演習

水力学

著者：国清行夫、木元知夫、長尾健

４． 日本機械学会基準
５． 機械振動学通論

JSME S 012-1998

配管内円柱状構造物の流力振動評価指針
著者：入江敏博

６． 機械工学便覧：A3
７． 応用熱力学

発行：森北出版力学

発行:朝倉書店

機械力学

著者：小野栄一郎

日本機械学会

日本機械学会編

校閲：菅原菅雄

発行：産業図書株式会社出版【参考資料】
以下に参考資料として、各種単位の換算表を掲げる。基本的に SI 単位系を主としているが、海外 JOB でフート･ポンド系単位もあるため、それらも入れた。末尾に、工学系で使用されるギリシャ文字の解説を付した。
１．長さの単位
m

km

in

ft

1

0.001

39.37

3.28084

1000

1

39370

3280.84

0.02540.0000254

1

0.083332

0.30480.0003048

12

1

lb/in3lb/ft3２．密度の単位
kg/m3

g/cm3

1

10

103

1

27680

-3

3.613×10-5t(英)/yd3

lb/gal(英)lb/gal (米)-40.010028.345×10-3

0.06243

7.52×100.03613

62.43

0.752

10.02

8.345

27.68

1

1728

20.83

277.4

231

16.01850.016020.0005787

10.01205

0.16050.1337

1329

1.329

0.04801

82.96

1

13.32

11.09

98.8

0.09980.003605

6.229

0.075081

0.8327119.8

0.11980.004329

7.481

0.090171.201

1３．圧力・応力の単位
MPa
N/mm2

Pa
N/m2

Mdyn/cm2bar

kgf/cm2
at

1

1×106

10

10.197

1×10

-6

1

0.1

1000009.80665×10
6.895×10

1×10

-2

-3

0.101325
1.33332×10-4
9.8065×10-39.80665×10

46895
101325(1013.25hPa)
133.329806.65-5lbf/in2psi

atm

145.041.0197×10-51.4504×10

mmHg
Torr

9.869
-4

mAq

7500.629.8693×10-67.50062×10

101.97-3

1.0197×10-4

11.019714.5040.9869750.062

10.1970.980665

1

14.2233

0.9678735.559100.06895

0.0703071

0.0680551.71

0.70311.01325

1.0332

14.6961760

10.3320.0013332

1.35951×10-3

0.019340.001316

10.0135951

1.4223-2

73.5559

1

0.0980665

0.1

-9-

9.678×10

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

４．力の単位（重力の加速度 g=9.807m/s2 とする）
ニュートン（SI 単位）
N

重量キログラムkgf

1 (kg)×1 (m/s2)

1 (kg･m/s2)/9.80665 (m/s2)=0.101972 (kgf)

1N

重量ポンド
lbf
2

0.2248 lbf

1kgf

9.80665 (N)1 (kg)×9.80665 (m/s )2.2046 lbf

1lbf

4.448 N0.4536 kgf

1 (lb) ×9.80665 (m/s2)

５．粘度の単位
５−１．粘度・絶対粘度
Pa･s

mPa･s

P

cP

kgf･s/m2

lbf･s/in2

1

1×103

1×10

1×103

1.0197×10-11.4504×10-4

1×10-31

1×10-21

1.0197×10-11.4504×10-7

1×10-1

1×102

1

1×102

1.0197×10-11.4504×10-5

1×10-31

1×10-21

1.0197×10-11.4504×10-7

9.80679.8067×10

9.8067×103

9.8067×10

11.4223×10-3

6.8948×103

6.8948×106

6.8948×104

6.8948×106

7.0307×10

1５−２．動粘度
m2/s

m2/s
St

mm2/s
cSt

1

1×104

1×106

1×10-41

1×102

1×10-6

1×10-2

1

６．ギリシャ文字
大文字

小文字

英語表記

A

α

alpha

B

β

Γ

読み方

大文字

小文字

アルファ

N

ν

Nu

ニュー

beta

ベータ

Ξ

ξ

xi

グザイ

γ

gamma

ガンマ

Ο

ο

omicron

オミクロン

Δ

δ

delta

デルタ

Π

π

pi

パイ

E

ε

epsilon

イプシロン

Ρ

ρ

rho

ロー

Z

ζ

zeta

ツエータ

Σ

σ

sigms

シグマ

H

η

eta

イータ

Τ

τ

tau

タウ

Θ

θ

theta

シータ

Υ

υ

upsilon

ウプシロン

I

ι

iota

イオタ

Φ

φ

phi

ファイ

Κ

κ

kappa

カッパ

X

χ

chi

カイ

Λ

λ

lambda

ラムダ

Ψ

ψ

psi

プサイ

M

μ

mu

ミュー

Ω

ω

omega

オメガ

-10-

英語表記

読み方

OKAZAKI MANUFACTURING COMPANY

TD-1116R4【計算例】
１． 片持ちはりの計算例
図−Ａに示す片持ちはりの根元に生じる応力 σ と先端のたわみ δ を求める。
ここで、はりは外径φ33 の SUS304 丸棒とし、材料密度 ρm：7980 kg/m3
縦弾性係数 E：195000 N/mm2

図−Aであるとする。

片持ちはり（１） はりに加わる荷重
はりに加わる荷重は自重分の当分布荷重で、w N/m とすれば、その値は下式で表される。

ｇ：重力の加速度 m/s2
w＝ρ × g × π/4 × D2
＝7980×9.80665×π/4×0.0332＝

66.93 N/m

断面 2 次モメント
I = π/64×D4 = π/64×0.0334 = 5.82 10-8 m4
先端のたわみ量 δ m

機械工学便覧より

!

𝛿=

𝑤∙𝐿
66.93×2!
=
8∙𝐸∙𝐼
8×195×10! ×5.82×10!!

= 0.011795 m= 11.8 mm（２） 根元の曲げモメント M
𝑀=

機械工学便覧より𝑤 ∙ 𝐿!
66.93×2!
=
2
2= 133.86 N･m
断面係数
𝑍=

Zは

𝜋 !
𝜋
𝐷 =
0.033!
32
32= 3.528×10-6 m4
根元の曲げ応力σ N/mm2

σ＝M/Z＝133.86 / 3.528×10-6 ＝37.18×106 N/m2
＝37.94 N/mm2
２． 流体中の円柱に働く力
図—B に示す円柱が、流速 V m/s の流対中に曝されていた場合に円柱に働く力を計算する。ここで、円柱は外径φ25 の SUS304 丸棒とし、材料密度 ρm：7980 kg/m3
縦弾性係数 E：183000 N/mm2
流体は水とし、
流速：3m/s
温度：200℃
圧力：5MPaであるとする。

図−Ｂ-11-

流体中の円柱OKAZAKI MANUFACTURING COMPANY

TD-1116R4（１） 流体による力
流体により、円柱が受ける力 F（N）は本文中の式（２）で表される。流速は全長にわたって
均一に加わる物とした。
𝐹 = 𝐶!

𝜌! !
𝑉 ∙𝐴
2

ここで、
CD：抗力係数0.6（以下の理由による）Re 数から、抗力係数を求める。
5MPa、200℃における水の粘度 µ と密度 ρ は蒸気表からµ = 0.135 mP･s（0.000135Pa･s）ρ = 868.4 kg/m3
動粘度 ν は ν = µ/ρ で表されるので、
Re 数は
𝑅𝑒 =

𝑉 ∙ 𝐷 𝑉 ∙ 𝐷 ∙ 𝜌 3×0.025×868.4==
𝜈
𝜇0.000135= 4.82×105
よって、図−４から CD=0.6 とした。
A：円柱の投影面積 m2
ρf：流体密度 868.4 kg/m3
従って、
𝐹 = 0.6×

868.4 !
3 ×0.025×0.2
2= 11.725 N（２） 流体により受ける応力
円柱が F（N）の力を受けた時に、円柱の根元に働く応力を評価。
通常、根元に働く応力 σ は下式で表される。
𝜎=𝑀𝑍

M：曲げモーメント kg･m
Z：断面係数 m3
ここで、流速による力は、等分布加重の合計であるため、円柱の中心部に集中して加わる力と考えることが出来る。曲げモーメントは M = F×L/2

であらわされる。

従って、根元に働く応力は𝐿
200
𝐹×11.725×
2
2
𝜎= 𝜋
=
𝜋
!
!
𝐷
25
32
32
= 0.76 N/mm2
この応力が、材料の許容応力以下であれば使用可能と判定できる。
通常、許容応力は「高圧ガス保安法

特定設備検査規則」に規定されている値を採用する。

３． 流対中の円柱の振動
図−Ｂに示す流対中の円柱の後流には、渦が出来る。流体条件によって出来る渦の種類は異なるが、その渦により円柱の前後又は左右に圧力差が生じることで、円柱に渦数に比例した繰り返しの力が加わ-12-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4

り、円柱の固有振動数と一致した場合に共振等の現象を起こす。以下に、円柱の固有振動数 fn を計算し、円柱の後流に発生する渦（カルマン渦）数 fs との共振が生じないか検証する。
（１） 固有振動数 fn と強制振動数 fs の計算
本文に示した式（６）に示す計算式で、円柱の固有振動数が計算できる。
𝑓! =

𝜆!
𝐸∙𝐼
∙ 10!"!
2𝜋𝐿 𝑤 ∙ 𝐴!

𝜋
183×10! ∙×25!
1.875!
64=∙ 10!"𝜋
2𝜋×200!
7980 ∙ ×25!
4
＝1.399×10-5×29.93×106 ＝418.7 Hz
強制振動数は式（３）で計算できる。
𝑓! =𝑆𝑡 ∙ 𝑉𝐷=

0.21×30.025= 25.2 Hz
振動比 fs/fn = 0.06 < 0.2 となり、使用可能な条件となる。
【注記１】
流体密度が小さいガス又は蒸気の場合は、固有振動数 fn と強制振動数 fs が一致しないように
fs/fn ＜ 0.8

で使用可能とするが、流体密度が大きな液体の場合は対称渦による共振を考慮して

Vr ＜ 1 の条件とほぼ同等である fs/fn ＜ 0.2 で使用可能とする。
【注記２】
式（６）中の単位を合わせる手法を紹介する。
𝑓! =

𝜆!
𝐸∙𝐼
∙ 10!"!
2𝜋𝐿 𝑤 ∙ 𝐴!

それぞれの記号を上式で使われる単位として、単位のみ計算し、fn Hz（1/s）となることを確認。
𝑁∙ 𝑚𝑚 !
1𝑚𝑚 !
𝑚𝑚 ! 𝑘𝑔∙ 𝑚𝑚 !
𝑚!

ここで、𝑁 =𝑘𝑔𝑚 ∙ 𝑠!

𝑚 = 10! 𝑚𝑚
1𝑚𝑚 !

𝑘𝑔 ∙ 𝑚 1
𝑘𝑔 ∙ 10! 𝑚𝑚 1
∙ 𝑚𝑚 !
∙ 𝑚𝑚 !
1𝑠 ! 𝑚𝑚 !
𝑠!
𝑚𝑚 !
=
𝑘𝑔
𝑘𝑔
𝑚𝑚 !
∙ 𝑚𝑚 !
∙ 𝑚𝑚 !
!
𝑚
10! 𝑚𝑚 !

1
=
𝑚𝑚 !

=

𝑘𝑔 ∙ 10! 𝑚𝑚 1
∙ 𝑚𝑚 !
1𝑘𝑔 ∙ 10! 𝑚𝑚 ∙ 𝑚𝑚 ! ∙ 10! 𝑚𝑚 !
110!" 𝑚𝑚 !
𝑠!
𝑚𝑚 !
=
=
𝑘𝑔
𝑚𝑚 !
𝑠 ! ∙ 𝑚𝑚 ! ∙ 𝑘𝑔 ∙ 𝑚𝑚 !
𝑚𝑚 ! 𝑠 ! ∙ 𝑚𝑚 !
∙ 𝑚𝑚 !
10! 𝑚𝑚 !

1
10!" 𝑚𝑚 ! 1
=
10!"
𝑚𝑚 !
𝑠!
𝑠

-13-

OKAZAKI MANUFACTURING COMPANY

TD-1116R4以上から、単位は 1/s となり、乗数として 1012 が平方内に残ることが確認できた。
４． 完全ガスの計算
ガスの流量は標準状態での流量である、Nm3/h で与えられることが多い。この場合、実使用条件での流量に換算して、流速を求める。下記の条件の場合、実使用条件での流量はボイル-シャールの法則から
流量 QN：100,000 Nm3/h
圧力 PO：4MPa

標準状態 PN：0MPa

P=PN+PO

温度 TO：300℃

標準状態 TN：0℃（273.15K）

T=TN+TO𝑄 = 𝑄! ×𝑃! ×𝑇
101325×(273.15 + 300)
= 1000000×
𝑃×𝑇!
(4×10! + 101325)×(273.15)= 1000000×

101325×573.15
4101325×273.15= 5183.94 m3/h
配管内径：φ574.6 とした場合の流速は下式
𝑄5183.94𝑣= =𝜋𝐴
×0.5746!
4= 19991.2 m/h

A:配管内径断面積（m2）= 5.55 m/s
容積流量を与えられた場合は、流体密度が不明でも流速は計算可能。

-14-