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​ = 25 kNP2​ = 45 kNABRA​RB​ = ?1 m0.5 m4.5 m

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

L = 6 mP1 = 25 kN 位于 1 mP2 = 45 kN 位于 1.5 m
RB = kN
2.梁的支座反力困难
w0​ = 15 kN/mP = 70 kNABMA​ = ?4 m6 m

悬臂梁在 A 处固定在墙上。梁承受线性变化的荷载,在墙处为零,到自由端 B 增大到 15 kN/m;另在距墙 4 m 处承受集中荷载 P。 求墙处反力矩 MA 的大小。

L = 10 mw0 = 15 kN/mP = 70 kN 位于 4 m
MA = kN·m
3.轴向应力与应变中等
Ø 30 mm铝0.5 mP = 15 kN

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

P = 15 kNd = 30 mmL = 0.5 mE = 70 GPa
δ = mm
4.轴向应力与应变困难
Ø 40 mm黄铜1.5 mØ 20 mm黄铜1.5 mP = 75 kN

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

P = 75 kN第 1 段:Ø40 mm,1.5 m,黄铜(E = 100 GPa)第 2 段:Ø20 mm,1.5 m,黄铜(E = 100 GPa)
δ = mm
5.轴的扭转中等
do​ = ?T = 2000 N·m

一根空心轴须承受 2000 N·m 的扭矩。 其内径为外径的 0.8 倍,切应力不得超过 80 MPa。求最小外径 do。

T = 2000 N·mdi/do = 0.8τallow = 80 MPa
do = mm
6.轴的扭转困难
do​ = ?75 kW, 1450 rpm

一根空心轴须在 1450 rpm 下传递 75 kW。 其内径为外径的 0.75 倍,切应力不得超过 60 MPa。求最小外径 do。

P = 75 kWN = 1450 rpmdi/do = 0.75τallow = 60 MPa
do = mm
7.剪力和弯矩中等
4 kN/m50 kNABCDEFG2 m2 m2 m2 m2 m2 m

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

50 kN @ 8 m4 kN/m, 4–8 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 = 8 kN/mP = 20 kNAB2.5 m2.5 m

简支梁全跨承受均布荷载 w,并承受一个集中荷载 P。求梁中的最大弯矩 Mmax。

L = 5 mP = 20 kN 位于 2.5 mw = 8 kN/m,从 0 到 5 m
Mmax = kN·m
9.莫尔圆与应力变换中等
σx​ = −60σy​ = −20τxy​ = 30MPa

薄板中一点处于平面应力状态:σₓ = −60 MPa,σᵧ = −20 MPa,τₓᵧ = 30 MPa,σz = 0。求该点的绝对最大切应力 τabs(考虑三维中所有平面)。

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

薄板中一点处于平面应力状态:σₓ = 30 MPa,σᵧ = −50 MPa,τₓᵧ = 40 MPa(σz = 0)。 材料的屈服强度为 420 MPa。按冯·米塞斯(畸变能)准则求防止屈服的安全系数 n。

σx = 30 MPaσy = −50 MPaτxy = 40 MPaσY = 420 MPa
n =
11.梁的挠度中等
P = 12 kNABδ0.5 m2 m

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

P = 12 kNa = 0.5 mL = 2.5 mE = 70 GPa (铝)I = 23 × 10⁶ mm⁴
δB = mm
12.梁的挠度困难
w = 10 kN/mP = 25 kNδL = 3 m

简支梁承受均布荷载 w,并在跨中受集中力 P。用叠加法求跨中挠度 δmax。

P = 25 kNw = 10 kN/mL = 3 mE = 70 GPa (铝)I = 84 × 10⁶ mm⁴
δmax = mm
13.柱的屈曲中等
PL = 6 m30 × 80 mm横截面

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

L = 6 mE = 70 GPa (铝)b × h = 30 × 80 mm
Pcr = kN
14.柱的屈曲困难
P = 10 kNL = 2 md = ?横截面

一根长为 L 的铝实心圆柱,底部固定、顶部自由。它需承受 10 kN,屈曲安全系数为 2。求最小直径 d。

P = 10 kNn = 2L = 2 mE = 70 GPa (铝)
d = mm
15.抛体运动中等
v0​ = ?35°R = 45 m

