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​ = 15 kNP2​ = 40 kNABRA​ = ?RB​2.5 m0.5 m1 m

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

L = 4 mP1 = 15 kN 位置 2.5 mP2 = 40 kN 位置 3 m
RA = kN
2.はりの支点反力むずかしい
w0​ = 12 kN/mP = 65 kNABMA​ = ?5 m3 m

片持ちばりは A で壁に固定されている。壁で 0 で自由端 B で 12 kN/m まで直線的に増加する分布荷重と、壁から 5 m の位置の集中荷重 P が作用する。 壁における反力モーメント MA の大きさを求めよ。

L = 8 mw0 = 12 kN/mP = 65 kN 位置 5 m
MA = kN·m
3.軸応力とひずみふつう
P = ?12600 mm400 mmØ30Ø60δ = 0.3 mm

棒が、鋼 の区間(直径 30 mm、長さ 600 mm、E = 200 GPa)と、青銅 でできた区間(直径 60 mm、長さ 400 mm、E = 110 GPa)を端で接合したものでできています。一端が固定され、もう一端が引かれます。全体で 0.3 mm 伸ばす軸荷重 P を求めなさい。

区間 1: 鋼, Ø30 mm, 600 mm, E = 200 GPa区間 2: 青銅, Ø60 mm, 400 mm, E = 110 GPaδ = 0.3 mm
P = kN
4.軸応力とひずみむずかしい
Ø 25 mm銅1 mØ 20 mm黄銅0.5 mP = 25 kN

段付き棒は壁に固定され、自由端に荷重を受ける。棒の最大垂直応力 σmax を求めよ。

P = 25 kN区間 1:Ø25 mm、1 m、銅(E = 117 GPa)区間 2:Ø20 mm、0.5 m、黄銅(E = 100 GPa)
σmax = MPa
5.軸のねじりふつう
d1​ = 40 mmd2​ = 30 mmT = 400 N·mABC2.5 m1.2 m

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

T = 400 N·md1 = 40 mm, L1 = 2.5 md2 = 30 mm, L2 = 1.2 m
τmax = MPa
6.軸のねじりむずかしい
Tdo​ = 50 mm, di​ = 25 mm

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

P = 75 kWN = 750 rpmdo = 50 mmdi = 25 mm
τmax = MPa
7.せん断力と曲げモーメントふつう
P1​ = 70 kNP2​ = 60 kNAB2 m1.5 m1.5 m

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

L = 5 mP1 = 70 kN 位置 2 mP2 = 60 kN 位置 3.5 m
Mmax = kN·m
8.せん断力と曲げモーメントむずかしい
w = 5 kN/mP = 40 kNAB2.5 m3.5 m

単純ばりの全スパンに等分布荷重 w が作用し、さらに集中荷重 P が作用する。はりの最大曲げモーメント Mmax を求めよ。

L = 6 mP = 40 kN 位置 2.5 mw = 5 kN/m(0 から 6 m まで)
Mmax = kN·m
9.モールの応力円と応力の座標変換ふつう
σx​ = 70σy​ = 40τxy​ = 20MPa

ある点で σₓ = 70 MPa、σᵧ = 40 MPa、τₓᵧ = 20 MPa である(引張を正とする)。せん断応力が最大の正の値 +R をとる面について、x 軸からその面の法線までの角度 θs を −90° から 90° の範囲で求めよ(せん断応力は τₓᵧ と同じ向きを正とする)。

σx = 70 MPaσy = 40 MPaτxy = 20 MPa
θs = °
10.モールの応力円と応力の座標変換むずかしい
σx​ = 50σy​ = 90τxy​ = 30MPa

薄板内の点は平面応力状態にある:σₓ = 50 MPa、σᵧ = 90 MPa、τₓᵧ = 30 MPa(σz = 0)。 材料の降伏強さは 275 MPa である。トレスカ(最大せん断応力)説による降伏に対する安全率 n を求めよ。

σx = 50 MPaσy = 90 MPaτxy = 30 MPaσY = 275 MPa
n =
11.はりのたわみふつう
w = 10 kN/mδL = 2.5 m

