FE Mechanical 試験対策:直前対策パック

公式集、全トピックの解法の要点、42 問の演習問題。完全な解答は巻末にあります。

新しい問題科目に戻る

印刷ウィンドウで「PDF として保存」を選ぶと、コピーを端末に保存できます。

1公式と解法

はりの支点反力

つり合い
ΣFx = 0, ΣFy = 0, ΣM = 0
等分布荷重
合力 w·L は荷重の中央に作用
モーメント
M = F·d

解き方

  1. 自由体図を描き、各支点を反力に置き換えます。
  2. 分布荷重はそれぞれ、その図心に作用する 1 つの力に置き換えます。
  3. 一方の支点まわりのモーメントをとり、もう一方の支点反力を求めます。
  4. ΣFy = 0 を使って残りの反力を求め、もう 1 つのモーメントの式で確かめます。

軸応力とひずみ

垂直応力
σ = F/A
ひずみ
ε = ΔL/L
フックの法則
σ = E·ε
伸び
ΔL = FL/(EA)
熱ひずみ
εT = α·ΔT

軸のねじり

せん断応力
τ = T·r/J
中実軸
J = πd⁴/32
中空軸
J = π(D⁴ − d⁴)/32
ねじれ角
φ = TL/(GJ)
電力
P = T·ω = 2πnT/60

解き方

  1. 断面二次極モーメントを求めます:中実軸では J = πd⁴/32。
  2. 表面のせん断応力:τ = T·r/J。r は外半径です。
  3. ねじれ角:φ = TL/(GJ)(ラジアン)。
  4. 動力が与えられている場合は、まず P = T·ω、ω = 2πn/60 からトルクを求めます。

せん断力と曲げモーメント

せん断力と曲げモーメント
V = 片側の ΣF、 M = 片側の ΣF·d
曲げ応力
σ = M·y/I
長方形
I = bh³/12
断面係数
W = I/ymax
単純支持、等分布荷重
Mmax = wL²/8

解き方

  1. はり全体のつり合いから支点反力を求めます。
  2. 必要な断面ではりを切断し、片側だけを考えます。
  3. せん断力:その側に働く鉛直方向の力の和。曲げモーメント:それらの力の切断面まわりのモーメントの和。
  4. 最大曲げモーメントは、せん断力がゼロになる位置を求めて調べます。

モールの応力円と応力の座標変換

中心
C = (σx + σy)/2
半径
R = √(((σx − σy)/2)² + τxy²)
主応力
σ1,2 = C ± R
主応力の方向
tan 2θp = 2τxy/(σx − σy)
斜め断面
σθ = C + (σx − σy)/2·cos 2θ + τxy·sin 2θ

解き方

  1. モールの応力円の中心を計算します:C = (σx + σy)/2。
  2. 半径を計算します:R = √(((σx − σy)/2)² + τxy²)。
  3. 主応力は C ± R、面内の最大せん断応力は R です。
  4. 円上の角度は要素上の角度の 2 倍です:tan 2θp = 2τxy/(σx − σy)。

はりのたわみ

支間中央の集中荷重
δ = PL³/(48EI)
等分布荷重、単純支持
δ = 5wL⁴/(384EI)
片持ちばり、先端荷重
δ = PL³/(3EI)
片持ちばり、等分布荷重
δ = wL⁴/(8EI)
重ね合わせの理
δ = δ1 + δ2

柱の座屈

オイラー座屈荷重
Pcr = π²EI/(KL)²
有効座屈長さ係数
K = 1 ピン–ピン、0.5 固定–固定、0.7 固定–ピン、2 固定–自由
座屈応力
σcr = Pcr/A
細長比
λ = KL/r, r = √(I/A)

放物運動

成分
vx = v0cos θ, vy = v0sin θ
位置
x = vxt, y = h + vyt − ½gt²
最高点の高さ
H = h + vy²/(2g)
水平到達距離(水平な地面)
R = v0²·sin 2θ/g

解き方

  1. 初速度を成分に分けます:vx = v₀·cos θ、vy = v₀·sin θ。
  2. 水平方向は等速運動、鉛直方向は加速度 −g の運動です。
  3. 鉛直方向の運動から飛行時間を求めます。
  4. その時間を水平方向の運動に使い、水平到達距離を求めます。

仕事とエネルギー

運動エネルギー
Ek = ½mv²
位置エネルギー
Ep = mgh
エネルギー収支
Ek1 + Ep1 + W = Ek2 + Ep2 + losses
摩擦力がする仕事
W = μ·N·d
電力
P = W/t

機械振動

固有振動数
ωn = √(k/m), fn = ωn/(2π)
静たわみから
ωn = √(g/δst)
臨界減衰
cc = 2√(km) = 2mωn
減衰比
ζ = c/cc
減衰固有振動数
ωd = ωn√(1 − ζ²)
対数減衰率
δ = (1/n)·ln(x0/xn), ζ = δ/√(4π² + δ²)
並列ばね
k = k1 + k2
直列ばね
1/k = 1/k1 + 1/k2

熱機関と冷凍機

効率
η = W/QH = 1 − QC/QH
カルノー効率
ηC = 1 − TC/TH (kelvin)
冷凍機の COP
COPR = QC/W
ヒートポンプの COP
COPHP = QH/W = COPR + 1
エネルギー収支
QH = W + QC

ガスサイクル

等エントロピー変化(体積比)
T2/T1 = rk−1
等エントロピー変化(圧力比)
T2/T1 = rp(k−1)/k
オットーサイクルの効率
η = 1 − 1/rk−1
ディーゼルサイクルの効率
η = 1 − (1/rk−1)·(rck − 1)/(k(rc − 1))
ブレイトンサイクルの効率
η = 1 − 1/rp(k−1)/k
加熱量
オットー:cv(T3 − T2) ディーゼル、ブレイトン:cp(T3 − T2)
ブレイトンサイクルの正味仕事
w = cp[(T3 − T4) − (T2 − T1)]
断熱効率
wc = wc,s/ηc, wt = ηt·wt,s
冷空気標準
k = 1.4, cp = 1.005, cv = 0.718 kJ/(kg·K)

連続の式とベルヌーイ

連続の式
A1v1 = A2v2
ベルヌーイの定理
p1 + ½ρv1² + ρgz1 = p2 + ½ρv2² + ρgz2
流量
Q = A·v
円の面積
A = πd²/4

管摩擦と損失水頭

レイノルズ数
Re = V·D/ν = ρ·V·D/μ
流れの状態
Re ≈ 2000 未満で層流、Re ≈ 4000 を超えると乱流
平均流速
V = Q/A, A = π·D²/4
ダルシー・ワイスバッハの式
hf = f·(L/D)·V²/(2g)
局所損失
hm = ΣK·V²/(2g)
代表的な K の値
入口 0.5、出口 1.0、90° 曲がり ≈ 0.3
圧力損失
Δp = ρ·g·hf
2 つの開放貯水槽
Δz = (f·L/D + ΣK)·V²/(2g)
自由放出
H = (1 + f·L/D + ΣK)·V²/(2g)

ポンプとポンプ動力

水動力
Ph = ρ·g·Q·H = Q·Δp
軸動力
Pshaft = Ph/ηpump
電力
Pel = Pshaft/ηmotor
ポンプ揚程
H = Δz + Δp/(ρ·g) + hL
管路系の損失
hL = (f·L/D + ΣK)·V²/(2g)

熱伝導

フーリエの法則
Q̇ = kAΔT/L
熱伝導抵抗
R = L/(kA)
対流熱抵抗
R = 1/(hA)
直列の熱抵抗
Q̇ = ΔTtotal/ΣR
U 値
U = 1/(ΣR·A)

歯車と歯車列

かみ合う歯車
n1·N1 = n2·N2
トルク
T2 = T1·(N2/N1)·η
電力
P = T·ω, ω = 2π·n/60
複合歯車列
全体の速比 = 各段の速比の積

定常偏差と PID 制御

