2026年《机械优化设计》试卷及答案

2026年《机械优化设计》试卷及答案一、单项选择题（每题2分，共10分）1.机械优化设计中，若目标函数为f(x)=x₁²+2x₂²-3x₁x₂+5，其梯度向量∇f(x)为（）。A.[2x₁-3x₂,4x₂-3x₁]B.[2x₁+3x₂,4x₂+3x₁]C.[x₁-3x₂,2x₂-3x₁]D.[2x₁-3x₂,2x₂-3x₁]2.以下不属于无约束优化算法的是（）。A.最速下降法B.牛顿法C.惩罚函数法D.共轭梯度法3.遗传算法中，若群体规模为50，交叉概率为0.8，变异概率为0.01，则每代预计产生的新个体数约为（）。A.40B.50C.80D.104.对于约束优化问题minf(x)s.t.g₁(x)≤0,g₂(x)≥0，其可行域是（）。A.g₁(x)≤0与g₂(x)≥0的交集B.g₁(x)≤0与g₂(x)≥0的并集C.g₁(x)≤0或g₂(x)≥0D.g₁(x)≥0与g₂(x)≤0的交集5.某机械结构优化中，设计变量为截面宽度b（mm）和高度h（mm），目标是最小化重量，则目标函数可表示为（）（材料密度为ρ）。A.f(b,h)=ρbhB.f(b,h)=ρbhL（L为长度）C.f(b,h)=bhD.f(b,h)=ρ(b+h)L二、填空题（每空2分，共20分）1.机械优化设计的三要素是________、________和________。2.目标函数的等值线（面）越密集，说明该方向上函数的________变化越剧烈。3.梯度法的搜索方向是目标函数的________方向，牛顿法的搜索方向由________矩阵的逆矩阵与梯度的乘积确定。4.遗传算法的基本操作包括________、________和________。5.对于约束优化问题，当设计点位于可行域内部时，起作用的约束数量为________；若位于边界，则至少有一个约束的________等于零。三、简答题（每题8分，共32分）1.简述机械优化设计中“可行解”与“最优解”的区别，并说明在实际工程中为何通常只能得到“满意解”。2.比较最速下降法与牛顿法的优缺点，说明各自适用的场景。3.解释遗传算法中“群体多样性”的重要性，并列举3种提高群体多样性的方法。4.某齿轮传动优化设计中，需将接触疲劳强度作为约束条件，试推导该约束的数学表达式（已知接触应力公式为σ_H=Z_HZ_E√(F_tu/(bd₁))，许用接触应力为[σ_H]）。四、计算题（每题12分，共36分）1.用黄金分割法（初始区间[a,b]=[1,5]，精度要求ε=0.2）求单变量函数f(x)=x²-4x+5的极小值点。要求列出前3次迭代的区间缩小区间和计算点函数值。2.已知二元函数f(x)=x₁²+2x₂²+2x₁x₂-4x₁-6x₂，初始点x⁰=[0,0]^T，用最速下降法求其极小值点（迭代2次，步长α采用精确线搜索）。3.求解约束优化问题：minf(x)=x₁²+x₂²，s.t.g(x)=x₁+x₂-1≥0。要求用外点惩罚函数法构造惩罚函数（惩罚因子M=100），并计算初始点x⁰=[0,0]^T处的惩罚函数值及梯度。五、综合分析题（22分）某带式输送机传动系统中，需优化设计一级圆柱齿轮减速器的高速级齿轮副。已知输入功率P=10kW，输入转速n₁=1450r/min，传动比i=4，齿轮材料为40Cr（调质，齿面硬度240HBS），许用接触应力[σ_H]=600MPa，许用弯曲应力[σ_F]=200MPa，模数m≥2mm，齿数z₁≥20，齿宽系数φ_d=0.8（φ_d=b/d₁，d₁为小齿轮分度圆直径）。设计目标为最小化齿轮副的体积（体积近似为πd₁²b/4+πd₂²b/4，d₂为大齿轮分度圆直径）。要求：（1）确定设计变量；（2）建立目标函数的数学表达式；（3）列出所有约束条件（包括不等式约束和等式约束）；（4）说明选择优化算法时需考虑的因素，并推荐一种适合的算法（需说明理由）。答案一、单项选择题1.A2.C3.A4.A5.B二、填空题1.设计变量、目标函数、约束条件2.梯度（或变化率）3.负梯度、海森（Hessian）4.选择、交叉、变异5.0、约束函数值三、简答题1.可行解是满足所有约束条件的设计变量组合；最优解是可行解中使目标函数达到极值的解。实际工程中，由于目标函数可能存在多个局部极值、约束条件复杂（如非线性、离散性）、计算资源限制等，通常难以找到全局最优解，因此退而求其次选择满足工程要求的“满意解”。2.最速下降法优点：计算简单，只需求梯度；缺点：收敛速度慢（尤其在极值点附近呈锯齿形搜索）。适用于低精度或初始阶段的优化。牛顿法优点：收敛速度快（二次函数可一步收敛）；缺点：需计算海森矩阵（计算复杂且可能不可逆）。适用于目标函数二阶可导且海森矩阵正定的情况。3.群体多样性是遗传算法避免早熟收敛的关键，多样性不足会导致群体陷入局部最优。提高方法：增大初始群体规模；采用自适应交叉/变异概率（如群体多样性低时提高变异概率）；引入移民操作（定期加入新个体）；使用非均匀变异算子。4.接触强度约束要求σ_H≤[σ_H]。由σ_H=Z_HZ_E√(F_tu/(bd₁))≤[σ_H]，平方后得：Z_H²Z_E²(F_tu)/(bd₁)≤[σ_H]²。其中F_t=2T₁/d₁（T₁为小齿轮转矩，T₁=9550×10³P/n₁），u=i=4，d₁=mz₁，b=φ_dd₁=φ_dmz₁。代入后约束可表示为：Z_H²Z_E²×(2×9550×10³P/n₁×4)/(φ_dmz₁×mz₁)≤[σ_H]²，即Z_H²Z_E²×(76400×10³P)/(φ_dm²z₁²n₁)≤[σ_H]²。四、计算题1.黄金分割法：初始区间[a₀,b₀]=[1,5]，λ=0.618，计算点x₁=a₀+0.382(b₀-a₀)=1+0.382×4=2.528，x₂=a₀+0.618(b₀-a₀)=1+0.618×4=3.472。f(x₁)=2.528²-4×2.528+5≈(6.39)-10.11+5≈1.28；f(x₂)=3.472²-4×3.472+5≈(12.06)-13.89+5≈3.17。因f(x₁)<f(x₂)，新区间[a₁,b₁]=[1,3.472]。第二次迭代：x₁'=1+0.382×(3.472-1)=1+0.382×2.472≈1.944，x₂'=x₁=2.528（利用对称性）。f(x₁')=1.944²-4×1.944+5≈(3.78)-7.78+5=1；f(x₂')=1.28（同前）。因f(x₁')<f(x₂')，新区间[a₂,b₂]=[1,2.528]。第三次迭代：区间长度=2.528-1=1.528>0.2，继续。x₁''=1+0.382×1.528≈1.584，x₂''=x₁'=1.944。f(x₁'')=1.584²-4×1.584+5≈(2.51)-6.34+5=1.17；f(x₂'')=1。因f(x₂'')<f(x₁'')，新区间[a₃,b₃]=[1.584,2.528]。最终极小值点近似为区间中点≈(1.584+2.528)/2≈2.056（实际极小值点x=2，因f(x)=x²-4x+5=(x-2)²+1，极小值点x=2）。2.最速下降法：f(x)=x₁²+2x₂²+2x₁x₂-4x₁-6x₂，梯度∇f(x)=[2x₁+2x₂-4,4x₂+2x₁-6]。初始点x⁰=[0,0]^T，∇f(x⁰)=[-4,-6]^T，搜索方向d⁰=-∇f(x⁰)=[4,6]^T。步长α⁰通过精确线搜索求解：f(x⁰+αd⁰)=(4α)²+2(6α)²+2(4α)(6α)-4(4α)-6(6α)=16α²+72α²+48α²-16α-36α=136α²-52α。求导得272α-52=0→α⁰=52/272=13/68≈0.191。第一次迭代后点x¹=x⁰+α⁰d⁰=[4×13/68,6×13/68]=[52/68,78/68]=[13/17,39/34]≈[0.765,1.147]。计算∇f(x¹)=[2×13/17+2×39/34-4,4×39/34+2×13/17-6]=[(26/17+78/34)-4,(156/34+26/17)-6]=[(52/34+78/34)-4,(156/34+52/34)-6]=[130/34-4,208/34-6]=[130/34-136/34,208/34-204/34]=[-6/34,4/34]=[-3/17,2/17]。搜索方向d¹=-∇f(x¹)=[3/17,-2/17]。步长α¹：f(x¹+αd¹)=(13/17+3α/17)²+2(39/34-2α/17)²+2(13/17+3α/17)(39/34-2α/17)-4(13/17+3α/17)-6(39/34-2α/17)。展开后求导并令导数为0，解得α¹=(计算过程略)，最终α¹≈0.5。第二次迭代后点x²≈[0.765+3/17×0.5,1.147-2/17×0.5]≈[0.852,1.088]（实际极小值点通过求导得∇f(x)=0→2x₁+2x₂=4，2x₁+4x₂=6，解得x₁=1，x₂=1）。3.外点惩罚函数法：原问题为minf(x)=x₁²+x₂²，s.t.g(x)=x₁+x₂-1≥0→-g(x)=1-x₁-x₂≤0。外点惩罚函数为P(x,M)=x₁²+x₂²+M[max(0,1-x₁-x₂)]²。初始点x⁰=[0,0]，g(x⁰)=0+0-1=-1<0，违反约束，故max(0,1-x⁰₁-x⁰₂)=1-0-0=1。惩罚函数值P(0,0,100)=0+0+100×1²=100。梯度∇P(x,M)=[2x₁+2M(1-x₁-x₂)(-1),2x₂+2M(1-x₁-x₂)(-1)]（当1-x₁-x₂>0时）。在x⁰处，1-x₁-x₂=1>0，故∇P(0,0,100)=[0+2×100×1×(-1),0+2×100×1×(-1)]=[-200,-200]。五、综合分析题（1）设计变量：小齿轮齿数z₁，模数m（均为连续变量，实际设计中需圆整）。（2）目标函数：体积V=πd₁²b/4+πd₂²b/4=πb(d₁²+d₂²)/4。因d₁=mz₁，d₂=mz₂=miz₁（z₂=iz₁），b=φ_dd₁=φ_dmz₁。代入得V=πφ_dmz₁[(mz₁)²+(miz₁)²]/4=πφ_dmz₁m²z₁²(1+i²)/4=πφ_dm³z₁³(1+i²)/4。故目标函数f(z₁,m)=km³z₁³（k=πφ_d(1+i²)/4为常数）。（3）约束条件：①接触强度约束：σ_H≤[σ_H]，由σ_H=Z_HZ_E√(F_tu/(bd₁))≤[σ_H]，代入F_t=2T₁/d₁=2×9550×10³P/(n₁d₁)=2×9550×10³×10/(1450×mz₁)=191000000/(1450mz₁)=131724.14/(mz₁)（N），u=i=4，b=φ_dd₁=0.8mz₁，d₁=mz₁，得σ_H=Z_HZ_E√[(131724.14×4)/(0.8mz₁×mz₁)]=Z_HZ_E√[658620.7/(0.8m²z₁²)]=Z_HZ_E√(823275.88)/(mz₁)≤600MPa。取Z_H=2.5，Z_E=189.8MPa^(1/2)，则约束为2.5×189.8×√823275.88/(mz₁)≤600→(2.5×189.8×907.3)/(mz₁)≤600→(436,500)/(mz₁)≤600→mz₁≥727.5。②弯曲强度约束：σ_F≤[σ_F]，σ_F=Y_FY_SF_t/(bm)≤200MPa。Y_F（齿形系数）、Y_S（应力修正系数）与z₁相关（如z₁=20时Y_F=2.8，Y_S=1.55），F_t=131724.14/(mz₁)，b=0.8mz₁，故σ_F=2.8×1.55×131724.14/(0.8mz₁×m)=(574,000)/(0.8m²z₁)=717,500/(m²z₁)≤200→m²z₁≥3