FE 机械考试练习:速成包

公式表、每个专题的关键方法步骤和 42 道练习题,完整解答附在最后。

新题目返回课程

在打印窗口中选择“另存为 PDF”,即可在设备上保存一份副本。

1公式与方法

梁的支座反力

平衡
ΣFx = 0, ΣFy = 0, ΣM = 0
均布荷载
合力 w·L 作用在荷载中点
力矩
M = F·d

方法

  1. 画受力图,并用支座反力代替每个支座。
  2. 用作用在形心处的单个集中力代替每个分布荷载。
  3. 对一个支座取矩,求另一个支座的反力。
  4. 用 ΣFy = 0 求剩下的反力,并用第二个力矩方程验算。

轴向应力与应变

正应力
σ = 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. 圆上的角度是单元体上角度的两倍: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
两个敞口水池
Δ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
二阶目标
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​ = 80 kNP2​ = 10 kNABCRB​ = ?2 m2.5 m1.5 m

梁在 A 处为铰支座,在距 A 为 4.5 m 的 B 处为滚动支座,并从 B 外伸至自由端 C。求竖向反力 RB,以向上为正。

AB = 4.5 mAC = 6 mP1 = 80 kN 位于 2 mP2 = 10 kN 位于 6 m
RB = kN
2.梁的支座反力困难
w = 12 kN/mP = 20 kNABRA​ = ?RB​7 m2 m1 m

简支梁在 A 处为铰支座,在 B 处为滚动支座。求竖向支座反力 RA。

L = 10 mP = 20 kN 位于 7 mw = 12 kN/m,从 7 到 9 m
RA = kN
3.轴向应力与应变中等
150 kNØ25110 GPaØ80/60500 mmδ = ?

一根钢杆(直径 25 mm,E = 200 GPa)放在一根青铜管(外径 80 mm,内径 60 mm,E = 110 GPa)内。二者长度均为 500 mm,在刚性板之间被 150 kN 的力一起压缩。 求缩短量 δ。

P = 150 kN杆: Ø25 mm, E = 200 GPa管: Ø80/Ø60 mm, E = 110 GPaL = 500 mm
δ = mm
4.轴向应力与应变困难
30 kN12400 mm200 mmØ20Ø30σ2​ = ?

一段钢(直径 20 mm,长 400 mm,E = 200 GPa)与一段铝(直径 30 mm,长 200 mm,E = 70 GPa)首尾相接,固定在两堵刚性墙之间。一个 30 kN 的轴向载荷作用在接合处,指向第 2 段。求第 2 段中的应力大小 σ2。

第 1 段: 钢, Ø20 mm, 400 mm, E = 200 GPa第 2 段: 铝, Ø30 mm, 200 mm, E = 70 GPaP = 30 kN
σ2 = MPa
5.轴的扭转中等
d1​ = 40 mmd2​ = 30 mmT = 750 N·mABC3 m0.8 m

一根阶梯形实心钢轴在 A 端固定。AB 段直径 40 mm,长 3 m;BC 段直径 30 mm,长 0.8 m。自由端 C 作用 750 N·m 的扭矩。 求端部 C 的扭转角 φC,以度表示。

T = 750 N·md1 = 40 mm, L1 = 3 md2 = 30 mm, L2 = 0.8 mG = 80 GPa
φC = °
6.轴的扭转困难
Td = ?

一根实心轴需在 1500 rpm 下传递 20 kW,且切应力不超过 60 MPa。求满足要求的最小直径 d。

P = 20 kWN = 1500 rpmτallow = 60 MPa
d = mm
7.剪力和弯矩中等
P1​ = 80 kNP2​ = 40 kNAB2 m2 m1 m

简支梁承受两个集中荷载。求梁中的最大弯矩 Mmax。

L = 5 mP1 = 80 kN 位于 2 mP2 = 40 kN 位于 4 m
Mmax = kN·m
8.剪力和弯矩困难
10 kN/m50 kN15 kNABCDEFG1 m1 m1 m1 m1 m1 m

标记 A 到 G 把梁分成六等份。在哪个标记点处弯矩的绝对值最大,不论是正弯矩还是负弯矩?请在图中点击它。

50 kN @ 6 m15 kN @ 2 m10 kN/m, 0–5 m
  1. 1A (x = 0 m)
  2. 2B (x = 1 m)
  3. 3C (x = 2 m)
  4. 4D (x = 3 m)
  5. 5E (x = 4 m)
  6. 6F (x = 5 m)
  7. 7G (x = 6 m)
9.莫尔圆与应力变换中等
σx​ = 150σy​ = 60τxy​ = 30MPa

受载零件中某点的应力为 σₓ = 150 MPa,σᵧ = 60 MPa,τₓᵧ = 30 MPa(拉为正)。 求从 x 轴到最大主应力方向的角度 θp,取值在 −90° 到 90° 之间。

σx = 150 MPaσy = 60 MPaτxy = 30 MPa
θp = °
10.莫尔圆与应力变换困难
σx​ = −50σy​ = −70τxy​ = −60MPa

受载零件中某点的应力为 σₓ = −50 MPa,σᵧ = −70 MPa,τₓᵧ = −60 MPa(拉为正)。 求法线由 x 轴逆时针转过 20° 的截面上的正应力 σθ。

σx = −50 MPaσy = −70 MPaτxy = −60 MPaθ = 20°
σθ = MPa
11.梁的挠度中等
P = 18 kNABδ4 m1 m

简支梁在距支座 A 为 a = 4 m 处承受集中荷载 P。求荷载作用处的挠度 δP。

P = 18 kNa = 4 mL = 5 mE = 200 GPa (钢)I = 8.5 × 10⁶ mm⁴
δP = mm
12.梁的挠度困难
P = 55 kNAB3 m3 m

一端固定一端简支的梁在 A 端固定,在 B 端简支(支撑)。梁在跨中承受集中荷载 P。EI 为常数。 求固定端弯矩 MA 的大小。

