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.はりの支点反力ふつう
P1​ = 65 kNP2​ = 50 kNABCRB​ = ?0.5 m3.5 m1 m

はりは A でピン支点、A から 4 m の B でローラー支点に支持され、B から自由端 C まで張り出している。上向きを正として鉛直反力 RB を求めよ。

AB = 4 mAC = 5 mP1 = 65 kN 位置 0.5 mP2 = 50 kN 位置 5 m
RB = kN
2.はりの支点反力むずかしい
P = 75 kNADCMA​ = ?3 m4 m3 m

複合ばりは A で壁に固定され、C(A から 3 m)に中間ヒンジをもち、D(A から 10 m)でローラーに支持されている。荷重は A から 7 m の位置の集中荷重 P である。 壁における反力モーメント MA の大きさを求めよ。

AC = 3 mAD = 10 mP = 75 kN 位置 7 m
MA = kN·m
3.軸応力とひずみふつう
250 kNØ20100 GPaØ100/50400 mmδ = ?

鋼 の棒(直径 20 mm、E = 200 GPa)が、黄銅 でできた管(外径 100 mm、内径 50 mm、E = 100 GPa)の中に入っています。どちらも長さ 400 mm で、剛な板の間で 250 kN によって一緒に圧縮されます。 縮み δ を求めなさい。

P = 250 kN棒: Ø20 mm, E = 200 GPa管: Ø100/Ø50 mm, E = 100 GPaL = 400 mm
δ = mm
4.軸応力とひずみむずかしい
Ø 35 mm銅1 mØ 16 mm銅0.5 mP = 95 kN

段付き棒は壁に固定され、自由端に荷重を受ける。全体の伸び δ を求めよ。

P = 95 kN区間 1:Ø35 mm、1 m、銅(E = 117 GPa)区間 2:Ø16 mm、0.5 m、銅(E = 117 GPa)
δ = mm
5.軸のねじりふつう
d = ?T = 2000 N·mL = 1.2 m

長さ 1.2 m の鋼製の中実軸がトルク 2000 N·m を受ける。両端間のねじれ角は 1° を超えてはならない。この制限を満たす最小の直径 d を求めよ。

T = 2000 N·mL = 1.2 mG = 80 GPaφallow = 1°
d = mm
6.軸のねじりむずかしい
Td = ?

中実軸が 500 rpm で 10 kW を伝達し、せん断応力が 80 MPa を超えないようにしたい。条件を満たす最小直径 d を求めよ。

P = 10 kWN = 500 rpmτallow = 80 MPa
d = mm
7.せん断力と曲げモーメントふつう
P1​ = 20 kNP2​ = 80 kNAB2 m3 m3 m

単純ばりに 2 つの集中荷重が作用する。はりの最大曲げモーメント Mmax を求めよ。

L = 8 mP1 = 20 kN 位置 2 mP2 = 80 kN 位置 5 m
Mmax = kN·m
8.せん断力と曲げモーメントむずかしい
8 kN/m25 kN45 kNABCDEFG1 m1 m1 m1 m1 m1 m

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

25 kN @ 6 m45 kN @ 2 m8 kN/m, 0–4 m
  1. 1A (x = 0 m)
  2. 2B (x = 1 m)
  3. 3C (x = 2 m)
  4. 4D (x = 3 m)
  5. 5E (x = 4 m)
  6. 6F (x = 5 m)
  7. 7G (x = 6 m)
9.モールの応力円と応力の座標変換ふつう
σx​ = −30σy​ = −30τxy​ = −50MPa

ある点で σₓ = −30 MPa、σᵧ = −30 MPa、τₓᵧ = −50 MPa である(引張を正とする)。要素を反時計回りに 40°回転させる。回転した要素の、回転後の y′ 軸に垂直な面に働く垂直応力 σy′ を求めよ。

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

直径 25 mm の中実軸が曲げモーメント 100 N·mとトルク 500 N·m を同時に受ける。 曲げ応力が引張で最大となる上表面で、 面内最大せん断応力 τmax を求めよ。

d = 25 mmM = 100 N·mT = 500 N·m
τmax = MPa
11.はりのたわみふつう
w = 4 kN/mP = 27 kNABδ1.5 m

片持ちばりの全長に等分布荷重 w が、自由端に集中荷重 P が作用する。重ね合わせを用いて自由端のたわみ δ を求めよ。

