机械自动化控制理论试题详解

机械自动化控制理论试题详解机械自动化控制理论是工科学生，特别是机械工程、自动化及相关专业的核心课程。它不仅是理解现代工业自动化设备工作原理的基石，也是进行系统设计与优化的理论依据。掌握这门学科，离不开对基本概念的深刻理解和对解题方法的熟练运用。本文将通过对若干典型试题的详细解析，帮助读者回顾核心知识点，梳理解题思路，提升分析与解决实际控制问题的能力。一、控制理论核心知识点回顾在深入试题之前，我们先来简要回顾一下控制理论的核心框架，这对于解题至关重要：1.控制系统的基本概念：包括被控对象、输入量、输出量、扰动量、反馈、开环与闭环控制等。2.数学模型：是分析控制系统的基础。主要包括微分方程、传递函数、方框图、信号流图等。建立准确的数学模型是解决控制问题的第一步。3.时域分析法：直接在时间域内分析系统的动态性能（如超调量、调节时间、峰值时间等）和稳态性能（稳态误差）。重点是一阶、二阶系统的响应分析。4.根轨迹法：根据系统开环传递函数的零点和极点，研究当某个参数变化时闭环极点在s平面上的运动轨迹，进而分析系统性能的变化。5.频域分析法：利用频率特性（幅频特性、相频特性）来分析系统的稳定性、动态性能和稳态性能。主要工具包括Nyquist图、Bode图。6.稳定性判据：判断系统是否稳定是控制理论的核心问题之一。主要有Routh-Hurwitz判据（代数判据）、Nyquist稳定判据（几何判据）等。7.控制系统的校正与设计：根据系统的性能指标要求，设计或选择合适的控制器（如PID控制器），以改善系统性能。二、典型试题详解试题一：数学模型的建立与化简题目：考虑如图所示的机械平移系统。质量块m的质量为m，弹簧的弹性系数为k，阻尼器的阻尼系数为c。外力F(t)为输入，质量块的位移x(t)为输出。(1)试建立该系统的微分方程。(2)若初始条件为零，试求该系统的传递函数G(s)=X(s)/F(s)。(3)若在质量块与阻尼器之间再串联一个弹性系数为k2的弹簧，试画出此时系统的力学模型，并简述如何修改传递函数的推导过程（无需详细计算）。解题思路：本题主要考察机械系统数学模型的建立，核心是应用牛顿第二定律。对于线性系统，微分方程是时域模型，传递函数是复频域模型，两者通过拉普拉斯变换相联系。详细解答：(1)建立微分方程：对质量块m进行受力分析。质量块受到的力包括：*外力F(t)(输入)*弹簧的弹性力：-kx(t)(与位移方向相反)*阻尼器的粘性阻尼力：-cdx(t)/dt(与速度方向相反)根据牛顿第二定律ΣF=ma，有：F(t)-kx(t)-cdx(t)/dt=md²x(t)/dt²整理可得系统的微分方程：md²x(t)/dt²+cdx(t)/dt+kx(t)=F(t)---(1)(2)求传递函数G(s)=X(s)/F(s)：对微分方程(1)两边进行拉普拉斯变换，注意到初始条件为零（x(0)=0,dx(0)/dt=0），则：ms²X(s)+csX(s)+kX(s)=F(s)将X(s)提取出来：X(s)(ms²+cs+k)=F(s)因此，传递函数为：G(s)=X(s)/F(s)=1/(ms²+cs+k)---(2)(3)串联弹簧后的力学模型及传递函数推导思路：*力学模型：此时，阻尼器两端分别与质量块m和新弹簧k2相连，新弹簧k2的另一端固定（或连接到原固定端）。或者说，原系统的阻尼器与一个新弹簧k2串联后，再与原弹簧k并联（取决于原系统阻尼器的安装方式，通常理解为阻尼器与k2串联后整体再与k并联作用于质量块）。更准确地说，若原系统是质量块连接弹簧k到固定端，连接阻尼器c到固定端；则修改后是质量块连接弹簧k到固定端，连接阻尼器c到一个弹簧k2的一端，弹簧k2的另一端连接到固定端。*传递函数推导修改：此时，阻尼器的位移不再是x(t)，而是质量块位移x(t)与弹簧k2的位移x2(t)之差。设流经阻尼器的力为fd(t)=cd(x(t)-x2(t))/dt。而弹簧k2的力fk2(t)=k2x2(t)。由于阻尼器和弹簧k2串联，其受力相等，即fd(t)=fk2(t)。因此，可以得到一个关于x(t)和x2(t)的方程组，消去中间变量x2(t)后，即可得到新的输入输出关系和传递函数。这会引入一个新的储能元件，系统阶次可能保持二阶，但传递函数的分母系数会发生变化。点评与拓展：本题是最基本的单自由度机械系统建模。关键在于正确的受力分析和牛顿定律的应用。传递函数的推导则依赖于拉普拉斯变换的线性性质和微分性质。对于更复杂的机械系统（多质量、齿轮系、旋转系统等），方法类似，都是从元件的力-位移（或力矩-角位移）关系入手，列写运动微分方程。试题二：时域分析与性能指标题目：已知某单位负反馈控制系统的开环传递函数为G(s)=10/[s(s+2)]。(1)求系统的闭环传递函数Φ(s)。(2)分析该闭环系统的稳定性。(3)求系统单位阶跃响应的调节时间ts(±5%误差带)、超调量σ%和峰值时间tp。解题思路：对于单位负反馈系统，闭环传递函数Φ(s)=G(s)/(1+G(s))。系统稳定性可通过闭环极点或开环Nyquist图、Routh判据等判断。对于二阶系统，其典型单位阶跃响应性能指标有解析公式，可直接套用。详细解答：(1)求闭环传递函数Φ(s)：单位负反馈系统的闭环传递函数为：Φ(s)=G(s)/[1+G(s)H(s)]，其中H(s)=1。代入G(s)：Φ(s)=[10/(s(s+2))]/[1+10/(s(s+2))]=10/[s(s+2)+10]分母展开：s²+2s+10所以Φ(s)=10/(s²+2s+10)---(3)(2)分析闭环系统的稳定性：方法一：根据闭环传递函数的极点。闭环特征方程为：s²+2s+10=0解此二次方程：s=[-2±√(4-40)]/2=[-2±√(-36)]/2=[-2±j6]/2=-1±j3两个闭环极点均具有负实部（-1），因此系统是稳定的。方法二：Routh-Hurwitz判据。特征方程：s²+2s+10=0Routh表：s²:110s¹:20s⁰:10(由(2*10-1*0)/2=10得到)第一列元素均为正（1,2,10），故系统稳定。(3)求单位阶跃响应的性能指标：单位负反馈系统，输入R(s)=1/s。输出C(s)=Φ(s)R(s)=10/[s(s²+2s+10)]将闭环传递函数分母与标准二阶系统传递函数分母对比：标准二阶系统传递函数：Φ(s)=ωn²/(s²+2ζωns+ωn²)其中，ωn为无阻尼固有频率，ζ为阻尼比。对比式(3)Φ(s)=10/(s²+2s+10)，可得：ωn²=10→ωn=√10≈3.162rad/s2ζωn=2→ζ=2/(2ωn)=1/ωn=1/√10≈0.316(欠阻尼，ζ<1)对于欠阻尼二阶系统，单位阶跃响应的性能指标公式：*超调量σ%：σ%=e^(-πζ/√(1-ζ²))×100%*峰值时间tp：tp=π/(ωn√(1-ζ²))*调节时间ts(±5%误差带)：ts≈3/(ζωn)(当ζ<0.9时)代入数值计算：√(1-ζ²)=√(1-1/10)=√(9/10)=3/√10≈0.9487σ%=e^(-π*(1/√10)/(3/√10))×100%=e^(-π/3)×100%≈e^(-1.047)×100%≈0.358×100%≈35.8%tp=π/(ωn*(3/√10)))=π√10/(3ωn)。由于ωn=√10，故tp=π√10/(3√10))=π/3≈1.047秒。ts(±5%)≈3/(ζωn)。已知2ζωn=2→ζωn=1。因此ts≈3/1=3秒。点评与拓展：本题考察了闭环传递函数的计算、稳定性判断以及二阶系统时域性能指标的计算。对于二阶系统，熟记标准形式和性能指标公式至关重要。注意，调节时间的近似公式有不同的误差带（5%或2%），需看清题目要求。当系统并非严格二阶时，可能需要通过劳斯判据判断主导极点，然后用主导极点法近似为二阶系统进行分析。试题三：频域分析与稳定性题目：某最小相位系统的开环对数幅频特性渐近线如图所示（此处省略图形，假设根据图形可以得到以下信息：低频段斜率为-20dB/dec，延长线交0dB线于ω=1rad/s；在ω=2rad/s处有一个转折频率，斜率变为-40dB/dec；在ω=10rad/s处有另一个转折频率，斜率变为-20dB/dec）。(1)试根据开环对数幅频特性渐近线确定系统的开环传递函数G(s)。(2)若该系统为单位负反馈，试用Nyquist稳定判据判断闭环系统的稳定性。(3)若要求系统的相角裕度γ≥45°，简述可采取哪些校正措施（至少列举两种）。解题思路：根据Bode图（对数幅频特性渐近线）确定最小相位系统的开环传递函数，是频域分析的基本技能。关键在于识别各频段的斜率变化，从而确定对应的典型环节（比例环节、积分环节、微分环节、惯性环节、一阶微分环节、振荡环节等）。Nyquist稳定判据则需要根据开环传递函数的极点分布（特别是右半s平面的极点数P）和Nyquist曲线包围(-1,j0)点的圈数N来判断闭环稳定性。详细解答：(1)确定开环传递函数G(s)：最小相位系统的开环传递函数可由对数幅频特性渐近线唯一确定。分析步骤：a)确定低频段斜率：低频段（ω→0）的斜率为-20νdB/dec，其中ν为积分环节的个数（型别）。题目中低频段斜率为-20dB/dec，故ν=1(I型系统)。b)确定开环增益K：对于I型系统，低频段渐近线（或其延长线）与0dB线的交点频率ωc0满足K=ωc0^ν。此处ν=1，延长线交0dB线于ω=1rad/s，故K=ωc0=1。c)分析转折频率及对应的典型环节：*第一个转折频率ω1=2rad/s：斜率从-20dB/dec→-40dB/dec，斜率变化了-20dB/dec，对应一个惯性环节（1/(Ts+1)，T=1/ω1=1/2）。*第二个转折频率ω2=10rad/s：斜率从-40dB/dec→-20dB/dec，斜率变化了+20dB/dec，对应一个一阶微分环节（Ts+1，T=1/ω2=1/10）。综合以上分析，开环传递函数的结构为：G(s)=K*(T2s+1)/[s(T1s+1)]其中：K=1(由步骤b得)T1=1/ω1=1/2→1/(T1s+1)=1/((s/2)+1)=2/(s+2)T2=1/ω2=1/10→(T2s+1)=(s/10+1)=(s+10)/10代入得：G(s)=1*[(s+10)/10]/[s*(s+2)/2]=[(s+10)/10]*[2/(s(s+2))]=[2(s+10)]/[10s(s+2)]=(s+10)/[5s(s+2)]检查增益：当ω→0时，G(jω)≈K/(jω)=1/(jω)。而我们推导出的G(s)在ω→0时，G(jω)≈(10)/(5jω*2))→这里似乎有计算错误。让我们重新计算K。更准确的求K方法：低频段渐近线方程。对于I型系统，低频段（ω远小于第一个转折频率ω1）的对数幅频特性为：L(ω)=20lgK-20νlgω=20lgK-20lgω(因为ν=1)已知该直线过点(ω=1,L(ω)=0dB)，代入得：0=20lgK-20lg1→20lgK=0→K=1。现在，我们推导出的G(s)=(s+10)/(5s(s+2))，当ω→0时，s=jω→0，G(jω)≈10/(5jω*2))=10/(10jω)=1/(jω)，这与L(ω)=20lg|G(jω)|=20lg(1/ω)=-20lgω，当ω=1时，L(1)=0dB，完全一致。之前的中间步骤化简可能造成了误解，但最终结果是正确的。因此，开环传递函数为：G(s)=(s+10)/[5s(s+2)](2)用Nyquist稳定判据判断闭环系统稳定性：Nyquist稳定判据：闭环系统稳定的充要条件是Nyquist曲线G(jω)包围(-1,j0)点的圈数N与开环传递函数在右半s平面的极点数P满足