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​ = 50 kNP2​ = 50 kNABRA​RB​ = ?3 m3 m2 m

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

L = 8 mP1 = 50 kN 位于 3 mP2 = 50 kN 位于 6 m
RB = kN
2.梁的支座反力困难
P = 80 kN, 30°Q = 50 kNABRA​ = ?0.5 m1 m2.5 m

简支梁在 A 处为铰支座,在 B 处为滚动支座。距 A 为 0.5 m 处作用一个荷载 P,与水平方向成 30°;距 A 为 1.5 m 处作用一个竖向荷载 Q。 求铰支座 A 处总反力 RA 的大小。

L = 4 mP = 80 kN 位于 0.5 m, θ = 30°Q = 50 kN 位于 1.5 m
RA = kN
3.轴向应力与应变中等
Ø 12 mm黄铜2.5 mP = 35 kN

求黄铜杆的轴向应变 ε,以微应变表示 (µε = 10⁻⁶)。

P = 35 kNd = 12 mmL = 2.5 mE = 100 GPa
ε = µε
4.轴向应力与应变困难
P = ?Ø4070 GPaØ120/60300 mm

一根钢杆(直径 40 mm,E = 200 GPa,许用应力 150 MPa)放在一根铝管(外径 120 mm,内径 60 mm,E = 70 GPa,许用应力 70 MPa)内。二者在刚性板之间被一起压缩。求最大载荷 Pmax。

杆: Ø40 mm, E = 200 GPa, σallow = 150 MPa管: Ø120/Ø60 mm, E = 70 GPa, σallow = 70 MPa
Pmax = kN
5.轴的扭转中等
d = ?T = 750 N·mL = 1.5 m

一根长 1.5 m 的实心钢轴承受 750 N·m 的扭矩。两端之间的扭转角不得超过 3°。求满足该限值的最小直径 d。

T = 750 N·mL = 1.5 mG = 80 GPaφallow = 3°
d = mm
6.轴的扭转困难
Td = ?

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

P = 75 kWN = 1450 rpmτallow = 80 MPa
d = mm
7.剪力和弯矩中等
15 kN/m55 kNABCDEFG1 m1 m1 m1 m1 m1 m

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

55 kN @ 3 m15 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)
8.剪力和弯矩困难
4 kN/m50 kN10 kN35 kNABCDEFG2 m2 m2 m2 m2 m2 m

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

50 kN @ 0 m10 kN @ 12 m35 kN @ 4 m4 kN/m, 2–10 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)
9.莫尔圆与应力变换中等
σx​ = −30σy​ = 70τxy​ = −50MPa

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

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

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

σx = 60 MPaσy = 120 MPaτxy = −50 MPaθ = 25°
σθ = MPa
11.梁的挠度中等
w = 10 kN/mP = 18 kNABδ1.5 m

悬臂梁全长承受均布荷载 w,并在自由端承受集中荷载 P。用叠加法求自由端的挠度 δ。

P = 18 kNw = 10 kN/mL = 1.5 mE = 70 GPa (铝)I = 57 × 10⁶ mm⁴
δ = mm
12.梁的挠度困难
w = 6 kN/mAB8 m

一端固定一端简支的梁在 A 端固定,在 B 端简支(支撑)。梁全长承受均布荷载 w。EI 为常数。 求支撑处的反力 RB。

L = 8 mw = 6 kN/m
RB = kN
13.柱的屈曲中等
PL = 2.5 m40 × 100 mm横截面

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

L = 2.5 mE = 70 GPa (铝)b × h = 40 × 100 mm
Pcr = kN
14.柱的屈曲困难
PL = 1 md = 20 mm横截面

一根长为 L 的铝柱,两端铰接。求其长细比,并用欧拉公式求临界屈曲应力 σcr。

L = 1 mE = 70 GPa (铝)d = 20 mm
σcr = MPa
15.抛体运动中等
v0​ = 34 m/s30°射程

一个球 从平地上以高于水平面 30° 的角度抛出。忽略空气阻力。求球落地前的飞行时间 t。

v0 = 34 m/sθ = 30°g = 9.81 m/s²
t = s
16.抛体运动困难
22 mv0​ = 32 m/s55°x = ?

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

