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​ = 35 kNP2​ = 15 kNABRA​RB​ = ?1 m6 m3 m

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

L = 10 mP1 = 35 kN 位于 1 mP2 = 15 kN 位于 7 m
RB = kN
2.梁的支座反力困难
w = 6 kN/mP = 45 kNABRA​ = ?RB​1 m2.5 m0.5 m1 m

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

L = 5 mP = 45 kN 位于 3.5 mw = 6 kN/m,从 1 到 4 m
RA = kN
3.轴向应力与应变中等
P = ?12400 mm500 mmØ30Ø60δ = 0.3 mm

一根杆由一段钢(直径 30 mm,长 400 mm,E = 200 GPa)与一段黄铜(直径 60 mm,长 500 mm,E = 100 GPa)首尾相接而成。它一端固定,另一端被拉。求使它总共伸长 0.3 mm 的轴向载荷 P。

第 1 段: 钢, Ø30 mm, 400 mm, E = 200 GPa第 2 段: 黄铜, Ø60 mm, 500 mm, E = 100 GPaδ = 0.3 mm
P = kN
4.轴向应力与应变困难
Ø 40 mm钢2 mØ 16 mm钢0.5 mP = 115 kN

阶梯杆一端固定在墙上,自由端受荷载。求杆中的最大正应力 σmax。

P = 115 kN第 1 段:Ø40 mm,2 m,钢(E = 200 GPa)第 2 段:Ø16 mm,0.5 m,钢(E = 200 GPa)
σmax = MPa
5.轴的扭转中等
T = 250 N·md = 25 mmL = 1 mφ = ?

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

T = 250 N·md = 25 mmL = 1 mG = 39 GPa
φ = °
6.轴的扭转困难
do​ = ?10 kW, 3000 rpm

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

P = 10 kWN = 3000 rpmdi/do = 0.7τallow = 60 MPa
do = mm
7.剪力和弯矩中等
50 kN40 kNABCDEFG2 m2 m2 m2 m2 m2 m

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

50 kN @ 10 m40 kN @ 6 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 = 12 kN/mAB4 m6 m

简支梁在部分跨度上承受均布荷载。 求梁中的最大弯矩 Mmax。

L = 10 mw = 12 kN/m,从 0 到 4 m
Mmax = kN·m
9.莫尔圆与应力变换中等
σx​ = 160σy​ = 40τxy​ = −20MPa

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

σx = 160 MPaσy = 40 MPaτxy = −20 MPaσz = 0
τabs = MPa
10.莫尔圆与应力变换困难
σx​τxy​轴的表面

一根直径 60 mm 的实心轴承受100 kN 的轴向压力,同时承受 1500 N·m 的扭矩。 在轴的表面处, 求面内最大切应力 τmax。

d = 60 mmF = −100 kNT = 1500 N·m
τmax = MPa
11.梁的挠度中等
w = 8 kN/mδL = 2.5 m

悬臂梁全长承受均布荷载 w。求自由端挠度 δ。

w = 8 kN/mL = 2.5 mE = 70 GPa (铝)I = 23 × 10⁶ mm⁴
δ = mm
12.梁的挠度困难
w = 8 kN/mP = 19 kNδL = 2.5 m

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

P = 19 kNw = 8 kN/mL = 2.5 mE = 11 GPa (结构木材)I = 160 × 10⁶ mm⁴
δmax = mm
13.柱的屈曲中等
PL = 3 md = 60 mm横截面

一根长为 L 的钢柱,底部固定、顶部铰接。用欧拉公式求其临界屈曲荷载 Pcr。

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

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

P = 20 kNn = 3L = 1 mE = 70 GPa (铝)
d = mm
15.抛体运动中等
v0​ = 34 m/s70°x = ?

一个球 从平地上以高于水平面 70° 的角度抛出。忽略空气阻力。求球落地前飞过的水平距离 x。

v0 = 34 m/sθ = 70°g = 9.81 m/s²
x = m
16.抛体运动困难
5 mv0​ = 31 m/s60°x = ?

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

