简介弹性力学徐芝纶课后习题及答案资料对应徐芝纶版《弹性力学》课程内容主要面向高校土木工程、工程力学、机械设计制造及其自动化等专业本科学生也适用于备战研究生入学考试的考生。资源将教材中的课后习题进行系统梳理并给出详细解答过程方便读者在课后自测、期末复习及考研强化阶段对照使用。资料包内共1个PDF文件压缩包大小约7MB文字公式清晰可按章节顺序查阅答案覆盖教材各章节代表性题目既有推导过程也有结果说明能帮助读者养成规范的分析习惯。针对弹性力学中常见的直角坐标解答、极坐标解答、能量原理等内容答案给出了必要的公式推导和关键步骤便于自查逻辑漏洞。目前已有222人学习浏览适合需要快速查找弹性力学各类典型题解法的学习者。借助其中的答案解析读者可以加深对平面应力、平面应变、应力函数、边界条件等知识点的理解有效提高计算能力与解题速度。无论平时作业还是期末冲刺都能从中获得明确参照。1. 弹性力学课后题不只是考研人的硬骨头《弹性力学徐芝纶》是国内工科院校固体力学、土木、机械类专业的标配教材但很多人翻开它的第一反应是“公式多、推导长、不知道算出来对不对”。这套课后习题与解答的价值恰好在于它把教材里偏理论的推导落成了可验证的数值结果覆盖平面问题、极坐标解、空间问题、变分法这几大模块每道题都给出了完整的应力函数、位移场和边界条件处理过程。对备考研究生入学考试的人它是刷题的主线对做有限元仿真的工程师它是验证仿真模型边界条件是否设对的解析基准。本文不按章节顺序念书而是按工程应用的方式来拆先建立全书的知识坐标再逐个攻克平面问题的应力函数解法和数值验证最后给出期末和考研场景下的刷题策略。2. 全书框架与习题分类先定位再动手2.1 教材章节结构与考查权重徐芝纶教材共分上下两册上册是平面问题的主战场下册进入空间问题和板壳理论。从考研真题和期末试题的分布看上册前三章平面应力与平面应变、直角坐标解、极坐标解占分通常超过60%变分法在考研中常以证明题或近似解形式出现下册的空间轴对称问题出题频率较低但一旦出现就是拉开差距的题。章节与题型对应关系可以按下表归类教材章节核心考点习题类型常见出题方式第一章 基本方程平衡微分方程、几何方程、物理方程推导与概念题写出平面应力问题的基本方程第二章 直角坐标解应力函数取法、边界条件校核计算大题三角形坝、悬臂梁受载第三章 极坐标解轴对称问题、孔口应力集中计算推导圆孔应力集中系数第四章 变分法最小势能原理、里兹法证明近似计算用里兹法求挠度近似解第五章 空间问题空间轴对称、扭转推导题圆杆扭转应力分量验证第六章 板壳理论薄板弯曲方程低频考点四边简支矩形板纳维解这个分类的意义在于不同目标的人精力分配完全不同。考研人要死磕前两章加变分法仿真工程师反而要把前两章当作标准算例库来用因为很多商业有限元软件的验证手册里用的就是这些题的解作为参考解。2.2 解题工具链从手推到数值验证的工作流做弹性力学习题最容易出现的状况是推导过程感觉对但最后算出一个奇怪的应力值却不知道哪里出问题。我的建议是不管是备考还是自学都配一套“手推 符号计算 数值验证”三件套手推理解应力函数怎么设、边界条件怎么列这是核心能力不能跳过符号计算用 Python 的 SymPy 或 Mathematica 快速展开偏导数和验证相容方程数值验证用有限差分或直接写个小脚本把应力分量画出来直观判断边界是否满足比如设应力函数为多项式时需要检验它是否满足双调和方程的问题手动展开四次多项式虽然可行但极慢。常见做法是用 SymPy 做一次符号校验确认应力函数满足相容方程后再把精力放在边界条件的计算上——因为那是考试真正的得分点。3. 平面问题的应力函数解法从设函数到边界条件校核3.1 应力函数选取的一般逻辑平面问题中如果体力为常量可以把应力分量用一个艾里应力函数关联起来。教材里给的是一套结论但习题要的是你会不会“反着用”。核心逻辑如下根据边界受力形式推测应力分量的量纲和变化规律由应力分量反推应力函数中应包含的多项式项代入双调和方程筛选出容许形式用力的边界条件定系数以教材第二章最经典的三角形重力坝问题为例坝体受自重作用两侧面为自由边界。教材的通用做法是将应力函数取为三次多项式组合。这时需要检验每一项是否满足双调和方程比较快的核对方式如下import sympy as sp # 定义应力函数为三次多项式组合 x, y, a, b, c, d sp.symbols(x y a b c d) phi a*x**3 b*x**2*y c*x*y**2 d*y**3 # 检验双调和方程 ∇^4 φ 0 d4phi sp.diff(phi, x, 4) 2*sp.diff(phi, x, 2, y, 2) sp.diff(phi, y, 4) print(sp.simplify(d4phi)) # 结果为0则满足这段代码先构造一般形式的三次多项式应力函数然后计算双调和算子。如果输出为 0说明该函数自动满足相容方程。注意三次多项式恒满足双调和方程次数不超过三时四阶偏导必为零所以重点不是验证它满足而是确认哪些边界条件约束着系数 a、b、c、d。考试中常见的丢分点恰恰是边界条件写不全而不是应力函数设错。边界条件校核的具体做法是将应力函数代入应力分量表达式再把边界面上的面力与已知外力对比建立方程组。例如三角形坝的左侧自由面 x y·tanβ需要把 x 替换为边界面方程联立求解系数。3.2 悬臂梁端部受集中力的真题拆解悬臂梁问题几乎每年考研都会换个数字考一遍但核心从未变过。题目要求矩形截面悬臂梁在自由端受集中力 P求应力分量。标准解法是取应力函数为包含 x·y、x·y³、x³·y 等项的组合。设应力函数为φ A·x·y B·x·y³代入双调和方程后发现需要约束 B 与 A 之间无耦合关系再由边界条件确定系数。这里最常见的困惑是为什么忽略 x³ 相关项——原因在于 x³ 项产生的正应力关于截面是常数分布与梁弯曲应力线性分布矛盾物理上不成立。应力分量的表达式为σx ∂²φ/∂y² 6B·x·yσy ∂²φ/∂x² 0τxy -∂²φ/∂x∂y -A - 3B·y²上端 y h/2 处剪应力为零、下端 y -h/2 处剪应力与外力平衡这两个条件能定出 A 和 B。实际计算中梁端部的应力分布与圣维南原理边界等效有关——端部细节不同只影响局部这正是有限元分析中验证“边界影响区”长度的好素材。import numpy as np import matplotlib.pyplot as plt # 几何和载荷参数 L, h, P 1.0, 0.2, 1000.0 I h**3 / 12 # 单位宽度截面惯性矩 # 网格采样 x_vals np.linspace(0, L, 50) y_vals np.linspace(-h/2, h/2, 50) X, Y np.meshgrid(x_vals, y_vals) # 应力分量教材解 sigma_x -P * (L - X) * Y / I tau_xy -P / (2 * I) * (h**2 / 4 - Y**2) # 绘制正应力云图 plt.figure(figsize(8, 4)) plt.contourf(X, Y, sigma_x, levels30, cmapRdBk_r) plt.colorbar(labelsigma_x (Pa)) plt.xlabel(x (m)) plt.ylabel(y (m)) plt.title(悬臂梁正应力分布) plt.axis(equal) plt.show()这段代码实现了悬臂梁端部受集中力时教材解的应力云图绘制。注意公式里用的是截面惯性矩 I对单位宽度梁取 h³/12。参数 X 和 Y 是网格坐标矩阵np.meshgrid 负责生成二维采样点sigma_x 表达式里 (L - X) 表示弯矩随 x 线性变化正应力的线性分布直观可见。运行后能明显看到固定端应力最大、自由端为零符合材料力学预期。3.3 极坐标系下的孔口应力集中问题极坐标解是第二道分水岭。圆孔应力集中问题的结论——无限大平板上圆孔边缘的最大应力是无孔时的三倍——几乎成了弹性力学的“品牌结论”。习题通常要求从拉梅解出发用边界条件定出系数。典型的题目无穷大平板受单向均匀拉伸半径为 a 的圆孔求孔边应力。解的过程是先写出应力函数的一般形式然后考虑无穷远处应力趋于均匀场孔边为自由边界。这里的关键步骤往往不是求偏导而是在极坐标下正确列出两个边界条件r a 处σr 0τrθ 0r → ∞ 处σx qσy 0τxy 0坐标转换时要把直角坐标的均匀应力场换算为极坐标分量这一步出错率极高。常见的技巧是先写出直角坐标下的应力张量再用转换公式变到极坐标而不是直接猜极坐标下的表达式。import sympy as sp r, a, q, theta sp.symbols(r a q theta, positiveTrue) # 教材解孔边环向应力表达式 sigma_theta q / 2 * (1 a**2 / r**2) - q / 2 * (1 3 * a**4 / r**4) * sp.cos(4 * theta) sigma_theta_at_r sp.simplify(sigma_theta.subs(r, a)) # 找到最大最小值 theta_vals [0, sp.pi/4, sp.pi/2, 3*sp.pi/4] for t in theta_vals: val sp.simplify(sigma_theta_at_r.subs(theta, t)) print(ftheta {sp.simplify(t)}: sigma_theta {val})这段代码把孔边 (r a) 处的环向应力计算出来取四个关键角度观察变化。输出结果会显示 θ 0 和 θ π 处应力为 -q压缩θ π/2 处为 3q拉伸θ 3π/4 处回到 -q。注意负值表示压应力三倍拉伸应力出现在与加载方向垂直的孔边位置与教材结论一致。如果你在有限元软件里建模时忘记细化孔边网格这个三倍应力是算不准的。4. 变分法与近似解当解析解推不动时怎么办4.1 最小势能原理与里兹法的实战操作教材后半部分的重要转折点是变分法很多实际问题找不到精确的应力函数转而用能量法求近似解。里兹法的思路是假设位移场为若干基函数的线性组合通过最小化总势能确定组合系数。以四边简支矩形薄板受均布载荷为例假设挠度场为二重三角级数形式。这个题目在教材中给了纳维解但题目经常要求用里兹法取少量项做近似。实操时的关键步骤设位移函数 w(x,y) ΣΣ Amn · sin(mπx/a) · sin(nπy/b)写出板的弯曲应变能和外力势能对每个 Amn 求偏导并令其为零解线性方程组得各系数import sympy as sp # 薄板几何与材料参数 a, b, q0, D sp.symbols(a b q0 D, positiveTrue) m, n sp.symbols(m n, positiveTrue, integerTrue) # 挠度假设单一项 x, y sp.symbols(x y) w sp.sin(m*sp.pi*x/a) * sp.sin(n*sp.pi*y/b) # 曲率表达式 w_xx sp.diff(w, x, 2) w_yy sp.diff(w, y, 2) w_xy sp.diff(sp.diff(w, x), y) # 应变能密度积分核 strain_energy D/2 * ((w_xx w_yy)**2 - 2*(1-sp.symbols(nu))*(w_xx*w_yy - w_xy**2)) U sp.integrate(sp.integrate(strain_energy, (x, 0, a)), (y, 0, b)) # 外力功 W sp.integrate(sp.integrate(q0*w, (x, 0, a)), (y, 0, b)) # 总势能 Pi U - W print(sp.simplify(Pi))这段代码完成了单一项里兹法的势能泛函推导。变量 a、b 为板边长q0 为均布载荷集度D 为弯曲刚度。代码逐项计算了薄板应变能中曲率项的组合再与外载荷做功相减得到总势能表达式。对 Amn 求导后就能得到挠度系数的显式解。注意代码中用了符号 nu 表示泊松比应变能公式包含了扭矩项这是薄板理论区别于简单梁理论的地方。4.2 里兹法与有限元法的边界条件差异不少做过有限元的同学会把里兹法当作“手算有限元”这个类比在思想上成立但有重要差异里兹法的基函数是全域定义的必须满足位移边界条件有限元的形函数是局部的逐片满足连续性要求。这意味着里兹法很难处理复杂几何但作为考试题几乎只需要套流程。学习这个章节时建议用同一个题目同时跑里兹法和有限元软件比较中心挠度。取一阶项时里兹法的相对误差通常在 10% 到 20% 之间取三阶以上可收敛到纳维解附近。这个对照过程能让你直观理解“收敛性”的含义而不是只把它当成考试公式背。5. 期末与考研冲刺按题型刷题的进阶技巧5.1 高频公式速查与易混点辨析进入冲刺阶段建议先做一次“公式整理”——把教材的公式用自己的逻辑重构一遍而不是照着目录背。下面是一份浓缩考点表可直接用于最后三天的快速翻阅考点模块核心公式/结论每章必背结论常见误区判断平面应力与应变σz区分平面应力 σz0平面应变 εz0混淆两者物理方程应力函数双调和方程多项式次数≤3恒满足忘记检查高阶项系数极坐标拉梅解孔边应力集中系数3θπ/2处取最大值变分法最小势能原理里兹法设函数满足位移边界忘记检查几何边界条件空间轴对称拉梅位移势圆筒受内外压解把轴对称当平面问题处理易错的关键点有三个一是平面应力和平面应变的弹性矩阵差一个因子导致结果完全不同二是应力函数法求出的应力分量必须满足所有的边界条件只满足平衡方程是不够的三是极坐标下的剪应力符号约定容易在应力转换时搞混建议每次都用张量转换公式验证一遍方向。5.2 真题-仿真对照验证法对已经在做仿真工作的读者我建议一个进阶玩法把徐芝纶教材里的经典题目当作 benchmark 用例反向验证有限元模型。悬臂梁端部受载模型是验证边界条件最灵敏的算例——如果固定端约束方式不对固定端附近应力分布就会偏离理论解孔口应力集中题则是验证网格收敛性的好工具。操作步骤很简单先用教材理论解算出关键位置的应力值然后在仿真软件中建同尺寸模型加密网格后对比。偏差超过 5% 时需要怀疑是网格问题、边界条件问题还是理论解的适用范围被超出。这个做法比盲目相信软件默认设置可靠得多。矩形板四边简支受均布载荷则是验证板壳单元的好题——注意理论解基于薄板假设长厚比低于 10 时误差会显著增大。5.3 最后阶段的刷题顺序冲刺期不建议从头到尾过题。有一个信息密度比较高的顺序先刷第三章极坐标题目因为这类题公式长但套路化程度最高属于最容易拿分的大题再刷第二章应力函数题这是计算量最大的部分最后用第四章变分法的推导题收尾因为变分法的题既考计算又考证明逻辑适合保持手感。每做完一道题在答案旁边标注“这题在考哪个边界条件”或“如果载荷方向变了哪里要改”。这种标注会让复习效率提高不少因为考试变题通常只改边界条件或载荷形式核心的应力函数结构和推导流程是不变的。本文还有配套的精品资源点击获取