L = 6 mP = 55 kN 位于 3 m
MA = kN·m
13.柱的屈曲中等
PL = 6 md = 50 mm横截面

一根长为 L 的钢柱,两端铰接。用欧拉公式求其临界屈曲荷载 Pcr。

L = 6 mE = 200 GPa (钢)d = 50 mm
Pcr = kN
14.柱的屈曲困难
P = 20 kNL = 1 md = ?横截面

一根长为 L 的铝实心圆柱,两端铰接。它需承受 20 kN,屈曲安全系数为 2。求最小直径 d。

P = 20 kNn = 2L = 1 mE = 70 GPa (铝)
d = mm
15.抛体运动中等
v0​ = 18 m/s35°射程H = ?

一个球 从平地上以高于水平面 35° 的角度抛出。忽略空气阻力。求球达到的最大高度 H。

v0 = 18 m/sθ = 35°g = 9.81 m/s²
H = m
16.抛体运动困难
11 mv0​ = 25 m/s45°x = ?

一个球从地面以 25 m/s、与水平方向成 45° 的角度抛出。它在下落途中落在一栋 11 m 高的建筑的平屋顶上。忽略空气阻力,取 g = 9.81 m/s²。求从抛出点到落点的水平距离 x。

v₀ = 25 m/sθ = 45°hroof = 11 mg = 9.81 m/s²
x = m
17.功和能中等
90 km/h8°800 kg

一辆 800 kg 的汽车以恒定的 90 km/h 驶上 8° 的斜坡。空气阻力和滚动阻力合计 300 N。求车轮须输出的功率 P,以 kW 表示。

m = 800 kgθ = 8°Fr = 300 Nv = 90 km/hg = 9.81 m/s²
P = kW
18.功和能困难
F = 75 N35°μ = 0.18 kg

一个 8 kg 的物块从 35° 斜坡的底部由静止开始,被一个平行于斜坡、大小为 75 N 的恒力沿斜坡向上拉。动摩擦因数为 0.1。求物块沿斜坡向上移动 2 m 后的速度 v。

m = 8 kgF = 75 Nθ = 35°μ = 0.1d = 2 mg = 9.81 m/s²
v = m/s
19.机械振动中等
k = 100 kN/mζ = 0.3m = 500 kg

一个 500 kg 的质量块挂在刚度为 100 kN/m 的弹簧上,黏性阻尼器使阻尼比为 0.3。求其自由振动的有阻尼固有频率 fd,以 Hz 计。

m = 500 kgk = 100 kN/mζ = 0.3
fd = Hz
20.机械振动困难
k1​ = 40k2​ = 80k3​ = 20m = 25 kg

一个 25 kg 的质量块悬挂在弹簧 k₃ = 20 kN/m 上。k₃ 的上端连接在一根刚性杆上,该杆由并排的两根弹簧 k₁ = 40 kN/m 和 k₂ = 80 kN/m 悬挂;杆保持水平。 忽略阻尼以及弹簧和杆的质量。求竖向振动的固有频率 fn,以 Hz 计。

m = 25 kgk1 = 40 kN/mk2 = 80 kN/mk3 = 20 kN/m
fn = Hz
21.热机和制冷机中等
TH​ = 550 °CTL​ = 25 °CEQH​QL​η = 26%

一台热机从 550 °C 的热源吸热,向 25 °C 的环境放热。其实际热效率为 26%。求其第二定律效率 ηII,即实际效率与卡诺效率之比。

TH = 550 °CTL = 25 °Cη = 26%
ηII = %
22.热机和制冷机困难
TH​ = 800 °CTL​ = 10 °CEQH​ = 800 kJQL​ = 384 kJW = 416 kJ

一台热机每个循环从 800 °C 的热源吸收 800 kJ 热量,产生 416 kJ 的功,其余热量排向 10 °C 的冷源。求每个循环产生的总熵 Sgen。

TH = 800 °CTL = 10 °CQH = 800 kJW = 416 kJ
Sgen = kJ/K
23.气体动力循环中等
vP1234qin​qout​

一个理想狄塞尔循环的压缩比 r = 19,预胀比 rc = v₃/v₂ = 1.5。 压缩开始时空气温度为 300 K。 采用冷空气标准假设:空气为比热容恒定的理想气体,k = 1.4,cp = 1.005 kJ/(kg·K),cv = 0.718 kJ/(kg·K)。 求每千克空气的吸热量 qin,以 kJ/kg 计。

r = 19rc = 1.5T1 = 300 Kk = 1.4cp = 1.005 kJ/(kg·K)
qin = kJ/kg
24.气体动力循环困难
vP1234qin​qout​

一座燃气轮机电站按布雷顿循环运行,压气机压比为 8。空气以 300 K 进入压气机,以 1100 K 进入涡轮。 压气机的等熵效率为 80%,涡轮为 88%。 采用冷空气标准假设:空气为比热容恒定的理想气体,k = 1.4,cp = 1.005 kJ/(kg·K),cv = 0.718 kJ/(kg·K)。 求每千克空气的净输出功 wnet,以 kJ/kg 计。

rp = 8T1 = 300 KT3 = 1100 Kηc = 80%ηt = 88%k = 1.4cp = 1.005 kJ/(kg·K)
wnet = kJ/kg
25.连续性和伯努利中等
12d1​ = 200 mm, p1​ = 140 kPad2​ = 50 mm, p2​ = ?

水流经一根水平管道,管径从 200 mm 缩小到 50 mm。宽段压强为 140 kPa,流速为 0.75 m/s。忽略损失,求窄段的压强 p2。

d1 = 200 mmd2 = 50 mmp1 = 140 kPav1 = 0.75 m/sρ = 1000 kg/m³
p2 = kPa
26.连续性和伯努利困难
12d1​ = 75 mm, p1​ = 370 kPad2​ = 40 mm, p2​ = ?Δz = 2 m

