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​ = 45 kNP2​ = 35 kNABCRA​ = ?2 m4 m2 m

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

AB = 6 mAC = 8 mP1 = 45 kN 位于 2 mP2 = 35 kN 位于 8 m
RA = kN
2.梁的支座反力困难
w = 10 kN/mP = 10 kNABRA​ = ?RB​5 m2 m3 m

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

L = 10 mP = 10 kN 位于 5 mw = 10 kN/m,从 7 到 10 m
RA = kN
3.轴向应力与应变中等
Ø 12 mm铜2 mP = 90 kN

求铜杆在轴向荷载下的伸长量 δ。

P = 90 kNd = 12 mmL = 2 mE = 117 GPa
δ = mm
4.轴向应力与应变困难
Ø 40 mm铝2 mØ 12 mm黄铜1.5 mP = 60 kN

阶梯杆一端固定在墙上,自由端受荷载。求总伸长量 δ。

P = 60 kN第 1 段:Ø40 mm,2 m,铝(E = 70 GPa)第 2 段:Ø12 mm,1.5 m,黄铜(E = 100 GPa)
δ = mm
5.轴的扭转中等
T = 500 N·md = 20 mmL = 0.8 mφ = ?

一根钢实心轴一端固定在墙上,自由端受扭矩 T。求自由端的扭转角 φ,以度为单位。

T = 500 N·md = 20 mmL = 0.8 mG = 80 GPa
φ = °
6.轴的扭转困难
Td = ?

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

P = 30 kWN = 300 rpmτallow = 50 MPa
d = mm
7.剪力和弯矩中等
4 kN/m55 kNABCDEFG1 m1 m1 m1 m1 m1 m

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

55 kN @ 3 m4 kN/m, 0–6 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)
8.剪力和弯矩困难
12 kN/m40 kN45 kNABCDEFG1 m1 m1 m1 m1 m1 m

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

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

受载零件中某点的应力为 σₓ = −20 MPa,σᵧ = −50 MPa,τₓᵧ = −60 MPa(拉为正)。 求面内最大切应力 τmax。

σx = −20 MPaσy = −50 MPaτxy = −60 MPa
τmax = MPa
10.莫尔圆与应力变换困难
σx​ = −50σy​ = 20τxy​ = 30MPa

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

σx = −50 MPaσy = 20 MPaτxy = 30 MPaθ = 30°
σθ = MPa
11.梁的挠度中等
P = 13 kNABδ2.5 m1.5 m

长为 L 的悬臂梁在 A 端固定,并在距 A 为 a = 2.5 m 处承受集中荷载 P。 求荷载作用处的挠度 δP。

P = 13 kNa = 2.5 mL = 4 mE = 11 GPa (结构木材)I = 160 × 10⁶ mm⁴
δP = mm
12.梁的挠度困难
w = 8 kN/mδL = 2 m

简支梁承受均布荷载 w,其挠度不得超过跨度/360。求满足该限值的最小惯性矩 Ireq。

w = 8 kN/mL = 2 mE = 70 GPa (铝)
Ireq = × 10⁶ mm⁴
13.柱的屈曲中等
PL = 5 md = 20 mm横截面

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

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

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

P = 30 kNn = 3L = 6 mE = 70 GPa (铝)
d = mm
15.抛体运动中等
48 my = ?v0​ = 26 m/s30°

一个球从水平地面以 26 m/s、与水平方向成 30° 的角度朝 48 m 外的一面墙踢出。忽略空气阻力,取 g = 9.81 m/s²。球经过墙时的高度 y 是多少?

v₀ = 26 m/sθ = 30°x = 48 mg = 9.81 m/s²
y = m
16.抛体运动困难
h = 40 mv0​ = 21 m/s40°H = ?

一个球从 40 m 高的悬崖顶部以 21 m/s、与水平方向成 40° 的角度抛出。忽略空气阻力,取 g = 9.81 m/s²。求它高出悬崖脚下地面的最大高度 H。

v₀ = 21 m/sθ = 40°h = 40 mg = 9.81 m/s²
H = m
17.功和能中等
72 km/h8°800 kg

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