システムの型
G(s) の s = 0 にある極の個数
偏差定数
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
ステップ A
ess = A/(1 + Kp)
ランプ A·t
ess = A/Kv
放物線 A·t²/2
ess = A/Ka
PID 制御器
C(s) = Kp + Ki/s + Kd·s
2 次系の目標
s² + 2ζωns + ωn²
オーバーシュートと時間
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

解き方

  1. 開ループ C(s)·G(s) の s = 0 にある極を数えます:それがシステムの型です。
  2. 入力に対応する偏差定数を求めます:ステップには Kp、ランプには Kv、放物線には Ka です。
  3. ess = A/(1 + Kp)、A/Kv、または A/Ka を使うか、それを必要なゲインについて解きます。
  4. 減衰や速度の仕様では、閉ループの特性多項式を書き、s² + 2ζω_n·s + ω_n² と係数を合わせます。

2演習問題

1.はりの支点反力ふつう
w = 12 kN/mP = 45 kNABMA​ = ?4 m6 m

片持ちばりは A で壁に固定されている。壁から 4 m まで等分布荷重 w が作用し、壁から 10 m の位置に集中荷重 P が作用する。 壁における反力モーメント MA の大きさを求めよ。

L = 10 mP = 45 kN 位置 10 mw = 12 kN/m(0 から 4 m まで)
MA = kN·m
2.はりの支点反力むずかしい
P = 10 kNQ = 40 kNADCRD​ = ?1.5 m0.5 m3.5 m0.5 m

複合ばりは A で壁に固定され、C(A から 2 m)に中間ヒンジをもち、D(A から 6 m)でローラーに支持されている。荷重は A から 5.5 m の位置の集中荷重 P および A から 1.5 m の位置の集中荷重 Q である。 ローラーの反力 RD を求めよ。

AC = 2 mAD = 6 mP = 10 kN 位置 5.5 mQ = 40 kN 位置 1.5 m
RD = kN
3.軸応力とひずみふつう
Ø 12 mm銅3 mP = 90 kN

軸荷重を受ける銅製の棒の伸び δ を求めよ。

P = 90 kNd = 12 mmL = 3 mE = 117 GPa
δ = mm
4.軸応力とひずみむずかしい
P = ?Ø25110 GPaØ80/60600 mm

鋼 の棒(直径 25 mm、E = 200 GPa、許容応力 150 MPa)が、青銅 でできた管(外径 80 mm、内径 60 mm、E = 110 GPa、許容応力 80 MPa)の中に入っています。両者は剛な板の間で一緒に圧縮されます。最大荷重 Pmax を求めなさい。

棒: Ø25 mm, E = 200 GPa, σallow = 150 MPa管: Ø80/Ø60 mm, E = 110 GPa, σallow = 80 MPa
Pmax = kN
5.軸のねじりふつう
d1​ = 75 mmd2​ = 60 mmT = 250 N·mABC2 m2 m

黄銅製の段付き中実軸が A で固定されている。区間 AB は直径 75 mm、長さ 2 m、区間 BC は直径 60 mm、長さ 2 m である。自由端 C にトルク 250 N·m が作用する。 軸の最大せん断応力 τmax を求めよ。

T = 250 N·md1 = 75 mm, L1 = 2 md2 = 60 mm, L2 = 2 m
τmax = MPa
6.軸のねじりむずかしい
Tdo​ = 50 mm, di​ = 35 mm

中空軸が 300 rpm で 15 kW を伝達する。軸の最大せん断応力 τmax を求めよ。

P = 15 kWN = 300 rpmdo = 50 mmdi = 35 mm
τmax = MPa
7.せん断力と曲げモーメントふつう
50 kN60 kNABCDEFG2 m2 m2 m2 m2 m2 m

点 A から G は、梁を 6 つの等しい部分に分けています。曲げモーメントの大きさが最大になるのは、正曲げでも負曲げでも、どの点ですか? 図でそれをクリックしなさい。

50 kN @ 10 m60 kN @ 6 m
  1. 1A (x = 0 m)
  2. 2B (x = 2 m)
  3. 3C (x = 4 m)
  4. 4D (x = 6 m)
  5. 5E (x = 8 m)
  6. 6F (x = 10 m)
  7. 7G (x = 12 m)
8.せん断力と曲げモーメントむずかしい
w = 10 kN/mP = 45 kNAB断面1 m1.5 m5.5 m

支点 A から x = 2.5 m の断面における曲げモーメント Mx を求めよ。下側が引張となる曲げモーメントを正、左側断面で上向きのせん断力を正とする。

L = 8 mP = 45 kN 位置 1 mw = 10 kN/m(0 から 8 m まで)x = 2.5 m
Mx = kN·m
9.モールの応力円と応力の座標変換ふつう
σx​ = 160σy​ = −40τxy​ = 30MPa

荷重を受ける部品内のある点で、σₓ = 160 MPa、σᵧ = −40 MPa、τₓᵧ = 30 MPa である(引張を正とする)。 x 軸から最大主応力の方向までの角度 θp を、−90° から 90° の範囲で求めよ。

σx = 160 MPaσy = −40 MPaτxy = 30 MPa
θp = °
10.モールの応力円と応力の座標変換むずかしい
σx​τxy​軸の表面

直径 30 mm の中実軸が軸圧縮力 140 kNとトルク 400 N·m を同時に受ける。 軸の表面で、 最大主応力 σ1 を求めよ。

d = 30 mmF = −140 kNT = 400 N·m
σ1 = MPa
11.はりのたわみふつう
w = 6 kN/mABθ3 m

単純ばりの全スパンに等分布荷重 w が作用する。支点 A におけるはりのたわみ角 θA をラジアンで求めよ。

w = 6 kN/mL = 3 mE = 11 GPa (構造用木材)I = 120 × 10⁶ mm⁴
θA = rad
12.はりのたわみむずかしい
w = 15 kN/mAB3 m

一端固定・他端支持ばりは A で固定され、B で単純支持されている。全長に等分布荷重 w が作用する。EI は一定である。 支点 B の反力 RB を求めよ。

L = 3 mw = 15 kN/m
RB = kN
13.柱の座屈ふつう
PL = 2.5 m40 × 50 mm断面

アルミニウム製の長さ L の柱が下端固定・上端自由である。オイラーの式を用いて座屈荷重 Pcr を求めよ。

L = 2.5 mE = 70 GPa (アルミニウム)b × h = 40 × 50 mm
Pcr = kN
14.柱の座屈むずかしい
P = 80 kNL = 5 md = ?断面

鋼製の長さ L の中実円柱が下端固定・上端ピンである。座屈に対する安全率 2 で 80 kN を支える必要がある。最小直径 d を求めよ。

P = 80 kNn = 2L = 5 mE = 200 GPa (鋼)
d = mm
15.放物運動ふつう
v0​ = 38 m/s65°到達距離

ボールを 平らな地面から水平から 65° 上向きに投げる。空気抵抗は無視する。ボールが地面に落ちるまでの飛行時間 t を求めよ。

v0 = 38 m/sθ = 65°g = 9.81 m/s²
t = s
16.放物運動むずかしい
h = 35 mv0​ = 35 m/s25°x = ?

ボールを 高さ 35 m の崖の上から水平から 25° 上向きに投げる。空気抵抗は無視する。ボールが地面に落ちるまでに進む水平距離 x を求めよ。

v0 = 35 m/sθ = 25°h = 35 mg = 9.81 m/s²
x = m
17.仕事とエネルギーふつう
25°d = 11 mμ = 0.15

質量 10 kg のブロックが水平と 25° をなす斜面を 11 m すべり下りる。動摩擦係数は 0.15 である。摩擦に逆らってした仕事 Wf を求めよ。

m = 10 kgθ = 25°μ = 0.15d = 11 mg = 9.81 m/s²
Wf = J
18.仕事とエネルギーむずかしい
k = 800 N/m15°μ = 0.35 kg