水在一根管道中向上流动,管径从 75 mm 缩小到 40 mm,同时升高 2 m。截面 1 处压强为 370 kPa,流速为 2 m/s。忽略损失,求截面 2 处的压强 p2。

d1 = 75 mmd2 = 40 mmp1 = 370 kPav1 = 2 m/sΔz = 2 mρ = 1000 kg/m³
p2 = kPa
27.管道摩擦与水头损失中等
Q = 20 L/sD = 150 mmL = 400 mf = 0.015

水以 20 L/s 的流量流过一根长 400 m、内径 150 mm 的管道。 取达西摩擦系数 f = 0.015,g = 9.81 m/s²。求沿程水头损失 hf。

Q = 20 L/sL = 400 mD = 150 mmf = 0.015g = 9.81 m/s²
hf = m
28.管道摩擦与水头损失困难
K = 2Δz = 24 mL = 150 m, D = 250 mm, f = 0.0254× 弯头 (K = 0.3)

两个敞口水池由一根长 150 m、内径 250 mm、达西摩擦系数 f = 0.025 的管道连接。上游水池水面比下游水池水面高 24 m。局部损失包括一个直角进口(K = 0.5)、4 个弯头(每个 K = 0.3)、一个阀门(K = 2)以及流入下游水池的出口(K = 1.0)。取 g = 9.81 m/s²。求流量 Q,以 L/s 计。

Δz = 24 mL = 150 mD = 250 mmf = 0.025Kentrance = 0.5Kbend = 0.3 (×4)Kvalve = 2Kexit = 1.0
Q = L/s
29.水泵与泵功率中等
PQ = 100 L/sH = 72 m

在一次试验中,水泵以 72 m 的扬程输送 100 L/s 的水(ρ = 1000 kg/m³),轴输入功率为 86.1 kW。取 g = 9.81 m/s²。求水泵效率 η,以百分比计。

Q = 100 L/sH = 72 mPshaft = 86.1 kWg = 9.81 m/s²
η = %
30.水泵与泵功率困难
PQ = 80 L/sΔz = 26 mL = 500 m, D = 200 mm, f = 0.018, ΣK = 8

一台水泵把 80 L/s 的水(ρ = 1000 kg/m³)从敞口下游水池输送到水面高 26 m 的上游水箱。管道长 500 m,内径 200 mm,达西摩擦系数 f = 0.018;局部损失系数之和为 ΣK = 8(含进口和出口)。取 g = 9.81 m/s²。 水泵效率为 72%。求轴功率 Pshaft,以 kW 计。

Q = 80 L/sΔz = 26 mL = 500 mD = 200 mmf = 0.018ΣK = 8η = 72%
Pshaft = kW
31.热传导中等
混凝土, k = 1.4250 mm木材, k = 0.1725 mm20 °C-8 °CTi​ = ?Q̇

一面墙由 250 mm 的混凝土和 25 mm 的木材组成,面积为 20 m²。外表面温度为 20 °C 和 -8 °C。 求两层交界面处的温度 Ti。

L1 = 250 mm, k1 = 1.4 W/(m·K)L2 = 25 mm, k2 = 0.17 W/(m·K)A = 20 m²T1 = 20 °CT3 = -8 °C
Ti = °C
32.热传导困难
矿棉, k = 0.0450 mm木材, k = 0.17100 mm室内 20 °C室外 -13 °CQ̇ = ?

由 50 mm 矿棉和 100 mm 木材组成的墙,将 20 °C 的房间与 -13 °C 的室外空气隔开。对流换热系数室内为 8 W/(m²·K),室外为 25 W/(m²·K),墙面积为 8 m²。 求通过墙体的热损失 Q̇。

L1 = 50 mm, k1 = 0.04 W/(m·K)L2 = 100 mm, k2 = 0.17 W/(m·K)A = 8 m²hi = 8 W/(m²·K)ho = 25 W/(m²·K)Troom = 20 °CTout = -13 °C
Q̇ = W
33.齿轮与齿轮系中等
N1​ = 30N2​ = 48

一台 5.5 kW 的电动机以 2900 rpm 运转,驱动一个齿数为 30 的小齿轮,小齿轮与一个齿数为 48 的大齿轮啮合。 该齿轮副的效率为 95%。 求齿轮轴上的输出扭矩 T2。

P = 5.5 kWn1 = 2900 rpmN1 = 30N2 = 48η = 95%
T2 = N·m
34.齿轮与齿轮系困难
A = 45B = 80C = 60D = 72

一台 3 kW、2900 rpm 的电动机驱动齿轮 A(45 齿),A 与齿轮 B(80 齿)啮合。与 B 同轴的齿轮 C(60 齿)驱动输出轴上的齿轮 D(72 齿)。每对齿轮啮合的效率为 95%。求输出转矩 Tout。

P = 3 kWnin = 2900 rpmA = 45B = 80C = 60D = 72ηmesh = 95%
Tout = N·m
35.拉普拉斯变换中等
ReIm−1−2−3−4−5−61

F(s) = (2s + 8)/(s² + 8s + 12)。它的逆变换形如 f(t) = A·e−2t + B·e−6t。用部分分式法求系数 A。

F(s) = (2s + 8)/(s² + 8s + 12)
A =
36.拉普拉斯变换困难
t (s)y(t)t = 0.3y(0.3) = ?

用拉普拉斯变换解 y″ + 10y′ + 24y = 0,初始条件 y(0) = 5,y′(0) = −9,并求 y(0.3)。

y″ + 10y′ + 24y = 0y(0) = 5y′(0) = −9t = 0.3 s
y(0.3) =
37.方框图中等
R(s)+18s(s + 7)G(s)C(s)−

单位负反馈回路的 G(s) = 18/(s(s + 7))。参考输入是斜坡 r(t) = 10·t。求稳态误差 ess。