m = 800 kgθ = 8°Fr = 300 Nv = 72 km/hg = 9.81 m/s²
P = kW
18.功和能困难
h = 2 mL = 4 mμ = 0.2k = 4000 N/m

一个 0.5 kg 的物块从光滑坡道上 2 m 高处由静止释放。到达底部后,它通过一段长 4 m 的粗糙水平段(μ = 0.2),然后在光滑表面上撞上一根弹簧(k = 4000 N/m)。 求弹簧的最大压缩量 x。

m = 0.5 kgh = 2 mL = 4 mμ = 0.2k = 4000 N/mg = 9.81 m/s²
x = cm
19.机械振动中等
k = 10 kN/mζ = 0.35m = 15 kg

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

m = 15 kgk = 10 kN/mζ = 0.35
fd = Hz
20.机械振动困难
k1​ = 5k2​ = 40m = 10 kg

一个 10 kg 的质量块悬挂在并排(并联)的两根弹簧上,刚度分别为 k₁ = 5 kN/m 和 k₂ = 40 kN/m。质量块保持水平。 忽略阻尼和弹簧的质量。求竖向振动的固有频率 fn,以 Hz 计。

m = 10 kgk1 = 5 kN/mk2 = 40 kN/m
fn = Hz
21.热机和制冷机中等
TH​ = 400 °CTL​ = 20 °CEQH​QL​W

一台热机在 400 °C 的高温热源和 20 °C 的低温热源之间工作。 求可能的最高热效率 ηmax。

TH = 400 °CTL = 20 °C
ηmax = %
22.热机和制冷机困难
TH​ = 900 °CT0​ = 30 °CTL​ = -5 °CERWQ̇H​ = 150 kWQ̇L​ = ?

一台卡诺热机从 900 °C 的热源吸收 150 kW 热量,并向 30 °C 的环境放热。它的全部功用于驱动一台卡诺制冷机,该制冷机从 -5 °C 的冷冻室排出热量,也向 30 °C 的环境放热。求从冷冻室排出热量的速率 Q̇L。

TH = 900 °CT0 = 30 °CTL = -5 °CQ̇H = 150 kW
Q̇L = kW
23.气体动力循环中等
vP1234qin​qout​

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

r = 15rc = 2.2T1 = 300 Kk = 1.4
T3 = K
24.气体动力循环困难
vP1234qin​qout​

一个理想布雷顿循环的压气机压比为 14。 采用冷空气标准假设:空气为比热容恒定的理想气体,k = 1.4,cp = 1.005 kJ/(kg·K),cv = 0.718 kJ/(kg·K)。 求其热效率 η,以百分比计。

rp = 14k = 1.4
η = %
25.连续性和伯努利中等
12d1​ = 150 mm, p1​ = ?d2​ = 75 mm, p2​ = 90 kPa

水以 26 L/s 流过一根水平管道,管径从 150 mm 缩小到 75 mm。窄段的压强为 90 kPa。忽略损失,求宽段的压强 p1。

Q = 26 L/sd1 = 150 mmd2 = 75 mmp2 = 90 kPaρ = 1000 kg/m³
p1 = kPa
26.连续性和伯努利困难
12d1​ = 150 mm, p1​ = 220 kPad2​ = 75 mm, p2​ = ?Δz = 6 m

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

d1 = 150 mmd2 = 75 mmp1 = 220 kPav1 = 1.75 m/sΔz = 6 mρ = 1000 kg/m³
p2 = kPa
27.管道摩擦与水头损失中等
V = 2.8 m/sD = 50 mmL = 50 mf = 0.03

水以 2.8 m/s 的平均流速流过一根长 50 m、内径 50 mm 的管道。 取达西摩擦系数 f = 0.03,g = 9.81 m/s²。求沿程水头损失 hf。

V = 2.8 m/sL = 50 mD = 50 mmf = 0.03g = 9.81 m/s²
hf = m
28.管道摩擦与水头损失困难
K = 2Δz = 14 mL = 400 m, D = 200 mm, f = 0.0154× 弯头 (K = 0.3)

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

Δz = 14 mL = 400 mD = 200 mmf = 0.015Kentrance = 0.5Kbend = 0.3 (×4)Kvalve = 2Kexit = 1.0
Q = L/s
29.水泵与泵功率中等
PQ = 15 L/sH = 14 m

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

