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​ = 20 kNP2​ = 40 kNABRA​ = ?RB​1 m0.5 m4.5 m

単純ばりは A がピン支点、B がローラー支点である。鉛直支点反力 RA を求めよ。

L = 6 mP1 = 20 kN 位置 1 mP2 = 40 kN 位置 1.5 m
RA = kN
2.はりの支点反力むずかしい
w = 8 kN/mP = 25 kNABRA​ = ?RB​1 m1 m2 m

単純ばりは A がピン支点、B がローラー支点である。鉛直支点反力 RA を求めよ。

L = 4 mP = 25 kN 位置 2 mw = 8 kN/m(1 から 2 m まで)
RA = kN
3.軸応力とひずみふつう
100 kNØ4070 GPaØ80/50800 mmδ = ?

鋼 の棒(直径 40 mm、E = 200 GPa)が、アルミニウム でできた管(外径 80 mm、内径 50 mm、E = 70 GPa)の中に入っています。どちらも長さ 800 mm で、剛な板の間で 100 kN によって一緒に圧縮されます。 縮み δ を求めなさい。

P = 100 kN棒: Ø40 mm, E = 200 GPa管: Ø80/Ø50 mm, E = 70 GPaL = 800 mm
δ = mm
4.軸応力とひずみむずかしい
P = ?Ø20100 GPaØ120/60500 mm

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

棒: Ø20 mm, E = 200 GPa, σallow = 160 MPa管: Ø120/Ø60 mm, E = 100 GPa, σallow = 60 MPa
Pmax = kN
5.軸のねじりふつう
T = 200 N·md = 25 mmL = 1.2 mφ = ?

鋼製の中実軸が壁に固定され、自由端でトルク T によりねじられる。自由端のねじれ角 φ を度で求めよ。

T = 200 N·md = 25 mmL = 1.2 mG = 80 GPa
φ = °
6.軸のねじりむずかしい
Tdo​ = 75 mm, di​ = 38 mm

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

P = 20 kWN = 500 rpmdo = 75 mmdi = 38 mm
τmax = MPa
7.せん断力と曲げモーメントふつう
P1​ = 70 kNP2​ = 40 kNABC1.5 m4.5 m1.5 m

はりは A でピン支点、B でローラー支点に支持され、B から自由端 C まで 1.5 m 張り出している。支点間に集中荷重 P1、C に集中荷重 P2 が作用する。 支点 B における曲げモーメント MB を求めよ。 下側が引張となる曲げモーメントを正、左側断面で上向きのせん断力を正とする。

L = 7.5 mP1 = 70 kN 位置 1.5 mP2 = 40 kN 位置 7.5 m支点の位置 x = 0 と x = 6 m
MB = kN·m
8.せん断力と曲げモーメントむずかしい
w = 10 kN/mABC5 m1 m

はりは A でピン支点、B でローラー支点に支持され、B から自由端 C まで 1 m 張り出している。全長 A–C に等分布荷重 w が作用する。 A から支点間の反曲点(曲げモーメントの符号が変わる点)までの距離 xc を求めよ。

L = 6 mw = 10 kN/m(0 から 6 m まで)支点の位置 x = 0 と x = 5 m
xc = m
9.モールの応力円と応力の座標変換ふつう
σx​ = 40σy​ = 100τxy​ = −50MPa

薄板内の点は平面応力状態にある:σₓ = 40 MPa、σᵧ = 100 MPa、τₓᵧ = −50 MPa、σz = 0。その点での絶対最大せん断応力 τabs を求めよ(3 次元のすべての面を考える)。

σx = 40 MPaσy = 100 MPaτxy = −50 MPaσz = 0
τabs = MPa
10.モールの応力円と応力の座標変換むずかしい
σx​ = −20σy​ = 70τxy​ = −20MPa

薄板内の点は平面応力状態にある:σₓ = −20 MPa、σᵧ = 70 MPa、τₓᵧ = −20 MPa(σz = 0)。 材料の降伏強さは 275 MPa である。ミーゼス(せん断ひずみエネルギー)説による降伏に対する安全率 n を求めよ。

σx = −20 MPaσy = 70 MPaτxy = −20 MPaσY = 275 MPa
n =
11.はりのたわみふつう
P = 20 kNABδ2 m0.5 m