v₀ = 32 m/sθ = 55°hroof = 22 mg = 9.81 m/s²
x = m
17.功和能中等
30°d = 15 mμ = 0.4

一个物块由静止沿倾角 30° 的斜面下滑。动摩擦因数为 0.4。求它沿斜面滑行 15 m 后的速度 v。

θ = 30°μ = 0.4d = 15 mg = 9.81 m/s²
v = m/s
18.功和能困难
h = 2.5 mL = 4 mμ = 0.3k = 4000 N/m

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

m = 2 kgh = 2.5 mL = 4 mμ = 0.3k = 4000 N/mg = 9.81 m/s²
x = cm
19.机械振动中等
tx30 mm3 个周期17 mm

在有阻尼质量–弹簧系统的自由振动试验中,位移峰值 30 mm 之后 3 个周期出现峰值 17 mm。 求阻尼比 ζ。

x0 = 30 mmx3 = 17 mmn = 3
ζ =
20.机械振动困难
k1​ = 80k2​ = 12k3​ = 12m = 10 kg

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

m = 10 kgk1 = 80 kN/mk2 = 12 kN/mk3 = 12 kN/m
fn = Hz
21.热机和制冷机中等
TH​ = 225 °CTL​ = 20 °CEQH​QL​W

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

TH = 225 °CTL = 20 °C
ηmax = %
22.热机和制冷机困难
TH​ = 600 °CT0​ = 27 °CTL​ = -5 °CERWQ̇H​ = 160 kWQ̇L​ = ?

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

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

一个理想狄塞尔循环的压缩比 r = 15,预胀比 rc = v₃/v₂ = 2.5。 采用冷空气标准假设:空气为比热容恒定的理想气体,k = 1.4,cp = 1.005 kJ/(kg·K),cv = 0.718 kJ/(kg·K)。 求其热效率 η,以百分比计。

r = 15rc = 2.5k = 1.4
η = %
24.气体动力循环困难
vP1234qin​qout​

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

rp = 15k = 1.4
η = %
25.连续性和伯努利中等
12d1​ = 75 mm, p1​ = 110 kPad2​ = 40 mm, p2​ = ?

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

d1 = 75 mmd2 = 40 mmp1 = 110 kPav1 = 2.75 m/sρ = 1000 kg/m³
p2 = kPa
26.连续性和伯努利困难
C213 m3.5 md = 40 mmQ = ?

一根内径 40 mm 的软管从敞口水箱中虹吸取水。软管的最高点 C 高出水面 3 m,出口(通大气)低于水面 3.5 m。忽略损失,取 g = 9.81 m/s²。 求流量 Q。

d = 40 mmzC = +3 mz2 = −3.5 mg = 9.81 m/s²
Q = L/s
27.管道摩擦与水头损失中等
V = 1.6 m/sD = 200 mmL = 500 mf = 0.015

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

V = 1.6 m/sL = 500 mD = 200 mmf = 0.015g = 9.81 m/s²
hf = m
28.管道摩擦与水头损失困难
K = 2Δz = ?L = 100 m, D = 250 mm, f = 0.0222× 弯头 (K = 0.3)

两个敞口水池由一根长 100 m、内径 250 mm、达西摩擦系数 f = 0.022 的管道连接。局部损失包括一个直角进口(K = 0.5)、2 个弯头(每个 K = 0.3)、一个阀门(K = 2)以及流入下游水池的出口(K = 1.0)。取 g = 9.81 m/s²。两水池水位差 Δz 为多少时流量为 75 L/s?

Q = 75 L/sL = 100 mD = 250 mmf = 0.022Kentrance = 0.5Kbend = 0.3 (×2)Kvalve = 2Kexit = 1.0
Δz = m
29.水泵与泵功率中等
PQ = 25 L/sH = 16 m

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

Q = 25 L/sH = 16 mPshaft = 4.6 kWg = 9.81 m/s²
η = %
30.水泵与泵功率困难
PQ = 20 L/sΔz = 5 mL = 200 m, D = 100 mm, f = 0.015, ΣK = 6