Q = 15 L/sH = 14 mPshaft = 3.2 kWg = 9.81 m/s²
η = %
30.水泵与泵功率困难
PQ = 60 L/sΔz = 11 mL = 400 m, D = 200 mm, f = 0.015, ΣK = 6

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

Q = 60 L/sΔz = 11 mL = 400 mD = 200 mmf = 0.015ΣK = 6η = 85%
Pshaft = kW
31.热传导中等
矿棉, k = 0.0450 mm木材, k = 0.1775 mm15 °C-1 °CQ̇ = ?

一面墙由 50 mm 的矿棉和 75 mm 的木材组成,面积为 20 m²。外表面温度为 15 °C 和 -1 °C。 求通过墙体的热流量 Q̇。

L1 = 50 mm, k1 = 0.04 W/(m·K)L2 = 75 mm, k2 = 0.17 W/(m·K)A = 20 m²T1 = 15 °CT3 = -1 °C
Q̇ = W
32.热传导困难
木材, k = 0.17150 mm发泡聚苯乙烯, k = 0.03550 mm室内 21 °C室外 -5 °CQ̇

由 150 mm 木材和 50 mm 发泡聚苯乙烯组成的墙,将 21 °C 的房间与 -5 °C 的室外空气隔开。对流换热系数室内为 10 W/(m²·K),室外为 30 W/(m²·K),墙面积为 15 m²。 求墙体的总传热系数 U。

L1 = 150 mm, k1 = 0.17 W/(m·K)L2 = 50 mm, k2 = 0.035 W/(m·K)A = 15 m²hi = 10 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = -5 °C
U = W/(m²·K)
33.齿轮与齿轮系中等
d1​ = 80 mmd2​ = 300 mm2880 rpmn2​ = ?

一个 80 mm 的带轮以 2880 rpm 通过皮带驱动一个 300 mm 的带轮,打滑率为 1 %。求从动轮的转速 n2。

n1 = 2880 rpmd1 = 80 mmd2 = 300 mms = 1 %
n2 = rpm
34.齿轮与齿轮系困难
A = 12B = 50C = 16D = ?

在一个复合轮系中,输入轴上的齿轮 A(12 齿)以 960 rpm 转动,驱动齿轮 B(50 齿)。齿轮 C(16 齿)与 B 固定在同一根轴上,驱动输出轴上的齿轮 D。输出轴须以 153.6 rpm 转动。齿轮 D 须有多少个齿 ND?

nin = 960 rpmA = 12B = 50C = 16nout = 153.6 rpm
ND =
35.拉普拉斯变换中等
ReIm−5−10−15−20

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

F(s) = (2s + 38)/(s² + 5s + 6)
A =
36.拉普拉斯变换困难
t (s)y(t)t = 1y(1) = ?

用拉普拉斯变换解 y″ + 4y′ + 3y = 0,初始条件 y(0) = 7,y′(0) = 5,并求 y(1)。

y″ + 4y′ + 3y = 0y(0) = 7y′(0) = 5t = 1 s
y(1) =
37.方框图中等
R(s)+6s + 9G(s)C(s)0.1H−

被控对象 G(s) = 6/(s + 9) 由一个负反馈回路控制,该回路的传感器增益为 H = 0.1。 求闭环的时间常数 τ。

G(s) = 6/(s + 9)H = 0.1
τ = s
38.方框图困难
R(s)+15s + 3G₁+1s + 4G₂C(s)0.5H₁−0.2H₂−

化简方框图,其中 G₁ = 15/(s + 3),G₂ = 1/(s + 4),G₂ 周围有内部负反馈 H₁ = 0.5,外部有负反馈 H₂ = 0.2。求 C(s)/R(s) 的直流增益 T(0)。

G₁ = 15/(s + 3)G₂ = 1/(s + 4)H₁ = 0.5H₂ = 0.2
T(0) =
39.阶跃响应中等
t (s)y(t)0.5

某系统的传递函数为 G(s) = 22.5/(s² + 8s + 45),输入为单位阶跃。 极点为 s = σ ± jωd。求实部 σ。