G(s) = 18/(s(s + 7))r(t) = 10·t
ess =
38.方框图困难
R(s)+Ks(s + 1)(s + 6)G(s)C(s)−

被控对象 G(s) = K/(s(s + 1)(s + 6)) 放在单位负反馈回路中。 在最大稳定增益时回路发生振荡。求振荡频率 ω,单位 rad/s。

G(s) = K/(s(s + 1)(s + 6))单位负反馈
ω = rad/s
39.阶跃响应中等
t (s)y(t)3?

某系统的传递函数为 G(s) = 189/(s² + 2s + 63),输入为单位阶跃。 求阶跃响应的百分比超调量 %OS。

G(s) = 189/(s² + 2s + 63)
%OS = %
40.阶跃响应困难
R(s)+Ks(s + 6)G(s)C(s)−

单位负反馈回路的 G(s) = K/(s(s + 6))。求使阶跃响应有 30 % 超调量的增益 K。

G(s) = K/(s(s + 6))%OS = 30 %
K =
41.稳态误差与 PID 控制中等
R(s)+K = 925s + 1C(s)−

过程 G(s) = 2/(5s + 1) 由比例增益 K = 9 在单位负反馈回路中控制。求阶跃输入的稳态误差 ess,以阶跃大小的百分数表示。

G(s) = 2/(5s + 1)K = 9
ess = %
42.稳态误差与 PID 控制困难
R(s)+Kp​ + Ki​/s4s + 4C(s)−

PI 控制器 C(s) = Kp + Ki/s 在单位负反馈回路中控制对象 G(s) = 4/(s + 4)。闭环必须满足 ωₙ = 10 rad/s 和 ζ = 1。求积分增益 Ki。

G(s) = 4/(s + 4)ωₙ = 10 rad/sζ = 1
Ki =

3溶液

1.梁的支座反力

答案: RB = 48.89 kN

静力平衡:静止的物体满足 ΣFx = 0、ΣFy = 0,且对任意点 ΣM = 0。把每个分布荷载替换为其合力(荷载图的面积),作用在荷载图的形心处;对某个支座取矩,可从方程中消去该支座的反力。

对 A 取矩

逆时针为正。

ΣMA = 0: RB·4.5 = 80·2 + 10·6
RB = 220/4.5 = 48.89 kN

竖向平衡

RA = 80 + 10 − 48.89 = 41.11 kN

答案

RB = 48.89 kN
2.梁的支座反力

答案: RA = 10.8 kN

静力平衡:静止物体对任意点满足 ΣFx = 0、ΣFy = 0 和 ΣM = 0。对一个支座取矩可从方程中消去该支座的反力。

w = 12 kN/mP = 20 kNABRA​RB​7 m2 m1 m

受力图

用反力代替支座。滚动支座提供一个竖向反力;只有竖向荷载时,铰支座也只提供一个竖向反力。

平衡

对 A 取矩,逆时针为正。然后列竖向力平衡。

均布荷载合力:W = w·ℓ = 12·2 = 24 kN,作用于 x̄ = 8 m
ΣMA = 0: RB·10 = 20·7 + 12·2·8
RB = 332/10 = 33.2 kN
ΣFy = 0: RA = 44 − 33.2 = 10.8 kN

答案

RA = 10.8 kN
3.轴向应力与应变

答案: δ = 0.2205 mm

两种材料:串联(一个接一个)时,每一段承受相同的力,伸长量相加:δ = P·Σ L/(E·A)。并联(在刚性板之间并排)时,两部分缩短相同的量,因此应变相等:刚度大的部分承担更大的份额,Pi = P·EiAi/ΣEA,σi = Ei·ε。夹在两堵刚性墙之间、载荷作用在接合处的杆是超静定的:两段按刚度 k = EA/L 的比例分担载荷。

面积

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

组合轴向刚度

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

缩短量

δ = 150000·500/3.401 × 10⁸ = 0.2205 mm

答案

δ = 0.2205 mm
4.轴向应力与应变

答案: σ2 = 25.96 MPa

两种材料:串联(一个接一个)时,每一段承受相同的力,伸长量相加:δ = P·Σ L/(E·A)。并联(在刚性板之间并排)时,两部分缩短相同的量,因此应变相等:刚度大的部分承担更大的份额,Pi = P·EiAi/ΣEA,σi = Ei·ε。夹在两堵刚性墙之间、载荷作用在接合处的杆是超静定的:两段按刚度 k = EA/L 的比例分担载荷。

刚度

k₁ = 200000·314.2/400 = 157080 N/mm
k₂ = 70000·706.9/200 = 247400 N/mm

第 2 段的份额

F2 = 30000·247400/404480 = 18350 N

应力

σ2 = 18350/706.9 = 25.96 MPa (压缩)

答案

σ2 = 25.96 MPa
5.轴的扭转

答案: φC = 11.82 °

圆轴扭转:切应力从圆心处的零线性增大到表面的 τmax = T·r/J,实心轴 J = πd⁴/32,空心轴 J = π(do⁴ − di⁴)/32。扭转角 φ = TL/(GJ),单位为弧度。功率与扭矩的关系为 P = T·ω。

极惯性矩

J1 = π·40⁴/32 = 251327 mm⁴
J2 = π·30⁴/32 = 79522 mm⁴

各段的扭转角

T = 750 N·m = 750000 N·mm
G = 80 GPa = 80000 MPa
φAB = 750000·3000/(80000·251327) = 0.1119 rad
φBC = 750000·800/(80000·79522) = 0.09431 rad

总扭转角

φC = 0.1119 + 0.09431 = 0.2062 rad = 11.82°

答案

φC = 11.82°
6.轴的扭转

答案: d = 22.11 mm

圆轴扭转:切应力从圆心处的零线性增大到表面的 τmax = T·r/J,实心轴 J = πd⁴/32,空心轴 J = π(do⁴ − di⁴)/32。扭转角 φ = TL/(GJ),单位为弧度。功率与扭矩的关系为 P = T·ω。