片持ちばりの全長に等分布荷重 w が作用する。自由端のたわみ δ を求めよ。

w = 10 kN/mL = 2.5 mE = 11 GPa (構造用木材)I = 120 × 10⁶ mm⁴
δ = mm
12.はりのたわみむずかしい
w = 8 kN/mP = 18 kNδL = 3 m

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

P = 18 kNw = 8 kN/mL = 3 mE = 70 GPa (アルミニウム)I = 120 × 10⁶ mm⁴
δmax = mm
13.柱の座屈ふつう
PL = 3 m25 × 50 mm断面

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

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

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

P = 70 kNn = 3L = 3 mE = 200 GPa (鋼)
d = mm
15.放物運動ふつう
v0​ = ?35°R = 135 m

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

θ = 35°R = 135 mg = 9.81 m/s²
v0 = m/s
16.放物運動むずかしい
h = 40 mv0​ = 32 m/s45°到達距離

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

v0 = 32 m/sθ = 45°h = 40 mg = 9.81 m/s²
t = s
17.仕事とエネルギーふつう
30°d = 4 mμ = 0.1

ブロックが水平と 30° をなす斜面を静止状態からすべり下りる。動摩擦係数は 0.1 である。斜面に沿って 4 m すべった後の速さ v を求めよ。

θ = 30°μ = 0.1d = 4 mg = 9.81 m/s²
v = m/s
18.仕事とエネルギーむずかしい
h = 6 mR = 1.5 mN = ?

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

m = 2 kgh = 6 mR = 1.5 mg = 9.81 m/s²
N = N
19.機械振動ふつう
k = 40 kN/mζ = 0.4m = 15 kg

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

m = 15 kgk = 40 kN/mζ = 0.4
fd = Hz
20.機械振動むずかしい
k1​ = 60k2​ = 60m = 5 kg

5 kg の質量が、横に並べた(並列の)2 本のばね k₁ = 60 kN/m、k₂ = 60 kN/m につり下げられている。質量は水平を保つ。 減衰とばねの質量は無視する。上下振動の固有振動数 fn を Hz で求めよ。

m = 5 kgk1 = 60 kN/mk2 = 60 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 350 °CTL​ = 30 °CEQH​QL​η = 26%

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

TH = 350 °CTL = 30 °Cη = 26%
ηII = %
22.熱機関と冷凍機むずかしい
TH​ = 25 °CTL​ = -10 °CRQ̇H​Q̇L​ = 2.5 kWẆ = ?

冷蔵庫が、25 °C のキッチンで冷蔵室を -10 °C に保ち、冷蔵室から 2.5 kW の熱を取り除く。 必要な最小入力動力 Ẇmin を求めよ。

TH = 25 °CTL = -10 °CQ̇L = 2.5 kW
Ẇmin = kW
23.ガスサイクルふつう
vP1234qin​qout​

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

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

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

rp = 12k = 1.4
η = %
25.連続の式とベルヌーイふつう
Q = ?d = 50 mm, Cd​ = 0.64h = 1.5 m

大きな開放水槽から、直径 50 mm の刃形オリフィスを通って水が流出する。水位はオリフィスより 1.5 m 上に保たれ、流量係数は 0.64 である。流量 Q を求めよ。g = 9.81 m/s² とする。

d = 50 mmh = 1.5 mCd = 0.64g = 9.81 m/s²
Q = L/s
26.連続の式とベルヌーイむずかしい
12d1​ = 100 mm, p1​ = 240 kPad2​ = 75 mm, p2​ = ?Δz = 1 m

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

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

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

V = 2.4 m/sL = 1000 mD = 50 mmf = 0.025g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 10自由放出口H = ?L = 100 m, D = 100 mm, f = 0.0222× 曲がり (K = 0.3)