G(s) = 22.5/(s² + 8s + 45)
σ =
40.阶跃响应困难
R(s)+Ks(s + 20)G(s)C(s)−

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

G(s) = K/(s(s + 20))%OS = 10 %
K =
41.稳态误差与 PID 控制中等
r(t) = 2 + 5t+5s(s + 2)(s + 4)C(s)−

输入 r(t) = 2 + 5t(阶跃加斜坡)加到 G(s) = 5/(s(s + 2)(s + 4)) 的单位负反馈回路上。求稳态误差 ess。

G(s) = 5/(s(s + 2)(s + 4))r(t) = 2 + 5t
ess =
42.稳态误差与 PID 控制困难
R(s)+K = 410s(s + 10)C(s)−

比例控制器 K = 4 在单位负反馈回路中控制对象 G(s) = 10/(s(s + 10))。 求阶跃响应的峰值时间 tp。

G(s) = 10/(s(s + 10))K = 4
tp = s

3溶液

1.梁的支座反力

答案: RA = 18.33 kN

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

对 A 取矩

逆时针为正。

ΣMA = 0: RB·6 = 45·2 + 35·8
RB = 370/6 = 61.67 kN

竖向平衡

RA = 45 + 35 − 61.67 = 18.33 kN

答案

RA = 18.33 kN
2.梁的支座反力

答案: RA = 9.5 kN

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

w = 10 kN/mP = 10 kNABRA​RB​5 m2 m3 m

受力图

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

平衡

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

均布荷载合力:W = w·ℓ = 10·3 = 30 kN,作用于 x̄ = 8.5 m
ΣMA = 0: RB·10 = 10·5 + 10·3·8.5
RB = 305/10 = 30.5 kN
ΣFy = 0: RA = 40 − 30.5 = 9.5 kN

答案

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

答案: δ = 13.6 mm

胡克定律:在弹性范围内 σ = E·ε。结合 σ = P/A 和 ε = δ/L 得 δ = PL/(AE)。

应力

A = π·12²/4 = 113.1 mm²
σ = P/A = 90·1000/113.1 = 795.8 MPa

应变(胡克定律)

E for 铜 = 117 GPa = 117000 MPa.

ε = σ/E = 795.8/117000 = 0.006801 = 6801 µε

伸长量

δ = ε·L = 0.006801·2000 mm = 13.6 mm

答案

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

答案: δ = 9.322 mm

轴向受力构件:每个等截面段满足 δ = PL/(AE)。串联构件的伸长量相加。

两段受力相同

自由端承受 P,所以杆的每个截面都承受相同的轴力 P。只有面积和材料不同。

A1 = π·40²/4 = 1257 mm², A2 = π·12²/4 = 113.1 mm²
σ1 = 47.75 MPa
σ2 = 530.5 MPa

叠加各段伸长量

δ1 = σ1L1/E1 = 47.75·2000/70000 = 1.364 mm
δ2 = σ2L2/E2 = 530.5·1500/100000 = 7.958 mm
δ = δ1 + δ2 = 9.322 mm

答案

δ = 9.322 mm
5.轴的扭转

答案: φ = 18.24 °

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

极惯性矩

J = π·20⁴/32 = 15708 mm⁴

扭转角

φ = TL/(GJ) = 500·1000·800/(80000·15708) = 0.3183 rad
φ = 0.3183·180/π = 18.24°

答案

φ = 18.24°
6.轴的扭转

答案: d = 45.99 mm

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

扭矩

ω = 2πN/60 = 2π·300/60 = 31.42 rad/s
T = P/ω = 30000/31.42 = 954.9 N·m

由许用应力求直径

d = ∛(16T/(π·τallow)) = ∛(16·954930/(π·50)) = 45.99 mm

答案

d = 45.99 mm

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

7.剪力和弯矩

答案: D (x = 3 m)

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

支座反力

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

R1 = 39.5 kN, R2 = 39.5 kN

每个标记点处的弯矩

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

MA = 0, MB = 37.5, MC = 71, MD = 100.5 kN·m, ME = 71, MF = 37.5, MG = 0

答案