15° の斜面の下端にあるばね定数 800 N/m のばねを、質量 5 kg のブロックに押し当てて 20 cm 縮めてから放す。ブロックと斜面の間の動摩擦係数は 0.3 である。ブロックは放した点から斜面をどれだけの距離 d 上がって止まるか。(ばねが押す短い距離は無視し、摩擦と重力は全距離 d にわたって作用するとする。)

m = 5 kgk = 800 N/mx = 20 cmθ = 15°μ = 0.3g = 9.81 m/s²
d = m
19.機械振動ふつう
tx40 mm1 周期29.2 mm

減衰のある質量–ばね系の自由振動試験で、変位のピーク 40 mm の 1 周期後に 29.2 mm のピークが現れた。 減衰比 ζ を求めよ。

x0 = 40 mmx1 = 29.2 mmn = 1
ζ =
20.機械振動むずかしい
k1​ = 12k2​ = 60m = 25 kg

25 kg の質量が、縦につないだ(直列の)2 本のばねにつり下げられている:上のばねは k₁ = 12 kN/m、下のばねは k₂ = 60 kN/m である。 減衰とばねの質量は無視する。上下振動の固有振動数 fn を Hz で求めよ。

m = 25 kgk1 = 12 kN/mk2 = 60 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 400 °CTL​ = 10 °CEQH​QL​η = 29%

熱機関が 400 °C の熱源から熱を受け取り、10 °C の周囲へ熱を放出する。実際の熱効率は 29% である。第二法則効率(実際の効率とカルノー効率の比)ηII を求めよ。

TH = 400 °CTL = 10 °Cη = 29%
ηII = %
22.熱機関と冷凍機むずかしい
TH​ = 550 °CTL​ = 35 °CEQH​ = 3200 kJQL​ = 1600 kJW = 1600 kJ

熱機関が 1 サイクルあたり 550 °C の熱源から 3200 kJ の熱を受け取り、1600 kJ の仕事をして、残りを 35 °C の熱源へ放出する。1 サイクルあたりの全エントロピー生成量 Sgen を求めよ。

TH = 550 °CTL = 35 °CQH = 3200 kJW = 1600 kJ
Sgen = kJ/K
23.ガスサイクルふつう
vP1234qin​qout​

理想ディーゼルサイクルの圧縮比は r = 18、締切比は rc = v₃/v₂ = 3 である。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 理論熱効率 η をパーセントで求めよ。

r = 18rc = 3k = 1.4
η = %
24.ガスサイクルむずかしい
vP1234qin​qout​

理想ブレイトンサイクルが圧縮機の圧力比 16 で運転されている。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 理論熱効率 η をパーセントで求めよ。

rp = 16k = 1.4
η = %
25.連続の式とベルヌーイふつう
v = ?p0​pp0​ − p = 37 kPaよどみ点

水道管内のピトー静圧管が、よどみ点圧と静圧の差として 37 kPa を示す。先端での水の流速 v を求めよ。

p0 − p = 37 kPaρ = 1000 kg/m³
v = m/s
26.連続の式とベルヌーイむずかしい
d1​ = 150 mmd2​ = 75 mm12Δh = 50 mm水銀Q = ?

水平なベンチュリ管(入口 150 mm、のど部 75 mm)を水が流れる。入口とのど部の間に接続した水銀 U 字管マノメータの液面差は 50 mm である。流量係数は 0.97 である。流量 Q を求めよ。水は ρ = 1000 kg/m³、水銀は 13 600 kg/m³、g = 9.81 m/s² とする。

d1 = 150 mmd2 = 75 mmΔh = 50 mmCd = 0.97ρHg = 13 600 kg/m³g = 9.81 m/s²
Q = L/s
27.管摩擦と損失水頭ふつう
V = 1.6 m/sD = 150 mmL = 50 mf = 0.03

水(ρ = 1000 kg/m³)が、長さ 50 m、内径 150 mm の水平管内を 1.6 m/s で流れている。ダルシーの管摩擦係数 f = 0.03、g = 9.81 m/s² とする。管での圧力損失 Δp を kPa で求めよ。

V = 1.6 m/sL = 50 mD = 150 mmf = 0.03ρ = 1000 kg/m³
Δp = kPa
28.管摩擦と損失水頭むずかしい
K = 5Δz = 28 mL = 100 m, D = 150 mm, f = 0.0254× 曲がり (K = 0.3)

2 つの開放貯水槽が、長さ 100 m、内径 150 mm、ダルシーの管摩擦係数 f = 0.025 の管でつながれている。上の貯水槽の水面は下の貯水槽の水面より 28 m 高い。局所損失は 角ばった入口(K = 0.5)、4 個の曲がり(各 K = 0.3)、弁(K = 5)、下の貯水槽への出口(K = 1.0) である。g = 9.81 m/s² とする。流量 Q を L/s で求めよ。

Δz = 28 mL = 100 mD = 150 mmf = 0.025Kentrance = 0.5Kbend = 0.3 (×4)Kvalve = 5Kexit = 1.0
Q = L/s
29.ポンプとポンプ動力ふつう
PQ = 10 L/sH = 12 m

ポンプが揚程 12 m に対して水(ρ = 1000 kg/m³)を 10 L/s 送っている。ポンプ効率は 75%、これを駆動する電動機の効率は 95% である。g = 9.81 m/s² とする。電動機が消費する電力 Pel を kW で求めよ。

Q = 10 L/sH = 12 mηpump = 75%ηmotor = 95%g = 9.81 m/s²
Pel = kW
30.ポンプとポンプ動力むずかしい
PQ = 40 L/sΔz = 17 mL = 500 m, D = 150 mm, f = 0.025, ΣK = 8

ポンプが 40 L/s の水(ρ = 1000 kg/m³)を、開放された下の貯水槽から、水面が 17 m 高い上のタンクへ送る。管は長さ 500 m、内径 150 mm、ダルシーの管摩擦係数 f = 0.025 で、局所損失係数の合計は ΣK = 8(入口と出口を含む)である。g = 9.81 m/s² とする。 ポンプ効率は 78% である。軸動力 Pshaft を kW で求めよ。

Q = 40 L/sΔz = 17 mL = 500 mD = 150 mmf = 0.025ΣK = 8η = 78%
Pshaft = kW
31.熱伝導ふつう
ロックウール, k = 0.04150 mmコンクリート, k = 1.475 mm15 °C-1 °CQ̇ = ?

壁はロックウール 150 mm とコンクリート 75 mm からなり、面積は 12 m² である。外側の両面は 15 °C と -1 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 150 mm, k1 = 0.04 W/(m·K)L2 = 75 mm, k2 = 1.4 W/(m·K)A = 12 m²T1 = 15 °CT3 = -1 °C
Q̇ = W
32.熱伝導むずかしい
ロックウール, k = 0.04200 mm木材, k = 0.1725 mm室内 22 °C屋外 -3 °CQ̇

200 mm のロックウールと 25 mm の木材からなる壁が、22 °C の部屋と -3 °C の外気を隔てている。熱伝達率は室内側 7.7 W/(m²·K)、屋外側 20 W/(m²·K)、壁の面積は 6 m² である。 壁の熱通過率 U を求めよ。

L1 = 200 mm, k1 = 0.04 W/(m·K)L2 = 25 mm, k2 = 0.17 W/(m·K)A = 6 m²hi = 7.7 W/(m²·K)ho = 20 W/(m²·K)Troom = 22 °CTout = -3 °C
U = W/(m²·K)
33.歯車と歯車列ふつう
N1​ = 32N2​ = 72

1450 rpm で回転する 11 kW のモータが歯数 32 のピニオンを駆動し、ピニオンは歯数 72 の歯車とかみ合っている。 この歯車対の効率は 97% である。 歯車軸の出力トルク T2 を求めよ。

P = 11 kWn1 = 1450 rpmN1 = 32N2 = 72η = 97%
T2 = N·m
34.歯車と歯車列むずかしい
d1​ = 80 mmd2​ = 450 mm960 rpmC = 1200 mmT1​ = 1000 NT2​P = ?

