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.はりの支点反力ふつう
P = 75 kNM = 15 kN·mABRA​ = ?2.5 m2.5 m1 m

単純ばりは A がピン支点、B がローラー支点である。A から 2.5 m の位置に集中荷重 P が、A から 5 m の位置に時計回りの偶力 M が作用する。上向きを正として鉛直反力 RA を求めよ。

L = 6 mP = 75 kN 位置 2.5 mM = 15 kN·m (時計回り) 位置 5 m
RA = kN
2.はりの支点反力むずかしい
w = 12 kN/mP = 50 kNABCRB​ = ?2.5 m0.5 m1 m2 m

はりは A でピン支点、A から 4 m の B でローラー支点に支持され、B から自由端 C まで張り出している。A から 2.5 m の位置に集中荷重 P が作用し、3 m から 6 m の自由端まで等分布荷重 w が作用する。上向きを正として鉛直反力 RB を求めよ。

AB = 4 mAC = 6 mP = 50 kN 位置 2.5 mw = 12 kN/m(3 から 6 m まで)
RB = kN
3.軸応力とひずみふつう
Ø 25 mmアルミニウム3 mP = 115 kN

アルミニウム製の棒の軸ひずみ ε をマイクロひずみ (µε = 10⁻⁶) で求めよ。

P = 115 kNd = 25 mmL = 3 mE = 70 GPa
ε = µε
4.軸応力とひずみむずかしい
100 kN12400 mm200 mmØ40Ø50σ2​ = ?

鋼 の区間(直径 40 mm、長さ 400 mm、E = 200 GPa)と、銅 でできた区間(直径 50 mm、長さ 200 mm、E = 120 GPa)が端で接合され、2 つの剛な壁の間に固定されています。接合部に、区間 2 に向かって軸荷重 100 kN が作用します。区間 2 の応力の大きさ σ2 を求めなさい。

区間 1: 鋼, Ø40 mm, 400 mm, E = 200 GPa区間 2: 銅, Ø50 mm, 200 mm, E = 120 GPaP = 100 kN
σ2 = MPa
5.軸のねじりふつう
d = 12 mmT = 30 N·m200 mmφ = 4.33°

ねじり試験で、直径 12 mm の中実試験片をトルク 30 N·m でねじる。標点距離 200 mm の両端は互いに 4.33° 回転する(材料はまだ弾性域にある)。材料の横弾性係数 G を求めよ。

T = 30 N·md = 12 mmL = 200 mmφ = 4.33°
G = GPa
6.軸のねじりむずかしい
Td = ?

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

P = 100 kWN = 3000 rpmτallow = 60 MPa
d = mm
7.せん断力と曲げモーメントふつう
10 kN35 kNABCDEFG2 m2 m2 m2 m2 m2 m

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

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

単純ばりのスパンの一部に等分布荷重が作用する。 はりの最大曲げモーメント Mmax を求めよ。

L = 8 mw = 20 kN/m(5 から 7 m まで)
Mmax = kN·m
9.モールの応力円と応力の座標変換ふつう
σx​ = 30σy​ = 20τxy​ = 50MPa

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

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

荷重を受ける部品内のある点で、σₓ = 70 MPa、σᵧ = 90 MPa、τₓᵧ = 30 MPa である(引張を正とする)。 法線が x 軸から反時計回りに 45° 回転した面上の垂直応力 σθ を求めよ。

σx = 70 MPaσy = 90 MPaτxy = 30 MPaθ = 45°
σθ = MPa
11.はりのたわみふつう
w = 10 kN/mδL = 6 m

単純ばりの全スパンに等分布荷重 w が作用する。スパン中央の最大たわみ δmax を求めよ。

w = 10 kN/mL = 6 mE = 200 GPa (鋼)I = 120 × 10⁶ mm⁴
δmax = mm
12.はりのたわみむずかしい
w = 8 kN/mδL = 2 m

単純ばりに等分布荷重 w が作用する。たわみはスパン/360 を超えてはならない。この制限を満たす最小の断面二次モーメント Ireq を求めよ。