标记点 D: MD = 100.5 kN·m
8.剪力和弯矩

答案: E (x = 4 m)

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

支座反力

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

R1 = 37.75 kN, R2 = 95.25 kN

每个标记点处的弯矩

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

MA = 0, MB = 31.75, MC = 6.5, MD = −30.75, ME = −80 kN·m, MF = −40, MG = 0

答案

标记点 E: ME = −80 kN·m
9.莫尔圆与应力变换

答案: τmax = 61.85 MPa

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

στxyCτmax​

莫尔圆半径

τmax = R = √(((−20 − −50)/2)² + −60²) = 61.85 MPa

答案

τmax = 61.85 MPa
10.莫尔圆与应力变换

答案: σθ = −6.519 MPa

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

变换公式

σθ = −15 + −35·0.5 + 30·0.866 = −6.519 MPa

答案

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

答案: δP = 38.47 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 = 13000 N, a = 2500 mm, L = 4000 mm
EI = 1.76 × 10¹² N·mm²

荷载作用处的挠度

δP = Pa³/(3EI) = 13000·2500³/(3·1.76 × 10¹²) = 38.47 mm

答案

δP = 38.47 mm
12.梁的挠度

答案: Ireq = 4.286 × 10⁶ mm⁴

梁的挠度:挠度由曲率 M/(EI) 积分两次得到。常见情况有表可查,例如端部受集中力的悬臂梁为 PL³/(3EI),全跨均布荷载的简支梁为 5wL⁴/(384EI)。各荷载产生的挠度可以叠加。

挠度限值

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

所需惯性矩

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

答案

Ireq = 4.286 × 10⁶ mm⁴
13.柱的屈曲

答案: Pcr = 2.481 kN

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

最小惯性矩

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

I = πd⁴/64 = π·20⁴/64 = 7854 mm⁴

计算长度

KL = 0.5·5000 = 2500 mm

临界荷载

Pcr = π²EI/(KL)² = π²·200000·7854/2500² = 2481 N = 2.481 kN

答案

Pcr = 2.481 kN
14.柱的屈曲

答案: d = 98.87 mm

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

所需临界荷载

Pcr = n·P = 3·30000 = 90000 N
KL = 1·6000 = 6000 mm

求解直径

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

答案

d = 98.87 mm

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

15.抛体运动

答案: y = 5.423 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 = 26·cos 30° = 22.52 m/s
t = x/vx = 48/22.52 = 2.132 s

此时的高度

y = 13·2.132 − ½·9.81·2.132² = 5.423 m

答案

y = 5.423 m

球到达墙时已经在下落。

16.抛体运动

答案: H = 49.29 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.

竖直方向的初速度

vy = 21·sin 40° = 13.5 m/s

高出悬崖的上升高度

Δh = vy²/(2g) = 13.5²/(2·9.81) = 9.287 m

离地面的高度

H = 40 + 9.287 = 49.29 m

答案

H = 49.29 m
17.功和能

答案: P = 27.84 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 = 72/3.6 = 20 m/s
P = F·v = 1392·20 = 27845 W = 27.84 kW

答案

P = 27.84 kW
18.功和能

答案: x = 5.425 cm

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

通过粗糙段后剩余的能量

坡道光滑,只有粗糙段消耗能量。

mgh − μmgL = 4.905·(2 − 0.2·4) = 5.886 J

最大压缩量

½kx² = 5.886 → x = √(2·5.886/4000) = 0.05425 m = 5.425 cm

答案

x = 5.425 cm
19.机械振动

答案: fd = 3.849 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 = √(10000/15) = 25.82 rad/s
fn = 25.82/(2π) = 4.109 Hz

有阻尼固有频率

fd = fn·√(1 − ζ²) = 4.109·√(1 − 0.35²) = 3.849 Hz

答案

fd = 3.849 Hz
20.机械振动

答案: fn = 10.68 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 = 5 + 40 = 45 kN/m

固有频率

ωn = √(45000/10) = 67.08 rad/s
fn = 67.08/(2π) = 10.68 Hz

答案

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

答案: ηmax = 56.45 %