平行掛けのベルトが、80 rpm の直径 960 mm の駆動プーリと直径 450 mm のプーリをつなぎ、中心は 1200 mm 離れています。摩擦係数は 0.3 で、張り側の張力は 1000 N を超えてはなりません。ベルトが滑らずに伝達できる最大動力 P を kW で求めなさい(遠心力の影響は無視します)。

d1 = 80 mmn1 = 960 rpmd2 = 450 mmC = 1200 mmμ = 0.3T1 = 1000 N
P = kW
35.ラプラス変換ふつう
ReIm−1−2−3−4−5−6−7−8−91

F(s) = 34/(s² + 10s + 16) である。その逆変換は f(t) = A·e−2t + B·e−8t の形になる。部分分数展開により係数 A を求めよ。

F(s) = 34/(s² + 10s + 16)
A =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.2y(0.2) = ?

ある系が y″ + 7y′ + 6y = 12·u(t) に従い、y(0) = 5、y′(0) = 0 である。t = 0 でステップ入力 u(t) = 1 を加える。ラプラス変換を用いて y(0.2) を求めよ。

y″ + 7y′ + 6y = 12·u(t)y(0) = 5y′(0) = 0u(t) = 1t = 0.2 s
y(0.2) =
37.ブロック線図ふつう
R(s)+Ks(s + 20)G(s)C(s)−

プラント G(s) = K/(s(s + 20)) を単位負帰還ループに入れる。閉ループが臨界減衰(等しい 2 つの実数極)になるゲイン K を求めよ。

G(s) = K/(s(s + 20))単位負帰還
K =
38.ブロック線図むずかしい
R(s)+12s + 5G₁+5s + 7G₂C(s)2H₁−1H₂−

G₁ = 12/(s + 5)、G₂ = 5/(s + 7)、G₂ の周りの内側の負帰還 H₁ = 2、外側の負帰還 H₂ = 1 のブロック線図を簡単にせよ。C(s)/R(s) の直流ゲイン T(0) を求めよ。

G₁ = 12/(s + 5)G₂ = 5/(s + 7)H₁ = 2H₂ = 1
T(0) =
39.ステップ応答ふつう
t (s)y(t)1

ある系の伝達関数は G(s) = 9/(s² + 2s + 9) で、単位ステップ入力を受ける。 減衰比 ζ を求めよ。

G(s) = 9/(s² + 2s + 9)
ζ =
40.ステップ応答むずかしい
R(s)+Ks(s + 5)G(s)C(s)−

単位負帰還ループの G(s) = K/(s(s + 5)) である。ステップ応答のオーバーシュートが 25 % になるゲイン K を求めよ。

G(s) = K/(s(s + 5))%OS = 25 %
K =
41.定常偏差と PID 制御ふつう
R(s)+K(s + 1)(s + 4)C(s)−

単位フィードバックのループで G(s) = K/((s + 1)(s + 4)) です。ステップ入力に対する定常偏差が 10 % となるゲイン K を求めなさい。

G(s) = K/((s + 1)(s + 4))ess = 10 %
K =
42.定常偏差と PID 制御むずかしい
r(t) = 5t+2 + 10/s4(s + 1)(s + 10)C(s)−

PI 制御器 C(s) = 2 + 10/s が、単位フィードバックのループで制御対象 G(s) = 4/((s + 1)(s + 10)) を制御しています。ランプ入力 r(t) = 5t に対する定常偏差 ess を求めなさい。

C(s) = 2 + 10/sG(s) = 4/((s + 1)(s + 10))r(t) = 5t
ess =

3溶液

1.はりの支点反力

答え: MA = 546 kN·m

片持ちばりのつり合い:固定支点は上向きに押すことも回転に抵抗することもできるので、鉛直反力 RA と反力モーメント MA を与えます。ΣFy = 0 から RA が、壁まわりのモーメントから MA = Σ(荷重 × 壁からの距離) が求まります。分布荷重は、荷重図の図心に作用する合力として働きます。

分布荷重の合力

W = 12·4 = 48 kN, x̄ = 2 m

壁まわりのモーメント

MA = P·a + W·x̄ = 45·10 + 48·2 = 546 kN·m

答え

MA = 546 kN·m
2.はりの支点反力

答え: RD = 8.75 kN

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。各分布荷重を、荷重図の図心に作用する合力(その面積)に置き換え、支点まわりのモーメントをとってその支点の反力を式から消去します。 中間ヒンジは曲げモーメントを伝えられないので、追加の式が 1 つ得られます:ヒンジのどちら側についても、その側の力のヒンジまわりのモーメントは 0 です。

部分 CD

ヒンジにはモーメントがないので、CD だけについての C まわりのモーメントがつり合わなければなりません。

ΣMC = 0: RD·4 = 10·3.5
RD = 8.75 kN
VC = 10 − 8.75 = 1.25 kN

答え

RD = 8.75 kN
3.軸応力とひずみ

答え: δ = 20.4 mm

フックの法則:弾性範囲では σ = E·ε。σ = P/A と ε = δ/L を組み合わせると δ = PL/(AE)。

応力

A = π·12²/4 = 113.1 mm²
σ = P/A = 90·1000/113.1 = 795.8 MPa

ひずみ(フックの法則)

E for 銅 = 117 GPa = 117000 MPa.

ε = σ/E = 795.8/117000 = 0.006801 = 6801 µε

伸び

δ = ε·L = 0.006801·3000 mm = 20.4 mm

答え

δ = 20.4 mm
4.軸応力とひずみ

答え: Pmax = 247.3 kN

2 種類の材料:直列(端から端へ)では、どの区間も同じ力を受け、伸びが足されます:δ = P·Σ L/(E·A)。並列(剛な板の間に並べた場合)では、両方が同じ量だけ縮むので、ひずみが等しくなります:剛性の高いほうがより大きな割合を受け持ち、Pi = P·EiAi/ΣEA、σi = Ei·ε です。2 つの剛な壁の間に固定され、接合部に荷重がかかる棒は不静定です:2 つの区間は、剛性 k = EA/L に比例して荷重を分担します。

面積

Ar = π·25²/4 = 490.9 mm²
At = π(80² − 60²)/4 = 2199 mm²

合成した軸剛性

ΣEA = 200·10³·490.9 + 110·10³·2199 = 3.401 × 10⁸ N

支配するひずみ

ε = min(0.00075, 0.0007273) = 0.0007273 (青銅)

最大荷重

Pmax = 0.0007273·3.401 × 10⁸ = 247329 N = 247.3 kN

答え

Pmax = 247.3 kN
5.軸のねじり

答え: τmax = 5.895 MPa

丸軸のねじり:せん断応力は中心で 0、表面で τmax = T·r/J まで直線的に増えます。中実軸は J = πd⁴/32、中空軸は J = π(do⁴ − di⁴)/32。ねじれ角は φ = TL/(GJ)(ラジアン)。動力とトルクは P = T·ω で結ばれます。

各区間の応力

T = 250 N·m = 250000 N·mm
τAB = 16·250000/(π·75³) = 3.018 MPa
τBC = 16·250000/(π·60³) = 5.895 MPa

答え

τmax = 5.895 MPa
6.軸のねじり

答え: τmax = 25.6 MPa

丸軸のねじり:せん断応力は中心で 0、表面で τmax = T·r/J まで直線的に増えます。中実軸は J = πd⁴/32、中空軸は J = π(do⁴ − di⁴)/32。ねじれ角は φ = TL/(GJ)(ラジアン)。動力とトルクは P = T·ω で結ばれます。

トルク

ω = 2πN/60 = 2π·300/60 = 31.42 rad/s
T = P/ω = 15000/31.42 = 477.5 N·m

断面二次極モーメント

J = π(50⁴ − 35⁴)/32 = 466269 mm⁴

最大せん断応力