扭矩

ω = 2πN/60 = 2π·1500/60 = 157.1 rad/s
T = P/ω = 20000/157.1 = 127.3 N·m

由许用应力求直径

d = ∛(16T/(π·τallow)) = ∛(16·127324/(π·60)) = 22.11 mm

答案

d = 22.11 mm

实际中会向上取到下一个标准尺寸。

7.剪力和弯矩

答案: Mmax = 112 kN·m

截面法:在需要求内力的位置截开梁,保留一侧。剪力 V 平衡该侧的竖向力,弯矩 M 平衡这些力对截面的矩。 最大弯矩:由于 dM/dx = V,弯矩的峰值出现在剪力过零的位置,可能在集中荷载处,也可能在分布荷载范围内。

支座反力

整根梁的平衡。

ΣMA = 0: RB·5 = 80·2 + 40·4
RB = 320/5 = 64 kN
ΣFy = 0: RA = 120 − 64 = 56 kN

各荷载处的弯矩

保留每个荷载左侧的部分。

荷载点 P1: M = RA·2 = 56·2 = 112 kN·m
荷载点 P2: M = RB·1 = 64·1 = 64 kN·m

最大值

Mmax = 112 kN·m (荷载点 P1)

答案

Mmax = 112 kN·m
8.剪力和弯矩

答案: F (x = 5 m)

最大弯矩:只要剪力变化,M 的斜率就随之变化,所以它的峰值出现在剪力过零的地方:集中荷载下、均布荷载内部,或外伸梁的支座上方。悬臂梁固定端的弯矩在墙处最大。比较这些位置上 M 的大小,不论是正弯矩还是负弯矩。

支座反力

对一个支座取矩,然后用竖向力平衡。

R1 = 24 kN, R2 = 91 kN

每个标记点处的弯矩

把截面一侧所有力的力矩相加。使梁下凹(正弯矩)为正,使梁上凸(负弯矩)为负。

MA = 0, MB = 19, MC = 28, MD = 12, ME = −14, MF = −50 kN·m, MG = 0

答案

标记点 F: MF = −50 kN·m
9.莫尔圆与应力变换

答案: θp = 16.85 °

平面应力变换:在转过 θ 的截面上,σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。莫尔圆可同时表示所有截面:圆心位于 C = (σₓ + σᵧ)/2,半径为 R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主应力为 C ± R,面内最大切应力为 R。

στxyCσ1​σ2​

圆上的角度

在莫尔圆上,点 x 位于距 σ 轴 2θp 的角度处;象限由 σₓ − σᵧ 和 τₓᵧ 的符号确定。

2θp = atan2(2·30, 150 − 60) = 33.69°

单元体上的角度

θp = 16.85°

答案

θp = 16.85°

逆时针方向的角度为正。

10.莫尔圆与应力变换

答案: σθ = −90.91 MPa

平面应力变换:在转过 θ 的截面上,σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。莫尔圆可同时表示所有截面:圆心位于 C = (σₓ + σᵧ)/2,半径为 R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主应力为 C ± R,面内最大切应力为 R。

变换公式

σθ = −60 + 10·0.766 + −60·0.6428 = −90.91 MPa

答案

σθ = −90.91 MPa
11.梁的挠度

答案: δP = 11.29 mm

梁的挠度:挠度由曲率 M/(EI) 积分两次得到。常见情况有表可查,例如端部受集中力的悬臂梁为 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, a = 4000 mm, b = 1000 mm, L = 5000 mm
EI = 1.7 × 10¹² N·mm²

荷载作用处的挠度

δP = Pa²b²/(3EIL) = 18000·4000²·1000²/(3·1.7 × 10¹²·5000) = 11.29 mm

答案

δP = 11.29 mm
12.梁的挠度

答案: MA = 61.88 kN·m

梁的挠度:挠度由曲率 M/(EI) 积分两次得到。常见情况有表可查,例如端部受集中力的悬臂梁为 PL³/(3EI),全跨均布荷载的简支梁为 5wL⁴/(384EI)。各荷载产生的挠度可以叠加。 对超静定梁,去掉一个多余支座,并令该处挠度为零(变形协调条件)。

去掉支撑

去掉支撑后,梁是一根悬臂梁。其自由端在荷载作用下会产生挠度。

δload = P(L/2)²(3L − L/2)/(6EI) = 5PL³/(48EI)

支撑处的变形协调

支撑力 RB 单独作用时使自由端上抬 RBL³/(3EI)。支撑使 B 点挠度为零,所以两者必须抵消。EI 被约去。

RB·L³/(3EI) = 5PL³/(48EI) → RB = 5P/16 = 5·55/16 = 17.19 kN

固定端的弯矩

荷载和支撑力对 A 取矩。

MA = P·L/2 − RB·L = 165 − 103.1 = 61.88 kN·m

答案

MA = 61.88 kN·m
13.柱的屈曲

答案: Pcr = 16.82 kN

欧拉屈曲:细长柱在 Pcr = π²EI/(KL)² 时失稳,KL 为计算长度:两端铰接 K = 1,一端固定一端自由 2,一端固定一端铰接 0.7,两端固定 0.5。柱绕 I 最小的轴屈曲。该公式适用于临界应力低于屈服应力的细长柱。

最小惯性矩

圆形截面对任意轴的 I 都相同。

I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴

计算长度

KL = 1·6000 = 6000 mm

临界荷载

Pcr = π²EI/(KL)² = π²·200000·306796/6000² = 16822 N = 16.82 kN

答案

Pcr = 16.82 kN
14.柱的屈曲

答案: d = 32.96 mm

欧拉屈曲:细长柱在 Pcr = π²EI/(KL)² 时失稳,KL 为计算长度:两端铰接 K = 1,一端固定一端自由 2,一端固定一端铰接 0.7,两端固定 0.5。柱绕 I 最小的轴屈曲。该公式适用于临界应力低于屈服应力的细长柱。