一台水泵把 20 L/s 的水(ρ = 1000 kg/m³)从敞口下游水池输送到水面高 5 m 的上游水箱。管道长 200 m,内径 100 mm,达西摩擦系数 f = 0.015;局部损失系数之和为 ΣK = 6(含进口和出口)。取 g = 9.81 m/s²。 求水泵须提供的扬程 H。

Q = 20 L/sΔz = 5 mL = 200 mD = 100 mmf = 0.015ΣK = 6
H = m
31.热传导中等
矿棉, k = 0.0450 mm混凝土, k = 1.4150 mm21 °C-2 °CQ̇ = ?

一面墙由 50 mm 的矿棉和 150 mm 的混凝土组成,面积为 8 m²。外表面温度为 21 °C 和 -2 °C。 求通过墙体的热流量 Q̇。

L1 = 50 mm, k1 = 0.04 W/(m·K)L2 = 150 mm, k2 = 1.4 W/(m·K)A = 8 m²T1 = 21 °CT3 = -2 °C
Q̇ = W
32.热传导困难
砖, k = 0.72250 mm发泡聚苯乙烯, k = 0.035150 mm室内 21 °C室外 -5 °CQ̇ = ?

由 250 mm 砖和 150 mm 发泡聚苯乙烯组成的墙,将 21 °C 的房间与 -5 °C 的室外空气隔开。对流换热系数室内为 8 W/(m²·K),室外为 20 W/(m²·K),墙面积为 15 m²。 求通过墙体的热损失 Q̇。

L1 = 250 mm, k1 = 0.72 W/(m·K)L2 = 150 mm, k2 = 0.035 W/(m·K)A = 15 m²hi = 8 W/(m²·K)ho = 20 W/(m²·K)Troom = 21 °CTout = -5 °C
Q̇ = W
33.齿轮与齿轮系中等
N1​N2​ = ?

两根平行轴相距 337.5 mm。须用一对模数为 5 mm、减速比为 3.5 : 1 的直齿圆柱齿轮把它们连接起来。 大齿轮须有多少个齿 N2?

a = 337.5 mmm = 5 mmi = 3.5
N2 =
34.齿轮与齿轮系困难
A = 24B = 25C = 15D = ?

在一个复合轮系中,输入轴上的齿轮 A(24 齿)以 1500 rpm 转动,驱动齿轮 B(25 齿)。齿轮 C(15 齿)与 B 固定在同一根轴上,驱动输出轴上的齿轮 D。输出轴须以 1080 rpm 转动。齿轮 D 须有多少个齿 ND?

nin = 1500 rpmA = 24B = 25C = 15nout = 1080 rpm
ND =
35.拉普拉斯变换中等
t (s)f(t)t = 0.8f(0.8) = ?

求 F(s) = (5s + 24)/(s² + 2s + 1) 对应的 f(t),并计算 f(0.8)。

F(s) = (5s + 24)/(s² + 2s + 1)t = 0.8 s
f(0.8) =
36.拉普拉斯变换困难
ReIm−0.5−1−1.5−20.51

某系统由微分方程 y″ + 3y′ + 2y = 4u′ + 3u 描述,输入为 u(t),输出为 y(t)。求它的传递函数 G(s) = Y(s)/U(s)(零初始条件),并在 t = 0.8 s 时计算冲激响应 g(t)。求 g(0.8)。

y″ + 3y′ + 2y = 4u′ + 3ut = 0.8 s
g(0.8) =
37.方框图中等
R(s)+Ks(s + 5)G(s)C(s)−

被控对象 G(s) = K/(s(s + 5)) 放在单位负反馈回路中。求使闭环处于临界阻尼(两个相等的实极点)的增益 K。

G(s) = K/(s(s + 5))单位负反馈
K =
38.方框图困难
R(s)+K控制器+D(s)+10s + 3被控对象C(s)−

比例控制器 K 在单位负反馈回路中驱动被控对象 G(s) = 10/(s + 3)。在 R = 0 时,大小为 1 的扰动阶跃在被控对象输入处进入。求使稳态输出不超过 0.3333 的最小增益 K。

G(s) = 10/(s + 3)D = 1css ≤ 0.3333
K =
39.阶跃响应中等
t (s)y(t)1

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