第二定律的限制:在绝对温度为 TH 和 TL 的热源之间工作的任何循环都不能超过可逆(卡诺)循环:ηmax = 1 − TL/TH,COPR,max = TL/(TH − TL),COPHP,max = TH/(TH − TL)。能量总是守恒:QH = W + QL。

绝对温度

TH = 400 + 273.15 = 673.1 K
TL = 20 + 273.15 = 293.1 K

卡诺效率

ηmax = 1 − TL/TH = 1 − 293.1/673.1 = 56.45 %

答案

ηmax = 56.45 %
22.热机和制冷机

答案: Q̇L = 852.2 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 = 900 + 273.15 = 1173 K
T0 = 30 + 273.15 = 303.1 K
TL = -5 + 273.15 = 268.1 K

热机的功

η = 1 − 303.1/1173 = 0.7416
Ẇ = 0.7416·150 = 111.2 kW

制冷机排出的热量

COPR = 268.1/(303.1 − 268.1) = 7.661
Q̇L = 7.661·111.2 = 852.2 kW

答案

Q̇L = 852.2 kW
23.气体动力循环

答案: T3 = 1950 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·150.4 = 886.3 K

定压加热

T3 = rc·T2 = 2.2·886.3 = 1950 K

答案

T3 = 1950 K
24.气体动力循环

答案: η = 52.95 %

布雷顿循环(燃气轮机):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/140.2857 = 1 − 1/2.126 = 52.95%

答案

η = 52.95%
25.连续性和伯努利

答案: p1 = 106.2 kPa

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.

流速

v1 = Q/A1 = 0.026/0.01767 = 1.471 m/s
v2 = Q/A2 = 0.026/0.004418 = 5.885 m/s

伯努利方程

p1 = p2 + ½ρ(v2² − v1²) = 90 + ½·1000·(5.885² − 1.471²)/1000 = 90 + 16.24 = 106.2 kPa

答案

p1 = 106.2 kPa
26.连续性和伯努利

答案: p2 = 138.2 kPa

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

连续性和伯努利

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

v2 = v1·(d1/d2)² = 1.75·(150/75)² = 7 m/s
½ρ(v1² − v2²) = ½·1000·(1.75² − 7²)/1000 = −22.97 kPa
ρgΔz = 1000·9.81·6/1000 = 58.86 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 220 + (−22.97) − 58.86 = 138.2 kPa

答案

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

答案: hf = 11.99 m

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

速度水头

V²/(2g) = 2.8²/(2·9.81) = 0.3996 m

达西–魏斯巴赫公式

hf = f·(L/D)·V²/(2g) = 0.03·(50/0.05)·2.8²/(2·9.81) = 11.99 m

答案

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

答案: Q = 88.39 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.015·400/0.2 = 30
ΣK = 0.5 + 4·0.3 + 2 + 1.0 = 4.7

能量方程

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

14 = (30 + 4.7)·V²/(2·9.81)
V = √(2·9.81·14/34.7) = 2.814 m/s

流量

Q = V·π·D²/4 = 2.814·0.03142 = 0.08839 m³/s = 88.39 L/s

答案

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

答案: η = 64.38 %

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

水功率

Ph = ρ·g·Q·H = 1000·9.81·0.015·14 = 2.06 kW

效率

η = Ph/Pshaft = 2.06/3.2 = 64.38%

答案

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

答案: Pshaft = 12.25 kW

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

速度

V = Q/A = 0.06/0.03142 = 1.91 m/s
V²/(2g) = 0.1859 m

水头损失

hL = (f·L/D + ΣK)·V²/(2g) = (30 + 6)·0.1859 = 6.693 m

泵扬程

H = Δz + hL = 11 + 6.693 = 17.69 m

轴功率

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.06·17.69/0.85 = 12.25 kW

答案

Pshaft = 12.25 kW
31.热传导

答案: Q̇ = 189.2 W

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

各层热阻

R1 = L1/(k1A) = (50/1000)/(0.04·20) = 0.0625 K/W
R2 = L2/(k2A) = (75/1000)/(0.17·20) = 0.02206 K/W

热流量

Q̇ = (T1 − T3)/(R1 + R2) = 16/0.08456 = 189.2 W

答案

Q̇ = 189.2 W
32.热传导