v₀ = 31 m/sθ = 60°hroof = 5 mg = 9.81 m/s²
x = m
17.功和能中等
45°d = 10 mμ = 0.15

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

θ = 45°μ = 0.15d = 10 mg = 9.81 m/s²
v = m/s
18.功和能困难
h = 2 mL = 5 mμ = 0.2k = 2000 N/m

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

m = 0.5 kgh = 2 mL = 5 mμ = 0.2k = 2000 N/mg = 9.81 m/s²
x = cm
19.机械振动中等
k = 2 kN/mζ = 0.7m = 12 kg

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

m = 12 kgk = 2 kN/mζ = 0.7
fd = Hz
20.机械振动困难
k1​ = 100k2​ = 25m = 10 kg

一个 10 kg 的质量块悬挂在并排(并联)的两根弹簧上,刚度分别为 k₁ = 100 kN/m 和 k₂ = 25 kN/m。质量块保持水平。 求弹簧系统的等效刚度 keq,以 kN/m 计。

m = 10 kgk1 = 100 kN/mk2 = 25 kN/m
keq = kN/m
21.热机和制冷机中等
TH​ = 750 °CTL​ = 35 °CEQH​ = 3900 kJQL​W = ?

一台热机在 750 °C 的高温热源和 35 °C 的低温热源之间工作。 它吸收 3900 kJ 热量。求它最多能输出的功 Wmax。

TH = 750 °CTL = 35 °CQH = 3900 kJ
Wmax = kJ
22.热机和制冷机困难
TH​ = 900 °CTL​ = 40 °CEQH​ = 5000 kJQL​ = 2050 kJW = 2950 kJ

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

TH = 900 °CTL = 40 °CQH = 5000 kJW = 2950 kJ
Sgen = kJ/K
23.气体动力循环中等
vP1234qin​qout​

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

r = 19rc = 2T1 = 320 Kk = 1.4
T3 = K
24.气体动力循环困难
vP1234qin​qout​

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

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

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

d1 = 100 mmd2 = 75 mmp1 = 160 kPav1 = 0.5 m/sρ = 1000 kg/m³
p2 = kPa
26.连续性和伯努利困难
12d1​ = 100 mm, p1​ = 280 kPad2​ = 75 mm, p2​ = ?Δz = 1 m

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

d1 = 100 mmd2 = 75 mmp1 = 280 kPav1 = 2.25 m/sΔz = 1 mρ = 1000 kg/m³
p2 = kPa
27.管道摩擦与水头损失中等
V = 2.4 m/sD = 300 mmL = 300 mf = 0.018

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

V = 2.4 m/sL = 300 mD = 300 mmf = 0.018ρ = 1000 kg/m³
Δp = kPa
28.管道摩擦与水头损失困难
K = 2自由出流H = ?L = 400 m, D = 300 mm, f = 0.0183× 弯头 (K = 0.3)

一根长 400 m、内径 300 mm(达西摩擦系数 f = 0.018)的管道从敞口水池引出,以自由射流排入大气。局部损失包括一个直角进口(K = 0.5)、3 个弯头(每个 K = 0.3)和一个阀门(K = 2)。取 g = 9.81 m/s²。要使流量为 150 L/s,水池水面须比管道出口高多少 H?

Q = 150 L/sL = 400 mD = 300 mmf = 0.018Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 2
H = m
29.水泵与泵功率中等
PQ = 10 L/sH = 70 m

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

Q = 10 L/sH = 70 mη = 60%g = 9.81 m/s²
Pshaft = kW
30.水泵与泵功率困难
150 kPa (表压)PQ = 20 L/sΔz = 17 mL = 150 m, D = 100 mm, f = 0.02, ΣK = 5

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

Q = 20 L/sΔz = 17 mL = 150 mD = 100 mmf = 0.02ΣK = 5ptank = 150 kPa
H = m
31.热传导中等
石膏板, k = 0.17200 mm木材, k = 0.1750 mm23 °C9 °CTi​ = ?Q̇