w = 8 kN/mL = 2 mE = 200 GPa (鋼)
Ireq = × 10⁶ mm⁴
13.柱の座屈ふつう
PL = 5 m40 × 60 mm断面

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

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

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

P = 5 kNn = 3L = 3 mE = 70 GPa (アルミニウム)
d = mm
15.放物運動ふつう
53 my = ?v0​ = 34 m/s20°

水平な地面からボールを 34 m/s、水平から 20° 上向きに、53 m 先の壁に向かって蹴る。空気抵抗を無視し、g = 9.81 m/s² とする。ボールはどの高さ y で壁を通過するか。

v₀ = 34 m/sθ = 20°x = 53 mg = 9.81 m/s²
y = m
16.放物運動むずかしい
h = 55 mv0​ = 27 m/s20°x = ?

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

v0 = 27 m/sθ = 20°h = 55 mg = 9.81 m/s²
x = m
17.仕事とエネルギーふつう
30°L = 0.5 m1 kg

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

m = 1 kgL = 0.5 mθ = 30°g = 9.81 m/s²
T = N
18.仕事とエネルギーむずかしい
k = 800 N/m15°μ = 0.22 kg

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

m = 2 kgk = 800 N/mx = 12 cmθ = 15°μ = 0.2g = 9.81 m/s²
d = m
19.機械振動ふつう
tx40 mm1 周期15.4 mm

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

x0 = 40 mmx1 = 15.4 mmn = 1
ζ =
20.機械振動むずかしい
k1​ = 80k2​ = 40m = 40 kg

40 kg の質量が、縦につないだ(直列の)2 本のばねにつり下げられている:上のばねは k₁ = 80 kN/m、下のばねは k₂ = 40 kN/m である。 ばね系の等価ばね定数 keq を kN/m で求めよ。

m = 40 kgk1 = 80 kN/mk2 = 40 kN/m
keq = kN/m
21.熱機関と冷凍機ふつう
TH​ = 30 °CTL​ = 0 °CRQ̇H​ = ?Q̇L​ = 1.5 kWẆ

可逆冷凍機が 1.5 kW の熱を取り除いて冷蔵室を 0 °C に保ち、30 °C の周囲へ熱を放出する。放出する熱量 Q̇H を求めよ。

TL = 0 °CTH = 30 °CQ̇L = 1.5 kW
Q̇H = kW
22.熱機関と冷凍機むずかしい
TH​ = 32 °CTL​ = -25 °CRQ̇H​Q̇L​ = 2 kWẆ = ?

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

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

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

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

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

rp = 12T1 = 290 KT3 = 1400 Kηc = 82%ηt = 85%k = 1.4cp = 1.005 kJ/(kg·K)
wnet = kJ/kg
25.連続の式とベルヌーイふつう
Q = ?d = 25 mm, Cd​ = 0.64h = 3.5 m

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

d = 25 mmh = 3.5 mCd = 0.64g = 9.81 m/s²
Q = L/s
26.連続の式とベルヌーイむずかしい
pair​ = 50 kPav = ?d = 25 mmh = 4 m

密閉水槽の水面上に、ゲージ圧 50 kPa の圧縮空気がある。水面より 4 m 下にある直径 25 mm の孔から、水が大気中に流出する。損失を無視して、次を求めよ 噴流速度 v。 g = 9.81 m/s² とする。

pair = 50 kPa (ゲージ圧)h = 4 md = 25 mmρ = 1000 kg/m³g = 9.81 m/s²
v = m/s
27.管摩擦と損失水頭ふつう
Q = 40 L/sD = 150 mmL = 800 mf = 0.03

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

Q = 40 L/sL = 800 mD = 150 mmf = 0.03g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 10Δz = 28 mL = 300 m, D = 250 mm, f = 0.022× 曲がり (K = 0.3)

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

Δz = 28 mL = 300 mD = 250 mmf = 0.02Kentrance = 0.5Kbend = 0.3 (×2)Kvalve = 10Kexit = 1.0
Q = L/s
29.ポンプとポンプ動力ふつう
PQ = 10 L/sH = 72 m

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