長さ 100 m、内径 100 mm(ダルシーの管摩擦係数 f = 0.022)の管が開放貯水槽から出て、自由噴流として大気中に放出している。局所損失は 角ばった入口(K = 0.5)、2 個の曲がり(各 K = 0.3)、弁(K = 10) である。g = 9.81 m/s² とする。流量 15 L/s を得るには、貯水槽の水面を管の出口より何 m H 高くしなければならないか。

Q = 15 L/sL = 100 mD = 100 mmf = 0.022Kentrance = 0.5Kbend = 0.3 (×2)Kvalve = 10
H = m
29.ポンプとポンプ動力ふつう
PQ = 60 L/sH = 30 m

試験で、ポンプが揚程 30 m に対して水(ρ = 1000 kg/m³)を 60 L/s 送り、そのとき軸動力は 27.2 kW である。g = 9.81 m/s² とする。ポンプ効率 η をパーセントで求めよ。

Q = 60 L/sH = 30 mPshaft = 27.2 kWg = 9.81 m/s²
η = %
30.ポンプとポンプ動力むずかしい
100 kPa (ゲージ圧)PQ = 40 L/sΔz = 7 mL = 500 m, D = 150 mm, f = 0.022, ΣK = 8

ポンプが 40 L/s の水(ρ = 1000 kg/m³)を、開放された下の貯水槽から、水面が 7 m 高い上のタンクへ送る。管は長さ 500 m、内径 150 mm、ダルシーの管摩擦係数 f = 0.022 で、局所損失係数の合計は ΣK = 8(入口と出口を含む)である。g = 9.81 m/s² とする。 上のタンクは密閉されており、水面上の空気のゲージ圧は 100 kPa である。ポンプが与えなければならない揚程 H を求めよ。

Q = 40 L/sΔz = 7 mL = 500 mD = 150 mmf = 0.022ΣK = 8ptank = 100 kPa
H = m
31.熱伝導ふつう
れんが, k = 0.7250 mm発泡ポリスチレン, k = 0.03525 mm19 °C3 °CTi​ = ?Q̇

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

L1 = 50 mm, k1 = 0.72 W/(m·K)L2 = 25 mm, k2 = 0.035 W/(m·K)A = 6 m²T1 = 19 °CT3 = 3 °C
Ti = °C
32.熱伝導むずかしい
発泡ポリスチレン, k = 0.035100 mmロックウール, k = 0.0475 mm室内 22 °C屋外 -11 °CQ̇

100 mm の発泡ポリスチレンと 75 mm のロックウールからなる壁が、22 °C の部屋と -11 °C の外気を隔てている。熱伝達率は室内側 10 W/(m²·K)、屋外側 25 W/(m²·K)、壁の面積は 8 m² である。 壁の熱通過率 U を求めよ。

L1 = 100 mm, k1 = 0.035 W/(m·K)L2 = 75 mm, k2 = 0.04 W/(m·K)A = 8 m²hi = 10 W/(m²·K)ho = 25 W/(m²·K)Troom = 22 °CTout = -11 °C
U = W/(m²·K)
33.歯車と歯車列ふつう
N1​ = 18N2​ = 32

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

P = 3 kWn1 = 2900 rpmN1 = 18N2 = 32η = 95%
T2 = N·m
34.歯車と歯車列むずかしい
d1​ = 100 mmd2​ = 450 mm1500 rpmC = 500 mmT1​ = 2000 NT2​P = ?

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

d1 = 100 mmn1 = 1500 rpmd2 = 450 mmC = 500 mmμ = 0.25T1 = 2000 N
P = kW
35.ラプラス変換ふつう
t (s)f(t)t = 0.3f(0.3) = ?

F(s) = 10/(s(s + 3)(s + 6)) の f(t) を部分分数展開で求め、f(0.3) を計算せよ。

F(s) = 10/(s(s + 3)(s + 6))t = 0.3 s
f(0.3) =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.6y(0.6) = ?

ラプラス変換を用いて、y(0) = 8、y′(0) = 3 の条件で y″ + 4y′ + 3y = 0 を解き、y(0.6) を求めよ。

y″ + 4y′ + 3y = 0y(0) = 8y′(0) = 3t = 0.6 s
y(0.6) =
37.ブロック線図ふつう
R(s)+K(s + 9)(s + 10)G(s)C(s)−

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