一个球从水平地面以与水平方向成 35° 的角度抛出,落在 45 m 远处。忽略空气阻力,取 g = 9.81 m/s²。求抛出速度 v0。

θ = 35°R = 45 mg = 9.81 m/s²
v0 = m/s
16.抛体运动困难
11 mv0​ = 36 m/s50°x = ?

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

v₀ = 36 m/sθ = 50°hroof = 11 mg = 9.81 m/s²
x = m
17.功和能中等
72 km/h2°2000 kg

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

m = 2000 kgθ = 2°Fr = 500 Nv = 72 km/hg = 9.81 m/s²
P = kW
18.功和能困难
2 kg2 kgμ = 0.15d = 1.5 m

一个 2 kg 的物块放在粗糙的水平桌面上(μ = 0.15)。一根轻绳从它出发,跨过桌边的光滑滑轮,连接一个悬挂的 2 kg 物块。系统由静止释放。求悬挂物块下降 1.5 m 后两物块的速度 v。

m1 = 2 kgm2 = 2 kgμ = 0.15d = 1.5 mg = 9.81 m/s²
v = m/s
19.机械振动中等
tx10 mm3 个周期2.2 mm

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

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

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

m = 50 kgk1 = 30 kN/mk2 = 80 kN/m
fn = Hz
21.热机和制冷机中等
TH​ = 40 °CTL​ = 5 °CRQ̇H​ = ?Q̇L​ = 1.5 kWẆ

一台可逆制冷机通过排出 1.5 kW 的热量使冷库保持在 5 °C,并向 40 °C 的环境放热。求它放热的速率 Q̇H。

TL = 5 °CTH = 40 °CQ̇L = 1.5 kW
Q̇H = kW
22.热机和制冷机困难
TH​ = 21 °CTL​ = -10 °CHPQ̇H​ = 12 kWQ̇L​Ẇ

一台热泵在室外 -10 °C 时将房屋保持在 21 °C,向房屋提供 12 kW 热量。 求可能的最高性能系数 COPmax。

TH = 21 °CTL = -10 °CQ̇H = 12 kW
COPmax =
23.气体动力循环中等
vP1234qin​qout​

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

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

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

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

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

d1 = 75 mmd2 = 50 mmp1 = 200 kPav1 = 1.5 m/sρ = 1000 kg/m³
p2 = kPa
26.连续性和伯努利困难
pair​ = 50 kPav = ?d = 30 mmh = 2 m

一个封闭水箱中水面上方是压缩空气,表压为 50 kPa。水从水面以下 2 m 处直径 30 mm 的孔口流入大气。忽略损失,求 射流速度 v。 取 g = 9.81 m/s²。

pair = 50 kPa (表压)h = 2 md = 30 mmρ = 1000 kg/m³g = 9.81 m/s²
v = m/s
27.管道摩擦与水头损失中等
V = 1 m/sD = 200 mmL = 200 mf = 0.022

水(ρ = 1000 kg/m³)以 1 m/s 的速度流过一根长 200 m、内径 200 mm 的水平管道。取达西摩擦系数 f = 0.022,g = 9.81 m/s²。求管道的压降 Δp,以 kPa 计。

V = 1 m/sL = 200 mD = 200 mmf = 0.022ρ = 1000 kg/m³
Δp = kPa
28.管道摩擦与水头损失困难
K = 0.2Δz = ?L = 100 m, D = 200 mm, f = 0.0152× 弯头 (K = 0.3)

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

Q = 100 L/sL = 100 mD = 200 mmf = 0.015Kentrance = 0.5Kbend = 0.3 (×2)Kvalve = 0.2Kexit = 1.0
Δz = m
29.水泵与泵功率中等
PQ = 30 L/sH = 28 m

一台效率为 78% 的水泵以 28 m 的扬程输送 30 L/s 的水(ρ = 1000 kg/m³)。取 g = 9.81 m/s²。求所需的轴功率 Pshaft,以 kW 计。

Q = 30 L/sH = 28 mη = 78%g = 9.81 m/s²
Pshaft = kW
30.水泵与泵功率困难
PQ = 80 L/sΔz = 11 mL = 500 m, D = 200 mm, f = 0.02, ΣK = 6

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