長さ L の片持ちばりは A で固定され、A から距離 a = 2 m の位置に集中荷重 P が作用する。 自由端 B のたわみ δB を求めよ。

P = 20 kNa = 2 mL = 2.5 mE = 70 GPa (アルミニウム)I = 84 × 10⁶ mm⁴
δB = mm
12.はりのたわみむずかしい
w = 15 kN/mAB5 m

一端固定・他端支持ばりは A で固定され、B で単純支持されている。全長に等分布荷重 w が作用する。EI は一定である。 固定端における曲げモーメントの大きさ MA を求めよ。

L = 5 mw = 15 kN/m
MA = kN·m
13.柱の座屈ふつう
PL = 4 md = 30 mm断面

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

L = 4 mE = 70 GPa (アルミニウム)d = 30 mm
Pcr = kN
14.柱の座屈むずかしい
P = 45 kNL = 5 md = ?断面

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

P = 45 kNn = 3L = 5 mE = 70 GPa (アルミニウム)
d = mm
15.放物運動ふつう
v0​ = 14 m/s30°到達距離H = ?

ボールを 平らな地面から水平から 30° 上向きに投げる。空気抵抗は無視する。ボールが達する最高高さ H を求めよ。

v0 = 14 m/sθ = 30°g = 9.81 m/s²
H = m
16.放物運動むずかしい
9 mv0​ = 21 m/s70°x = ?

地面の高さからボールを 21 m/s、水平から 70° 上向きに投げる。ボールは下降中に、高さ 9 m の建物の平らな屋上に着地する。空気抵抗を無視し、g = 9.81 m/s² とする。投げた点から着地点までの水平距離 x を求めよ。

v₀ = 21 m/sθ = 70°hroof = 9 mg = 9.81 m/s²
x = m
17.仕事とエネルギーふつう
90 km/h8°1000 kg

質量 1000 kg の自動車が 8° の坂を一定の 90 km/h で上る。空気抵抗と転がり抵抗の合計は 800 N である。車輪が出さなければならない動力 P を kW で求めよ。

m = 1000 kgθ = 8°Fr = 800 Nv = 90 km/hg = 9.81 m/s²
P = kW
18.仕事とエネルギーむずかしい
8 kg5 kgμ = 0.25d = 2.5 m

質量 8 kg のブロックがあらい水平な机の上にある(μ = 0.25)。軽い糸がそこから机の端のなめらかな滑車を通って、つり下げられた質量 5 kg のブロックにつながっている。系を静かに放す。つり下げられたブロックが 2.5 m 下がった後のブロックの速さ v を求めよ。

m1 = 8 kgm2 = 5 kgμ = 0.25d = 2.5 mg = 9.81 m/s²
v = m/s
19.機械振動ふつう
k = 80 kN/mζ = 0.4m = 12 kg

ばね定数 80 kN/m のばねにつながれた 12 kg の質量に、減衰比 0.4 となる粘性ダンパがついている。自由振動の減衰固有振動数 fd を Hz で求めよ。

m = 12 kgk = 80 kN/mζ = 0.4
fd = Hz
20.機械振動むずかしい
k1​ = 5k2​ = 10k3​ = 10m = 25 kg

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

m = 25 kgk1 = 5 kN/mk2 = 10 kN/mk3 = 10 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 750 °CTL​ = 15 °CEQH​ = 1500 kJQL​W = ?

熱機関が 750 °C の高温熱源と 15 °C の低温熱源の間で働く。 この機関は 1500 kJ の熱を受け取る。得られる最大の仕事 Wmax を求めよ。

TH = 750 °CTL = 15 °CQH = 1500 kJ
Wmax = kJ
22.熱機関と冷凍機むずかしい
TH​ = 21 °CTL​ = 5 °CHPQ̇H​ = 9 kWQ̇L​Ẇ = ?

ヒートポンプが、外気 5 °C のときに家を 21 °C に保ち、家に 9 kW の熱を供給する。 必要な最小入力動力 Ẇmin を求めよ。

TH = 21 °CTL = 5 °CQ̇H = 9 kW
Ẇmin = kW
23.ガスサイクルふつう
vP1234qin​qout​

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

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

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

rp = 20k = 1.4
η = %
25.連続の式とベルヌーイふつう
12d1​ = 75 mm, p1​ = 190 kPad2​ = 50 mm, p2​ = ?

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