τmax = T·(do/2)/J = 477465·25/466269 = 25.6 MPa

答え

τmax = 25.6 MPa
7.せん断力と曲げモーメント

答え: D (x = 6 m)

最大曲げモーメント:せん断力が変わるところでは M の傾きが変わるので、そのピークはせん断力がゼロを横切る位置にあります:集中荷重の下、分布荷重の内側、または張り出し梁の支点の上です。片持ち梁の固定端では、モーメントは壁で最大です。それらの位置で M の大きさを、正曲げでも負曲げでも比べます。

支点反力

一方の支点まわりのモーメントをとり、次に鉛直方向の力の釣り合いを使います。

R1 = 38.33 kN, R2 = 71.67 kN

各点での曲げモーメント

切断面の一方の側にあるすべての力のモーメントを合計します。正曲げ(下に凸)を正、負曲げ(上に凸)を負とします。

MA = 0, MB = 76.67, MC = 153.3, MD = 230 kN·m, ME = 186.7, MF = 143.3, MG = 0

答え

点 D: MD = 230 kN·m
8.せん断力と曲げモーメント

答え: Mx = 99.69 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 10·8 = 80 kN、作用点 x̄ = 4 m
ΣMA = 0: RB·8 = 45·1 + 10·8·4
RB = 365/8 = 45.63 kN
ΣFy = 0: RA = 125 − 45.63 = 79.38 kN
w = 10 kN/mP = 45 kNABRA​ = 79.38RB​ = 45.63断面1 m1.5 m5.5 m

x = 2.5 m で切断し、左側を残す

切断面より左の力だけを考えます。分布荷重の一部は、その部分の中央に働く 1 つの力として扱います。

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 79.38·2.5 − 45·(2.5 − 1) − 10·2.5·1.25 = 99.69 kN·m

答え

Mx = 99.69 kN·m
9.モールの応力円と応力の座標変換

答え: θp = 8.35 °

平面応力の座標変換:θ だけ回転した面では σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。モールの応力円はすべての面を一度に表します:中心は C = (σₓ + σᵧ)/2、半径は R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主応力は C ± R、面内の最大せん断応力は R です。

στxyCσ1​σ2​

円上の角度

モールの応力円上で、点 x は σ 軸から角度 2θp の位置にあります。象限は σₓ − σᵧ と τₓᵧ の符号で決まります。

2θp = atan2(2·30, 160 − −40) = 16.7°

要素上の角度

θp = 8.35°

答え

θp = 8.35°

反時計回りの角度を正とします。

10.モールの応力円と応力の座標変換

答え: σ1 = 25.47 MPa

平面応力の座標変換:θ だけ回転した面では σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。モールの応力円はすべての面を一度に表します:中心は C = (σₓ + σᵧ)/2、半径は R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主応力は C ± R、面内の最大せん断応力は R です。

表面の応力

σ = F/A = −140000/(π·30²/4) = −140000/706.9 = −198.1 MPa
τ = 16T/(πd³) = 16·400000/(π·30³) = 75.45 MPa
στxyCσ1​σ2​

モールの応力円

C = σ/2 = −99.03 MPa
R = √((−99.03)² + 75.45²) = 124.5 MPa
σ1 = C + R = 25.47 MPa
σ2 = C − R = −223.5 MPa

答え

σ1 = 25.47 MPa
11.はりのたわみ

答え: θA = 0.005114 rad

はりのたわみ:たわみは曲率 M/(EI) を 2 回積分して得られます。代表的な場合は表にまとめられており、たとえば先端集中荷重の片持ちばりは PL³/(3EI)、全長等分布荷重の単純ばりは 5wL⁴/(384EI) です。別々の荷重によるたわみは足し合わせられます(重ね合わせ)。 たわみ角 θ = dv/dx は M/(EI) の 1 回目の積分です。単位をそろえればラジアンで得られます。

単位をそろえる

N と mm で計算します:P は N、w は N/mm (kN/m = N/mm)、L は mm、E は N/mm² (GPa × 1000)、I は mm⁴。

w = 6 N/mm, L = 3000 mm
EI = 1.32 × 10¹² N·mm²

端部のたわみ角

θA = wL³/(24EI) = 6·3000³/(24·1.32 × 10¹²) = 0.005114 rad

答え

θA = 0.005114 rad

これは 0.293° です。

12.はりのたわみ

答え: RB = 16.88 kN

はりのたわみ:たわみは曲率 M/(EI) を 2 回積分して得られます。代表的な場合は表にまとめられており、たとえば先端集中荷重の片持ちばりは PL³/(3EI)、全長等分布荷重の単純ばりは 5wL⁴/(384EI) です。別々の荷重によるたわみは足し合わせられます(重ね合わせ)。 不静定ばりでは、余分な支点を 1 つ取り除き、そこでのたわみが 0 になる条件(適合条件)を課します。

支点 B を取り除く

支点 B がなければ、はりは片持ちばりです。その自由端は荷重によってたわみます。

δload = wL⁴/(8EI)

支点 B での適合条件

支点反力 RB だけで自由端は RBL³/(3EI) 持ち上がります。支点 B はたわみを 0 に保つので、両者は打ち消し合わなければなりません。EI は消去されます。

RB·L³/(3EI) = wL⁴/(8EI) → RB = 3wL/8 = 3·15·3/8 = 16.88 kN

答え

RB = 16.88 kN
13.柱の座屈

答え: Pcr = 7.369 kN

オイラー座屈:細長い柱は Pcr = π²EI/(KL)² で不安定になります。KL は有効座屈長さで、K = 1(ピン–ピン)、2(固定–自由)、0.7(固定–ピン)、0.5(固定–固定)です。柱は I が最小の軸まわりに座屈します。この式は、座屈応力が降伏応力より小さい細長い柱に適用できます。

最小の断面二次モーメント

長方形は弱軸まわりに座屈するので、小さい方の寸法の 3 乗を使います。

Imin = h·b³/12 = 50·40³/12 = 266667 mm⁴

有効座屈長さ

KL = 2·2500 = 5000 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·266667/5000² = 7369 N = 7.369 kN

答え

Pcr = 7.369 kN
14.柱の座屈

答え: d = 67.06 mm

オイラー座屈:細長い柱は Pcr = π²EI/(KL)² で不安定になります。KL は有効座屈長さで、K = 1(ピン–ピン)、2(固定–自由)、0.7(固定–ピン)、0.5(固定–固定)です。柱は I が最小の軸まわりに座屈します。この式は、座屈応力が降伏応力より小さい細長い柱に適用できます。

必要な座屈荷重

Pcr = n·P = 2·80000 = 160000 N
KL = 0.7·5000 = 3500 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·160000·3500²/(π³·200000))^(1/4) = 67.06 mm

答え

d = 67.06 mm

確認:細長比 KL/r は 208.8 なので柱は細長く、オイラーの式が使えます。

15.放物運動

答え: t = 7.021 s

放物運動:空気抵抗がなければ x = vxt、y = h + vyt − ½gt²。2 つの方向は独立で、共通なのは時間だけです。

初速度を分解する

水平速度は一定で、鉛直運動は下向きの一定加速度 g = 9.81 m/s² をもちます。

vx = v0·cos θ = 38·cos 65° = 16.06 m/s
vy = v0·sin θ = 38·sin 65° = 34.44 m/s

滞空時間

発射点と着地点が同じ高さ:t = 2vy/g。

t = 2·34.44/9.81 = 7.021 s

答え

t = 7.021 s
16.放物運動

答え: x = 145.1 m

放物運動:空気抵抗がなければ x = vxt、y = h + vyt − ½gt²。2 つの方向は独立で、共通なのは時間だけです。

初速度を分解する

水平速度は一定で、鉛直運動は下向きの一定加速度 g = 9.81 m/s² をもちます。

vx = v0·cos θ = 35·cos 25° = 31.72 m/s
vy = v0·sin θ = 35·sin 25° = 14.79 m/s