Q = 80 L/sΔz = 11 mL = 500 mD = 200 mmf = 0.02ΣK = 6η = 75%
Pshaft = kW
31.热传导中等
砖, k = 0.72250 mm矿棉, k = 0.0475 mm24 °C0 °CTi​ = ?Q̇

一面墙由 250 mm 的砖和 75 mm 的矿棉组成,面积为 10 m²。外表面温度为 24 °C 和 0 °C。 求两层交界面处的温度 Ti。

L1 = 250 mm, k1 = 0.72 W/(m·K)L2 = 75 mm, k2 = 0.04 W/(m·K)A = 10 m²T1 = 24 °CT3 = 0 °C
Ti = °C
32.热传导困难
混凝土, k = 1.4100 mm玻璃, k = 0.96100 mm室内 21 °C室外 1 °CQ̇

由 100 mm 混凝土和 100 mm 玻璃组成的墙,将 21 °C 的房间与 1 °C 的室外空气隔开。对流换热系数室内为 10 W/(m²·K),室外为 30 W/(m²·K),墙面积为 10 m²。 求墙体的总传热系数 U。

L1 = 100 mm, k1 = 1.4 W/(m·K)L2 = 100 mm, k2 = 0.96 W/(m·K)A = 10 m²hi = 10 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = 1 °C
U = W/(m²·K)
33.齿轮与齿轮系中等
N1​ = 30N2​ = 60

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

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

在一个复合轮系中,输入轴上的齿轮 A(32 齿)以 1440 rpm 转动,驱动齿轮 B(48 齿)。齿轮 C(72 齿)与 B 固定在同一根轴上,驱动输出轴上的齿轮 D(80 齿)。求输出转速 nout。

nin = 1440 rpmA = 32B = 48C = 72D = 80
nout = rpm
35.拉普拉斯变换中等
ReIm−5−10−15−20−25

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

F(s) = (s + 25)/(s² + 11s + 24)
A =
36.拉普拉斯变换困难
ReIm−2−4−6−8−10−12−14−16−1824

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

y″ + 4y′ + 20y = u′ + 16ut = 0.25 s
g(0.25) =
37.方框图中等
R(s)+4s + 7G(s)C(s)0.5H−

被控对象 G(s) = 4/(s + 7) 由一个负反馈回路控制,该回路的传感器增益为 H = 0.5。 在 R 处施加大小为 5 的阶跃。求稳态输出 css。

G(s) = 4/(s + 7)H = 0.5R = 5
css =
38.方框图困难
R(s)+30s(s + 3)G(s)C(s)1H−

被控对象 G(s) = 30/(s(s + 3)) 放在传感器增益为 H = 1 的负反馈回路中。闭环极点是复数。 求它们的实部 σ。

G(s) = 30/(s(s + 3))H = 1
σ =
39.阶跃响应中等
t (s)y(t)2ts​ = ?

某系统的传递函数为 G(s) = 92/(s² + 6s + 46),输入为单位阶跃。 用 ts ≈ 4/(ζωₙ),求阶跃响应的 2 % 调节时间 ts。

G(s) = 92/(s² + 6s + 46)
ts = s
40.阶跃响应困难
t (s)y(t)112 %

一个没有零点的二阶系统,其阶跃响应必须有 12 % 的超调量。求阻尼比 ζ。

%OS = 12 %
ζ =
41.稳态误差与 PID 控制中等
r(t) = 1 + 0.5t+20s(s + 2)(s + 4)C(s)−

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

G(s) = 20/(s(s + 2)(s + 4))r(t) = 1 + 0.5t
ess =
42.稳态误差与 PID 控制困难
R(s)+K = 2510s(s + 8)C(s)−

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

G(s) = 10/(s(s + 8))K = 25
tp = s

3溶液

1.梁的支座反力

答案: RB = 15.42 kN

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

P1​ = 25 kNP2​ = 45 kNABRA​RB​1 m0.5 m4.5 m

受力图

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

平衡

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

ΣMA = 0: RB·6 = 25·1 + 45·1.5
RB = 92.5/6 = 15.42 kN
ΣFy = 0: RA = 70 − 15.42 = 54.58 kN