Q = 10 L/sH = 72 mPshaft = 9.8 kWg = 9.81 m/s²
η = %
30.ポンプとポンプ動力むずかしい
PQ = 20 L/sΔz = 36 mL = 150 m, D = 100 mm, f = 0.018, ΣK = 8

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

Q = 20 L/sΔz = 36 mL = 150 mD = 100 mmf = 0.018ΣK = 8η = 75%
Pshaft = kW
31.熱伝導ふつう
ガラス, k = 0.96150 mmれんが, k = 0.7225 mm24 °C-13 °CQ̇ = ?

壁はガラス 150 mm とれんが 25 mm からなり、面積は 10 m² である。外側の両面は 24 °C と -13 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 150 mm, k1 = 0.96 W/(m·K)L2 = 25 mm, k2 = 0.72 W/(m·K)A = 10 m²T1 = 24 °CT3 = -13 °C
Q̇ = W
32.熱伝導むずかしい
ガラス, k = 0.96250 mm発泡ポリスチレン, k = 0.03575 mm室内 22 °C屋外 2 °CTs​ = ?Q̇

250 mm のガラスと 75 mm の発泡ポリスチレンからなる壁が、22 °C の部屋と 2 °C の外気を隔てている。熱伝達率は室内側 10 W/(m²·K)、屋外側 30 W/(m²·K)、壁の面積は 12 m² である。 壁の内表面温度 Ts を求めよ。

L1 = 250 mm, k1 = 0.96 W/(m·K)L2 = 75 mm, k2 = 0.035 W/(m·K)A = 12 m²hi = 10 W/(m²·K)ho = 30 W/(m²·K)Troom = 22 °CTout = 2 °C
Ts = °C
33.歯車と歯車列ふつう
Tw​ = ?z = 3, 2900 rpm

2900 rpm の 1.1 kW のモータが条数 3 のウォームを駆動する。ウォームホイールの歯数は 45、伝動効率は 55% である。ホイール軸の出力トルク Tw を求めよ。

z = 3Nw = 45n = 2900 rpmP = 1.1 kWη = 55%
Tw = N·m
34.歯車と歯車列むずかしい
Ns​ = 15リング固定Np​ = 12nc​ = ?

遊星歯車装置で、太陽歯車の歯数は 15、各遊星歯車の歯数は 12 である。リングギヤを固定し、太陽歯車を 900 rpm で駆動する。遊星キャリヤの回転速度 nc を求めよ。

Ns = 15Np = 12ns = 900 rpm
nc = rpm
35.ラプラス変換ふつう
t (s)f(t)t = 0.3f(0.3) = ?

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

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

ある系が y″ + 5y′ + 4y = 6·u(t) に従う。系は静止状態(y(0) = 0、y′(0) = 0)から始まり、入力は t ≥ 0 で大きさ 10 のステップ u(t) である。ラプラス変換を用いて出力 y(0.6) を求めよ。

y″ + 5y′ + 4y = 6·u(t)u(t) = 10y(0) = 0, y′(0) = 0t = 0.6 s
y(0.6) =
37.ブロック線図ふつう
R(s)+2(s + 3)(s + 4)(s + 7)G(s)C(s)−

単位負帰還ループの G(s) = 2(s + 3)/((s + 4)(s + 7)) である。基準入力に大きさ 5 のステップを加える。定常偏差 ess を求めよ。

G(s) = 2(s + 3)/((s + 4)(s + 7))R = 5
ess =
38.ブロック線図むずかしい
R(s)+K(s + 2)(s + 4)(s + 8)G(s)C(s)−

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

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

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

G(s) = 249/(s² + 11s + 83)
ts = s
40.ステップ応答むずかしい
+−10 VR = 100 ΩL = 10 mHC = 1 µFvC​

R = 100 Ω、L = 10 mH、C = 1 µF の直列 RLC 回路に、t = 0 で 10 V のステップを加える。コンデンサは充電されておらず、ステップの前には電流は流れていない。 コンデンサ電圧のオーバーシュート率 %OS を求めよ。