G(s) = 29/(s² + 2s + 29)
σ =
40.阶跃响应困难
R(s)+135s(s + 10)G(s)C(s)−

单位负反馈回路的 G(s) = 135/(s(s + 10))。用 ts ≈ 4/(ζωₙ),求闭环阶跃响应的 2 % 调节时间 ts。

G(s) = 135/(s(s + 10))
ts = s
41.稳态误差与 PID 控制中等
r(t) = 5 + 5t+4s(s + 4)(s + 5)C(s)−

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

G(s) = 4/(s(s + 4)(s + 5))r(t) = 5 + 5t
ess =
42.稳态误差与 PID 控制困难
R(s)+Kp​ + Ki​/s10s + 5C(s)−

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

G(s) = 10/(s + 5)ωₙ = 4 rad/sζ = 0.7
Kp =

3溶液

1.梁的支座反力

答案: RB = 56.25 kN

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

P1​ = 50 kNP2​ = 50 kNABRA​RB​3 m3 m2 m

受力图

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

平衡

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

ΣMA = 0: RB·8 = 50·3 + 50·6
RB = 450/8 = 56.25 kN
ΣFy = 0: RA = 100 − 56.25 = 43.75 kN

答案

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

答案: RA = 95.86 kN

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

P 的分量

Pv = 80·sin 30° = 40 kN
Ph = 80·cos 30° = 69.28 kN

反力

ΣMA = 0: RB·4 = 40·0.5 + 50·1.5 → RB = 23.75 kN
Ay = 40 + 50 − 23.75 = 66.25 kN
Ax = Ph = 69.28 kN

A 处的合力

RA = √(69.28² + 66.25²) = 95.86 kN

答案

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

答案: ε = 3095 µε

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

应力

A = π·12²/4 = 113.1 mm²
σ = P/A = 35·1000/113.1 = 309.5 MPa

应变(胡克定律)

E for 黄铜 = 100 GPa = 100000 MPa.

ε = σ/E = 309.5/100000 = 0.003095 = 3095 µε

答案

ε = 3095 µε
4.轴向应力与应变

答案: Pmax = 633.8 kN

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

面积

Ar = π·40²/4 = 1257 mm²
At = π(120² − 60²)/4 = 8482 mm²

组合轴向刚度

ΣEA = 200·10³·1257 + 70·10³·8482 = 8.451 × 10⁸ N

控制应变

ε = min(0.00075, 0.001) = 0.00075 (钢)

最大载荷

Pmax = 0.00075·8.451 × 10⁸ = 633816 N = 633.8 kN

答案

Pmax = 633.8 kN
5.轴的扭转

答案: d = 40.67 mm

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

所需的极惯性矩

T = 750 N·m = 750000 N·mm
L = 1.5 m = 1500 mm
G = 80 GPa = 80000 MPa
φ = 3·π/180 = 0.05236 rad
J = TL/(Gφ) = 750000·1500/(80000·0.05236) = 268574 mm⁴

直径

d = (32J/π)1/4 = 40.67 mm

答案

d = 40.67 mm

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

6.轴的扭转

答案: d = 31.56 mm

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

扭矩

ω = 2πN/60 = 2π·1450/60 = 151.8 rad/s
T = P/ω = 75000/151.8 = 493.9 N·m

由许用应力求直径

d = ∛(16T/(π·τallow)) = ∛(16·493929/(π·80)) = 31.56 mm

答案

d = 31.56 mm

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

7.剪力和弯矩

答案: D (x = 3 m)

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

支座反力

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

R1 = 71.25 kN, R2 = 58.75 kN

每个标记点处的弯矩

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

MA = 0, MB = 63.75, MC = 112.5, MD = 146.3 kN·m, ME = 110, MF = 58.75, MG = 0

答案

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

答案: B (x = 2 m)

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

支座反力

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

R1 = 102.3 kN, R2 = 24.75 kN

每个标记点处的弯矩

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

MA = 0, MB = −100 kN·m, MC = −3.5, MD = 7, ME = 1.5, MF = −20, MG = 0

答案