一面墙由 200 mm 的石膏板和 50 mm 的木材组成,面积为 12 m²。外表面温度为 23 °C 和 9 °C。 求两层交界面处的温度 Ti。

L1 = 200 mm, k1 = 0.17 W/(m·K)L2 = 50 mm, k2 = 0.17 W/(m·K)A = 12 m²T1 = 23 °CT3 = 9 °C
Ti = °C
32.热传导困难
木材, k = 0.17200 mm矿棉, k = 0.0425 mm室内 22 °C室外 -15 °CTs​ = ?Q̇

由 200 mm 木材和 25 mm 矿棉组成的墙,将 22 °C 的房间与 -15 °C 的室外空气隔开。对流换热系数室内为 6 W/(m²·K),室外为 25 W/(m²·K),墙面积为 4 m²。 求墙内表面的温度 Ts。

L1 = 200 mm, k1 = 0.17 W/(m·K)L2 = 25 mm, k2 = 0.04 W/(m·K)A = 4 m²hi = 6 W/(m²·K)ho = 25 W/(m²·K)Troom = 22 °CTout = -15 °C
Ts = °C
33.齿轮与齿轮系中等
N1​ = 12N2​ = 24

一台 15 kW 的电动机以 2900 rpm 运转,驱动一个齿数为 12 的小齿轮,小齿轮与一个齿数为 24 的大齿轮啮合。 忽略损失。 求电动机轴上的扭矩 T1。

P = 15 kWn1 = 2900 rpmN1 = 12N2 = 24η = 100%
T1 = N·m
34.齿轮与齿轮系困难
A = 36B = 60C = 24D = ?

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

nin = 1440 rpmA = 36B = 60C = 24nout = 576 rpm
ND =
35.拉普拉斯变换中等
t (s)f(t)t = 0.5f(0.5) = ?

求 F(s) = (2s + 34)/(s² + 8s + 16) 对应的 f(t),并计算 f(0.5)。

F(s) = (2s + 34)/(s² + 8s + 16)t = 0.5 s
f(0.5) =
36.拉普拉斯变换困难
t (s)y(t)t = 0.1y(0.1) = ?

用拉普拉斯变换解 y″ + 9y′ + 18y = 0,初始条件 y(0) = 10,y′(0) = 3,并求 y(0.1)。

y″ + 9y′ + 18y = 0y(0) = 10y′(0) = 3t = 0.1 s
y(0.1) =
37.方框图中等
R(s)+31s + 9G(s)C(s)0.25H−

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

G(s) = 31/(s + 9)H = 0.25R = 10
css =
38.方框图困难
R(s)+68s(s + 12)G(s)C(s)1H−

被控对象 G(s) = 68/(s(s + 12)) 放在传感器增益为 H = 1 的负反馈回路中。闭环极点是复数。 求它们的虚部 ωd(取正值),单位 rad/s。

G(s) = 68/(s(s + 12))H = 1
ωd = rad/s
39.阶跃响应中等
t (s)y(t)3?

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

G(s) = 150/(s² + 3s + 50)
%OS = %
40.阶跃响应困难
k = 2000 N/mc = 51 N·s/mm = 8 kgF = 100 Nx(t)

质量 m = 8 kg 的物体挂在弹簧 k = 2000 N/m 上,带有黏性阻尼器 c = 51 N·s/m(m·x″ + c·x′ + k·x = F)。在静止平衡状态下,突然施加恒力 F = 100 N。 求阻尼比 ζ。

m = 8 kgk = 2000 N/mc = 51 N·s/mF = 100 N
ζ =
41.稳态误差与 PID 控制中等
r(t) = 2 + 5t+4s(s + 4)C(s)−

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

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

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

G(s) = 1/(s(s + 4))K = 40
tp = s

3溶液

1.梁的支座反力

答案: RB = 14 kN

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

P1​ = 35 kNP2​ = 15 kNABRA​RB​1 m6 m3 m

受力图

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

平衡

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

ΣMA = 0: RB·10 = 35·1 + 15·7
RB = 140/10 = 14 kN
ΣFy = 0: RA = 50 − 14 = 36 kN