R = 100 ΩL = 10 mHC = 1 µFV = 10 V
%OS = %
41.定常偏差と PID 制御ふつう
r(t) = 5t+20(s + 2)s(s + 4)(s + 5)C(s)−

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

G(s) = 20(s + 2)/(s(s + 4)(s + 5))r(t) = 5t
ess =
42.定常偏差と PID 制御むずかしい
R(s)+Kp​ + Ki​/s1s + 2C(s)−

PI 制御器 C(s) = Kp + Ki/s が、単位フィードバックのループで制御対象 G(s) = 1/(s + 2) を制御しています。閉ループは ωₙ = 5 rad/s、ζ = 0.5 でなければなりません。積分ゲイン Ki を求めなさい。

G(s) = 1/(s + 2)ωₙ = 5 rad/sζ = 0.5
Ki =

3溶液

1.はりの支点反力

答え: RA = 41.25 kN

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

A まわりのモーメント

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

ΣMA = 0: RB·6 − 75·2.5 − 15 = 0
RB = 202.5/6 = 33.75 kN

鉛直方向のつり合い

RA = 75 − 33.75 = 41.25 kN

答え

RA = 41.25 kN

負の反力は下向きに作用します。

2.はりの支点反力

答え: RB = 71.75 kN

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

分布荷重の合力

W = 12·3 = 36 kN, x̄ = (3 + 6)/2 = 4.5 m

つり合い

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

ΣMA = 0: RB·4 = 50·2.5 + 36·4.5
RB = 287/4 = 71.75 kN
RA = 86 − 71.75 = 14.25 kN

答え

RB = 71.75 kN
3.軸応力とひずみ

答え: ε = 3347 µε

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

応力

A = π·25²/4 = 490.9 mm²
σ = P/A = 115·1000/490.9 = 234.3 MPa

ひずみ(フックの法則)

E for アルミニウム = 70 GPa = 70000 MPa.

ε = σ/E = 234.3/70000 = 0.003347 = 3347 µε

答え

ε = 3347 µε
4.軸応力とひずみ

答え: σ2 = 33.21 MPa

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

剛性

k₁ = 200000·1257/400 = 628319 N/mm
k₂ = 120000·1963/200 = 1.178 × 10⁶ N/mm

区間 2 の分担

F2 = 100000·1.178 × 10⁶/1.806 × 10⁶ = 65217 N

応力

σ2 = 65217/1963 = 33.21 MPa (圧縮)

答え

σ2 = 33.21 MPa
5.軸のねじり

答え: G = 39 GPa

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

断面二次極モーメントと角度

J = π·12⁴/32 = 2036 mm⁴
φ = 4.33·π/180 = 0.07557 rad

横弾性係数

T = 30 N·m = 30000 N·mm
G = TL/(Jφ) = 30000·200/(2036·0.07557) = 39000 MPa
G = 39 GPa

答え

G = 39 GPa

比較:鋼 ≈ 80 GPa、黄銅 ≈ 39 GPa、アルミニウム ≈ 26 GPa。

6.軸のねじり

答え: d = 30.01 mm

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

トルク

ω = 2πN/60 = 2π·3000/60 = 314.2 rad/s
T = P/ω = 100000/314.2 = 318.3 N·m

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

d = ∛(16T/(π·τallow)) = ∛(16·318310/(π·60)) = 30.01 mm

答え

d = 30.01 mm

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

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

答え: B (x = 2 m)

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

支点反力

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

R1 = 32.5 kN, R2 = 12.5 kN

各点での曲げモーメント

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

MA = 0, MB = 65 kN·m, MC = 60, MD = 55, ME = 50, MF = 25, MG = 0

答え

点 B: MB = 65 kN·m
8.せん断力と曲げモーメント

答え: Mmax = 52.5 kN·m

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

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 20·2 = 40 kN、作用点 x̄ = 6 m
ΣMA = 0: RB·8 = 20·2·6
RB = 240/8 = 30 kN
ΣFy = 0: RA = 40 − 30 = 10 kN

せん断力 0 の位置

せん断力は RA から始まり、荷重のない区間では一定のまま、荷重区間内では勾配 w で減少します。