P = 27 kNw = 4 kN/mL = 1.5 mE = 70 GPa (アルミニウム)I = 57 × 10⁶ mm⁴
δ = mm
12.はりのたわみむずかしい
w = 5 kN/mP = 28 kNδL = 2.5 m

単純ばりに等分布荷重 w とスパン中央の集中荷重 P が作用する。重ね合わせを用いて中央たわみ δmax を求めよ。

P = 28 kNw = 5 kN/mL = 2.5 mE = 70 GPa (アルミニウム)I = 160 × 10⁶ mm⁴
δmax = mm
13.柱の座屈ふつう
PL = 3 m30 × 60 mm断面

鋼製の長さ L の柱が両端固定である。オイラーの式を用いて座屈荷重 Pcr を求めよ。

L = 3 mE = 200 GPa (鋼)b × h = 30 × 60 mm
Pcr = kN
14.柱の座屈むずかしい
P = 30 kNL = 4 md = ?断面

アルミニウム製の長さ L の中実円柱が両端固定である。座屈に対する安全率 3 で 30 kN を支える必要がある。最小直径 d を求めよ。

P = 30 kNn = 3L = 4 mE = 70 GPa (アルミニウム)
d = mm
15.放物運動ふつう
v0​ = ?20°R = 85 m

水平な地面からボールを水平から 20° 上向きに打ち出すと、85 m 先に着地する。空気抵抗を無視し、g = 9.81 m/s² とする。打ち出した速さ v0 を求めよ。

θ = 20°R = 85 mg = 9.81 m/s²
v0 = m/s
16.放物運動むずかしい
h = 55 mv0​ = 33 m/s25°x = ?

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

v0 = 33 m/sθ = 25°h = 55 mg = 9.81 m/s²
x = m
17.仕事とエネルギーふつう
70°L = 1 m10 kg

質量 10 kg の小球が長さ 1 m の軽い糸につり下げられている。糸が鉛直と 70° をなすまで横に引き、静かに放す。小球が最下点を通過するときの糸の張力 T を求めよ。

m = 10 kgL = 1 mθ = 70°g = 9.81 m/s²
T = N
18.仕事とエネルギーむずかしい
h = 10.4 mR = 4 mN = ?

質量 1 kg の台車が、半径 4 m の鉛直な円形ループにつながるなめらかな軌道上の高さ 10.4 m から静かに放される。 ループの最高点で軌道が台車に及ぼす垂直抗力 N を求めよ。

m = 1 kgh = 10.4 mR = 4 mg = 9.81 m/s²
N = N
19.機械振動ふつう
tx50 mm5 周期7.6 mm

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

x0 = 50 mmx5 = 7.6 mmn = 5
ζ =
20.機械振動むずかしい
k1​ = 25k2​ = 20k3​ = 12m = 40 kg

40 kg の質量がばね k₃ = 12 kN/m につり下げられている。k₃ の上端は剛体の棒に取り付けられ、その棒は横に並べた 2 本のばね k₁ = 25 kN/m と k₂ = 20 kN/m からつり下げられている。棒は水平を保つ。 減衰とばね・棒の質量は無視する。上下振動の固有振動数 fn を Hz で求めよ。

m = 40 kgk1 = 25 kN/mk2 = 20 kN/mk3 = 12 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 475 °CTL​ = 35 °CEQH​QL​W

熱機関が 475 °C の高温熱源と 35 °C の低温熱源の間で働く。 可能な最大の熱効率 ηmax を求めよ。

TH = 475 °CTL = 35 °C
ηmax = %
22.熱機関と冷凍機むずかしい
TH​ = 650 °CTL​ = 35 °CEQH​ = 1600 kJQL​ = 848 kJW = 752 kJ

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

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

理想ディーゼルサイクルの圧縮比は r = 20、締切比は rc = v₃/v₂ = 2 である。 圧縮開始時の空気の温度は 295 K である。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 加熱終わりの温度 T3 を K で求めよ。

r = 20rc = 2T1 = 295 Kk = 1.4
T3 = K
24.ガスサイクルむずかしい
vP1234qin​qout​

ガスタービン発電設備が、圧縮機の圧力比 16 のブレイトンサイクルで運転されている。空気は 288 K で圧縮機に入り、1500 K でタービンに入る。 圧縮機の断熱効率は 82%、タービンの断熱効率は 85% である。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 空気 1 kg あたりの正味仕事 wnet を kJ/kg で求めよ。