d1 = 75 mmd2 = 50 mmp1 = 190 kPav1 = 2.25 m/sρ = 1000 kg/m³
p2 = kPa
26.連続の式とベルヌーイむずかしい
C210.5 m1 md = 25 mmQ = ?

内径 25 mm のホースで開放水槽から水をサイホンで吸い出す。ホースの最高点 C は水面より 0.5 m 上にあり、出口(大気に開放)は水面より 1 m 下にある。損失を無視し、g = 9.81 m/s² とする。 流量 Q を求めよ。

d = 25 mmzC = +0.5 mz2 = −1 mg = 9.81 m/s²
Q = L/s
27.管摩擦と損失水頭ふつう
Q = 20 L/sD = 150 mmL = 1000 mf = 0.025

水が、長さ 1000 m、内径 150 mm の管内を 20 L/s で流れている。 ダルシーの管摩擦係数 f = 0.025、g = 9.81 m/s² とする。摩擦損失水頭 hf を求めよ。

Q = 20 L/sL = 1000 mD = 150 mmf = 0.025g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 5Δz = ?L = 150 m, D = 200 mm, f = 0.0153× 曲がり (K = 0.3)

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

Q = 40 L/sL = 150 mD = 200 mmf = 0.015Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 5Kexit = 1.0
Δz = m
29.ポンプとポンプ動力ふつう
PQ = 50 L/sH = 64 m

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

Q = 50 L/sH = 64 mη = 70%g = 9.81 m/s²
Pshaft = kW
30.ポンプとポンプ動力むずかしい
PQ = 75 L/sΔz = 20 mL = 150 m, D = 200 mm, f = 0.025, ΣK = 5

ポンプが 75 L/s の水(ρ = 1000 kg/m³)を、開放された下の貯水槽から、水面が 20 m 高い上のタンクへ送る。管は長さ 150 m、内径 200 mm、ダルシーの管摩擦係数 f = 0.025 で、局所損失係数の合計は ΣK = 5(入口と出口を含む)である。g = 9.81 m/s² とする。 ポンプが与えなければならない揚程 H を求めよ。

Q = 75 L/sΔz = 20 mL = 150 mD = 200 mmf = 0.025ΣK = 5
H = m
31.熱伝導ふつう
れんが, k = 0.72100 mm発泡ポリスチレン, k = 0.035150 mm20 °C10 °CTi​ = ?Q̇

壁はれんが 100 mm と発泡ポリスチレン 150 mm からなり、面積は 8 m² である。外側の両面は 20 °C と 10 °C である。 2 つの層の境界面の温度 Ti を求めよ。

L1 = 100 mm, k1 = 0.72 W/(m·K)L2 = 150 mm, k2 = 0.035 W/(m·K)A = 8 m²T1 = 20 °CT3 = 10 °C
Ti = °C
32.熱伝導むずかしい
ロックウール, k = 0.04250 mmコンクリート, k = 1.4100 mm室内 22 °C屋外 -4 °CTs​ = ?Q̇

250 mm のロックウールと 100 mm のコンクリートからなる壁が、22 °C の部屋と -4 °C の外気を隔てている。熱伝達率は室内側 8 W/(m²·K)、屋外側 25 W/(m²·K)、壁の面積は 10 m² である。 壁の内表面温度 Ts を求めよ。

L1 = 250 mm, k1 = 0.04 W/(m·K)L2 = 100 mm, k2 = 1.4 W/(m·K)A = 10 m²hi = 8 W/(m²·K)ho = 25 W/(m²·K)Troom = 22 °CTout = -4 °C
Ts = °C
33.歯車と歯車列ふつう
N1​ = ?N2​

2 本の平行軸の間隔は 70 mm である。これらをモジュール 2 mm、減速比 2.5 : 1 の平歯車対で連結しなければならない。 ピニオンの歯数 N1 はいくつ必要か。

a = 70 mmm = 2 mmi = 2.5
N1 =
34.歯車と歯車列むずかしい
A = 45B = 80C = 50D = 60

2900 rpm の 4 kW のモータが歯車 A(歯数 45)を駆動し、A は歯車 B(歯数 80)とかみ合う。B と同じ軸にある歯車 C(歯数 50)が出力軸の歯車 D(歯数 60)を駆動する。各かみ合いの効率は 96% である。出力トルク Tout を求めよ。