答案: U = 0.4091 W/(m²·K)

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

串联热阻

Ri = 1/(hiA) = 1/(10·15) = 0.006667 K/W
R1 = L1/(k1A) = (150/1000)/(0.17·15) = 0.05882 K/W
R2 = L2/(k2A) = (50/1000)/(0.035·15) = 0.09524 K/W
Ro = 1/(hoA) = 1/(30·15) = 0.002222 K/W
ΣR = 0.163 K/W

总传热系数

U = 1/(A·ΣR) = 1/(15·0.163) = 0.4091 W/(m²·K)

答案

U = 0.4091 W/(m²·K)
33.齿轮与齿轮系

答案: n2 = 760.3 rpm

皮带传动:无打滑时,两个带轮的轮缘线速度相同,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μθ,其中 θ 是小带轮上的包角(弧度);对开口皮带,θ = π − 2·asin((d₂ − d₁)/(2C)),C 为中心距。

无打滑时的转速

2880·80/300 = 768 rpm

有打滑

n2 = 768·0.99 = 760.3 rpm

答案

n2 = 760.3 rpm
34.齿轮与齿轮系

答案: ND = 24

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.

中间轴转速

nB = nC = 960·12/50 = 230.4 rpm

齿轮 D 的齿数

ND = C·nC/nout = 16·230.4/153.6 = 24

答案

ND = 24
35.拉普拉斯变换

答案: A = 34

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

遮盖法

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

A = (2s + 38)/(s + 3) |s = −2 = 34/1 = 34

答案

A = 34
(B = −32,所以 f(t) = 34·e−2t − 32·e−3t)
36.拉普拉斯变换

答案: y(1) = 4.484

拉普拉斯变换: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 − 7s − 5 + 4(sY − 7) + 3Y = 0
Y(s) = (7s + 33)/((s + 1)(s + 3))

部分分式

A = 26/2 = 13
B = 12/(−2) = −6

结果

y(t) = 13·e−t − 6·e−3t
y(1) = 4.484

答案

y(1) = 4.484

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

37.方框图

答案: τ = 0.1042 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)。

闭环传递函数

T(s) = G/(1 + G·H) = 6/(s + 9 + 0.6) = 6/(s + 9.6)

结果

τ = 1/9.6 = 0.1042 s

答案

τ = 0.1042 s
38.方框图

答案: T(0) = 0.9091

方框图代数:串联的方框相乘(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)。

内环

Gᵢ = G₂/(1 + G₂H₁) = 1/(s + 4.5)

外环

T(s) = G₁Gᵢ/(1 + G₁GᵢH₂) = 15/((s + 3)(s + 4.5) + 3)

直流增益

T(0) = 15/(3·4.5 + 3) = 15/16.5 = 0.9091

答案

T(0) = 0.9091
39.阶跃响应

答案: σ = −4

阶跃响应:一阶系统 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 + ωₙ² 对应。

ωₙ² = 45 → ωₙ = 6.708 rad/s
2ζωₙ = 8 → ζ = 8/(2·6.708) = 0.5963

结果

σ = −ζωₙ = −8/2 = −4

答案

σ = −4
40.阶跃响应

答案: K = 286.2

阶跃响应:一阶系统 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² + 20s + K)

由超调量求阻尼比

ζ = −ln(0.1)/√(π² + ln²(0.1)) = 0.5912

增益

ωₙ = 20/(2ζ) = 16.92 rad/s
K = ωₙ² = 286.2

答案

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

答案: ess = 8

稳态误差:在单位负反馈回路中,设开环传递函数为 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) = 0.625

阶跃部分

Kp = ∞ → estep = 0

斜坡误差

ess = 5/0.625 = 8

答案

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

答案: tp = 0.8112 s

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)。

特征多项式

s(s + 10) + 4·10 = s² + 10s + 40

ωₙ 和 ζ

ωₙ = √40 = 6.325 rad/s
ζ = 10/(2·6.325) = 0.7906

峰值时间

ωd = 6.325·√(1 − 0.7906²) = 3.873 rad/s
tp = π/3.873 = 0.8112 s

答案

tp = 0.8112 s

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

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