rp = 16T1 = 288 KT3 = 1500 Kηc = 82%ηt = 85%k = 1.4cp = 1.005 kJ/(kg·K)
wnet = kJ/kg
25.連続の式とベルヌーイふつう
12d1​ = 75 mm, p1​ = 280 kPad2​ = 40 mm, p2​ = ?

水が、75 mm から 40 mm に細くなる水平管を流れる。太い部分の圧力は 280 kPa、流速は 1 m/s である。損失を無視して、細い部分の圧力 p2 を求めよ。

d1 = 75 mmd2 = 40 mmp1 = 280 kPav1 = 1 m/sρ = 1000 kg/m³
p2 = kPa
26.連続の式とベルヌーイむずかしい
12d1​ = 75 mm, p1​ = 230 kPad2​ = 50 mm, p2​ = ?Δz = 1 m

水が、75 mm から 50 mm に細くなりながら 1 m 上昇する管を上向きに流れる。断面 1 での圧力は 230 kPa、流速は 1 m/s である。損失を無視して、断面 2 での圧力 p2 を求めよ。

d1 = 75 mmd2 = 50 mmp1 = 230 kPav1 = 1 m/sΔz = 1 mρ = 1000 kg/m³
p2 = kPa
27.管摩擦と損失水頭ふつう
V = 2.6 m/sD = 75 mmL = 500 mf = 0.022

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

V = 2.6 m/sL = 500 mD = 75 mmf = 0.022ρ = 1000 kg/m³
Δp = kPa
28.管摩擦と損失水頭むずかしい
K = 0.2Δz = ?L = 150 m, D = 300 mm, f = 0.0153× 曲がり (K = 0.3)

2 つの開放貯水槽が、長さ 150 m、内径 300 mm、ダルシーの管摩擦係数 f = 0.015 の管でつながれている。局所損失は 角ばった入口(K = 0.5)、3 個の曲がり(各 K = 0.3)、弁(K = 0.2)、下の貯水槽への出口(K = 1.0) である。g = 9.81 m/s² とする。流量 80 L/s となる貯水槽間の水位差 Δz を求めよ。

Q = 80 L/sL = 150 mD = 300 mmf = 0.015Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 0.2Kexit = 1.0
Δz = m
29.ポンプとポンプ動力ふつう
PQ = 60 L/sH = 46 m

効率 65% のポンプが揚程 46 m に対して水(ρ = 1000 kg/m³)を 60 L/s 送っている。g = 9.81 m/s² とする。必要な軸動力 Pshaft を kW で求めよ。

Q = 60 L/sH = 46 mη = 65%g = 9.81 m/s²
Pshaft = kW
30.ポンプとポンプ動力むずかしい
PQ = 20 L/sΔz = 8 mL = 100 m, D = 100 mm, f = 0.025, ΣK = 3

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

Q = 20 L/sΔz = 8 mL = 100 mD = 100 mmf = 0.025ΣK = 3η = 80%
Pshaft = kW
31.熱伝導ふつう
ロックウール, k = 0.04250 mmガラス, k = 0.96150 mm21 °C1 °CQ̇ = ?

壁はロックウール 250 mm とガラス 150 mm からなり、面積は 6 m² である。外側の両面は 21 °C と 1 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 250 mm, k1 = 0.04 W/(m·K)L2 = 150 mm, k2 = 0.96 W/(m·K)A = 6 m²T1 = 21 °CT3 = 1 °C
Q̇ = W
32.熱伝導むずかしい
石こうボード, k = 0.17250 mm木材, k = 0.17150 mm室内 21 °C屋外 -7 °CTs​ = ?Q̇

250 mm の石こうボードと 150 mm の木材からなる壁が、21 °C の部屋と -7 °C の外気を隔てている。熱伝達率は室内側 6 W/(m²·K)、屋外側 20 W/(m²·K)、壁の面積は 4 m² である。 壁の内表面温度 Ts を求めよ。

L1 = 250 mm, k1 = 0.17 W/(m·K)L2 = 150 mm, k2 = 0.17 W/(m·K)A = 4 m²hi = 6 W/(m²·K)ho = 20 W/(m²·K)Troom = 21 °CTout = -7 °C
Ts = °C
33.歯車と歯車列ふつう
Nw​ = 80z = 2, 960 rpm

条数 2 のウォームが 960 rpm で回転し、歯数 80 のウォームホイールを駆動する。ホイールの回転速度 nw を求めよ。

z = 2Nw = 80n = 960 rpm
nw = rpm
34.歯車と歯車列むずかしい
A = 30B = 48C = 15D = 72