答案

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

答案: RA = 22.5 kN

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

w = 6 kN/mP = 45 kNABRA​RB​1 m2.5 m0.5 m1 m

受力图

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

平衡

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

均布荷载合力:W = w·ℓ = 6·3 = 18 kN,作用于 x̄ = 2.5 m
ΣMA = 0: RB·5 = 45·3.5 + 6·3·2.5
RB = 202.5/5 = 40.5 kN
ΣFy = 0: RA = 63 − 40.5 = 22.5 kN

答案

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

答案: P = 65.25 kN

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

面积

A₁ = 706.9 mm², A₂ = 2827 mm²

柔度

L₁/(E₁A₁) = 400/(200000·706.9) = 2.829 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 500/(100000·2827) = 1.768 × 10⁻⁶ mm/N

载荷

P = 0.3/4.598 × 10⁻⁶ = 65248 N = 65.25 kN

答案

P = 65.25 kN
4.轴向应力与应变

答案: σmax = 572 MPa

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

两段受力相同

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

A1 = π·40²/4 = 1257 mm², A2 = π·16²/4 = 201.1 mm²
σ1 = 91.51 MPa
σ2 = 572 MPa

最大应力

较细的一段面积较小,因此应力较大。

σmax = 572 MPa

答案

σmax = 572 MPa
5.轴的扭转

答案: φ = 9.577 °

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

极惯性矩

J = π·25⁴/32 = 38350 mm⁴

扭转角

φ = TL/(GJ) = 250·1000·1000/(39000·38350) = 0.1672 rad
φ = 0.1672·180/π = 9.577°

答案

φ = 9.577°
6.轴的扭转

答案: do = 15.26 mm

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

扭矩

ω = 2πN/60 = 2π·3000/60 = 314.2 rad/s
T = P/ω = 10000/314.2 = 31.83 N·m

由容许应力求外径

T = 31.83 N·m = 31831 N·mm
1 − k⁴ = 1 − 0.7⁴ = 0.7599
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·31831/(π·60·0.7599)) = 15.26 mm
di = 0.7·15.26 = 10.68 mm

答案

do = 15.26 mm

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

7.剪力和弯矩

答案: D (x = 6 m)

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

支座反力

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

R1 = 28.33 kN, R2 = 61.67 kN

每个标记点处的弯矩

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

MA = 0, MB = 56.67, MC = 113.3, MD = 170 kN·m, ME = 146.7, MF = 123.3, MG = 0

答案

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

答案: Mmax = 61.44 kN·m

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

支座反力

整根梁的平衡。

均布荷载合力:W = w·ℓ = 12·4 = 48 kN,作用于 x̄ = 2 m
ΣMA = 0: RB·10 = 12·4·2
RB = 96/10 = 9.6 kN
ΣFy = 0: RA = 48 − 9.6 = 38.4 kN

剪力为零

剪力从 RA 开始,在无荷载段保持不变,在荷载范围内以斜率 w 下降。

RA − w·(x0 − 0) = 0
x0 = 0 + 38.4/12 = 3.2 m

最大弯矩

Mmax = RA·x0 − w·(x0 − 0)²/2 = 38.4·3.2 − 12·3.2²/2 = 61.44 kN·m

答案

Mmax = 61.44 kN·m
9.莫尔圆与应力变换

答案: τabs = 81.62 MPa

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

στxyCσ1​σ2​

面内主应力

C = 100 MPa, R = √(60² + (−20)²) = 63.25 MPa
σ1 = 163.2 MPa, σ2 = 36.75 MPa, σz = 0

三个圆中最大的

τ = max(|σ1 − σ2|, |σ1|, |σ2|)/2 = max(126.5, 163.2, 36.75)/2 = 81.62 MPa

答案

τabs = 81.62 MPa

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

10.莫尔圆与应力变换

答案: τmax = 39.54 MPa

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

表面处的应力