P = 4 kWnin = 2900 rpmA = 45B = 80C = 50D = 60ηmesh = 96%
Tout = N·m
35.ラプラス変換ふつう
ReIm−1−2−3−4−5−61

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

F(s) = (2s + 7)/(s² + 10s + 24)
A =
36.ラプラス変換むずかしい
ReIm−2−4−6−82

ある系が微分方程式 y″ + 2y′ + 10y = 2u′ + 16u で表される(入力 u(t)、出力 y(t))。その伝達関数 G(s) = Y(s)/U(s)(初期条件はゼロ)を求め、t = 0.8 s におけるインパルス応答 g(t) を計算せよ。g(0.8) を求めよ。

y″ + 2y′ + 10y = 2u′ + 16ut = 0.8 s
g(0.8) =
37.ブロック線図ふつう
R(s)+39s + 1G(s)C(s)0.2H−

プラント G(s) = 39/(s + 1) を、センサのゲインが H = 0.2 の負帰還ループで制御する。 R に大きさ 5 のステップを加える。定常出力 css を求めよ。

G(s) = 39/(s + 1)H = 0.2R = 5
css =
38.ブロック線図むずかしい
R(s)+K(s + 2)(s + 3)(s + 6)G(s)C(s)−

プラント G(s) = K/((s + 2)(s + 3)(s + 6)) を単位負帰還ループに入れる。 安定な最大ゲインでは、ループは振動します。振動の周波数 ω を rad/s で求めよ。

G(s) = K/((s + 2)(s + 3)(s + 6))単位負帰還
ω = rad/s
39.ステップ応答ふつう
t (s)y(t)2?

ある系の伝達関数は G(s) = 146/(s² + 13s + 73) で、単位ステップ入力を受ける。 ステップ応答のオーバーシュート率 %OS を求めよ。

G(s) = 146/(s² + 13s + 73)
%OS = %
40.ステップ応答むずかしい
+−12 VR = 28 ΩL = 50 mHC = 10 µFvC​

R = 28 Ω、L = 50 mH、C = 10 µF の直列 RLC 回路に、t = 0 で 12 V のステップを加える。コンデンサは充電されておらず、ステップの前には電流は流れていない。 電圧の最初のピークの時刻 tp を ms で求めよ。

R = 28 ΩL = 50 mHC = 10 µFV = 12 V
tp = ms
41.定常偏差と PID 制御ふつう
r(t) = 2 + 1t+2(s + 4)s(s + 1)(s + 4)C(s)−

入力 r(t) = 2 + 1t(ステップとランプの和)が、G(s) = 2(s + 4)/(s(s + 1)(s + 4)) の単位フィードバックのループに加えられます。定常偏差 ess を求めなさい。

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

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

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

3溶液

1.はりの支点反力

答え: RA = 46.67 kN

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。一方の支点まわりのモーメントをとると、その支点反力が式から消えます。

P1​ = 20 kNP2​ = 40 kNABRA​RB​1 m0.5 m4.5 m

自由体図

支点を反力に置き換えます。ローラーは鉛直反力を 1 つ与えます。鉛直荷重のみなら、ピンも鉛直反力を 1 つ与えます。

つり合い

A まわりのモーメント(反時計回りを正)。次に鉛直方向の力のつり合い。

ΣMA = 0: RB·6 = 20·1 + 40·1.5
RB = 80/6 = 13.33 kN
ΣFy = 0: RA = 60 − 13.33 = 46.67 kN

答え

RA = 46.67 kN
2.はりの支点反力

答え: RA = 17.5 kN

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。一方の支点まわりのモーメントをとると、その支点反力が式から消えます。

w = 8 kN/mP = 25 kNABRA​RB​1 m1 m2 m

自由体図

支点を反力に置き換えます。ローラーは鉛直反力を 1 つ与えます。鉛直荷重のみなら、ピンも鉛直反力を 1 つ与えます。

つり合い

A まわりのモーメント(反時計回りを正)。次に鉛直方向の力のつり合い。

分布荷重の合力:W = w·ℓ = 8·1 = 8 kN、作用点 x̄ = 1.5 m
ΣMA = 0: RB·4 = 25·2 + 8·1·1.5
RB = 62/4 = 15.5 kN
ΣFy = 0: RA = 33 − 15.5 = 17.5 kN