複合歯車列において、1500 rpm で回転する入力軸上の歯車 A(歯数 30)が歯車 B(歯数 48)を駆動する。歯車 C(歯数 15)は B と同じ軸に固定されており、出力軸上の歯車 D(歯数 72)を駆動する。出力回転速度 nout を求めよ。

nin = 1500 rpmA = 30B = 48C = 15D = 72
nout = rpm
35.ラプラス変換ふつう
t (s)f(t)t = 1f(1) = ?

F(s) = (5s + 21)/(s² + 4s + 4) の f(t) を求め、f(1) を計算せよ。

F(s) = (5s + 21)/(s² + 4s + 4)t = 1 s
f(1) =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.1y(0.1) = ?

ラプラス変換を用いて、y(0) = −5、y′(0) = −1 の条件で y″ + 11y′ + 30y = 0 を解き、y(0.1) を求めよ。

y″ + 11y′ + 30y = 0y(0) = −5y′(0) = −1t = 0.1 s
y(0.1) =
37.ブロック線図ふつう
R(s)+54s(s + 2)G(s)C(s)−

単位負帰還ループの G(s) = 54/(s(s + 2)) である。基準入力はランプ r(t) = 5·t である。定常偏差 ess を求めよ。

G(s) = 54/(s(s + 2))r(t) = 5·t
ess =
38.ブロック線図むずかしい
R(s)+12コントローラ+D(s)+2s + 9プラントC(s)−

比例コントローラ K = 12 が、プラント G(s) = 2/(s + 9) を単位負帰還ループで駆動する。外乱 D はプラントの入力に加わる(コントローラの出力に加算される)。 大きさ 10 の基準ステップと大きさ 1 の外乱ステップを同時に加える。定常出力 css を求めよ。

K = 12G(s) = 2/(s + 9)D = 1R = 10
css =
39.ステップ応答ふつう
t (s)y(t)2

ある系の伝達関数は G(s) = 102/(s² + 5s + 51) で、単位ステップ入力を受ける。 極は s = σ ± jωd である。その実部 σ を求めよ。

G(s) = 102/(s² + 5s + 51)
σ =
40.ステップ応答むずかしい
R(s)+60s(s + 2)G(s)C(s)−

単位負帰還ループの G(s) = 60/(s(s + 2)) である。ts ≈ 4/(ζωₙ) を用いて、閉ループのステップ応答の 2 % 整定時間 ts を求めよ。

G(s) = 60/(s(s + 2))
ts = s
41.定常偏差と PID 制御ふつう
R(s)+K = 19210s + 1C(s)−

プロセス G(s) = 2/(10s + 1) を、比例ゲイン K = 19 で単位フィードバックのループにより制御します。ステップ入力に対する定常偏差 ess を、ステップの大きさに対する百分率で求めなさい。

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

比例制御器 K = 25 が、単位フィードバックのループで制御対象 G(s) = 4/(s(s + 5)) を制御しています。 ステップ応答のピーク時間 tp を求めよ。

G(s) = 4/(s(s + 5))K = 25
tp = s

3溶液

1.はりの支点反力

答え: RB = 70.63 kN

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

A まわりのモーメント

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

ΣMA = 0: RB·4 = 65·0.5 + 50·5
RB = 282.5/4 = 70.63 kN

鉛直方向のつり合い

RA = 65 + 50 − 70.63 = 44.38 kN

答え

RB = 70.63 kN
2.はりの支点反力

答え: MA = 96.43 kN·m

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

部分 CD

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

ΣMC = 0: RD·7 = 75·4
RD = 42.86 kN
VC = 75 − 42.86 = 32.14 kN

部分 AC

ヒンジ力 VC は AC に下向きに作用します(CD への作用と大きさが等しく向きが逆)。

RA = 32.14 = 32.14 kN
MA = 32.14·3 = 96.43 kN·m

答え

MA = 96.43 kN·m
3.軸応力とひずみ

答え: δ = 0.1534 mm

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

面積

Ar = π·20²/4 = 314.2 mm²
At = π(100² − 50²)/4 = 5890 mm²

合成した軸剛性

ΣEA = 200·10³·314.2 + 100·10³·5890 = 6.519 × 10⁸ N

縮み

δ = 250000·400/6.519 × 10⁸ = 0.1534 mm

答え

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

答え: δ = 2.863 mm

軸力部材:各等断面区間は δ = PL/(AE) に従います。直列の部材では伸びが足し合わされます。