答案

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

答案: MA = 780 kN·m

悬臂梁的平衡:固定端既能向上支承,也能阻止转动,所以它提供竖向反力 RA 和反力矩 MA。由 ΣFy = 0 得 RA,对墙取矩得 MA = Σ(荷载 × 到墙的距离)。分布荷载以其合力作用在荷载图的形心处。

三角形荷载的合力

W = ½·15·10 = 75 kN
x̄ = 2L/3 = 6.667 m

对墙的力矩

MA = W·x̄ + P·a = 75·6.667 + 70·4 = 780 kN·m

答案

MA = 780 kN·m
3.轴向应力与应变

答案: δ = 0.1516 mm

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

应力

A = π·30²/4 = 706.9 mm²
σ = P/A = 15·1000/706.9 = 21.22 MPa

应变(胡克定律)

E for 铝 = 70 GPa = 70000 MPa.

ε = σ/E = 21.22/70000 = 0.0003032 = 303.2 µε

伸长量

δ = ε·L = 0.0003032·500 mm = 0.1516 mm

答案

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

答案: δ = 4.476 mm

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

两段受力相同

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

A1 = π·40²/4 = 1257 mm², A2 = π·20²/4 = 314.2 mm²
σ1 = 59.68 MPa
σ2 = 238.7 MPa

叠加各段伸长量

δ1 = σ1L1/E1 = 59.68·1500/100000 = 0.8952 mm
δ2 = σ2L2/E2 = 238.7·1500/100000 = 3.581 mm
δ = δ1 + δ2 = 4.476 mm

答案

δ = 4.476 mm
5.轴的扭转

答案: do = 59.97 mm

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

由容许应力求外径

T = 2000 N·m = 2 × 10⁶ N·mm
1 − k⁴ = 1 − 0.8⁴ = 0.5904
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·2 × 10⁶/(π·80·0.5904)) = 59.97 mm
di = 0.8·59.97 = 47.97 mm

答案

do = 59.97 mm

实心轴需要 d = 50.31 mm;空心轴稍大一些,但轻得多。

6.轴的扭转

答案: do = 39.44 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

由容许应力求外径

T = 493.9 N·m = 493929 N·mm
1 − k⁴ = 1 − 0.75⁴ = 0.6836
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·493929/(π·60·0.6836)) = 39.44 mm
di = 0.75·39.44 = 29.58 mm

答案

do = 39.44 mm

实心轴需要 d = 34.74 mm;空心轴稍大一些,但轻得多。

7.剪力和弯矩

答案: E (x = 8 m)

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

支座反力

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

R1 = 24.67 kN, R2 = 41.33 kN

每个标记点处的弯矩

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

MA = 0, MB = 49.33, MC = 98.67, MD = 140, ME = 165.3 kN·m, MF = 82.67, MG = 0

答案

标记点 E: ME = 165.3 kN·m
8.剪力和弯矩

答案: Mmax = 50 kN·m

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

支座反力

整根梁的平衡。

均布荷载合力:W = w·ℓ = 8·5 = 40 kN,作用于 x̄ = 2.5 m
ΣMA = 0: RB·5 = 20·2.5 + 8·5·2.5
RB = 150/5 = 30 kN
ΣFy = 0: RA = 60 − 30 = 30 kN

V = 0 在哪里?

V(2.5⁻) = RA − w·2.5 = 10 kN
V(2.5⁺) = 10 − 20 = −10 kN
剪力在荷载处变号,所以最大值在荷载处。

最大弯矩

Mmax = RA·2.5 − w·2.5²/2 = 30·2.5 − 8·2.5²/2 = 50 kN·m

答案

Mmax = 50 kN·m

它出现在距 A 为 x = 2.5 m 处。

9.莫尔圆与应力变换

答案: τabs = 38.03 MPa

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

στxyCσ1​σ2​

面内主应力

C = −40 MPa, R = √((−20)² + 30²) = 36.06 MPa
σ1 = −3.944 MPa, σ2 = −76.06 MPa, σz = 0

三个圆中最大的

τ = max(|σ1 − σ2|, |σ1|, |σ2|)/2 = max(72.11, 3.944, 76.06)/2 = 38.03 MPa