答え

RA = 17.5 kN
3.軸応力とひずみ

答え: δ = 0.1718 mm

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

面積

Ar = π·40²/4 = 1257 mm²
At = π(80² − 50²)/4 = 3063 mm²

合成した軸剛性

ΣEA = 200·10³·1257 + 70·10³·3063 = 4.657 × 10⁸ N

縮み

δ = 100000·800/4.657 × 10⁸ = 0.1718 mm

答え

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

答え: Pmax = 546.6 kN

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

面積

Ar = π·20²/4 = 314.2 mm²
At = π(120² − 60²)/4 = 8482 mm²

合成した軸剛性

ΣEA = 200·10³·314.2 + 100·10³·8482 = 9.111 × 10⁸ N

支配するひずみ

ε = min(0.0008, 0.0006) = 0.0006 (黄銅)

最大荷重

Pmax = 0.0006·9.111 × 10⁸ = 546637 N = 546.6 kN

答え

Pmax = 546.6 kN
5.軸のねじり

答え: φ = 4.482 °

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

断面二次極モーメント

J = π·25⁴/32 = 38350 mm⁴

ねじれ角

φ = TL/(GJ) = 200·1000·1200/(80000·38350) = 0.07823 rad
φ = 0.07823·180/π = 4.482°

答え

φ = 4.482°
6.軸のねじり

答え: τmax = 4.937 MPa

丸軸のねじり:せん断応力は中心で 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/ω = 20000/52.36 = 382 N·m

断面二次極モーメント

J = π(75⁴ − 38⁴)/32 = 2.902 × 10⁶ mm⁴

最大せん断応力

τmax = T·(do/2)/J = 381972·37.5/2.902 × 10⁶ = 4.937 MPa

答え

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

答え: MB = −60 kN·m

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

支点反力

はり全体のつり合いを考えます。張出し部の荷重も A まわりのモーメントに寄与します。

ΣMA = 0: RB·6 = 70·1.5 + 40·7.5
RB = 405/6 = 67.5 kN
ΣFy = 0: RA = 110 − 67.5 = 42.5 kN

B で切断し、張出し部を残す

張出し部には P2 だけが腕の長さ a で作用します。これにより、はりは支点上で上に凸に曲がります(負の曲げ)。

MB = −P2·a = −40·1.5 = −60 kN·m

答え

MB = −60 kN·m
8.せん断力と曲げモーメント

答え: xc = 4.8 m

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

支点反力

張出し部の荷重を含めた、はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 10·6 = 60 kN、作用点 x̄ = 3 m
ΣMA = 0: RB·5 = 10·6·3
RB = 180/5 = 36 kN
ΣFy = 0: RA = 60 − 36 = 24 kN

モーメント 0 の位置

支点間では M(x) = RA·x − w·x²/2。x = 0 を除くと、x = 2RA/w で 0 になります。

RA·xc − w·xc²/2 = 0
xc = 2·24/10 = 4.8 m

答え

xc = 4.8 m
9.モールの応力円と応力の座標変換

答え: τabs = 64.15 MPa

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

στxyCσ1​σ2​

面内の主応力

C = 70 MPa, R = √((−30)² + (−50)²) = 58.31 MPa
σ1 = 128.3 MPa, σ2 = 11.69 MPa, σz = 0

3 つの円のうち最大のもの

τ = max(|σ1 − σ2|, |σ1|, |σ2|)/2 = max(116.6, 128.3, 11.69)/2 = 64.15 MPa

答え

τabs = 64.15 MPa

ここでは面外の円で決まります:絶対最大値は面内最大値 R より大きくなります。

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

答え: n = 3.094

平面応力の座標変換:θ だけ回転した面では σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。モールの応力円はすべての面を一度に表します:中心は C = (σₓ + σᵧ)/2、半径は R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主応力は C ± R、面内の最大せん断応力は R です。 降伏条件:ミーゼス説では σvM = √(σ₁² − σ₁σ₂ + σ₂²) が σY に達したとき、トレスカ説では絶対最大せん断応力が σY/2 に達したときに降伏すると予測します。

στxyCσ1​σ2​

主応力