両区間で同じ力

自由端が P を受けるので、棒のどの断面も同じ軸力 P を受けます。変わるのは断面積と材料だけです。

A1 = π·35²/4 = 962.1 mm², A2 = π·16²/4 = 201.1 mm²
σ1 = 98.74 MPa
σ2 = 472.5 MPa

各区間の伸びを足す

δ1 = σ1L1/E1 = 98.74·1000/117000 = 0.8439 mm
δ2 = σ2L2/E2 = 472.5·500/117000 = 2.019 mm
δ = δ1 + δ2 = 2.863 mm

答え

δ = 2.863 mm
5.軸のねじり

答え: d = 64.69 mm

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

必要な断面二次極モーメント

T = 2000 N·m = 2 × 10⁶ N·mm
L = 1.2 m = 1200 mm
G = 80 GPa = 80000 MPa
φ = 1·π/180 = 0.01745 rad
J = TL/(Gφ) = 2 × 10⁶·1200/(80000·0.01745) = 1.719 × 10⁶ mm⁴

直径

d = (32J/π)1/4 = 64.69 mm

答え

d = 64.69 mm

実際には、次の規格寸法に切り上げます。

6.軸のねじり

答え: d = 22.99 mm

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

トルク

ω = 2πN/60 = 2π·500/60 = 52.36 rad/s
T = P/ω = 10000/52.36 = 191 N·m

許容応力から直径を求める

d = ∛(16T/(π·τallow)) = ∛(16·190986/(π·80)) = 22.99 mm

答え

d = 22.99 mm

実際には、次の規格寸法に切り上げます。

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

答え: Mmax = 165 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。 最大曲げモーメント:dM/dx = V なので、曲げモーメントはせん断力が 0 を通過する位置、つまり集中荷重の位置か分布荷重の区間内で最大になります。

支点反力

はり全体のつり合い。

ΣMA = 0: RB·8 = 20·2 + 80·5
RB = 440/8 = 55 kN
ΣFy = 0: RA = 100 − 55 = 45 kN

各荷重位置のモーメント

各荷重より左側の部分を残します。

荷重点 P1: M = RA·2 = 45·2 = 90 kN·m
荷重点 P2: M = RB·3 = 55·3 = 165 kN·m

最大値

Mmax = 165 kN·m (荷重点 P2)

答え

Mmax = 165 kN·m
8.せん断力と曲げモーメント

答え: E (x = 4 m)

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

支点反力

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

R1 = 26 kN, R2 = 76 kN

各点での曲げモーメント

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

MA = 0, MB = 22, MC = 36, MD = −3, ME = −50 kN·m, MF = −25, MG = 0

答え

点 E: ME = −50 kN·m
9.モールの応力円と応力の座標変換

答え: σy′ = 19.24 MPa

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

y′ 面の座標変換の式

σy′ = −30 − 0·0.1736 − (−50)·0.9848 = 19.24 MPa

確認

垂直応力の和は回転しても変わりません。

σx′ = −79.24 MPa, σx′ + σy′ = −60 = σx + σy

答え

σy′ = 19.24 MPa
10.モールの応力円と応力の座標変換

答え: τmax = 166.2 MPa

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

表面の応力

σ = 32M/(πd³) = 32·100000/(π·25³) = 65.19 MPa
τ = 16T/(πd³) = 16·500000/(π·25³) = 163 MPa
στxyCτmax​

モールの応力円

C = σ/2 = 32.59 MPa
R = √(32.59² + 163²) = 166.2 MPa

答え

τmax = 166.2 MPa
11.はりのたわみ

答え: δ = 8.247 mm

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

単位をそろえる

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

P = 27000 N, w = 4 N/mm, L = 1500 mm
EI = 3.99 × 10¹² N·mm²

各荷重を単独で

δP = PL³/(3EI) = 27000·1500³/(3·3.99 × 10¹²) = 7.613 mm
δw = wL⁴/(8EI) = 4·1500⁴/(8·3.99 × 10¹²) = 0.6344 mm

足し合わせる

δ = 7.613 + 0.6344 = 8.247 mm

答え

δ = 8.247 mm
12.はりのたわみ

答え: δmax = 1.041 mm

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

単位をそろえる

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

P = 28000 N, w = 5 N/mm, L = 2500 mm, EI = 1.12 × 10¹³ N·mm²

各荷重を単独で