RA − w·(x0 − 5) = 0
x0 = 5 + 10/20 = 5.5 m

最大モーメント

Mmax = RA·x0 − w·(x0 − 5)²/2 = 10·5.5 − 20·0.5²/2 = 52.5 kN·m

答え

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

答え: θp = 42.14 °

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

στxyCσ1​σ2​

円上の角度

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

2θp = atan2(2·50, 30 − 20) = 84.29°

要素上の角度

θp = 42.14°

答え

θp = 42.14°

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

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

答え: σθ = 110 MPa

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

変換式

σθ = 80 + −10·0 + 30·1 = 110 MPa

答え

σθ = 110 MPa
11.はりのたわみ

答え: δmax = 7.031 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 = 6000 mm
EI = 2.4 × 10¹³ N·mm²

たわみ

δmax = 5wL⁴/(384EI) = 5·10·6000⁴/(384·2.4 × 10¹³) = 7.031 mm

答え

δmax = 7.031 mm
12.はりのたわみ

答え: Ireq = 1.5 × 10⁶ mm⁴

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

たわみの制限

δallow = L/360 = 2000/360 = 5.556 mm

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

I = 5wL⁴/(384·E·δallow) = 5·8·2000⁴/(384·200000·5.556) = 1.5 × 10⁶ mm⁴
Ireq = 1.5 × 10⁶ mm⁴

答え

Ireq = 1.5 × 10⁶ mm⁴
13.柱の座屈

答え: Pcr = 8.843 kN

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

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

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

Imin = h·b³/12 = 60·40³/12 = 320000 mm⁴

有効座屈長さ

KL = 1·5000 = 5000 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·320000/5000² = 8843 N = 8.843 kN

答え

Pcr = 8.843 kN
14.柱の座屈

答え: d = 31.58 mm

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

必要な座屈荷重

Pcr = n·P = 3·5000 = 15000 N
KL = 0.5·3000 = 1500 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·15000·1500²/(π³·70000))^(1/4) = 31.58 mm

答え

d = 31.58 mm

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

15.放物運動

答え: y = 5.793 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 = 34·cos 20° = 31.95 m/s
t = x/vx = 53/31.95 = 1.659 s

その時刻での高さ

y = 11.63·1.659 − ½·9.81·1.659² = 5.793 m

答え

y = 5.793 m

壁の位置では、ボールはすでに下降中です。

16.放物運動

答え: x = 112.1 m

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

初速度を分解する

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

vx = v0·cos θ = 27·cos 20° = 25.37 m/s
vy = v0·sin θ = 27·sin 20° = 9.235 m/s

滞空時間

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

0 = 55 + 9.235·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (9.235 + √(85.28 + 1079))/9.81 = 4.42 s

水平距離

x = vx·t = 25.37·4.42 = 112.1 m

答え

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

答え: T = 12.44 N

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

最下点での速さ

v² = 2gL(1 − cos θ) = 2·9.81·0.5·(1 − cos 30°) = 1.314 m²/s²

最下点での力

m·g = 1·9.81 = 9.81 N
T = mg + mv²/L = 9.81 + 1·1.314/0.5 = 12.44 N

答え

T = 12.44 N

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

18.仕事とエネルギー

答え: d = 0.6495 m

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

ばねのエネルギー

x = 12 cm = 0.12 m
½kx² = ½·800·0.12² = 5.76 J

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

m·g = 2·9.81 = 19.62 N
½kx² = mg·d·(sin θ + μ·cos θ)
d = 5.76/(19.62·(0.2588 + 0.2·0.9659)) = 0.6495 m

答え

d = 0.6495 m
19.機械振動

答え: ζ = 0.1502

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

対数減衰率

δ = (1/1)·ln(40/15.4) = 0.9545

減衰比

ζ = δ/√(4π² + δ²) = 0.9545/√(39.48 + 0.9111) = 0.1502

答え

ζ = 0.1502
20.機械振動

答え: keq = 26.67 kN/m

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

等価ばね定数

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

1/keq = 1/80 + 1/40 → keq = 80·40/(80 + 40) = 26.67 kN/m