C = 25 MPa, R = 49.24 MPa
σ1 = 74.24 MPa, σ2 = −24.24 MPa, σ3 = 0

ミーゼス

σvM = √(74.24² − 74.24·(−24.24) + (−24.24)²) = 88.88 MPa
n = 275/88.88 = 3.094

答え

n = 3.094

n > 1 はその点が降伏しないことを意味します。トレスカはミーゼスより少し安全側です。

11.はりのたわみ

答え: δB = 12.47 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 = 20000 N, a = 2000 mm, L = 2500 mm
EI = 5.88 × 10¹² N·mm²

自由端のたわみ

荷重位置のたわみに、荷重位置のたわみ角 × 残りの長さを加えます。

δP = Pa³/(3EI) = 9.07 mm
θP = Pa²/(2EI) = 0.006803 rad
δB = δP + θP·(L − a) = Pa²(3L − a)/(6EI) = 12.47 mm

答え

δB = 12.47 mm
12.はりのたわみ

答え: MA = 46.88 kN·m

はりのたわみ:たわみは曲率 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·5/8 = 28.13 kN

固定端のモーメント

荷重と支点 B の反力による A まわりのモーメント。

MA = wL²/2 − RB·L = 187.5 − 140.6 = 46.88 kN·m

答え

MA = 46.88 kN·m
13.柱の座屈

答え: Pcr = 3.504 kN

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

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

円形断面はどの軸まわりでも I が同じです。

I = πd⁴/64 = π·30⁴/64 = 39761 mm⁴

有効座屈長さ

KL = 0.7·4000 = 2800 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·39761/2800² = 3504 N = 3.504 kN

答え

Pcr = 3.504 kN
14.柱の座屈

答え: d = 141.3 mm

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

必要な座屈荷重

Pcr = n·P = 3·45000 = 135000 N
KL = 2·5000 = 10000 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·135000·10000²/(π³·70000))^(1/4) = 141.3 mm

答え

d = 141.3 mm

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

15.放物運動

答え: H = 2.497 m

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

初速度を分解する

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

vx = v0·cos θ = 14·cos 30° = 12.12 m/s
vy = v0·sin θ = 14·sin 30° = 7 m/s

最高点

最高点では鉛直速度が 0:vy² = 2gH。

H = vy²/(2g) = 7²/(2·9.81) = 2.497 m

答え

H = 2.497 m
16.放物運動

答え: x = 25.13 m

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.

初速度を分解する

vx = 7.182 m/s, vy = 19.73 m/s

屋上に着地するまでの時間

大きい方の解:ボールは下降中です。

9 = 19.73·t − 4.905·t²
t = (19.73 + √(19.73² − 2·9.81·9))/9.81 = 3.499 s

水平距離

x = vx·t = 7.182·3.499 = 25.13 m

答え

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

答え: P = 54.13 kW

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

駆動力

m·g = 1000·9.81 = 9810 N
mg·sin θ = 9810·sin 8° = 1365 N
F = 1365 + 800 = 2165 N

電力

v = 90/3.6 = 25 m/s
P = F·v = 2165·25 = 54132 W = 54.13 kW

答え

P = 54.13 kW
18.仕事とエネルギー

答え: v = 3.364 m/s

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

放出されるエネルギーと失われるエネルギー

m₂g·d = 49.05·2.5 = 122.6 J
μ·m₁g·d = 0.25·78.48·2.5 = 49.05 J

両方のブロックの運動エネルギー

½(m₁ + m₂)v² = 73.58 → v = √(2·73.58/13) = 3.364 m/s

答え

v = 3.364 m/s
19.機械振動

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

減衰のない固有振動数

ωn = √(80000/12) = 81.65 rad/s
fn = 81.65/(2π) = 12.99 Hz

減衰固有振動数

fd = fn·√(1 − ζ²) = 12.99·√(1 − 0.4²) = 11.91 Hz

答え

fd = 11.91 Hz
20.機械振動

答え: fn = 2.466 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 = 5 + 10 = 15 kN/m
keq = 15·10/(15 + 10) = 6 kN/m

固有振動数

ωn = √(6000/25) = 15.49 rad/s
fn = 15.49/(2π) = 2.466 Hz

答え

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

答え: Wmax = 1078 kJ

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

絶対温度