σ = F/A = −100000/(π·60²/4) = −100000/2827 = −35.37 MPa
τ = 16T/(πd³) = 16·1500000/(π·60³) = 35.37 MPa
στxyCτmax​

莫尔圆

C = σ/2 = −17.68 MPa
R = √((−17.68)² + 35.37²) = 39.54 MPa

答案

τmax = 39.54 MPa
11.梁的挠度

答案: δ = 24.26 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⁴。

w = 8 N/mm, L = 2500 mm
EI = 1.61 × 10¹² N·mm²

挠度

δ = wL⁴/(8EI) = 8·2500⁴/(8·1.61 × 10¹²) = 24.26 mm

答案

δ = 24.26 mm
12.梁的挠度

答案: δmax = 5.826 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 = 19000 N, w = 8 N/mm, L = 2500 mm, EI = 1.76 × 10¹² N·mm²

各荷载单独作用

δP = PL³/(48EI) = 19000·2500³/(48·1.76 × 10¹²) = 3.514 mm
δw = 5wL⁴/(384EI) = 5·8·2500⁴/(384·1.76 × 10¹²) = 2.312 mm

相加

δmax = 3.514 + 2.312 = 5.826 mm

答案

δmax = 5.826 mm
13.柱的屈曲

答案: Pcr = 284.8 kN

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

最小惯性矩

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

I = πd⁴/64 = π·60⁴/64 = 636173 mm⁴

计算长度

KL = 0.7·3000 = 2100 mm

临界荷载

Pcr = π²EI/(KL)² = π²·200000·636173/2100² = 284752 N = 284.8 kN

答案

Pcr = 284.8 kN
14.柱的屈曲

答案: d = 36.47 mm

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

所需临界荷载

Pcr = n·P = 3·20000 = 60000 N
KL = 1·1000 = 1000 mm

求解直径

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

答案

d = 36.47 mm

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

15.抛体运动

答案: x = 75.75 m

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

分解初速度

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

vx = v0·cos θ = 34·cos 70° = 11.63 m/s
vy = v0·sin θ = 34·sin 70° = 31.95 m/s

飞行时间

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

t = 2·31.95/9.81 = 6.514 s

水平距离

x = vx·t = 11.63·6.514 = 75.75 m

答案

x = 75.75 m
16.抛体运动

答案: x = 81.84 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 = 15.5 m/s, vy = 26.85 m/s

落到屋顶的时间

较大的根:球正在下降。

5 = 26.85·t − 4.905·t²
t = (26.85 + √(26.85² − 2·9.81·5))/9.81 = 5.28 s

水平距离

x = vx·t = 15.5·5.28 = 81.84 m

答案

x = 81.84 m
17.功和能

答案: v = 10.86 m/s

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

沿斜面的能量衡算

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

½v² = g·d·(sin θ − μ·cos θ)
v = √(2·9.81·10·(0.7071 − 0.15·0.7071)) = 10.86 m/s

答案

v = 10.86 m/s
18.功和能

答案: x = 7.004 cm

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

通过粗糙段后剩余的能量

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

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

最大压缩量

½kx² = 4.905 → x = √(2·4.905/2000) = 0.07004 m = 7.004 cm

答案

x = 7.004 cm
19.机械振动

答案: fd = 1.467 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 = √(2000/12) = 12.91 rad/s
fn = 12.91/(2π) = 2.055 Hz

有阻尼固有频率

fd = fn·√(1 − ζ²) = 2.055·√(1 − 0.7²) = 1.467 Hz

答案

fd = 1.467 Hz
20.机械振动

答案: keq = 125 kN/m

质量–弹簧系统的自由振动:无阻尼固有角频率为 ωn = √(k/m)(rad/s),fn = ωn/(2π)(Hz)。黏性阻尼器 c 给出阻尼比 ζ = c/cc,其中临界阻尼 cc = 2·√(k·m) = 2·m·ωn。欠阻尼系统(ζ < 1)以 ωd = ωn·√(1 − ζ²) 振荡。并联弹簧刚度相加,k = k₁ + k₂;串联弹簧按倒数相加,1/k = 1/k₁ + 1/k₂。