所需临界荷载

Pcr = n·P = 2·20000 = 40000 N
KL = 1·1000 = 1000 mm

求解直径

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·40000·1000²/(π³·70000))^(1/4) = 32.96 mm

答案

d = 32.96 mm

校核:长细比 KL/r 为 121.4,柱为细长柱,适用欧拉公式。

15.抛体运动

答案: H = 5.433 m

抛体运动:无空气阻力时,x = vxt,y = h + vyt − ½gt²。两个方向相互独立,只共享同一时间。

分解初速度

水平速度保持不变;竖直运动具有向下的恒定加速度 g = 9.81 m/s²。

vx = v0·cos θ = 18·cos 35° = 14.74 m/s
vy = v0·sin θ = 18·sin 35° = 10.32 m/s

最高点

在最高点竖直速度为零:vy² = 2gH。

H = vy²/(2g) = 10.32²/(2·9.81) = 5.433 m

答案

H = 5.433 m
16.抛体运动

答案: x = 49.57 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 = 17.68 m/s, vy = 17.68 m/s

落到屋顶的时间

较大的根:球正在下降。

11 = 17.68·t − 4.905·t²
t = (17.68 + √(17.68² − 2·9.81·11))/9.81 = 2.804 s

水平距离

x = vx·t = 17.68·2.804 = 49.57 m

答案

x = 49.57 m
17.功和能

答案: P = 34.81 kW

动能定理:动能的变化等于所有力做的功。把重力和弹簧写成势能(mgh 和 ½kx²)后得到能量衡算:初始能量 = 终了能量 + 摩擦损失的能量 (μN·d)。 功率等于力乘以速度:P = F·v。

驱动力

m·g = 800·9.81 = 7848 N
mg·sin θ = 7848·sin 8° = 1092 N
F = 1092 + 300 = 1392 N

功率

v = 90/3.6 = 25 m/s
P = F·v = 1392·25 = 34806 W = 34.81 kW

答案

P = 34.81 kW
18.功和能

答案: v = 3.432 m/s

动能定理:动能的变化等于所有力做的功。把重力和弹簧写成势能(mgh 和 ½kx²)后得到能量衡算:初始能量 = 终了能量 + 摩擦损失的能量 (μN·d)。

沿斜面方向的力

m·g = 8·9.81 = 78.48 N
mg·sin θ = 45.01 N
μ·mg·cos θ = 6.429 N

功–能关系

½mv² = (75 − 51.44)·2 = 47.11 J
v = √(2·47.11/8) = 3.432 m/s

答案

v = 3.432 m/s
19.机械振动

答案: fd = 2.147 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 = √(100000/500) = 14.14 rad/s
fn = 14.14/(2π) = 2.251 Hz

有阻尼固有频率

fd = fn·√(1 − ζ²) = 2.251·√(1 − 0.3²) = 2.147 Hz

答案

fd = 2.147 Hz
20.机械振动

答案: fn = 4.168 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 = 40 + 80 = 120 kN/m
keq = 120·20/(120 + 20) = 17.14 kN/m

固有频率

ωn = √(17143/25) = 26.19 rad/s
fn = 26.19/(2π) = 4.168 Hz

答案

fn = 4.168 Hz
21.热机和制冷机

答案: ηII = 40.77 %

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 = 550 + 273.15 = 823.1 K
TL = 25 + 273.15 = 298.1 K

卡诺效率

ηCarnot = 1 − 298.1/823.1 = 63.78 %

第二定律效率

ηII = η/ηCarnot = 26/63.78 = 40.77 %

答案

ηII = 40.77 %
22.热机和制冷机

答案: Sgen = 0.6107 kJ/K

熵产:在一个完整循环中,工质没有净熵变。高温热源失去 QH/TH,低温热源得到 QL/TL。两者之和就是熵产,可逆循环为零,所有实际循环都为正。

排放的热量

QL = QH − W = 800 − 416 = 384 kJ

绝对温度

TH = 800 + 273.15 = 1073 K
TL = 10 + 273.15 = 283.1 K

热源的熵变

ΔSH = −800/1073 = −0.7455 kJ/K
ΔSL = 384/283.1 = 1.356 kJ/K
Sgen = 1.356 − 0.7455 = 0.6107 kJ/K

答案

Sgen = 0.6107 kJ/K

Sgen > 0,所以该热机是可能的,但不可逆。

23.气体动力循环

答案: qin = 489.5 kJ/kg

狄塞尔循环(压燃式发动机):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·190.4 = 974.1 K
T3 = 1.5·974.1 = 1461 K

定压吸热量

qin = 1.005·(1461 − 974.1) = 489.5 kJ/kg

答案

qin = 489.5 kJ/kg
24.气体动力循环

答案: wnet = 130 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 = 300·80.2857 = 543.4 K
T4s = 1100/80.2857 = 607.2 K

压气机功和涡轮功

wc = 1.005·(543.4 − 300)/0.8 = 305.8 kJ/kg
wt = 0.88·1.005·(1100 − 607.2) = 435.8 kJ/kg

净功

wnet = 435.8 − 305.8 = 130 kJ/kg

答案

wnet = 130 kJ/kg

回功比:wc/wt = 70.17%。

25.连续性和伯努利

答案: p2 = 68.28 kPa

连续性方程和伯努利方程:对不可压缩流体,各截面流量相同,Q = A₁v₁ = A₂v₂。沿无损失的流线,p + ½ρv² + ρgz 为常数,所以流体越快或越高的地方压强越低。

连续性和伯努利

用 Pa 计算,再除以 1000 得到 kPa。

v2 = v1·(d1/d2)² = 0.75·(200/50)² = 12 m/s
½ρ(v1² − v2²) = ½·1000·(0.75² − 12²)/1000 = −71.72 kPa
p2 = p1 + ½ρ(v1² − v2²) = 140 + (−71.72) = 68.28 kPa