TH = 750 + 273.15 = 1023 K
TL = 15 + 273.15 = 288.1 K

カルノー効率

ηmax = 1 − TL/TH = 1 − 288.1/1023 = 71.84 %

最大仕事

Wmax = ηmax·QH = 0.7184·1500 = 1078 kJ

答え

Wmax = 1078 kJ
22.熱機関と冷凍機

答え: Ẇmin = 0.4895 kW

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

絶対温度

TH = 21 + 273.15 = 294.1 K
TL = 5 + 273.15 = 278.1 K

可能な最良の COP

COPHP,max = TH/(TH − TL) = 294.1/(294.1 − 278.1) = 18.38

最小動力

必要な仕事が最も少ないのは、COP が最大の可逆機械です。

Ẇmin = Q̇H/COP = 9/18.38 = 0.4895 kW

答え

Ẇmin = 0.4895 kW
23.ガスサイクル

答え: T3 = 2582 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 = 300·220.4 = 1033 K

等圧加熱

T3 = rc·T2 = 2.5·1033 = 2582 K

答え

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

答え: η = 57.51 %

ブレイトンサイクル(ガスタービン):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/200.2857 = 1 − 1/2.354 = 57.51%

答え

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

答え: p2 = 179.7 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 2.25·(75/50)² = 5.063 m/s
½ρ(v1² − v2²) = ½·1000·(2.25² − 5.063²)/1000 = −10.28 kPa
p2 = p1 + ½ρ(v1² − v2²) = 190 + (−10.28) = 179.7 kPa

答え

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

答え: Q = 2.174 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.

ホース内の流速

水面と出口はどちらも大気に開放されているので、圧力の項は打ち消し合います。

v = √(2g·1) = √(2·9.81·1) = 4.429 m/s

流量

Q = A·v = 0.0004909·4.429 = 0.002174 m³/s = 2.174 L/s

答え

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

答え: hf = 10.88 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 を含む)が得られます。

平均流速

A = π·0.15²/4 = 0.01767 m²
V = Q/A = 0.02/0.01767 = 1.132 m/s

速度水頭

V²/(2g) = 1.132²/(2·9.81) = 0.06529 m

ダルシー・ワイスバッハの式

hf = f·(L/D)·V²/(2g) = 0.025·(1000/0.15)·1.132²/(2·9.81) = 10.88 m

答え

hf = 10.88 m
28.管摩擦と損失水頭

答え: Δz = 1.541 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.04/0.03142 = 1.273 m/s
V²/(2g) = 0.08263 m

損失係数

f·L/D = 0.015·150/0.2 = 11.25
ΣK = 0.5 + 3·0.3 + 5 + 1.0 = 7.4

全損失水頭

hf = 11.25·0.08263 = 0.9296 m
hm = 7.4·0.08263 = 0.6114 m
Δz = hf + hm = 1.541 m

答え

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

答え: Pshaft = 44.85 kW

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.05·64 = 31.39 kW

軸動力

Pshaft = Ph/η = 31.39/0.7 = 44.85 kW

答え

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

答え: H = 26.9 m

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

速度

V = Q/A = 0.075/0.03142 = 2.387 m/s
V²/(2g) = 0.2905 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (18.75 + 5)·0.2905 = 6.899 m

ポンプ揚程

どちらの水面も開放されていてほぼ静止しているので、圧力と速度の項は消えます。

H = Δz + hL = 20 + 6.899 = 26.9 m

答え

H = 26.9 m
31.熱伝導

答え: Ti = 19.69 °C

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

各層の熱抵抗

R1 = L1/(k1A) = (100/1000)/(0.72·8) = 0.01736 K/W
R2 = L2/(k2A) = (150/1000)/(0.035·8) = 0.5357 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 10/0.5531 = 18.08 W

境界面温度

最初の層で温度は Q̇·R₁ だけ下がります。

Ti = T1 − Q̇·R1 = 20 − 18.08·0.01736 = 19.69 °C

答え

Ti = 19.69 °C
32.熱伝導

答え: Ts = 21.5 °C

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(8·10) = 0.0125 K/W
R1 = L1/(k1A) = (250/1000)/(0.04·10) = 0.625 K/W
R2 = L2/(k2A) = (100/1000)/(1.4·10) = 0.007143 K/W
Ro = 1/(hoA) = 1/(25·10) = 0.004 K/W
ΣR = 0.6486 K/W