答案

τabs = 38.03 MPa

这里由面外的圆控制:绝对最大值大于面内最大值 R。

10.莫尔圆与应力变换

答案: n = 4.264

平面应力变换:在转过 θ 的截面上,σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。莫尔圆可同时表示所有截面:圆心位于 C = (σₓ + σᵧ)/2,半径为 R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主应力为 C ± R,面内最大切应力为 R。 屈服准则:冯·米塞斯准则预测当 σvM = √(σ₁² − σ₁σ₂ + σ₂²) 达到 σY 时屈服;特雷斯卡准则预测当绝对最大切应力达到 σY/2 时屈服。

στxyCσ1​σ2​

主应力

C = −10 MPa, R = 56.57 MPa
σ1 = 46.57 MPa, σ2 = −66.57 MPa, σ3 = 0

冯·米塞斯

σvM = √(46.57² − 46.57·(−66.57) + (−66.57)²) = 98.49 MPa
n = 420/98.49 = 4.264

答案

n = 4.264

n > 1 表示该点不屈服;特雷斯卡准则比冯·米塞斯准则稍保守。

11.梁的挠度

答案: δB = 2.174 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 = 12000 N, a = 500 mm, L = 2500 mm
EI = 1.61 × 10¹² N·mm²

自由端的挠度

荷载处的挠度加上荷载处的转角乘以剩余长度。

δP = Pa³/(3EI) = 0.3106 mm
θP = Pa²/(2EI) = 0.0009317 rad
δB = δP + θP·(L − a) = Pa²(3L − a)/(6EI) = 2.174 mm

答案

δB = 2.174 mm
12.梁的挠度

答案: δmax = 4.185 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 = 25000 N, w = 10 N/mm, L = 3000 mm, EI = 5.88 × 10¹² N·mm²

各荷载单独作用

δP = PL³/(48EI) = 25000·3000³/(48·5.88 × 10¹²) = 2.392 mm
δw = 5wL⁴/(384EI) = 5·10·3000⁴/(384·5.88 × 10¹²) = 1.794 mm

相加

δmax = 2.392 + 1.794 = 4.185 mm

答案

δmax = 4.185 mm
13.柱的屈曲

答案: Pcr = 3.454 kN

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

最小惯性矩

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

Imin = h·b³/12 = 80·30³/12 = 180000 mm⁴

计算长度

KL = 1·6000 = 6000 mm

临界荷载

Pcr = π²EI/(KL)² = π²·70000·180000/6000² = 3454 N = 3.454 kN

答案

Pcr = 3.454 kN
14.柱的屈曲

答案: d = 55.42 mm

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

所需临界荷载

Pcr = n·P = 2·10000 = 20000 N
KL = 2·2000 = 4000 mm

求解直径

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·20000·4000²/(π³·70000))^(1/4) = 55.42 mm

答案

d = 55.42 mm

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

15.抛体运动

答案: v0 = 21.67 m/s

Projectile motion: with no air resistance, x = vxt and y = h + vyt − ½gt². The two directions are independent except that they share the same time.

从射程公式求 v₀

v₀ = √(g·R/sin 2θ) = √(9.81·45/sin 70°) = 21.67 m/s

答案

v0 = 21.67 m/s
16.抛体运动

答案: x = 120.1 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 = 23.14 m/s, vy = 27.58 m/s

落到屋顶的时间

较大的根:球正在下降。

11 = 27.58·t − 4.905·t²
t = (27.58 + √(27.58² − 2·9.81·11))/9.81 = 5.19 s

水平距离

x = vx·t = 23.14·5.19 = 120.1 m

答案

x = 120.1 m
17.功和能

答案: P = 23.69 kW

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

驱动力

m·g = 2000·9.81 = 19620 N
mg·sin θ = 19620·sin 2° = 684.7 N
F = 684.7 + 500 = 1185 N

功率

v = 72/3.6 = 20 m/s
P = F·v = 1185·20 = 23695 W = 23.69 kW

答案

P = 23.69 kW
18.功和能

答案: v = 3.537 m/s

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

释放和损失的能量

m₂g·d = 19.62·1.5 = 29.43 J
μ·m₁g·d = 0.15·19.62·1.5 = 4.415 J