标记点 B: MB = −100 kN·m
9.莫尔圆与应力变换

答案: θp = −67.5 °

平面应力变换:在转过 θ 的截面上,σ = (σₓ + σᵧ)/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 − 70) = −135°

单元体上的角度

θp = −67.5°

答案

θp = −67.5°

逆时针方向的角度为正。

10.莫尔圆与应力变换

答案: σθ = 32.41 MPa

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

变换公式

σθ = 90 + −30·0.6428 + −50·0.766 = 32.41 MPa

答案

σθ = 32.41 MPa
11.梁的挠度

答案: δ = 6.661 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, w = 10 N/mm, L = 1500 mm
EI = 3.99 × 10¹² N·mm²

各荷载单独作用

δP = PL³/(3EI) = 18000·1500³/(3·3.99 × 10¹²) = 5.075 mm
δw = wL⁴/(8EI) = 10·1500⁴/(8·3.99 × 10¹²) = 1.586 mm

相加

δ = 5.075 + 1.586 = 6.661 mm

答案

δ = 6.661 mm
12.梁的挠度

答案: RB = 18 kN

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

去掉支撑

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

δload = wL⁴/(8EI)

支撑处的变形协调

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

RB·L³/(3EI) = wL⁴/(8EI) → RB = 3wL/8 = 3·6·8/8 = 18 kN

答案

RB = 18 kN
13.柱的屈曲

答案: Pcr = 58.95 kN

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

最小惯性矩

矩形绕弱轴屈曲,使用较小边长的三次方。

Imin = h·b³/12 = 100·40³/12 = 533333 mm⁴

计算长度

KL = 1·2500 = 2500 mm

临界荷载

Pcr = π²EI/(KL)² = π²·70000·533333/2500² = 58954 N = 58.95 kN

答案

Pcr = 58.95 kN
14.柱的屈曲

答案: σcr = 17.27 MPa

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

最小惯性矩

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

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

回转半径和长细比

A = 314.2 mm²
r = √(I/A) = √(7854/314.2) = 5 mm
KL = 1·1000 = 1000 mm
KL/r = 1000/5 = 200

临界应力

σcr = π²E/(KL/r)² = π²·70000/200² = 17.27 MPa

答案

σcr = 17.27 MPa
15.抛体运动

答案: t = 3.466 s

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

分解初速度

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

vx = v0·cos θ = 34·cos 30° = 29.44 m/s
vy = v0·sin θ = 34·sin 30° = 17 m/s

飞行时间

起点与落点等高:t = 2vy/g。

t = 2·17/9.81 = 3.466 s

答案

t = 3.466 s
16.抛体运动

答案: x = 78.95 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 = 18.35 m/s, vy = 26.21 m/s

落到屋顶的时间

较大的根:球正在下降。

22 = 26.21·t − 4.905·t²
t = (26.21 + √(26.21² − 2·9.81·22))/9.81 = 4.301 s

水平距离

x = vx·t = 18.35·4.301 = 78.95 m

答案

x = 78.95 m
17.功和能

答案: v = 6.723 m/s

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

沿斜面的能量衡算

物块下降高度 d·sin θ,摩擦消耗 μmg·cos θ·d。质量约去。

½v² = g·d·(sin θ − μ·cos θ)
v = √(2·9.81·15·(0.5 − 0.4·0.866)) = 6.723 m/s

答案

v = 6.723 m/s
18.功和能

答案: x = 11.29 cm

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

通过粗糙段后剩余的能量

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

mgh − μmgL = 19.62·(2.5 − 0.3·4) = 25.51 J

最大压缩量

½kx² = 25.51 → x = √(2·25.51/4000) = 0.1129 m = 11.29 cm

答案

x = 11.29 cm
19.机械振动

答案: ζ = 0.03012

质量–弹簧系统的自由振动:无阻尼固有角频率为 ω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/3)·ln(30/17) = 0.1893

阻尼比

ζ = δ/√(4π² + δ²) = 0.1893/√(39.48 + 0.03585) = 0.03012

答案

ζ = 0.03012
20.机械振动

答案: fn = 5.185 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 = 80 + 12 = 92 kN/m
keq = 92·12/(92 + 12) = 10.62 kN/m