答案

p2 = 68.28 kPa
26.连续性和伯努利

答案: p2 = 327.7 kPa

连续性方程和伯努利方程:对不可压缩流体,各截面流量相同,Q = A₁v₁ = A₂v₂。沿无损失的流线,p + ½ρv² + ρgz 为常数,所以流体越快或越高的地方压强越低。

连续性和伯努利

用 Pa 计算,再除以 1000 得到 kPa。

v2 = v1·(d1/d2)² = 2·(75/40)² = 7.031 m/s
½ρ(v1² − v2²) = ½·1000·(2² − 7.031²)/1000 = −22.72 kPa
ρgΔz = 1000·9.81·2/1000 = 19.62 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 370 + (−22.72) − 19.62 = 327.7 kPa

答案

p2 = 327.7 kPa
27.管道摩擦与水头损失

答案: hf = 2.611 m

管流:雷诺数 Re = V·D/ν = ρ·V·D/μ 用于区分层流(Re 约低于 2000)和湍流(约高于 4000)。沿管道的摩擦造成达西–魏斯巴赫水头损失 hf = f·(L/D)·V²/(2g)。管件、阀门、进口和出口带来局部损失 hm = ΣK·V²/(2g)。在两个敞口水池液面之间,能量方程给出 Δ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.015·(400/0.15)·1.132²/(2·9.81) = 2.611 m

答案

hf = 2.611 m
28.管道摩擦与水头损失

答案: Q = 240 L/s

管流:雷诺数 Re = V·D/ν = ρ·V·D/μ 用于区分层流(Re 约低于 2000)和湍流(约高于 4000)。沿管道的摩擦造成达西–魏斯巴赫水头损失 hf = f·(L/D)·V²/(2g)。管件、阀门、进口和出口带来局部损失 hm = ΣK·V²/(2g)。在两个敞口水池液面之间,能量方程给出 Δz = hf + hm,其中包含出口损失 K = 1。

损失系数

f·L/D = 0.025·150/0.25 = 15
ΣK = 0.5 + 4·0.3 + 2 + 1.0 = 4.7

能量方程

两个水面都与大气相通且基本静止,所以全部 Δz 都损失在管道中。

24 = (15 + 4.7)·V²/(2·9.81)
V = √(2·9.81·24/19.7) = 4.889 m/s

流量

Q = V·π·D²/4 = 4.889·0.04909 = 0.24 m³/s = 240 L/s

答案

Q = 240 L/s
29.水泵与泵功率

答案: η = 82.03 %

水泵功率:传给水的功率是水功率 Ph = ρ·g·Q·H = Q·Δp。由于水泵有损失,轴必须输入更多功率:Pshaft = Ph/ηpump,电动机取用 Pel = Pshaft/ηmotor。泵扬程 H 由能量方程求得:静扬程 + 压力水头 + 全部沿程和局部损失。

水功率

Ph = ρ·g·Q·H = 1000·9.81·0.1·72 = 70.63 kW

效率

η = Ph/Pshaft = 70.63/86.1 = 82.03%

答案

η = 82.03%
30.水泵与泵功率

答案: Pshaft = 47.43 kW

水泵功率:传给水的功率是水功率 Ph = ρ·g·Q·H = Q·Δp。由于水泵有损失,轴必须输入更多功率:Pshaft = Ph/ηpump,电动机取用 Pel = Pshaft/ηmotor。泵扬程 H 由能量方程求得:静扬程 + 压力水头 + 全部沿程和局部损失。

速度

V = Q/A = 0.08/0.03142 = 2.546 m/s
V²/(2g) = 0.3305 m

水头损失

hL = (f·L/D + ΣK)·V²/(2g) = (45 + 8)·0.3305 = 17.52 m

泵扬程

H = Δz + hL = 26 + 17.52 = 43.52 m

轴功率

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.08·43.52/0.72 = 47.43 kW

答案

Pshaft = 47.43 kW
31.热传导

答案: Ti = 4.645 °C

热阻:通过串联各层的稳态热流类似于通过串联电阻的电流:Q̇ = ΔT/ΣR,导热层 R = L/(kA),对流表面 R = 1/(hA)。各层的温降与其热阻成正比。

各层热阻

R1 = L1/(k1A) = (250/1000)/(1.4·20) = 0.008929 K/W
R2 = L2/(k2A) = (25/1000)/(0.17·20) = 0.007353 K/W

热流量

Q̇ = (T1 − T3)/(R1 + R2) = 28/0.01628 = 1720 W

交界面温度

第一层中温度下降 Q̇·R₁。

Ti = T1 − Q̇·R1 = 20 − 1720·0.008929 = 4.645 °C

答案

Ti = 4.645 °C
32.热传导

答案: Q̇ = 131.8 W

热阻:通过串联各层的稳态热流类似于通过串联电阻的电流:Q̇ = ΔT/ΣR,导热层 R = L/(kA),对流表面 R = 1/(hA)。各层的温降与其热阻成正比。

串联热阻

Ri = 1/(hiA) = 1/(8·8) = 0.01563 K/W
R1 = L1/(k1A) = (50/1000)/(0.04·8) = 0.1563 K/W
R2 = L2/(k2A) = (100/1000)/(0.17·8) = 0.07353 K/W
Ro = 1/(hoA) = 1/(25·8) = 0.005 K/W
ΣR = 0.2504 K/W

热损失

Q̇ = (20 − (-13))/0.2504 = 131.8 W

答案

Q̇ = 131.8 W
33.齿轮与齿轮系

答案: T2 = 27.53 N·m

齿轮系:相互啮合的齿轮在轮齿处的线速度相同,所以 n₁·N₁ = n₂·N₂:齿数较少的齿轮转得更快。装在同一根轴上的齿轮转速相同。功率被传递下去(扣除损失),所以转速降低时扭矩增大:T₂ = T₁·(N₂/N₁)·η,且 P = T·ω,其中 ω = 2π·n/60。