两个物块的动能

½(m₁ + m₂)v² = 25.02 → v = √(2·25.02/4) = 3.537 m/s

答案

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

答案: ζ = 0.08007

质量–弹簧系统的自由振动:无阻尼固有角频率为 ω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(10/2.2) = 0.5047

阻尼比

ζ = δ/√(4π² + δ²) = 0.5047/√(39.48 + 0.2547) = 0.08007

答案

ζ = 0.08007
20.机械振动

答案: fn = 7.465 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 = 30 + 80 = 110 kN/m

固有频率

ωn = √(110000/50) = 46.9 rad/s
fn = 46.9/(2π) = 7.465 Hz

答案

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

答案: Q̇H = 1.689 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 = 40 + 273.15 = 313.1 K
TL = 5 + 273.15 = 278.1 K

可逆循环的热量比

Q̇H = Q̇L·TH/TL = 1.5·313.1/278.1 = 1.689 kW

答案

Q̇H = 1.689 kW

输入功率:Ẇ = 0.1887 kW。

22.热机和制冷机

答案: COPmax = 9.489

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

绝对温度

TH = 21 + 273.15 = 294.1 K
TL = -10 + 273.15 = 263.1 K

可能的最佳 COP

COPHP,max = TH/(TH − TL) = 294.1/(294.1 − 263.1) = 9.489

答案

COPmax = 9.489
23.气体动力循环

答案: qin = 1329 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 = 320·220.4 = 1102 K
T3 = 2.2·1102 = 2424 K

定压吸热量

qin = 1.005·(2424 − 1102) = 1329 kJ/kg

答案

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

答案: wnet = 139.4 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 = 295·100.2857 = 569.6 K
T4s = 1100/100.2857 = 569.7 K

压气机功和涡轮功

wc = 1.005·(569.6 − 295)/0.88 = 313.6 kJ/kg
wt = 0.85·1.005·(1100 − 569.7) = 453 kJ/kg

净功

wnet = 453 − 313.6 = 139.4 kJ/kg

答案

wnet = 139.4 kJ/kg

回功比:wc/wt = 69.22%。

25.连续性和伯努利

答案: p2 = 195.4 kPa

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

连续性和伯努利

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

v2 = v1·(d1/d2)² = 1.5·(75/50)² = 3.375 m/s
½ρ(v1² − v2²) = ½·1000·(1.5² − 3.375²)/1000 = −4.57 kPa
p2 = p1 + ½ρ(v1² − v2²) = 200 + (−4.57) = 195.4 kPa

答案

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

答案: v = 11.8 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·2)) = 11.8 m/s

答案

v = 11.8 m/s
27.管道摩擦与水头损失

答案: Δp = 11 kPa

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

水头损失

hf = f·(L/D)·V²/(2g) = 0.022·(200/0.2)·1²/(2·9.81) = 1.121 m

压降

管道水平且流速不变,所以能量方程简化为 Δp = ρ·g·hf。

Δp = ρ·g·hf = 1000·9.81·1.121 = 11000 Pa = 11 kPa

答案

Δp = 11 kPa
28.管道摩擦与水头损失

答案: Δz = 5.061 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.1/0.03142 = 3.183 m/s
V²/(2g) = 0.5164 m

损失系数

f·L/D = 0.015·100/0.2 = 7.5
ΣK = 0.5 + 2·0.3 + 0.2 + 1.0 = 2.3

总水头损失

hf = 7.5·0.5164 = 3.873 m
hm = 2.3·0.5164 = 1.188 m
Δz = hf + hm = 5.061 m

答案

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

答案: Pshaft = 10.56 kW

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

水功率

Ph = ρ·g·Q·H = 1000·9.81·0.03·28 = 8.24 kW

轴功率

Pshaft = Ph/η = 8.24/0.78 = 10.56 kW

答案

Pshaft = 10.56 kW
30.水泵与泵功率

答案: Pshaft = 30.88 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) = (50 + 6)·0.3305 = 18.51 m

泵扬程

H = Δz + hL = 11 + 18.51 = 29.51 m

轴功率

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.08·29.51/0.75 = 30.88 kW

答案

Pshaft = 30.88 kW
31.热传导