G(s) = K/((s + 9)(s + 10))単位負帰還
K =
38.ブロック線図むずかしい
R(s)+32s(s + 9)G(s)C(s)2H−

プラント G(s) = 32/(s(s + 9)) を、センサのゲイン H = 2 の負帰還ループに入れる。閉ループの極は複素数である。 その虚部 ωd(正の値)を rad/s で求めよ。

G(s) = 32/(s(s + 9))H = 2
ωd = rad/s
39.ステップ応答ふつう
t (s)y(t)2ts​ = ?

ある系の伝達関数は G(s) = 132/(s² + 4s + 66) で、単位ステップ入力を受ける。 ts ≈ 4/(ζωₙ) を用いて、ステップ応答の 2 % 整定時間 ts を求めよ。

G(s) = 132/(s² + 4s + 66)
ts = s
40.ステップ応答むずかしい
R(s)+Ks(s + 2)G(s)C(s)−

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

G(s) = K/(s(s + 2))%OS = 25 %
K =
41.定常偏差と PID 制御ふつう
r(t) = 0.5·t²/2+20(s + 2)s²(s + 2)C(s)−

放物線入力 r(t) = 0.5·t²/2 が、G(s) = 20(s + 2)/(s²(s + 2)) の単位フィードバックのループに加えられます。定常偏差 ess を求めなさい。

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

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

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

3溶液

1.はりの支点反力

答え: RA = 15.63 kN

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

P1​ = 15 kNP2​ = 40 kNABRA​RB​2.5 m0.5 m1 m

自由体図

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

つり合い

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

ΣMA = 0: RB·4 = 15·2.5 + 40·3
RB = 157.5/4 = 39.38 kN
ΣFy = 0: RA = 55 − 39.38 = 15.63 kN

答え

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

答え: MA = 581 kN·m

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

三角形分布荷重の合力

W = ½·12·8 = 48 kN
x̄ = 2L/3 = 5.333 m

壁まわりのモーメント

MA = W·x̄ + P·a = 48·5.333 + 65·5 = 581 kN·m

答え

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

答え: P = 54.25 kN

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

面積

A₁ = 706.9 mm², A₂ = 2827 mm²

柔軟性

L₁/(E₁A₁) = 600/(200000·706.9) = 4.244 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 400/(110000·2827) = 1.286 × 10⁻⁶ mm/N

荷重

P = 0.3/5.53 × 10⁻⁶ = 54247 N = 54.25 kN

答え

P = 54.25 kN
4.軸応力とひずみ

答え: σmax = 79.58 MPa

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

両区間で同じ力

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

A1 = π·25²/4 = 490.9 mm², A2 = π·20²/4 = 314.2 mm²
σ1 = 50.93 MPa
σ2 = 79.58 MPa

最大応力

細い区間は断面積が小さいので、応力が大きくなります。

σmax = 79.58 MPa

答え

σmax = 79.58 MPa
5.軸のねじり

答え: τmax = 75.45 MPa

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

各区間の応力

T = 400 N·m = 400000 N·mm
τAB = 16·400000/(π·40³) = 31.83 MPa
τBC = 16·400000/(π·30³) = 75.45 MPa

答え

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

答え: τmax = 41.5 MPa

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

トルク

ω = 2πN/60 = 2π·750/60 = 78.54 rad/s
T = P/ω = 75000/78.54 = 954.9 N·m

断面二次極モーメント

J = π(50⁴ − 25⁴)/32 = 575243 mm⁴

最大せん断応力

τmax = T·(do/2)/J = 954930·25/575243 = 41.5 MPa

答え

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

答え: Mmax = 120 kN·m

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

支点反力

はり全体のつり合い。

ΣMA = 0: RB·5 = 70·2 + 60·3.5
RB = 350/5 = 70 kN
ΣFy = 0: RA = 130 − 70 = 60 kN

各荷重位置のモーメント

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