δP = PL³/(48EI) = 28000·2500³/(48·1.12 × 10¹³) = 0.8138 mm
δw = 5wL⁴/(384EI) = 5·5·2500⁴/(384·1.12 × 10¹³) = 0.2271 mm

足し合わせる

δmax = 0.8138 + 0.2271 = 1.041 mm

答え

δmax = 1.041 mm
13.柱の座屈

答え: Pcr = 118.4 kN

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

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

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

Imin = h·b³/12 = 60·30³/12 = 135000 mm⁴

有効座屈長さ

KL = 0.5·3000 = 1500 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·200000·135000/1500² = 118435 N = 118.4 kN

答え

Pcr = 118.4 kN
14.柱の座屈

答え: d = 57.08 mm

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

必要な座屈荷重

Pcr = n·P = 3·30000 = 90000 N
KL = 0.5·4000 = 2000 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·90000·2000²/(π³·70000))^(1/4) = 57.08 mm

答え

d = 57.08 mm

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

15.放物運動

答え: v0 = 36.02 m/s

Projectile motion: with no air resistance, x = vxt and y = h + vyt − ½gt². The two directions are independent except that they share the same time.

水平到達距離の式を v₀ について解く

v₀ = √(g·R/sin 2θ) = √(9.81·85/sin 40°) = 36.02 m/s

答え

v0 = 36.02 m/s
16.放物運動

答え: x = 151.3 m

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

初速度を分解する

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

vx = v0·cos θ = 33·cos 25° = 29.91 m/s
vy = v0·sin θ = 33·sin 25° = 13.95 m/s

滞空時間

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

0 = 55 + 13.95·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (13.95 + √(194.5 + 1079))/9.81 = 5.06 s

水平距離

x = vx·t = 29.91·5.06 = 151.3 m

答え

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

答え: T = 227.2 N

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。 円軌道上では、中心に向かう合力は mv²/r に等しくなければなりません。

最下点での速さ

v² = 2gL(1 − cos θ) = 2·9.81·1·(1 − cos 70°) = 12.91 m²/s²

最下点での力

m·g = 10·9.81 = 98.1 N
T = mg + mv²/L = 98.1 + 10·12.91/1 = 227.2 N

答え

T = 227.2 N

確認:T = mg(3 − 2cos θ) = 227.2 N。

18.仕事とエネルギー

答え: N = 1.962 N

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。 円軌道上では、中心に向かう合力は mv²/r に等しくなければなりません。

最高点での速さ

v² = 2g(h − 2R) = 2·9.81·(10.4 − 8) = 47.09 m²/s²

最高点での力

m·g = 1·9.81 = 9.81 N
N = mv²/R − mg = 1·47.09/4 − 9.81 = 1.962 N

答え

N = 1.962 N

h = 2.5R のとき台車はぎりぎりで一周でき、最高点で N = 0 になります。

19.機械振動

答え: ζ = 0.05986

質量–ばね系の自由振動:減衰のない固有角振動数は ω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/5)·ln(50/7.6) = 0.3768

減衰比

ζ = δ/√(4π² + δ²) = 0.3768/√(39.48 + 0.142) = 0.05986

答え

ζ = 0.05986
20.機械振動

答え: fn = 2.449 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₂ と足し合わせます。

等価ばね定数

並列なのは k1 と k2 だけです。その組が k3 と直列なので、逆数で合成しましょう。

k12 = 25 + 20 = 45 kN/m
keq = 45·12/(45 + 12) = 9.474 kN/m

固有振動数

ωn = √(9474/40) = 15.39 rad/s
fn = 15.39/(2π) = 2.449 Hz

答え

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

答え: ηmax = 58.81 %

第 2 法則による限界:絶対温度 TH と TL の熱源の間で働くどんなサイクルも、可逆(カルノー)サイクルを上回れません:ηmax = 1 − TL/TH、COPR,max = TL/(TH − TL)、COPHP,max = TH/(TH − TL)。エネルギーは常に保存されます:QH = W + QL。

絶対温度

TH = 475 + 273.15 = 748.1 K
TL = 35 + 273.15 = 308.1 K

カルノー効率

ηmax = 1 − TL/TH = 1 − 308.1/748.1 = 58.81 %

答え

ηmax = 58.81 %
22.熱機関と冷凍機

答え: Sgen = 1.019 kJ/K

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

放出熱

QL = QH − W = 1600 − 752 = 848 kJ

絶対温度

TH = 650 + 273.15 = 923.1 K
TL = 35 + 273.15 = 308.1 K

熱源のエントロピー変化