熱損失

Q̇ = (22 − (-4))/0.6486 = 40.08 W
Ts = Troom − Q̇·Ri = 22 − 40.08·0.0125 = 21.5 °C

答え

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

答え: N1 = 20

Gear geometry: the module m sets the tooth size. A gear with N teeth has a pitch diameter d = m·N, two meshing gears have the same module, and their centre distance is a = (d₁ + d₂)/2 = m·(N₁ + N₂)/2. The pitch-line speed is v = π·d·n/60, and the tangential tooth force follows from the torque: Ft = 2T/d.

歯数の合計

N1 + N2 = 2a/m = 2·70/2 = 70

比で分ける

N1 = 70/(1 + 2.5) = 20

答え

N1 = 20
34.歯車と歯車列

答え: Tout = 25.9 N·m

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.

出力回転数

nout = 2900·(45/80)·(50/60) = 1359 rpm
ωout = 2π·1359/60 = 142.4 rad/s

出力動力とトルク

Pout = 4000·0.96² = 3686 W
Tout = Pout/ωout = 3686/142.4 = 25.9 N·m

答え

Tout = 25.9 N·m
35.ラプラス変換

答え: A = −0.5

ラプラス変換: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 + 24 = (s + 4)(s + 6)
F(s) = A/(s + 4) + B/(s + 6)

隠蔽法

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

A = (2s + 7)/(s + 6) |s = −4 = −1/2 = −0.5

答え

A = −0.5
(B = 2.5 なので、f(t) = −0.5·e−4t + 2.5·e−6t)
36.ラプラス変換

答え: g(0.8) = 0.7537

ラプラス変換: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)。

伝達関数

G(s) = Y(s)/U(s) = (2s + 16)/(s² + 2s + 10)

部分分数展開

s² + 2s + 10 = (s + 1)² + 3²
G(s) = 2·(s + 1)/((s + 1)² + 3²) + (14/3)·3/((s + 1)² + 3²)

インパルス応答

g(t) = e−t·(2·cos 3t + 4.667·sin 3t)
g(0.8) = 0.7537

答え

g(0.8) = 0.7537
37.ブロック線図

答え: css = 22.16

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

閉ループ伝達関数

T(s) = G/(1 + G·H) = 39/(s + 1 + 7.8) = 39/(s + 8.8)

結果

css = R·T(0) = 5·39/8.8 = 22.16

答え

css = 22.16
38.ブロック線図

答え: ω = 6 rad/s

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

固有方程式

s³ + 11s² + 36s + 36 + K = 0

ラウス表

s³ | 1 36
s² | 11 36 + K
s¹ | (11·36 − 36 − K)/11
s⁰ | 36 + K

安定限界

第 1 列のすべての要素が正でなければなりません。

Kmax = 11·36 − 36 = 360

振動の周波数

11·(jω)² + 36 + 360 = 0 → ω² = 396/11 = 36
ω = 6 rad/s

答え

ω = 6 rad/s
39.ステップ応答

答え: %OS = 2.516 %

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

ωₙ² = 73 → ωₙ = 8.544 rad/s
2ζωₙ = 13 → ζ = 13/(2·8.544) = 0.7608

結果

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−3.682 = 2.516 %

答え

%OS = 2.516 %
40.ステップ応答

答え: tp = 2.266 ms

ステップ応答: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) です。

固有角周波数と減衰比

ωₙ = 1/√(LC) = 1/√(0.05·1 × 10⁻⁵) = 1414 rad/s
ζ = (R/2)·√(C/L) = (28/2)·√(1 × 10⁻⁵/0.05) = 0.198

結果

ωd = 1386 rad/s
tp = π/ωd = 2.266 ms

答え

tp = 2.266 ms
41.定常偏差と PID 制御

答え: ess = 0.5

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

速度偏差定数

Kv = lim s·G(s) = 2

ステップ成分

Kp = ∞ → estep = 0

ランプ偏差

ess = 1/2 = 0.5

答え

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

答え: ess = 10

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 + 0.5 → s → 0 のとき 0.5。

Kv = 0.5·4/10 = 0.2

ランプ偏差

ess = 2/0.2 = 10

答え

ess = 10

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

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

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