滞空時間

y(t) = h + vy·t − ½gt² = 0 を解き、正の解をとります。

0 = 35 + 14.79·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (14.79 + √(218.8 + 686.7))/9.81 = 4.575 s

水平距離

x = vx·t = 31.72·4.575 = 145.1 m

答え

x = 145.1 m
17.仕事とエネルギー

答え: Wf = 146.7 J

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。

摩擦力

N = mg·cos θ = 98.1·0.9063 = 88.91 N
F = μN = 0.15·88.91 = 13.34 N

摩擦に逆らう仕事

Wf = F·d = 13.34·11 = 146.7 J

答え

Wf = 146.7 J
18.仕事とエネルギー

答え: d = 0.5946 m

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。

ばねのエネルギー

x = 20 cm = 0.2 m
½kx² = ½·800·0.2² = 16 J

斜面に沿ったエネルギー収支

m·g = 5·9.81 = 49.05 N
½kx² = mg·d·(sin θ + μ·cos θ)
d = 16/(49.05·(0.2588 + 0.3·0.9659)) = 0.5946 m

答え

d = 0.5946 m
19.機械振動

答え: ζ = 0.05003

質量–ばね系の自由振動:減衰のない固有角振動数は ωn = √(k/m)(rad/s)、固有振動数は fn = ωn/(2π)(Hz)です。粘性ダンパ c による減衰比は ζ = c/cc で、臨界減衰係数は cc = 2·√(k·m) = 2·m·ωn です。不足減衰系(ζ < 1)は ωd = ωn·√(1 − ζ²) で振動します。並列ばねは k = k₁ + k₂ と足し合わせ、直列ばねは逆数で 1/k = 1/k₁ + 1/k₂ と足し合わせます。

対数減衰率

δ = (1/1)·ln(40/29.2) = 0.3147

減衰比

ζ = δ/√(4π² + δ²) = 0.3147/√(39.48 + 0.09904) = 0.05003

答え

ζ = 0.05003
20.機械振動

答え: fn = 3.183 Hz

質量–ばね系の自由振動:減衰のない固有角振動数は ωn = √(k/m)(rad/s)、固有振動数は fn = ωn/(2π)(Hz)です。粘性ダンパ c による減衰比は ζ = c/cc で、臨界減衰係数は cc = 2·√(k·m) = 2·m·ωn です。不足減衰系(ζ < 1)は ωd = ωn·√(1 − ζ²) で振動します。並列ばねは k = k₁ + k₂ と足し合わせ、直列ばねは逆数で 1/k = 1/k₁ + 1/k₂ と足し合わせます。

等価ばね定数

直列のばねは同じ力を受け、たわみが足し合わされるので、逆数で合成します:1/k = 1/k1 + 1/k2。

1/keq = 1/12 + 1/60 → keq = 12·60/(12 + 60) = 10 kN/m

固有振動数

ωn = √(10000/25) = 20 rad/s
fn = 20/(2π) = 3.183 Hz

答え

fn = 3.183 Hz
21.熱機関と冷凍機

答え: ηII = 50.05 %

Second-law limits: no cycle working between reservoirs at absolute temperatures TH and TL can beat the reversible (Carnot) cycle: ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) and COPHP,max = TH/(TH − TL). Energy is always conserved: QH = W + QL.

絶対温度

TH = 400 + 273.15 = 673.1 K
TL = 10 + 273.15 = 283.1 K

カルノー効率

ηCarnot = 1 − 283.1/673.1 = 57.94 %

第二法則効率

ηII = η/ηCarnot = 29/57.94 = 50.05 %

答え

ηII = 50.05 %
22.熱機関と冷凍機

答え: Sgen = 1.305 kJ/K

エントロピー生成:1 サイクル全体で作動流体の正味のエントロピー変化はありません。高温熱源は QH/TH を失い、低温熱源は QL/TL を得ます。その和が生成されたエントロピーで、可逆サイクルでは 0、実際のサイクルではすべて正になります。

放出熱

QL = QH − W = 3200 − 1600 = 1600 kJ

絶対温度

TH = 550 + 273.15 = 823.1 K
TL = 35 + 273.15 = 308.1 K

熱源のエントロピー変化

ΔSH = −3200/823.1 = −3.888 kJ/K
ΔSL = 1600/308.1 = 5.192 kJ/K
Sgen = 5.192 − 3.888 = 1.305 kJ/K

答え

Sgen = 1.305 kJ/K

Sgen > 0 なので、この機関は実現可能ですが不可逆です。

23.ガスサイクル

答え: η = 58.91 %

ディーゼルサイクル(圧縮着火機関):1→2 等エントロピー圧縮、2→3 等圧加熱、3→4 等エントロピー膨張、4→1 等容放熱。圧縮比 r = v₁/v₂ と締切比 rc = v₃/v₂ を用いると、T₂ = T₁·rk−1、T₃ = rc·T₂、qin = cp·(T₃ − T₂)、η = 1 − (1/rk−1)·(rck − 1)/(k·(rc − 1)) です。

締切比の項

(31.4 − 1)/(1.4·(3 − 1)) = 1.306

熱効率

η = 1 − 0.3147·1.306 = 58.91%

答え

η = 58.91%
24.ガスサイクル

答え: η = 54.71 %

ブレイトンサイクル(ガスタービン):1→2 等エントロピー圧縮、2→3 等圧加熱、3→4 タービンでの等エントロピー膨張、4→1 等圧放熱。圧力比 rp = P₂/P₁ を用いると、T₂ = T₁·rp(k−1)/k、T₄ = T₃/rp(k−1)/k、η = 1 − 1/rp(k−1)/k です。正味仕事 w = cp·[(T₃ − T₄) − (T₂ − T₁)]。断熱効率を考えると、圧縮機には ws/ηc が必要で、タービンは ηt·ws を出します。

熱効率

η = 1 − 1/160.2857 = 1 − 1/2.208 = 54.71%

答え

η = 54.71%
25.連続の式とベルヌーイ

答え: v = 8.602 m/s

Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.

よどみ点へのベルヌーイの式

流れは先端で静止するので、その運動エネルギーが余分な圧力として現れます。

p0 − p = ½ρv²
v = √(2·37000/1000) = 8.602 m/s

答え

v = 8.602 m/s
26.連続の式とベルヌーイ

答え: Q = 15.56 L/s

Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.

マノメータによる圧力差

Δp = (13600 − 1000)·9.81·0.05 = 6180 Pa

のど部の流速

連続の式とベルヌーイの式を組み合わせると、理想的なのど部の流速が得られます。

β = 75/150 = 0.5
v2 = √(2·6180/(1000·(1 − 0.0625))) = 3.631 m/s

実際の流量

Q = Cd·A2·v2 = 0.97·0.004418·3.631 = 0.01556 m³/s = 15.56 L/s

答え

Q = 15.56 L/s
27.管摩擦と損失水頭

答え: Δp = 12.8 kPa

管内流れ:レイノルズ数 Re = V·D/ν = ρ·V·D/μ によって、層流(Re が約 2000 未満)と乱流(約 4000 超)を区別します。管に沿った摩擦によってダルシー・ワイスバッハの損失水頭 hf = f·(L/D)·V²/(2g) が生じます。継手、弁、入口、出口は局所損失 hm = ΣK·V²/(2g) を加えます。2 つの開放貯水槽の水面間では、エネルギー式から Δz = hf + hm(出口損失 K = 1 を含む)が得られます。

損失水頭

hf = f·(L/D)·V²/(2g) = 0.03·(50/0.15)·1.6²/(2·9.81) = 1.305 m

圧力損失

管は水平で流速も変わらないので、エネルギー式は Δp = ρ·g·hf になります。

Δp = ρ·g·hf = 1000·9.81·1.305 = 12800 Pa = 12.8 kPa

答え

Δp = 12.8 kPa
28.管摩擦と損失水頭

答え: Q = 83.91 L/s

