2026年弹性力学试题及答案

2026年弹性力学试题及答案一、填空题（每空2分，共20分）1.弹性力学中，应变协调方程的物理意义是保证物体变形的连续性。2.平面问题中，Airy应力函数φ需满足的基本条件是∇⁴φ=0（双调和方程）。3.圣维南原理的核心思想是分布在物体局部区域的力系，可用静力等效的力系替代，仅对附近区域的应力分布有显著影响，远处影响可忽略。4.极坐标下，轴对称问题的平衡微分方程为dσᵣ/dr+(σᵣσθ)/r=0。5.各向同性线弹性材料在平面应力状态下的广义胡克定律表达式为εₓ=(σₓ-νσᵧ)/E，εᵧ=(σᵧ-νσₓ)/E，γₓᵧ=τₓᵧ/G（其中G=E/[2(1+ν)]）。6.小变形假设下，位移分量u、v与应变分量的几何关系中，切应变γₓᵧ的表达式为γₓᵧ=∂u/∂y+∂v/∂x。7.弹性力学中，应力边界条件的数学表达式本质是面力分量与应力分量的积分关系（微元体平衡）。8.对于纯扭转的圆截面杆，其横截面上的切应力分布满足τ=Tr/J（T为扭矩，r为点到圆心距离，J为极惯性矩）。9.半平面体受法向集中力P作用时，其应力分布的极坐标表达式中，σᵣ与P的关系为σᵣ=-2Pcosθ/(πr)。10.弹性力学问题的基本解法分为位移法和应力法，其中应力法需引入应变协调方程。二、简答题（每题10分，共30分）1.简述平面应力问题与平面应变问题的主要区别。答：平面应力问题假设物体为等厚度薄板（厚度远小于板面尺寸），仅在板平面内有应力（σₓ、σᵧ、τₓᵧ），厚度方向应力σ_z≈0；应变在厚度方向（ε_z）不为零（由泊松效应引起）。物理方程中，应变与应力关系为εₓ=(σₓ-νσᵧ)/E，εᵧ=(σᵧ-νσₓ)/E，γₓᵧ=τₓᵧ/G，ε_z=-ν(σₓ+σᵧ)/E。平面应变问题假设物体为长柱体（长度远大于横截面尺寸），所有横截面变形相同，轴向位移w=0，故轴向应变为零（ε_z=0），但轴向应力σ_z=ν(σₓ+σᵧ)不为零；横截面上的应力（σₓ、σᵧ、τₓᵧ）与z无关。物理方程中，应变与应力关系需将E替换为E/(1-ν²)，ν替换为ν/(1-ν)，即εₓ=(σₓ-νσᵧ)(1+ν)/E，εᵧ=(σᵧ-νσₓ)(1+ν)/E，γₓᵧ=τₓᵧ/G。2.推导小变形条件下的几何方程（应变与位移的关系），并说明其物理意义。答：考虑微元体在x、y方向的位移u、v，沿x方向的正应变εₓ定义为单位长度的伸长量，即εₓ=∂u/∂x；同理，沿y方向的正应变εᵧ=∂v/∂y。切应变γₓᵧ定义为直角的改变量，微元体在x方向位移u引起y方向的偏转角为∂u/∂y，y方向位移v引起x方向的偏转角为∂v/∂x，故γₓᵧ=∂u/∂y+∂v/∂x（小变形下忽略高阶小量）。几何方程的物理意义是：将连续介质的变形（应变）与质点的相对位移（位移梯度）联系起来，确保变形的连续性（即相邻微元体变形后无重叠或分离）。3.虚功原理在弹性力学中的表述是什么？它与最小势能原理的关系如何？答：虚功原理表述为：在弹性体的平衡状态下，外力在任意虚位移上做的虚功等于弹性体内力在虚应变上做的虚功，即∫_Vσ_ijδε_ijdV=∫_ST_iδu_idS+∫_Vf_iδu_idV（σ_ij为应力，δε_ij为虚应变，T_i为面力，f_i为体力，δu_i为虚位移）。最小势能原理是虚功原理的能量表达形式，其核心是：在所有满足位移边界条件的可能位移中，真实位移使弹性体的总势能（应变能减去外力势能）取最小值。两者本质等价，虚功原理是能量守恒的微分形式，最小势能原理是其积分形式的变分表述。三、计算题（每题15分，共30分）1.无限长圆柱体，内半径a=100mm，外半径b=200mm，受内压p_i=50MPa、外压p_o=10MPa作用，材料弹性模量E=200GPa，泊松比ν=0.3。求圆柱体内的应力分布（用极坐标表示），并计算r=a、r=b处的环向应力σθ。解：由于圆柱体无限长且受轴对称载荷，属平面应变问题（ε_z=0）。极坐标下轴对称问题的应力分量仅与r有关，平衡微分方程为dσᵣ/dr+(σᵣσθ)/r=0（1）。几何方程：εᵣ=dε/dr，εθ=u/r（u为径向位移），γᵣθ=0（轴对称无剪切）。物理方程（平面应变）：εᵣ=(σᵣνσθ)/[(1+ν)E/(1-ν)]，εθ=(σθνσᵣ)/[(1+ν)E/(1-ν)]（化简后为εᵣ=(1+ν)(σᵣνσθ)/E，εθ=(1+ν)(σθνσᵣ)/E）。由几何方程得d(ru)/dr=u+rdu/dr=rεᵣ+u=rεθ（因εᵣ=du/dr，εθ=u/r），代入物理方程消去u，可得σθ=σᵣ+rdσᵣ/dr（由平衡方程推导）。假设应力函数为φ(r)=Alnr+Br²（轴对称问题的双调和方程解），则σᵣ=φ’/rφ''=(A/r²)+2B，σθ=-φ''=(-A/r²)+2B（验证满足平衡方程）。边界条件：r=a时，σᵣ=-p_i；r=b时，σᵣ=-p_o。代入得：A/a²+2B=-p_iA/b²+2B=-p_o解得：A=(p_op_i)a²b²/(b²a²)，B=(p_ia²p_ob²)/[2(b²a²)]。因此，应力分布为：σᵣ=(p_op_i)a²b²/[r²(b²a²)]+(p_ia²p_ob²)/(b²a²)σθ=-(p_op_i)a²b²/[r²(b²a²)]+(p_ia²p_ob²)/(b²a²)τᵣθ=0（轴对称无剪切）。计算r=a处σθ：σθ(a)=-(p_op_i)b²/(b²a²)+(p_ia²p_ob²)/(b²a²)=[-p_ob²+p_ib²+p_ia²p_ob²]/(b²a²)=[p_i(a²+b²)2p_ob²]/(b²a²)代入数值：p_i=50MPa，p_o=10MPa，a=100mm，b=200mm，σθ(a)=[50×(100²+200²)2×10×200²]/(200²-100²)=[50×500002×10×40000]/30000=(2,500,000800,000)/30,000=1,700,000/30,000≈56.67MPa。r=b处σθ：σθ(b)=-(p_op_i)a²/(b²a²)+(p_ia²p_ob²)/(b²a²)=[-p_oa²+p_ia²+p_ia²p_ob²]/(b²a²)=[2p_ia²p_o(a²+b²)]/(b²a²)代入数值：σθ(b)=[2×50×100²10×(100²+200²)]/(200²-100²)=[1,000,00010×50,000]/30,000=(1,000,000500,000)/30,000=500,000/30,000≈16.67MPa。2.矩形板（0≤x≤a，0≤y≤b）上边界y=b受线性载荷q(x)=q₀x/a（q₀为常数），下边界y=0自由（τₓᵧ=0，σᵧ=0），左右边界x=0、x=a无面力（τₓᵧ=0，σₓ=0）。试采用多项式应力函数法求解板内应力分量（σₓ、σᵧ、τₓᵧ）。解：假设应力函数为φ(x,y)=Ay³+By²x+Cyx³+Dx⁴（多项式形式，满足双调和方程∇⁴φ=0）。应力分量：σₓ=∂²φ/∂y²=6Ay+2Bx+6Cyxσᵧ=∂²φ/∂x²=6Cy+12Dxτₓᵧ=-∂²φ/∂x∂y=-2By3Cx²4Dx³边界条件：（1）y=0时，σᵧ=0（自由边界）：代入σᵧ=6C×0+12D×0=0（自动满足）；τₓᵧ=0（自由边界）：-2B×03Cx²4Dx³=0，对任意x成立，故C=0，D=0。（2）y=b时，σᵧ=-q(x)=-q₀x/a（载荷向下为负）：σᵧ=6C×b+12D×b=0（因C=D=0），矛盾，需调整应力函数。改用φ=Ay³x+Byx（含x的一次项）。重新假设φ=Ay³x+Byx，计算应力分量：σₓ=∂²φ/∂y²=6Ayx+0=6Ayxσᵧ=∂²φ/∂x²=0+0=0（不满足y=b处σᵧ=-q₀x/a，需增加x的一次项）。修正为φ=Ay³x+Byx+Cx⁴（含x⁴项），则σᵧ=∂²φ/∂x²=12Cx²。要求y=b时σᵧ=-q₀x/a，即12Cb²x²=-q₀x/a（不匹配x的幂次），改用φ=Ay³x+Byx+Cx²y³（含x²y³项），则σᵧ=∂²φ/∂x²=2Cy³。y=b时σᵧ=2Cb³=-q₀x/a（仍不匹配x的线性项，说明需引入x的一次项与y的三次项的乘积）。正确假设：φ(x,y)=y³(Ax+B)+y(Cx+D)（线性于x，三次于y）。计算应力分量：σₓ=∂²φ/∂y²=6y(Ax+B)σᵧ=∂²φ/∂x²=0（不满足y=b处σᵧ=-q₀x/a），需增加φ=Exy⁴（四次项），则σᵧ=∂²φ/∂x²=0（仍不行）。换用Airy应力函数满足双调和方程的一般解，对于线性载荷q(x)=kx（k=q₀/a），应力函数可设为φ=Ay⁴+By³x+Cy²x²+Dyx³+Ex⁴（四次多项式）。双调和方程∇⁴φ=24A+24C+24E=0→A+C+E=0。应力分量：σₓ=∂²φ/∂y²=12Ay²+6Byx+2Cx²σᵧ=∂²φ/∂x²=2Cy²+6Dyx+12Ex²τₓᵧ=-∂²φ/∂x∂y=-6By²4Cxy3Dx²边界条件：（1）y=0时，σᵧ=0→2C×0+6D×0+12E×0=0（成立）；τₓᵧ=0→-6B×04C×03D×0=0（成立）。（2）y=b时，σᵧ=-kx（k=q₀/a）→2Cb²+6Dbx+12Ex²=-kx。因右边为x的一次项，故12E=0（x²项系数为0），6Db=-k（x项系数），2Cb²=0（常数项为0）→E=0，D=-k/(6b)，C=0。（3）x=0时，σₓ=0（无面力）→12A×y²+6B×0×y+2×0×0=12Ay²=0→A=0（对任意y成立）。（4）x=a时，σₓ=0→12×0×y²+6Bay+2×0×a²=6Bay=0→B=0（对任意y成立）。由双调和方程A+C+E=0，A=C=E=0，故剩余D=-k/(6b)。此时应力分量：σₓ=0+0+0=0（与x=0、x=a边界条件一致）σᵧ=0+6×(-k/(6b))×yx+0=-kxy/b=-q₀xy/(ab)=-q₀xy/(ab)（y=b时σᵧ=-q₀x/a，符合）τₓᵧ=-003×(-k/(6b))x²=(kx²)/(2b)=q₀x²/(2ab)验证y=0时τₓᵧ=0（成立）；左右边界x=0、x=a时，τₓᵧ=0（x=0时成立，x=a时τₓᵧ=q₀a²/(2ab)=q₀a/(2b)≠0，矛盾）。需修正应力函数，添加φ=Fxy（线性项），则τₓᵧ=-∂²φ/∂x∂y=-F，令F=-q₀a/(2b)，则x=a时τₓᵧ=q₀a/(2b)F=0。最终应力分量：σₓ=0σᵧ=-q₀xy/(ab)τₓᵧ=q₀x²/(2ab)q₀a/(2b)四、综合题（20分）曲梁（内半径r₁=150mm，外半径r₂=250mm，圆心角α=π/2）受纯弯曲作用（两端受弯矩M=10kN·m），材料E=210GPa，ν=0.28。试通过选择合适的应力函数，推导其应力分布，并讨论当曲梁厚度（r₂-r₁）远小于平均曲率半径R=(r₁+r₂)/2时，结果如何简化为材料力学中的弯曲正应力公式σ=My/I。解：曲梁纯弯曲问题属平面应力（或平面应变）问题，采用极坐标，假设应力仅与r、θ有关，且对称于θ=0（α=π/2，两端θ=±α/2）。选择应力函数φ(r,θ)=f(r)cosθ（满足双调和方程∇⁴φ=0，因cosθ为调和函数，f(r)需满足r⁴f''''+2r³f'''r²f''+rf'f=0，其通解为f(r)=Ar+B/r+Cr³+Drlnr）。取φ(r,θ)=(Ar+B/r+Cr³)cosθ（忽略Drlnr项，因纯弯曲无体力），则应力分量：σᵣ=∂²φ/(r²∂θ²)+(1/r)∂φ/∂r=-(Ar+B/r+Cr³)cosθ/r²+(AB/r²+3Cr²)cosθ/r=[-A/rB/r³Cr+A/rB/r³+3Cr]cosθ=(-2B/r³+2Cr)cosθσθ=∂²φ/∂r²=(A+2B/r³+6Cr)cosθτᵣθ=-∂/∂r(∂φ/(r∂θ))=-∂/∂r[(Ar+B/r+Cr³)(-sinθ)/r]=(AB/r²+3Cr²)(sinθ)/r边界条件：（1）内边界r=r₁，外边界r=r₂无面力：σᵣ=0，τᵣθ=0→r=r₁时，-2B/r₁³+2Cr₁=0→B=Cr₁⁴r=r₂时，-2B/r₂³+2Cr₂=0→B=Cr₂⁴联立得C(r₁⁴r₂⁴)=0→C=0（因r₁≠r₂），故B=0。（2）两端θ=±α/2（θ=±π/4）受弯矩M，面力合力为零，合力矩为M：∫_{r₁}^{r₂}σθ(r,±π/4)rdθ=0（合力为零），∫_{r₁}^{r₂}σθ(r,±π/4)r²dθ=M×2（两侧力矩叠加）。由σθ=Acosθ（C=B=0时），代入θ=π/4，cosθ=√2/2，∫_{r₁}^{r₂}A(√2/2)rdr=0→A∫rdr=0（积分非零，故A=0，矛盾）。需修正应力函数为φ(r,θ)=f(r)cosθ，其中f(r)=Ar+B/r（C=0，D=0），则σθ=(A+2B/r³)cosθ。两端面力合力矩：M=∫_{r₁}^{r₂}σθ(r,θ)r²dr（θ方向单位长度，取θ=π/2，cosθ=0，错误，应取θ=0，曲梁两端θ=±α/2，正确积分应为M=∫_{r₁}^{r₂}σθ(r,θ)rsinθdr，θ=α/2时sinθ=sin(π/4)=√2/2）。重新推导：纯弯曲时，曲梁两端面力的合力矩为M，即∫_{r₁}^{r₂}σθ(r,θ)r(rdθ)=M（θ方向长度为rdθ）。对于θ=±α/2，σθ(r,θ)=σθ(r)（与θ无关，因对称），故M=∫_{r₁}^{r₂}σθ(r)r×αdr（α为圆心角弧度）。假设σθ(r)=(A+B/r²)cosθ（修正后），代入边界条件r=r₁、r=r₂时σᵣ=0，得：σᵣ=(-B/r³+Cr)c