固有频率

ωn = √(10615/10) = 32.58 rad/s
fn = 32.58/(2π) = 5.185 Hz

答案

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

答案: ηmax = 41.15 %

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

绝对温度

TH = 225 + 273.15 = 498.1 K
TL = 20 + 273.15 = 293.1 K

卡诺效率

ηmax = 1 − TL/TH = 1 − 293.1/498.1 = 41.15 %

答案

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

答案: Q̇L = 879.9 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 = 600 + 273.15 = 873.1 K
T0 = 27 + 273.15 = 300.1 K
TL = -5 + 273.15 = 268.1 K

热机的功

η = 1 − 300.1/873.1 = 0.6562
Ẇ = 0.6562·160 = 105 kW

制冷机排出的热量

COPR = 268.1/(300.1 − 268.1) = 8.38
Q̇L = 8.38·105 = 879.9 kW

答案

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

答案: η = 57.98 %

狄塞尔循环(压燃式发动机):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.3385·1.241 = 57.98%

答案

η = 57.98%
24.气体动力循环

答案: η = 53.87 %

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

答案

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

答案: p2 = 67.05 kPa

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

连续性和伯努利

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

v2 = v1·(d1/d2)² = 2.75·(75/40)² = 9.668 m/s
½ρ(v1² − v2²) = ½·1000·(2.75² − 9.668²)/1000 = −42.95 kPa
p2 = p1 + ½ρ(v1² − v2²) = 110 + (−42.95) = 67.05 kPa

答案

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

答案: Q = 10.41 L/s

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

软管中的流速

液面和出口都通大气,所以压强项相互抵消。

v = √(2g·3.5) = √(2·9.81·3.5) = 8.287 m/s

流量

Q = A·v = 0.001257·8.287 = 0.01041 m³/s = 10.41 L/s

答案

Q = 10.41 L/s
27.管道摩擦与水头损失

答案: hf = 4.893 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) = 1.6²/(2·9.81) = 0.1305 m

达西–魏斯巴赫公式

hf = f·(L/D)·V²/(2g) = 0.015·(500/0.2)·1.6²/(2·9.81) = 4.893 m

答案

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

答案: Δz = 1.535 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 = Q/A = 0.075/0.04909 = 1.528 m/s
V²/(2g) = 0.119 m

损失系数

f·L/D = 0.022·100/0.25 = 8.8
ΣK = 0.5 + 2·0.3 + 2 + 1.0 = 4.1

总水头损失

hf = 8.8·0.119 = 1.047 m
hm = 4.1·0.119 = 0.4878 m
Δz = hf + hm = 1.535 m

答案

Δz = 1.535 m
29.水泵与泵功率

答案: η = 85.3 %

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

水功率

Ph = ρ·g·Q·H = 1000·9.81·0.025·16 = 3.924 kW

效率

η = Ph/Pshaft = 3.924/4.6 = 85.3%

答案

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

答案: H = 16.9 m

水泵功率:传给水的功率是水功率 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) = (30 + 6)·0.3305 = 11.9 m

泵扬程

两个水面都是敞口的且基本静止,所以压强项和速度项消去。

H = Δz + hL = 5 + 11.9 = 16.9 m

答案

H = 16.9 m
31.热传导

答案: Q̇ = 135.6 W

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

各层热阻

R1 = L1/(k1A) = (50/1000)/(0.04·8) = 0.1563 K/W
R2 = L2/(k2A) = (150/1000)/(1.4·8) = 0.01339 K/W

热流量

Q̇ = (T1 − T3)/(R1 + R2) = 23/0.1696 = 135.6 W

答案

Q̇ = 135.6 W
32.热传导

答案: Q̇ = 81.12 W

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

串联热阻

Ri = 1/(hiA) = 1/(8·15) = 0.008333 K/W
R1 = L1/(k1A) = (250/1000)/(0.72·15) = 0.02315 K/W
R2 = L2/(k2A) = (150/1000)/(0.035·15) = 0.2857 K/W
Ro = 1/(hoA) = 1/(20·15) = 0.003333 K/W
ΣR = 0.3205 K/W

热损失