荷重点 P1: M = RA·2 = 60·2 = 120 kN·m
荷重点 P2: M = RB·1.5 = 70·1.5 = 105 kN·m

最大値

Mmax = 120 kN·m (荷重点 P1)

答え

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

答え: Mmax = 80.21 kN·m

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

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 5·6 = 30 kN、作用点 x̄ = 3 m
ΣMA = 0: RB·6 = 40·2.5 + 5·6·3
RB = 190/6 = 31.67 kN
ΣFy = 0: RA = 70 − 31.67 = 38.33 kN

V = 0 となるのはどこですか?

V(2.5⁻) = RA − w·2.5 = 25.83 kN
V(2.5⁺) = 25.83 − 40 = −14.17 kN
荷重の位置でせん断力の符号が変わるので、最大値は荷重の位置にあります。

最大モーメント

Mmax = RA·2.5 − w·2.5²/2 = 38.33·2.5 − 5·2.5²/2 = 80.21 kN·m

答え

Mmax = 80.21 kN·m

それは A から x = 2.5 m の位置に生じます。

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

答え: θs = −18.43 °

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

στxyCτmax​

円上の角度

モールの応力円上では頂点(+R)は主応力点 σ₁ から 90° 離れており、要素上では 45° です。

cos 2θs = 20/25 = 0.8, sin 2θs = −15/25 = −0.6
2θs = atan2(−15, 20) = −36.87°

要素上の角度

θs = −18.43°
確認: θs − θp = −18.43 − 26.57 = −45°

答え

θs = −18.43°

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

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

答え: n = 2.593

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

στxyCσ1​σ2​

主応力

C = 70 MPa, R = 36.06 MPa
σ1 = 106.1 MPa, σ2 = 33.94 MPa, σ3 = 0

トレスカ

σmax − σmin = 106.1 MPa
n = 275/106.1 = 2.593

答え

n = 2.593

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

11.はりのたわみ

答え: δ = 36.99 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⁴。

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

たわみ

δ = wL⁴/(8EI) = 10·2500⁴/(8·1.32 × 10¹²) = 36.99 mm

答え

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

答え: δmax = 2.21 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 = 18000 N, w = 8 N/mm, L = 3000 mm, EI = 8.4 × 10¹² N·mm²

各荷重を単独で

δP = PL³/(48EI) = 18000·3000³/(48·8.4 × 10¹²) = 1.205 mm
δw = 5wL⁴/(384EI) = 5·8·3000⁴/(384·8.4 × 10¹²) = 1.004 mm

足し合わせる

δmax = 1.205 + 1.004 = 2.21 mm

答え

δmax = 2.21 mm
13.柱の座屈

答え: Pcr = 57.12 kN

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

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

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

Imin = h·b³/12 = 50·25³/12 = 65104 mm⁴

有効座屈長さ

KL = 0.5·3000 = 1500 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·200000·65104/1500² = 57116 N = 57.12 kN

答え

Pcr = 57.12 kN
14.柱の座屈

答え: d = 46.99 mm

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

必要な座屈荷重

Pcr = n·P = 3·70000 = 210000 N
KL = 0.5·3000 = 1500 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·210000·1500²/(π³·200000))^(1/4) = 46.99 mm

答え

d = 46.99 mm

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

15.放物運動

答え: v0 = 37.54 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·135/sin 70°) = 37.54 m/s

答え

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

答え: t = 5.977 s

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

初速度を分解する

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

vx = v0·cos θ = 32·cos 45° = 22.63 m/s
vy = v0·sin θ = 32·sin 45° = 22.63 m/s

滞空時間

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

0 = 40 + 22.63·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (22.63 + √(512 + 784.8))/9.81 = 5.977 s

答え

t = 5.977 s
17.仕事とエネルギー

答え: v = 5.696 m/s

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

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

ブロックは高さ d·sin θ だけ下がり、摩擦が μmg·cos θ·d を奪います。質量は消えます。

½v² = g·d·(sin θ − μ·cos θ)
v = √(2·9.81·4·(0.5 − 0.1·0.866)) = 5.696 m/s