电动机扭矩

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

输出扭矩

T2 = 18.11·48/30·0.95 = 27.53 N·m

答案

T2 = 27.53 N·m
34.齿轮与齿轮系

答案: Tout = 19.02 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)·(60/72) = 1359 rpm
ωout = 2π·1359/60 = 142.4 rad/s

输出功率和转矩

Pout = 3000·0.95² = 2708 W
Tout = Pout/ωout = 2708/142.4 = 19.02 N·m

答案

Tout = 19.02 N·m
35.拉普拉斯变换

答案: A = 1

拉普拉斯变换: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² + 8s + 12 = (s + 2)(s + 6)
F(s) = A/(s + 2) + B/(s + 6)

遮盖法

乘以 (s + 2) 并令 s = −2:B 项消失。

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

答案

A = 1
(B = 1,所以 f(t) = e−2t + e−6t)
36.拉普拉斯变换

答案: y(0.3) = 2.253

拉普拉斯变换: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 − 5s − (−9) + 10(sY − 5) + 24Y = 0
Y(s) = (5s + 41)/((s + 4)(s + 6))

部分分式

A = 21/2 = 10.5
B = 11/(−2) = −5.5

结果

y(t) = 10.5·e−4t − 5.5·e−6t
y(0.3) = 2.253

答案

y(0.3) = 2.253

检验:y(0) = A + B = 5。

37.方框图

答案: ess = 3.889

方框图代数:串联的方框相乘(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);单位反馈时,对阶跃的稳态误差为 A/(1 + Kp),其中 Kp = G(0)。

速度误差系数

Kᵥ = lims→0 s·G(s) = 18/7 = 2.571

稳态误差

ess = A/Kᵥ = 10/2.571 = 3.889

答案

ess = 3.889
38.方框图

答案: ω = 2.449 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);单位反馈时,对阶跃的稳态误差为 A/(1 + Kp),其中 Kp = G(0)。

特征方程

s³ + 7s² + 6s + K = 0

劳斯阵列

s³ | 1 6
s² | 7 K
s¹ | (7·6 − K)/7
s⁰ | K

稳定性极限

第一列的所有元素都必须为正。

Kmax = 7·6 = 42

振荡频率

7·(jω)² + 42 = 0 → ω² = 42/7 = 6
ω = 2.449 rad/s

答案

ω = 2.449 rad/s
39.阶跃响应

答案: %OS = 67.1 %

阶跃响应:一阶系统 G(s) = K/(τs + 1) 对大小为 U 的阶跃的响应是 y(t) = K·U·(1 − e−t/τ):经过一个时间常数 τ 达到终值 K·U 的 63.2 %,经过 4τ 达到 98 %。二阶系统 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 + ωₙ² 对应。

ωₙ² = 63 → ωₙ = 7.937 rad/s
2ζωₙ = 2 → ζ = 2/(2·7.937) = 0.126

结果

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

答案

%OS = 67.1 %
40.阶跃响应

答案: K = 70.28

阶跃响应:一阶系统 G(s) = K/(τs + 1) 对大小为 U 的阶跃的响应是 y(t) = K·U·(1 − e−t/τ):经过一个时间常数 τ 达到终值 K·U 的 63.2 %,经过 4τ 达到 98 %。二阶系统 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² + 6s + K)

由超调量求阻尼比

ζ = −ln(0.3)/√(π² + ln²(0.3)) = 0.3579

增益

ωₙ = 6/(2ζ) = 8.383 rad/s
K = ωₙ² = 70.28

答案

K = 70.28
41.稳态误差与 PID 控制

答案: ess = 5.263 %

稳态误差:在单位负反馈回路中,设开环传递函数为 G(s),则误差为 E(s) = R(s)/(1 + G(s)),由终值定理得 ess = lims→0 s·E(s)。系统型别是 G(s) 在 s = 0 处的极点个数。误差常数为 Kp = lim G(s)、Kv = lim s·G(s) 和 Ka = lim s²·G(s),均取 s → 0。大小为 A 的阶跃输入给出 ess = A/(1 + Kp),斜坡 A·t 给出 A/Kv,抛物线 A·t²/2 给出 A/Ka。0 型回路跟踪阶跃有有限误差,1 型回路能无误差跟踪阶跃,跟踪斜坡则有有限误差。

位置误差常数

Kp = 9·2 = 18

阶跃误差

ess = 100/(1 + 18) = 5.263 %

答案

ess = 5.263 %
42.稳态误差与 PID 控制

答案: Ki = 25

P、PI 和 PID 控制:控制器 C(s) = Kp + Ki/s + Kd·s 作用于误差。比例增益使回路更快,并减小 0 型对象的误差,但无法将其消除。积分项增加一个 s = 0 处的极点:系统型别升高一级,阶跃输入可以被无误差地跟踪。微分项增加阻尼。选取增益时,写出闭环特征多项式,并与 s² + 2ζωns + ωn² 比较。然后 %OS = 100·e−ζπ/√(1−ζ²),tp = π/ωd,其中 ωd = ωn√(1 − ζ²),ts ≈ 4/(ζωn)。

特征多项式

1 + C(s)G(s) = 0,两边乘以 s(s + 4)。

s² + (4 + 4Kp)s + 4Ki

目标多项式

s² + 2·1·10s + 10² = s² + 20s + 100

匹配各项系数

4Ki = 100 → Ki = 25

答案

Ki = 25

积分项使回路成为 1 型,所以阶跃输入可以被无稳态误差地跟踪。

本内容为独立练习材料。ExamProblems 与该考试的主办机构没有隶属关系,也未获其认可,考试名称仅用于描述内容。 FE exam 由 NCEES 举办,该机构与 ExamProblems 无关,也不为其背书。

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