Q̇ = (21 − (-5))/0.3205 = 81.12 W

答案

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

答案: N2 = 105

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

总齿数

N1 + N2 = 2a/m = 2·337.5/5 = 135

按传动比分配

N1 = 135/(1 + 3.5) = 30
N2 = 3.5·30 = 105

答案

N2 = 105
34.齿轮与齿轮系

答案: ND = 20

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 = 1500·24/25 = 1440 rpm

齿轮 D 的齿数

ND = C·nC/nout = 15·1440/1080 = 20

答案

ND = 20
35.拉普拉斯变换

答案: f(0.8) = 9.076

拉普拉斯变换: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² + 2s + 1 = (s + 1)²

部分分式

5s + 24 = 5(s + 1) + 19
F(s) = 5/(s + 1) + 19/(s + 1)²

逆变换

f(t) = 5·e−t + 19·t·e−t
f(0.8) = (5 + 19·0.8)·e−0.8 = 9.076

答案

f(0.8) = 9.076
36.拉普拉斯变换

答案: g(0.8) = 0.5602

拉普拉斯变换:F(s) = ∫₀∞ f(t)·e−st dt 把导数变成乘以 s:ℒ{f′} = sF(s) − f(0),ℒ{f″} = s²F(s) − s·f(0) − f′(0)。常用变换对:1 ↔ 1/s,e−at ↔ 1/(s + a),t·e−at ↔ 1/(s + a)²,sin ωt ↔ ω/(s² + ω²),cos ωt ↔ s/(s² + ω²);平移 s → s + σ 相当于把时间函数乘以 e−σt。要对有理函数 F(s) 求逆变换,先分解成部分分式:单极点 p 的留数为 A = [(s − p)·F(s)] 在 s = p 处的值。终值定理:f(∞) = lims→0 s·F(s),仅当 s·F(s) 的所有极点都在左半平面时成立。初值定理:f(0⁺) = lims→∞ s·F(s)。

传递函数

G(s) = Y(s)/U(s) = (4s + 3)/(s² + 3s + 2)

部分分式

s² + 3s + 2 = (s + 1)(s + 2)
A = −1, B = 5

冲激响应

g(t) = −e−t + 5·e−2t
g(0.8) = 0.5602

答案

g(0.8) = 0.5602
37.方框图

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

特征方程

1 + K/(s(s + 5)) = 0 → s² + 5s + K = 0

重根

二次式 s² + βs + γ 在 β² = 4γ 时有二重根。

5² = 4K → K = 25/4 = 6.25

答案

K = 6.25
此时两个闭环极点都在 s = −2.5。
38.方框图

答案: K = 2.7

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

扰动响应

css = D·G(0)/(1 + K·G(0)) = 1·10/(3 + 10K)

解出 K

3 + 10K = 1·10/0.3333 = 30
K = (30 − 3)/10 = 2.7

答案

K = 2.7
39.阶跃响应

答案: σ = −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 + ωₙ² 对应。

ωₙ² = 29 → ωₙ = 5.385 rad/s
2ζωₙ = 2 → ζ = 2/(2·5.385) = 0.1857

结果

σ = −ζωₙ = −2/2 = −1

答案

σ = −1
40.阶跃响应

答案: ts = 0.8 s

阶跃响应:一阶系统 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) = 135/(s² + 10s + 135): ωₙ = 11.62 rad/s, ζ = 0.4303

调节时间

ζωₙ = 10/2 = 5
ts ≈ 4/(ζωₙ) = 0.8 s

答案

ts = 0.8 s

对这个回路,无论 K 取何值,ζωₙ = a/2:增大 K 会增加超调量,但不会缩短调节时间。

41.稳态误差与 PID 控制

答案: ess = 25

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

阶跃部分

Kp = ∞ → estep = 0

斜坡误差

ess = 5/0.2 = 25

答案

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

答案: Kp = 0.06

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

s² + (5 + 10Kp)s + 10Ki

目标多项式

s² + 2·0.7·4s + 4² = s² + 5.6s + 16

匹配各项系数

5 + 10Kp = 5.6 → Kp = 0.06

答案

Kp = 0.06

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

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

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