管内流れ:レイノルズ数 Re = V·D/ν = ρ·V·D/μ によって、層流(Re が約 2000 未満)と乱流(約 4000 超)を区別します。管に沿った摩擦によってダルシー・ワイスバッハの損失水頭 hf = f·(L/D)·V²/(2g) が生じます。継手、弁、入口、出口は局所損失 hm = ΣK·V²/(2g) を加えます。2 つの開放貯水槽の水面間では、エネルギー式から Δz = hf + hm(出口損失 K = 1 を含む)が得られます。

損失係数

f·L/D = 0.025·100/0.15 = 16.67
ΣK = 0.5 + 4·0.3 + 5 + 1.0 = 7.7

エネルギー式

どちらの水面も大気に開放され、ほぼ静止しているので、Δz はすべて管内で失われます。

28 = (16.67 + 7.7)·V²/(2·9.81)
V = √(2·9.81·28/24.37) = 4.748 m/s

流量

Q = V·π·D²/4 = 4.748·0.01767 = 0.08391 m³/s = 83.91 L/s

答え

Q = 83.91 L/s
29.ポンプとポンプ動力

答え: Pel = 1.652 kW

ポンプ動力:水に与えられる動力は水動力 Ph = ρ·g·Q·H = Q·Δp です。ポンプには損失があるので、軸はそれより大きな動力を伝えなければなりません:Pshaft = Ph/ηpump。モーターが消費する電力は Pel = Pshaft/ηmotor です。ポンプ揚程 H はエネルギー式から求めます:実揚程 + 圧力水頭 + すべての摩擦損失と局所損失。

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.01·12 = 1.177 kW

軸動力

Pshaft = 1.177/0.75 = 1.57 kW

電力

Pel = 1.57/0.95 = 1.652 kW

答え

Pel = 1.652 kW
30.ポンプとポンプ動力

答え: Pshaft = 20.55 kW

ポンプ動力:水に与えられる動力は水動力 Ph = ρ·g·Q·H = Q·Δp です。ポンプには損失があるので、軸はそれより大きな動力を伝えなければなりません:Pshaft = Ph/ηpump。モーターが消費する電力は Pel = Pshaft/ηmotor です。ポンプ揚程 H はエネルギー式から求めます:実揚程 + 圧力水頭 + すべての摩擦損失と局所損失。

速度

V = Q/A = 0.04/0.01767 = 2.264 m/s
V²/(2g) = 0.2611 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (83.33 + 8)·0.2611 = 23.85 m

ポンプ揚程

H = Δz + hL = 17 + 23.85 = 40.85 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.04·40.85/0.78 = 20.55 kW

答え

Pshaft = 20.55 kW
31.熱伝導

答え: Q̇ = 50.48 W

熱抵抗:直列の層を通る定常の熱流は、直列抵抗を流れる電流と同じように振る舞います:Q̇ = ΔT/ΣR。伝導層では R = L/(kA)、対流面では R = 1/(hA)。各層の温度降下はその熱抵抗に比例します。

各層の熱抵抗

R1 = L1/(k1A) = (150/1000)/(0.04·12) = 0.3125 K/W
R2 = L2/(k2A) = (75/1000)/(1.4·12) = 0.004464 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 16/0.317 = 50.48 W

答え

Q̇ = 50.48 W
32.熱伝導

答え: U = 0.1877 W/(m²·K)

熱抵抗:直列の層を通る定常の熱流は、直列抵抗を流れる電流と同じように振る舞います:Q̇ = ΔT/ΣR。伝導層では R = L/(kA)、対流面では R = 1/(hA)。各層の温度降下はその熱抵抗に比例します。

直列の熱抵抗

Ri = 1/(hiA) = 1/(7.7·6) = 0.02165 K/W
R1 = L1/(k1A) = (200/1000)/(0.04·6) = 0.8333 K/W
R2 = L2/(k2A) = (25/1000)/(0.17·6) = 0.02451 K/W
Ro = 1/(hoA) = 1/(20·6) = 0.008333 K/W
ΣR = 0.8878 K/W

熱通過率

U = 1/(A·ΣR) = 1/(6·0.8878) = 0.1877 W/(m²·K)

答え

U = 0.1877 W/(m²·K)
33.歯車と歯車列

答え: T2 = 158.1 N·m

歯車列:かみ合う歯車は歯の部分で同じ速さで動くので、n₁·N₁ = n₂·N₂ となり、歯数の少ない歯車ほど速く回ります。同じ軸に取り付けた歯車は同じ回転速度で回ります。動力は(損失を除いて)そのまま伝わるので、回転速度が下がるとトルクは上がります:T₂ = T₁·(N₂/N₁)·η。また P = T·ω、ただし ω = 2π·n/60。

モータのトルク

ω₁ = 2π·1450/60 = 151.8 rad/s
T1 = P/ω₁ = 11000/151.8 = 72.44 N·m

出力トルク

T2 = 72.44·72/32·0.97 = 158.1 N·m

答え

T2 = 158.1 N·m
34.歯車と歯車列

答え: P = 2.302 kW

ベルト駆動:滑りがなければ、両方のプーリの周速は同じ v = π·d·n/60 なので、n₁·d₁ = n₂·d₂ となり、速度比は d₂/d₁ です。滑り s があると、従動プーリは少し遅くなります:n₂ = n₁·(d₁/d₂)·(1 − s)。伝達される動力は P = (T₁ − T₂)·v で、T₁ は張り側の張力、T₂ は緩み側の張力です。摩擦が張力の比を制限します:T₁/T₂ = eμθ、θ は小さいほうのプーリへの巻き付け角(ラジアン)です。平行掛けのベルトでは、軸間距離 C を使って θ = π − 2·asin((d₂ − d₁)/(2C)) です。

巻き付け角

θ = 2.832 rad

緩み側の張力

eμθ = e0.8496 = 2.339, T₂ = 1000/2.339 = 427.6 N

電力

P = (1000 − 427.6)·4.021 = 2.302 kW

答え

P = 2.302 kW
35.ラプラス変換

答え: A = 5.667

ラプラス変換:F(s) = ∫₀∞ f(t)·e−st dt は、微分を s の掛け算に変えます:ℒ{f′} = sF(s) − f(0)、ℒ{f″} = s²F(s) − s·f(0) − f′(0)。基本的な変換対:1 ↔ 1/s、e−at ↔ 1/(s + a)、t·e−at ↔ 1/(s + a)²、sin ωt ↔ ω/(s² + ω²)、cos ωt ↔ s/(s² + ω²)。また、s → s + σ のシフトは、時間関数に e−σt を掛けることに対応します。有理式 F(s) を逆変換するには部分分数に分解します。単極 p の留数は A = [(s − p)·F(s)] の s = p での値です。最終値の定理:f(∞) = lims→0 s·F(s)。これは s·F(s) のすべての極が左半平面にあるときに限り成り立ちます。初期値の定理:f(0⁺) = lims→∞ s·F(s)。

分母を因数分解する

s² + 10s + 16 = (s + 2)(s + 8)
F(s) = A/(s + 2) + B/(s + 8)

隠蔽法

(s + 2) を掛けて s = −2 とおくと、B の項は消えます。

A = 34/(s + 8) |s = −2 = 34/6 = 5.667

答え

A = 5.667
(B = −5.667 なので、f(t) = 5.667·e−2t − 5.667·e−8t)
36.ラプラス変換

答え: y(0.2) = 4.767

ラプラス変換:F(s) = ∫₀∞ f(t)·e−st dt は、微分を s の掛け算に変えます:ℒ{f′} = sF(s) − f(0)、ℒ{f″} = s²F(s) − s·f(0) − f′(0)。基本的な変換対:1 ↔ 1/s、e−at ↔ 1/(s + a)、t·e−at ↔ 1/(s + a)²、sin ωt ↔ ω/(s² + ω²)、cos ωt ↔ s/(s² + ω²)。また、s → s + σ のシフトは、時間関数に e−σt を掛けることに対応します。有理式 F(s) を逆変換するには部分分数に分解します。単極 p の留数は A = [(s − p)·F(s)] の s = p での値です。最終値の定理:f(∞) = lims→0 s·F(s)。これは s·F(s) のすべての極が左半平面にあるときに限り成り立ちます。初期値の定理:f(0⁺) = lims→∞ s·F(s)。