答え

keq = 26.67 kN/m
21.熱機関と冷凍機

答え: Q̇H = 1.665 kW

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 = 30 + 273.15 = 303.1 K
TL = 0 + 273.15 = 273.1 K

可逆サイクルの熱量比

Q̇H = Q̇L·TH/TL = 1.5·303.1/273.1 = 1.665 kW

答え

Q̇H = 1.665 kW

入力動力:Ẇ = 0.1647 kW。

22.熱機関と冷凍機

答え: Ẇmin = 0.4594 kW

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

絶対温度

TH = 32 + 273.15 = 305.1 K
TL = -25 + 273.15 = 248.1 K

可能な最良の COP

COPR,max = TL/(TH − TL) = 248.1/(305.1 − 248.1) = 4.354

最小動力

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

Ẇmin = Q̇L/COP = 2/4.354 = 0.4594 kW

答え

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

答え: η = 61.77 %

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

締切比の項

(2.51.4 − 1)/(1.4·(2.5 − 1)) = 1.241

熱効率

η = 1 − 0.308·1.241 = 61.77%

答え

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

答え: wnet = 240.5 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 = 290·120.2857 = 589.8 K
T4s = 1400/120.2857 = 688.3 K

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

wc = 1.005·(589.8 − 290)/0.82 = 367.5 kJ/kg
wt = 0.85·1.005·(1400 − 688.3) = 608 kJ/kg

正味仕事

wnet = 608 − 367.5 = 240.5 kJ/kg

答え

wnet = 240.5 kJ/kg

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

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

答え: Q = 2.603 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·3.5) = 8.287 m/s

実際の流量

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

Q = Cd·A·√(2gh) = 0.64·0.0004909·8.287 = 0.002603 m³/s = 2.603 L/s

答え

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

答え: v = 13.36 m/s

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

水面から噴流へのベルヌーイの式

v = √(2·(50000/1000 + 9.81·4)) = 13.36 m/s

答え

v = 13.36 m/s
27.管摩擦と損失水頭

答え: hf = 41.78 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.04/0.01767 = 2.264 m/s

速度水頭

V²/(2g) = 2.264²/(2·9.81) = 0.2611 m

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

hf = f·(L/D)·V²/(2g) = 0.03·(800/0.15)·2.264²/(2·9.81) = 41.78 m

答え

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

答え: Q = 191.5 L/s

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

損失係数

f·L/D = 0.02·300/0.25 = 24
ΣK = 0.5 + 2·0.3 + 10 + 1.0 = 12.1

エネルギー式

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

28 = (24 + 12.1)·V²/(2·9.81)
V = √(2·9.81·28/36.1) = 3.901 m/s

流量

Q = V·π·D²/4 = 3.901·0.04909 = 0.1915 m³/s = 191.5 L/s

答え

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

答え: η = 72.07 %

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.01·72 = 7.063 kW

効率

η = Ph/Pshaft = 7.063/9.8 = 72.07%

答え

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

答え: Pshaft = 12.44 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) = (27 + 8)·0.3305 = 11.57 m

ポンプ揚程

H = Δz + hL = 36 + 11.57 = 47.57 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.02·47.57/0.75 = 12.44 kW

答え

Pshaft = 12.44 kW
31.熱伝導

答え: Q̇ = 1937 W

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

各層の熱抵抗

R1 = L1/(k1A) = (150/1000)/(0.96·10) = 0.01563 K/W
R2 = L2/(k2A) = (25/1000)/(0.72·10) = 0.003472 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 37/0.0191 = 1937 W

答え

Q̇ = 1937 W
32.熱伝導

答え: Ts = 21.21 °C

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(10·12) = 0.008333 K/W
R1 = L1/(k1A) = (250/1000)/(0.96·12) = 0.0217 K/W
R2 = L2/(k2A) = (75/1000)/(0.035·12) = 0.1786 K/W
Ro = 1/(hoA) = 1/(30·12) = 0.002778 K/W
ΣR = 0.2114 K/W

熱損失