答え

v = 5.696 m/s
18.仕事とエネルギー

答え: N = 58.86 N

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

最高点での速さ

v² = 2g(h − 2R) = 2·9.81·(6 − 3) = 58.86 m²/s²

最高点での力

m·g = 2·9.81 = 19.62 N
N = mv²/R − mg = 2·58.86/1.5 − 19.62 = 58.86 N

答え

N = 58.86 N

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

19.機械振動

答え: fd = 7.533 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 = √(40000/15) = 51.64 rad/s
fn = 51.64/(2π) = 8.219 Hz

減衰固有振動数

fd = fn·√(1 − ζ²) = 8.219·√(1 − 0.4²) = 7.533 Hz

答え

fd = 7.533 Hz
20.機械振動

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

等価ばね定数

並べたばねはたわみが共通で、力が足し合わされるので、ばね定数はそのまま足します:k = k1 + k2。

keq = 60 + 60 = 120 kN/m

固有振動数

ωn = √(120000/5) = 154.9 rad/s
fn = 154.9/(2π) = 24.66 Hz

答え

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

答え: ηII = 50.63 %

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 = 350 + 273.15 = 623.1 K
TL = 30 + 273.15 = 303.1 K

カルノー効率

ηCarnot = 1 − 303.1/623.1 = 51.35 %

第二法則効率

ηII = η/ηCarnot = 26/51.35 = 50.63 %

答え

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

答え: Ẇmin = 0.3325 kW

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

絶対温度

TH = 25 + 273.15 = 298.1 K
TL = -10 + 273.15 = 263.1 K

可能な最良の COP

COPR,max = TL/(TH − TL) = 263.1/(298.1 − 263.1) = 7.519

最小動力

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

Ẇmin = Q̇L/COP = 2.5/7.519 = 0.3325 kW

答え

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

答え: η = 57.96 %

ディーゼルサイクル(圧縮着火機関):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.322·1.306 = 57.96%

答え

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

答え: η = 50.83 %

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

答え

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

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

理想的な噴流速度

√(2gh) = √(2·9.81·1.5) = 5.425 m/s

実際の流量

流量係数は、噴流の縮流と縁での摩擦を考慮したものです。

Q = Cd·A·√(2gh) = 0.64·0.001963·5.425 = 0.006817 m³/s = 6.817 L/s

答え

Q = 6.817 L/s
26.連続の式とベルヌーイ

答え: p2 = 222 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 2.75·(100/75)² = 4.889 m/s
½ρ(v1² − v2²) = ½·1000·(2.75² − 4.889²)/1000 = −8.169 kPa
ρgΔz = 1000·9.81·1/1000 = 9.81 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 240 + (−8.169) − 9.81 = 222 kPa

答え

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

答え: hf = 146.8 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²/(2g) = 2.4²/(2·9.81) = 0.2936 m

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

hf = f·(L/D)·V²/(2g) = 0.025·(1000/0.05)·2.4²/(2·9.81) = 146.8 m

答え

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

答え: H = 6.34 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.015/0.007854 = 1.91 m/s
V²/(2g) = 0.1859 m

損失係数

f·L/D = 0.022·100/0.1 = 22
ΣK = 0.5 + 2·0.3 + 10 = 11.1

エネルギー式

出口の噴流はまだ速度水頭をもっています。これは損失ではありませんが、まかなう必要があります。

H = (1 + 22 + 11.1)·0.1859 = 6.34 m

答え

H = 6.34 m
29.ポンプとポンプ動力

答え: η = 64.92 %

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.06·30 = 17.66 kW

効率

η = Ph/Pshaft = 17.66/27.2 = 64.92%

答え

η = 64.92%
30.ポンプとポンプ動力

答え: H = 38.43 m

ポンプ動力:水に与えられる動力は水動力 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) = (73.33 + 8)·0.2611 = 21.24 m

圧力水頭

p/(ρ·g) = 100000/(1000·9.81) = 10.19 m

ポンプ揚程

どちらの水面もほぼ静止しているので、速度水頭の項はありません。