等效刚度

并排的弹簧变形相同、力相加,所以刚度直接相加:k = k1 + k2。

keq = 100 + 25 = 125 kN/m

答案

keq = 125 kN/m
21.热机和制冷机

答案: Wmax = 2725 kJ

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

绝对温度

TH = 750 + 273.15 = 1023 K
TL = 35 + 273.15 = 308.1 K

卡诺效率

ηmax = 1 − TL/TH = 1 − 308.1/1023 = 69.88 %

最大功

Wmax = ηmax·QH = 0.6988·3900 = 2725 kJ

答案

Wmax = 2725 kJ
22.热机和制冷机

答案: Sgen = 2.284 kJ/K

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

排放的热量

QL = QH − W = 5000 − 2950 = 2050 kJ

绝对温度

TH = 900 + 273.15 = 1173 K
TL = 40 + 273.15 = 313.1 K

热源的熵变

ΔSH = −5000/1173 = −4.262 kJ/K
ΔSL = 2050/313.1 = 6.546 kJ/K
Sgen = 6.546 − 4.262 = 2.284 kJ/K

答案

Sgen = 2.284 kJ/K

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

23.气体动力循环

答案: T3 = 2078 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 = 320·190.4 = 1039 K

定压加热

T3 = rc·T2 = 2·1039 = 2078 K

答案

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

答案: wnet = 171 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 = 288·60.2857 = 480.5 K
T4s = 1100/60.2857 = 659.3 K

压气机功和涡轮功

wc = 1.005·(480.5 − 288)/0.85 = 227.6 kJ/kg
wt = 0.9·1.005·(1100 − 659.3) = 398.6 kJ/kg

净功

wnet = 398.6 − 227.6 = 171 kJ/kg

答案

wnet = 171 kJ/kg

回功比:wc/wt = 57.1%。

25.连续性和伯努利

答案: p2 = 159.7 kPa

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

连续性和伯努利

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

v2 = v1·(d1/d2)² = 0.5·(100/75)² = 0.8889 m/s
½ρ(v1² − v2²) = ½·1000·(0.5² − 0.8889²)/1000 = −0.2701 kPa
p2 = p1 + ½ρ(v1² − v2²) = 160 + (−0.2701) = 159.7 kPa

答案

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

答案: p2 = 264.7 kPa

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

连续性和伯努利

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

v2 = v1·(d1/d2)² = 2.25·(100/75)² = 4 m/s
½ρ(v1² − v2²) = ½·1000·(2.25² − 4²)/1000 = −5.469 kPa
ρgΔz = 1000·9.81·1/1000 = 9.81 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 280 + (−5.469) − 9.81 = 264.7 kPa

答案

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

答案: Δp = 51.84 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.018·(300/0.3)·2.4²/(2·9.81) = 5.284 m

压降

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

Δp = ρ·g·hf = 1000·9.81·5.284 = 51840 Pa = 51.84 kPa

答案

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

答案: H = 6.518 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.15/0.07069 = 2.122 m/s
V²/(2g) = 0.2295 m

损失系数

f·L/D = 0.018·400/0.3 = 24
ΣK = 0.5 + 3·0.3 + 2 = 3.4

能量方程

出口射流仍带有速度水头;它不是损失,但必须由水头提供。

H = (1 + 24 + 3.4)·0.2295 = 6.518 m

答案

H = 6.518 m
29.水泵与泵功率

答案: Pshaft = 11.45 kW

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

水功率

Ph = ρ·g·Q·H = 1000·9.81·0.01·70 = 6.867 kW

轴功率

Pshaft = Ph/η = 6.867/0.6 = 11.45 kW

答案

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

答案: H = 43.86 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 + 5)·0.3305 = 11.57 m

压力水头

p/(ρ·g) = 150000/(1000·9.81) = 15.29 m

泵扬程

两个水面都基本静止,所以没有速度水头项。

H = Δz + p/(ρ·g) + hL = 17 + 15.29 + 11.57 = 43.86 m