ΔSH = −1600/923.1 = −1.733 kJ/K
ΔSL = 848/308.1 = 2.752 kJ/K
Sgen = 2.752 − 1.733 = 1.019 kJ/K

答え

Sgen = 1.019 kJ/K

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

23.ガスサイクル

答え: T3 = 1956 K

ディーゼルサイクル(圧縮着火機関):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)) です。

圧縮の終わり

T2 = 295·200.4 = 977.8 K

等圧加熱

T3 = rc·T2 = 2·977.8 = 1956 K

答え

T3 = 1956 K
24.ガスサイクル

答え: wnet = 274.6 kJ/kg

ブレイトンサイクル(ガスタービン):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 を出します。

等エントロピー変化後の温度

T2s = 288·160.2857 = 636 K
T4s = 1500/160.2857 = 679.3 K

圧縮機仕事とタービン仕事

wc = 1.005·(636 − 288)/0.82 = 426.5 kJ/kg
wt = 0.85·1.005·(1500 − 679.3) = 701.1 kJ/kg

正味仕事

wnet = 701.1 − 426.5 = 274.6 kJ/kg

答え

wnet = 274.6 kJ/kg

後退仕事比(バックワーク比):wc/wt = 60.83%。

25.連続の式とベルヌーイ

答え: p2 = 274.3 kPa

連続の式とベルヌーイの定理:非圧縮性流体では、どの断面でも流量は同じです:Q = A₁v₁ = A₂v₂。損失のない流線に沿って p + ½ρv² + ρgz は一定なので、流れが速いところや高いところほど圧力は低くなります。

連続の式とベルヌーイ

Pa で計算し、1000 で割って kPa にします。

v2 = v1·(d1/d2)² = 1·(75/40)² = 3.516 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 3.516²)/1000 = −5.68 kPa
p2 = p1 + ½ρ(v1² − v2²) = 280 + (−5.68) = 274.3 kPa

答え

p2 = 274.3 kPa
26.連続の式とベルヌーイ

答え: p2 = 218.2 kPa

連続の式とベルヌーイの定理:非圧縮性流体では、どの断面でも流量は同じです:Q = A₁v₁ = A₂v₂。損失のない流線に沿って p + ½ρv² + ρgz は一定なので、流れが速いところや高いところほど圧力は低くなります。

連続の式とベルヌーイ

Pa で計算し、1000 で割って kPa にします。

v2 = v1·(d1/d2)² = 1·(75/50)² = 2.25 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 2.25²)/1000 = −2.031 kPa
ρgΔz = 1000·9.81·1/1000 = 9.81 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 230 + (−2.031) − 9.81 = 218.2 kPa

答え

p2 = 218.2 kPa
27.管摩擦と損失水頭

答え: Δp = 495.7 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.022·(500/0.075)·2.6²/(2·9.81) = 50.53 m

圧力損失

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

Δp = ρ·g·hf = 1000·9.81·50.53 = 495733 Pa = 495.7 kPa

答え

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

答え: Δz = 0.6594 m

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

速度

V = Q/A = 0.08/0.07069 = 1.132 m/s
V²/(2g) = 0.06529 m

損失係数

f·L/D = 0.015·150/0.3 = 7.5
ΣK = 0.5 + 3·0.3 + 0.2 + 1.0 = 2.6

全損失水頭

hf = 7.5·0.06529 = 0.4896 m
hm = 2.6·0.06529 = 0.1697 m
Δz = hf + hm = 0.6594 m

答え

Δz = 0.6594 m
29.ポンプとポンプ動力

答え: Pshaft = 41.65 kW

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.06·46 = 27.08 kW

軸動力

Pshaft = Ph/η = 27.08/0.65 = 41.65 kW

答え

Pshaft = 41.65 kW
30.ポンプとポンプ動力

答え: Pshaft = 4.232 kW

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

速度

V = Q/A = 0.02/0.007854 = 2.546 m/s
V²/(2g) = 0.3305 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (25 + 3)·0.3305 = 9.254 m

ポンプ揚程

H = Δz + hL = 8 + 9.254 = 17.25 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.02·17.25/0.8 = 4.232 kW

答え

Pshaft = 4.232 kW
31.熱伝導

答え: Q̇ = 18.73 W

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

各層の熱抵抗

R1 = L1/(k1A) = (250/1000)/(0.04·6) = 1.042 K/W
R2 = L2/(k2A) = (150/1000)/(0.96·6) = 0.02604 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 20/1.068 = 18.73 W