H = Δz + p/(ρ·g) + hL = 7 + 10.19 + 21.24 = 38.43 m

答え

H = 38.43 m
31.熱伝導

答え: Ti = 17.58 °C

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

各層の熱抵抗

R1 = L1/(k1A) = (50/1000)/(0.72·6) = 0.01157 K/W
R2 = L2/(k2A) = (25/1000)/(0.035·6) = 0.119 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 16/0.1306 = 122.5 W

境界面温度

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

Ti = T1 − Q̇·R1 = 19 − 122.5·0.01157 = 17.58 °C

答え

Ti = 17.58 °C
32.熱伝導

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

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(10·8) = 0.0125 K/W
R1 = L1/(k1A) = (100/1000)/(0.035·8) = 0.3571 K/W
R2 = L2/(k2A) = (75/1000)/(0.04·8) = 0.2344 K/W
Ro = 1/(hoA) = 1/(25·8) = 0.005 K/W
ΣR = 0.609 K/W

熱通過率

U = 1/(A·ΣR) = 1/(8·0.609) = 0.2052 W/(m²·K)

答え

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

答え: T2 = 16.68 N·m

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

モータのトルク

ω₁ = 2π·2900/60 = 303.7 rad/s
T1 = P/ω₁ = 3000/303.7 = 9.879 N·m

出力トルク

T2 = 9.879·32/18·0.95 = 16.68 N·m

答え

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

答え: P = 7.144 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.426 rad

緩み側の張力

eμθ = e0.6066 = 1.834, T₂ = 2000/1.834 = 1090 N

電力

P = (2000 − 1090)·7.854 = 7.144 kW

答え

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

答え: f(0.3) = 0.1956

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

部分分数展開

F(s) = A/s + B/(s + 3) + C/(s + 6)
A = 10/(3·6) = 0.5556
B = 10/((−3)·3) = −1.111
C = 10/((−6)·(−3)) = 0.5556

逆変換

f(t) = 0.5556 − 1.111·e−3t + 0.5556·e−6t

計算する

f(0.3) = 0.5556 − 1.111·e−3·0.3 + 0.5556·e−6·0.3 = 0.1956

答え

f(0.3) = 0.1956
36.ラプラス変換

答え: y(0.6) = 6.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²Y − 8s − 3 + 4(sY − 8) + 3Y = 0
Y(s) = (8s + 35)/((s + 1)(s + 3))

部分分数展開

A = 27/2 = 13.5
B = 11/(−2) = −5.5

結果

y(t) = 13.5·e−t − 5.5·e−3t
y(0.6) = 6.5

答え

y(0.6) = 6.5

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

37.ブロック線図

答え: K = 0.25

ブロック線図の代数:直列のブロックは掛け算(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 + 9)(s + 10)) = 0 → s² + 19s + 90 + K = 0

重根

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

19² = 4(90 + K) → K = 361/4 − 90 = 0.25

答え

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

答え: ωd = 6.614 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(s + 9) + 32·2 = 0 → s² + 9s + 64 = 0

根

s = −4.5 ± j·√(64 − 20.25) = −4.5 ± j6.614
ωd = 6.614 rad/s

答え

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

答え: ts = 2 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) です。

標準形と比べる

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

ωₙ² = 66 → ωₙ = 8.124 rad/s
2ζωₙ = 4 → ζ = 4/(2·8.124) = 0.2462

結果

ζωₙ = 2 → ts ≈ 4/2 = 2 s

答え

ts = 2 s
40.ステップ応答

答え: K = 6.136

ステップ応答: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² + 2s + K)

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

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

利得

ωₙ = 2/(2ζ) = 2.477 rad/s
K = ωₙ² = 6.136

答え

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

答え: ess = 0.025

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

加速度偏差定数

Ka = 20·2/2 = 20

放物線偏差

ess = 0.5/20 = 0.025

答え

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

答え: ess = 2

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

Kv = 0.5·20/10 = 1

ランプ偏差

ess = 2/1 = 2

答え

ess = 2

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

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

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