方程式を変換する

(s² + 7s + 6)·Y(s) = 5s + 35 + 12/s
Y(s) = (5s² + 35s + 12)/(s(s + 1)(s + 6))

部分分数展開

A = 12/(1·6) = 2
B = −18/((−1)·5) = 3.6
C = −18/((−6)·(−5)) = −0.6

結果

y(t) = 2 + 3.6·e−t − 0.6·e−6t
y(0.2) = 4.767

答え

y(0.2) = 4.767

検算:y(0) = A + B + C = 5。

37.ブロック線図

答え: K = 100

ブロック線図の代数:直列のブロックは掛け算(G₁·G₂)、並列のブロックは足し算(G₁ + G₂)になります。前向き経路 G とフィードバック経路 H のループは、負帰還なら T = G/(1 + G·H)、正帰還なら G/(1 − G·H) にまとめられます。内側のループから先に簡単にし、外側へ進みます。閉ループの極は特性方程式 1 + G(s)H(s) = 0 の根です。ループが安定なら、大きさ A のステップ入力に対する定常出力は A·T(0) です。単位フィードバックでは、ステップに対する定常偏差は Kp = G(0) として A/(1 + Kp) です。

固有方程式

1 + K/(s(s + 20)) = 0 → s² + 20s + K = 0

重根

2 次式 s² + βs + γ は、β² = 4γ のとき重根をもちます。

20² = 4K → K = 400/4 = 100

答え

K = 100
このとき、閉ループの極は 2 つとも s = −10 にあります。
38.ブロック線図

答え: T(0) = 0.4138

ブロック線図の代数:直列のブロックは掛け算(G₁·G₂)、並列のブロックは足し算(G₁ + G₂)になります。前向き経路 G とフィードバック経路 H のループは、負帰還なら T = G/(1 + G·H)、正帰還なら G/(1 − G·H) にまとめられます。内側のループから先に簡単にし、外側へ進みます。閉ループの極は特性方程式 1 + G(s)H(s) = 0 の根です。ループが安定なら、大きさ A のステップ入力に対する定常出力は A·T(0) です。単位フィードバックでは、ステップに対する定常偏差は Kp = G(0) として A/(1 + Kp) です。

内側のループ

Gᵢ = G₂/(1 + G₂H₁) = 5/(s + 17)

外側のループ

T(s) = G₁Gᵢ/(1 + G₁GᵢH₂) = 60/((s + 5)(s + 17) + 60)

直流ゲイン

T(0) = 60/(5·17 + 60) = 60/145 = 0.4138

答え

T(0) = 0.4138
39.ステップ応答

答え: ζ = 0.3333

ステップ応答:1 次系 G(s) = K/(τs + 1) は、大きさ U のステップに対して y(t) = K·U·(1 − e−t/τ) と応答します。時定数 τ の後に最終値 K·U の 63.2 %、4τ の後に 98 % に達します。2 次系 G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) の極は s = −ζωₙ ± jωd(ωd = ωₙ√(1 − ζ²))です。0 < ζ < 1 のとき、ステップ応答にはオーバーシュートがあります:オーバーシュート率 %OS = 100·e−ζπ/√(1−ζ²)、ピーク時間 tp = π/ωd、2 % 整定時間 ts ≈ 4/(ζωₙ)。分母を s² + 2ζωₙs + ωₙ² と比べると、ωₙ と ζ が読み取れます。質量・ばね・ダンパでは ωₙ = √(k/m)、ζ = c/(2√(km))、直列 RLC 回路では ωₙ = 1/√(LC)、ζ = (R/2)·√(C/L) です。

標準形と比べる

分母は s² + 2ζωₙs + ωₙ² と一致しなければなりません。

ωₙ² = 9 → ωₙ = 3 rad/s
2ζωₙ = 2 → ζ = 2/(2·3) = 0.3333

結果

ζ = 0.3333

答え

ζ = 0.3333
40.ステップ応答

答え: K = 38.35

ステップ応答:1 次系 G(s) = K/(τs + 1) は、大きさ U のステップに対して y(t) = K·U·(1 − e−t/τ) と応答します。時定数 τ の後に最終値 K·U の 63.2 %、4τ の後に 98 % に達します。2 次系 G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) の極は s = −ζωₙ ± jωd(ωd = ωₙ√(1 − ζ²))です。0 < ζ < 1 のとき、ステップ応答にはオーバーシュートがあります:オーバーシュート率 %OS = 100·e−ζπ/√(1−ζ²)、ピーク時間 tp = π/ωd、2 % 整定時間 ts ≈ 4/(ζωₙ)。分母を s² + 2ζωₙs + ωₙ² と比べると、ωₙ と ζ が読み取れます。質量・ばね・ダンパでは ωₙ = √(k/m)、ζ = c/(2√(km))、直列 RLC 回路では ωₙ = 1/√(LC)、ζ = (R/2)·√(C/L) です。

閉ループ伝達関数

T(s) = K/(s² + 5s + K)

オーバーシュートから求める減衰比

ζ = −ln(0.25)/√(π² + ln²(0.25)) = 0.4037

利得

ωₙ = 5/(2ζ) = 6.193 rad/s
K = ωₙ² = 38.35

答え

K = 38.35
41.定常偏差と PID 制御

答え: K = 36

定常偏差:開ループ伝達関数 G(s) をもつ単位フィードバックのループでは、偏差は E(s) = R(s)/(1 + G(s)) で、最終値の定理により ess = lims→0 s·E(s) です。システムの型は、G(s) の s = 0 にある極の個数です。偏差定数は Kp = lim G(s)、Kv = lim s·G(s)、Ka = lim s²·G(s) で、いずれも s → 0 の極限です。大きさ A のステップ入力では ess = A/(1 + Kp)、ランプ入力 A·t では A/Kv、放物線入力 A·t²/2 では A/Ka になります。0 型のループはステップ入力に有限の偏差で追従し、1 型のループはステップ入力に偏差なく追従し、ランプ入力には有限の偏差で追従します。

必要な偏差定数

1/(1 + Kp) = 0.1 → Kp = 9

利得

K/4 = 9 → K = 36

答え

K = 36

閉ループの安定性も確認しましょう:この 2 次のループでは、K > 0 ならどんな値でも安定です。

42.定常偏差と PID 制御

答え: ess = 1.25

P 制御、PI 制御、PID 制御:制御器 C(s) = Kp + Ki/s + Kd·s は偏差に作用します。比例ゲインはループを速くし、0 型の制御対象の偏差を小さくしますが、なくすことはできません。積分項は s = 0 に極を加えます:システムの型が 1 つ上がり、ステップ入力に偏差なく追従します。微分項は減衰を加えます。ゲインを決めるには、閉ループの特性多項式を書き、s² + 2ζωns + ωn² と比較します。すると %OS = 100·e−ζπ/√(1−ζ²)、tp = π/ωd(ωd = ωn√(1 − ζ²))、ts ≈ 4/(ζωn) となります。

速度偏差定数

s·C(s) = 2s + 10 → s → 0 のとき 10。

Kv = 10·4/10 = 4

ランプ偏差

ess = 5/4 = 1.25

答え

ess = 1.25

PI 制御器はステップ偏差を完全になくします。ランプ入力には一定の遅れで追従します。

これは独自の練習教材です。ExamProblems はこの試験を実施する団体とは提携しておらず、その承認も受けていません。試験名は内容を説明するためにのみ使用しています。 FE examはNCEESが実施する試験です。同団体はExamProblemsと提携しておらず、推奨もしていません。

ExamProblems · examproblems.com/courses/fe-exam-mechanical