答え

Q̇ = 18.73 W
32.熱伝導

答え: Ts = 19.18 °C

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(6·4) = 0.04167 K/W
R1 = L1/(k1A) = (250/1000)/(0.17·4) = 0.3676 K/W
R2 = L2/(k2A) = (150/1000)/(0.17·4) = 0.2206 K/W
Ro = 1/(hoA) = 1/(20·4) = 0.0125 K/W
ΣR = 0.6424 K/W

熱損失

Q̇ = (21 − (-7))/0.6424 = 43.59 W
Ts = Troom − Q̇·Ri = 21 − 43.59·0.04167 = 19.18 °C

答え

Ts = 19.18 °C
33.歯車と歯車列

答え: nw = 24 rpm

Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.

ホイールの回転速度

nw = n·z/Nw = 960·2/80 = 24 rpm

答え

nw = 24 rpm
34.歯車と歯車列

答え: nout = 195.3 rpm

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

第 1 段

nB = 1500·30/48 = 937.5 rpm = nC

第 2 段

nout = 937.5·15/72 = 195.3 rpm

答え

nout = 195.3 rpm

全体の速比 = 7.68 : 1

35.ラプラス変換

答え: f(1) = 2.165

ラプラス変換: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² + 4s + 4 = (s + 2)²

部分分数展開

5s + 21 = 5(s + 2) + 11
F(s) = 5/(s + 2) + 11/(s + 2)²

逆変換

f(t) = 5·e−2t + 11·t·e−2t
f(1) = (5 + 11·1)·e−2·1 = 2.165

答え

f(1) = 2.165
36.ラプラス変換

答え: y(0.1) = −4.533

ラプラス変換: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²Y − (−5)s − (−1) + 11(sY − (−5)) + 30Y = 0
Y(s) = (−5s − 56)/((s + 5)(s + 6))

部分分数展開

A = −31/1 = −31
B = −26/(−1) = 26

結果

y(t) = −31·e−5t + 26·e−6t
y(0.1) = −4.533

答え

y(0.1) = −4.533

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

37.ブロック線図

答え: ess = 0.1852

ブロック線図の代数:直列のブロックは掛け算(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) です。

速度偏差定数

Kᵥ = lims→0 s·G(s) = 54/2 = 27

定常偏差

ess = A/Kᵥ = 5/27 = 0.1852

答え

ess = 0.1852
38.ブロック線図

答え: css = 7.333

ブロック線図の代数:直列のブロックは掛け算(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) です。

伝達関数

C/D = G/(1 + K·G) = 2/(s + 9 + 24)
C/R = K·G/(1 + K·G) = 24/(s + 9 + 24)

定常状態

css = 10·24/33 + 1·2/33 = 7.333

答え

css = 7.333

K を大きくすると外乱の影響は小さくなりますが、比例コントローラでは 0 にはなりません。

39.ステップ応答

答え: σ = −2.5

ステップ応答: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 + ωₙ² と一致しなければなりません。

ωₙ² = 51 → ωₙ = 7.141 rad/s
2ζωₙ = 5 → ζ = 5/(2·7.141) = 0.3501

結果

σ = −ζωₙ = −5/2 = −2.5

答え

σ = −2.5
40.ステップ応答

答え: ts = 4 s

ステップ応答: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) = 60/(s² + 2s + 60): ωₙ = 7.746 rad/s, ζ = 0.1291

整定時間

ζωₙ = 2/2 = 1
ts ≈ 4/(ζωₙ) = 4 s

答え

ts = 4 s

このループでは K によらず ζωₙ = a/2 です:K を大きくするとオーバーシュートは増えますが、整定時間は短くなりません。

41.定常偏差と PID 制御

答え: ess = 2.564 %

定常偏差:開ループ伝達関数 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 型のループはステップ入力に偏差なく追従し、ランプ入力には有限の偏差で追従します。

位置偏差定数

Kp = 19·2 = 38

ステップ偏差

ess = 100/(1 + 38) = 2.564 %

答え

ess = 2.564 %
42.定常偏差と PID 制御

答え: tp = 0.3245 s

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(s + 5) + 25·4 = s² + 5s + 100

ωₙ と ζ

ωₙ = √100 = 10 rad/s
ζ = 5/(2·10) = 0.25

ピーク時間

ωd = 10·√(1 − 0.25²) = 9.682 rad/s
tp = π/9.682 = 0.3245 s

答え

tp = 0.3245 s

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

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