Q̇ = (22 − (2))/0.2114 = 94.61 W
Ts = Troom − Q̇·Ri = 22 − 94.61·0.008333 = 21.21 °C

答え

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

答え: Tw = 29.88 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.

ホイールの回転速度

nw = n·z/Nw = 2900·3/45 = 193.3 rpm

出力トルク

ωw = 2π·193.3/60 = 20.25 rad/s
Tw = η·P/ωw = 0.55·1100/20.25 = 29.88 N·m

答え

Tw = 29.88 N·m
34.歯車と歯車列

答え: nc = 250 rpm

遊星歯車装置:太陽歯車、リング、キャリヤの回転速度はウィリスの式 (ns − nc)/(nr − nc) = −Nr/Ns で結ばれます:キャリヤから見ると、遊星歯車をアイドラとする普通の歯車列です。1 つの要素を固定すると、単一の減速比をもつ減速機になります。

リングギヤの歯数

Nr = Ns + 2·Np = 15 + 2·12 = 39

リング固定時のウィリスの式

nr = 0 とおき、キャリヤの回転速度について解きます。

(900 − nc)/(0 − nc) = −39/15
nc = ns·Ns/(Ns + Nr) = 900·15/54 = 250 rpm

答え

nc = 250 rpm

キャリヤは太陽歯車と同じ向きに回転します。

35.ラプラス変換

答え: f(0.3) = 0.519

ラプラス変換: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 + 1) + C/(s + 5)
A = 20/(1·5) = 4
B = 20/((−1)·4) = −5
C = 20/((−5)·(−4)) = 1

逆変換

f(t) = 4 − 5·e−t + e−5t

計算する

f(0.3) = 4 − 5·e−0.3 + e−5·0.3 = 0.519

答え

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

答え: y(0.6) = 4.477

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

方程式を変換する

Y(s) = 60/(s(s + 1)(s + 4))

部分分数展開

A = 60/(1·4) = 15
B = 60/((−1)·3) = −20
C = 60/((−4)·(−3)) = 5

排泄量

y(t) = 15 − 20·e−t + 5·e−4t
y(0.6) = 4.477

答え

y(0.6) = 4.477
37.ブロック線図

答え: ess = 4.118

ブロック線図の代数:直列のブロックは掛け算(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ₚ = G(0) = 2·3/(4·7) = 0.2143

定常偏差

ess = A/(1 + Kₚ) = 5/(1 + 0.2143) = 4.118

答え

ess = 4.118
38.ブロック線図

答え: ω = 7.483 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³ + 14s² + 56s + 64 + K = 0

ラウス表

s³ | 1 56
s² | 14 64 + K
s¹ | (14·56 − 64 − K)/14
s⁰ | 64 + K

安定限界

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

Kmax = 14·56 − 64 = 720

振動の周波数

14·(jω)² + 64 + 720 = 0 → ω² = 784/14 = 56
ω = 7.483 rad/s

答え

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

答え: ts = 0.7273 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 + ωₙ² と一致しなければなりません。

ωₙ² = 83 → ωₙ = 9.11 rad/s
2ζωₙ = 11 → ζ = 11/(2·9.11) = 0.6037

結果

ζωₙ = 5.5 → ts ≈ 4/5.5 = 0.7273 s

答え

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

答え: %OS = 16.3 %

ステップ応答: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.01·1 × 10⁻⁶) = 10000 rad/s
ζ = (R/2)·√(C/L) = (100/2)·√(1 × 10⁻⁶/0.01) = 0.5

結果

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

答え

%OS = 16.3 %
41.定常偏差と PID 制御

答え: ess = 2.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

ランプ偏差

ess = 5/2 = 2.5

答え

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

答え: Ki = 25

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

特性多項式

1 + C(s)G(s) = 0 に s(s + 2) を掛けます。

s² + (2 + 1Kp)s + 1Ki

目標の多項式

s² + 2·0.5·5s + 5² = s² + 5s + 25

係数を合わせる

1Ki = 25 → Ki = 25

答え

Ki = 25

積分項によりループは 1 型になるので、ステップ入力に定常偏差なく追従します。

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

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