答案

H = 43.86 m
31.热传导

答案: Ti = 11.8 °C

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

各层热阻

R1 = L1/(k1A) = (200/1000)/(0.17·12) = 0.09804 K/W
R2 = L2/(k2A) = (50/1000)/(0.17·12) = 0.02451 K/W

热流量

Q̇ = (T1 − T3)/(R1 + R2) = 14/0.1225 = 114.2 W

交界面温度

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

Ti = T1 − Q̇·R1 = 23 − 114.2·0.09804 = 11.8 °C

答案

Ti = 11.8 °C
32.热传导

答案: Ts = 18.93 °C

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

串联热阻

Ri = 1/(hiA) = 1/(6·4) = 0.04167 K/W
R1 = L1/(k1A) = (200/1000)/(0.17·4) = 0.2941 K/W
R2 = L2/(k2A) = (25/1000)/(0.04·4) = 0.1563 K/W
Ro = 1/(hoA) = 1/(25·4) = 0.01 K/W
ΣR = 0.502 K/W

热损失

Q̇ = (22 − (-15))/0.502 = 73.7 W
Ts = Troom − Q̇·Ri = 22 − 73.7·0.04167 = 18.93 °C

答案

Ts = 18.93 °C
33.齿轮与齿轮系

答案: T1 = 49.39 N·m

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

电动机扭矩

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

答案

T1 = 49.39 N·m
34.齿轮与齿轮系

答案: ND = 36

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 = 1440·36/60 = 864 rpm

齿轮 D 的齿数

ND = C·nC/nout = 24·864/576 = 36

答案

ND = 36
35.拉普拉斯变换

答案: f(0.5) = 2.03

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

部分分式

2s + 34 = 2(s + 4) + 26
F(s) = 2/(s + 4) + 26/(s + 4)²

逆变换

f(t) = 2·e−4t + 26·t·e−4t
f(0.5) = (2 + 26·0.5)·e−4·0.5 = 2.03

答案

f(0.5) = 2.03
36.拉普拉斯变换

答案: y(0.1) = 9.52

拉普拉斯变换: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 − 10s − 3 + 9(sY − 10) + 18Y = 0
Y(s) = (10s + 93)/((s + 3)(s + 6))

部分分式

A = 63/3 = 21
B = 33/(−3) = −11

结果

y(t) = 21·e−3t − 11·e−6t
y(0.1) = 9.52

答案

y(0.1) = 9.52

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

37.方框图

答案: css = 18.51

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

结果

css = R·T(0) = 10·31/16.75 = 18.51

答案

css = 18.51
38.方框图

答案: ωd = 5.657 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(s + 12) + 68·1 = 0 → s² + 12s + 68 = 0

根

s = −6 ± j·√(68 − 36) = −6 ± j5.657
ωd = 5.657 rad/s

答案

ωd = 5.657 rad/s
39.阶跃响应

答案: %OS = 50.56 %

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

ωₙ² = 50 → ωₙ = 7.071 rad/s
2ζωₙ = 3 → ζ = 3/(2·7.071) = 0.2121

结果

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

答案

%OS = 50.56 %
40.阶跃响应

答案: ζ = 0.2016

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

固有频率和阻尼比

ωₙ = √(k/m) = √(2000/8) = 15.81 rad/s
ζ = c/(2√(km)) = 51/(2·√(2000·8)) = 0.2016

结果

ζ = 0.2016

答案

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

答案: ess = 5

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

速度误差常数

Kv = lim s·G(s) = 1

阶跃部分

Kp = ∞ → estep = 0

斜坡误差

ess = 5/1 = 5

答案

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

答案: tp = 0.5236 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 + 4) + 40·1 = s² + 4s + 40

ωₙ 和 ζ

ωₙ = √40 = 6.325 rad/s
ζ = 4/(2·6.325) = 0.3162

峰值时间

ωd = 6.325·√(1 − 0.3162²) = 6 rad/s
tp = π/6 = 0.5236 s

答案

tp = 0.5236 s

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

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