答案: Ti = 20.25 °C

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

各层热阻

R1 = L1/(k1A) = (250/1000)/(0.72·10) = 0.03472 K/W
R2 = L2/(k2A) = (75/1000)/(0.04·10) = 0.1875 K/W

热流量

Q̇ = (T1 − T3)/(R1 + R2) = 24/0.2222 = 108 W

交界面温度

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

Ti = T1 − Q̇·R1 = 24 − 108·0.03472 = 20.25 °C

答案

Ti = 20.25 °C
32.热传导

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

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

串联热阻

Ri = 1/(hiA) = 1/(10·10) = 0.01 K/W
R1 = L1/(k1A) = (100/1000)/(1.4·10) = 0.007143 K/W
R2 = L2/(k2A) = (100/1000)/(0.96·10) = 0.01042 K/W
Ro = 1/(hoA) = 1/(30·10) = 0.003333 K/W
ΣR = 0.03089 K/W

总传热系数

U = 1/(A·ΣR) = 1/(10·0.03089) = 3.237 W/(m²·K)

答案

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

答案: T2 = 283.5 N·m

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

电动机扭矩

ω₁ = 2π·960/60 = 100.5 rad/s
T1 = P/ω₁ = 15000/100.5 = 149.2 N·m

输出扭矩

T2 = 149.2·60/30·0.95 = 283.5 N·m

答案

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

答案: nout = 864 rpm

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

第一级

nB = 1440·32/48 = 960 rpm = nC

第二级

nout = 960·72/80 = 864 rpm

答案

nout = 864 rpm

总传动比 = 1.667 : 1

35.拉普拉斯变换

答案: A = 4.4

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

遮盖法

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

A = (s + 25)/(s + 8) |s = −3 = 22/5 = 4.4

答案

A = 4.4
(B = −3.4,所以 f(t) = 4.4·e−3t − 3.4·e−8t)
36.拉普拉斯变换

答案: g(0.25) = 2.114

拉普拉斯变换: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) = (s + 16)/(s² + 4s + 20)

部分分式

s² + 4s + 20 = (s + 2)² + 4²
G(s) = 1·(s + 2)/((s + 2)² + 4²) + (14/4)·4/((s + 2)² + 4²)

冲激响应

g(t) = e−2t·(cos 4t + 3.5·sin 4t)
g(0.25) = 2.114

答案

g(0.25) = 2.114
37.方框图

答案: css = 2.222

方框图代数:串联的方框相乘(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) = 4/(s + 7 + 2) = 4/(s + 9)

结果

css = R·T(0) = 5·4/9 = 2.222

答案

css = 2.222
38.方框图

答案: σ = −1.5

方框图代数:串联的方框相乘(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(s + 3) + 30·1 = 0 → s² + 3s + 30 = 0

根

s = −1.5 ± j·√(30 − 2.25) = −1.5 ± j5.268
σ = −1.5

答案

σ = −1.5
39.阶跃响应

答案: ts = 1.333 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)。

与标准形式比较

分母必须与 s² + 2ζωₙs + ωₙ² 对应。

ωₙ² = 46 → ωₙ = 6.782 rad/s
2ζωₙ = 6 → ζ = 6/(2·6.782) = 0.4423

结果

ζωₙ = 3 → ts ≈ 4/3 = 1.333 s

答案

ts = 1.333 s
40.阶跃响应

答案: ζ = 0.5594

阶跃响应:一阶系统 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)。

阻尼比

ζ = −ln(0.12)/√(π² + ln²(0.12)) = 2.12/√(9.87 + 4.496) = 0.5594

答案

ζ = 0.5594
41.稳态误差与 PID 控制

答案: ess = 0.2

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

阶跃部分

Kp = ∞ → estep = 0

斜坡误差

ess = 0.5/2.5 = 0.2

答案

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

答案: tp = 0.2054 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 + 8) + 25·10 = s² + 8s + 250

ωₙ 和 ζ

ωₙ = √250 = 15.81 rad/s
ζ = 8/(2·15.81) = 0.253

峰值时间

ωd = 15.81·√(1 − 0.253²) = 15.3 rad/s
tp = π/15.3 = 0.2054 s

答案

tp = 0.2054 s

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

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