hzf第2章-平面问题的基本理论

弹性力学第2章平面问题的基本理论Contents本章内容框架弹性力学第二章——平面问题的基本理论01平面问题的基本概念与分类02平面问题的基本方程体系03边界条件与圣维南原理04平面问题的基本解法05常体力下的应力函数解法CHAPTER01平面问题的基本概念与分类从几何特征与受力特征理解平面应力问题和平面应变问题的本质区别Chapter02·Overview平面问题的研究意义与课程定位弹性力学平面问题是二维问题，能够以最简洁的形式阐明弹性力学的核心概念与理论框架。作为从三维到二维的关键简化，平面问题既是理论学习的起点，也是工程中应用最广泛的弹性力学模型之一。二维简化模型平面问题是弹性力学中最基本的二维简化模型，通过降维处理将复杂三维问题转化为可解析求解的形式。降维求解三大方程体系本章系统建立的平衡方程、几何方程和物理方程三大体系，构成弹性力学理论的核心骨架。核心骨架工程应用广泛薄板弯曲、重力坝截面分析、隧道围岩应力等大量实际问题均可简化为平面问题求解。薄板·坝体·隧道方法论基础逆解法、半逆解法、应力函数法等求解方法为后续空间问题和特殊问题奠定方法论基础。逆解法·应力函数Chapter2·平面问题的基本理论平面应力问题：几何特征与力学条件平面应力问题对应z方向几何尺寸远小于x、y方向的等厚度薄板类构件。由于板面为自由表面，沿厚度方向的应力分量σz、τzx、τzy均为零，只剩下σx、σy、τxy三个平面内应力分量。01几何特征：物体沿z方向的厚度尺寸远小于x、y方向的平面尺寸，典型构件为等厚度薄板02受力特征：板的两表面为自由面，不承受任何面力，体力与面力均平行于板面且沿厚度均匀分布03应力简化：σz=τzx=τzy=0，仅剩σx、σy、τxy三个非零应力分量，且它们仅为x、y的函数04注意：虽然σz=0，但εz一般不为零——薄板在面内受力后厚度方向会产生泊松效应引起的变形等厚度薄板构件·平面应力典型模型PLANESTRAINPROBLEM平面应变问题：几何特征与力学条件平面应变问题对应z方向尺寸远大于x、y方向的常截面长柱体。由于纵向约束，所有横截面保持相同的变形状态，εz=γzx=γzy=0，应力和应变仅为x、y的函数，三维问题退化为二维分析。典型平面应变构件——重力坝截面几何特征：物体沿z方向的尺寸远大于x、y方向的截面尺寸，典型构件为等截面长柱体如水坝、隧道受力特征：体力与面力均垂直于z轴且分布规律不随z变化，位移约束和支承条件沿z方向完全相同应变简化：εz=γzx=γzy=0，仅剩εx、εy、γxy三个非零应变分量，且仅为x、y的函数注意：虽然εz=0，但σz一般不为零——σz=μ(σx+σy)，这是纵向约束引起的被动应力ELASTICITY·CH.2两类平面问题的对比与弹性常数等价替换两类平面问题的平衡微分方程、几何方程和边界条件完全相同，仅物理方程系数不同。通过将平面应力问题物理方程中的弹性常数进行等价替换（E→E/(1-μ²)，μ→μ/(1-μ)），即可直接获得平面应变问题的解答。共同特征01平衡微分方程、几何方程完全相同，均不含材料参数，反映的是力学和几何的普遍规律02应力边界条件和位移边界条件的表达形式完全一致，仅边界上的已知量取值不同差异与转换01物理方程是唯一差异所在：平面应力中σz=0而εz≠0，平面应变中εz=0而σz≠002弹性常数替换关系：将E换为E/(1-μ²)、μ换为μ/(1-μ)，平面应力解答即转化为平面应变解答CHAPTER02·PLANEPROBLEM平面问题的未知函数体系平面问题共有8个未知函数需要求解：3个应力分量、3个应变分量和2个位移分量，通过三大基本方程和边界条件相互联系，构成完整数学模型。01应力分量σx、σy、τxy共3个独立分量，描述弹性体任意一点的受力状态3个02应变分量εx、εy、γxy共3个独立分量，描述弹性体任意一点的变形状态3个03位移分量u、v共2个独立分量，表示各点沿x、y方向的位移，最基本的运动学变量2个04方程体系2个平衡微分方程+3个几何方程+3个物理方程=8个独立方程8=8CHAPTER02平面问题的基本方程体系从静力学、几何学和物理学三个维度建立弹性力学平面问题的数学模型EQUILIBRIUMEQUATIONS平衡微分方程的推导与物理含义平衡微分方程描述弹性体内任意微元体在体力作用下必须满足的静力平衡条件。通过在dxdy微元体上列写x、y两方向的合力为零方程，并令高阶小量趋于零，得到联系应力分量与体力的两个偏微分方程。01取微元体dxdy×1（单位厚度），各面应力按Taylor展开保留一阶项σx+(∂σx/∂x)·dx02x方向合力为零，化简得第一平衡方程∂σx/∂x+∂τyx/∂y+fx=003y方向合力为零，化简得第二平衡方程∂σy/∂y+∂τxy/∂x+fy=004力矩平衡给出剪应力互等，独立应力分量由4个减为3个τxy=τyxChapter2·平面问题基本理论几何方程的推导与变形分析几何方程从纯运动学角度描述位移与应变的关系，不涉及材料属性。通过分析微元体变形前后各边的伸长和角度变化，建立3个应变分量与2个位移分量之间的微分关系，是连接位移场与应变场的桥梁。线应变定义εx=∂u/∂x（x方向位移沿x方向的变化率），εy=∂v/∂y（y方向位移沿y方向的变化率）εx=∂u/∂xεy=∂v/∂y剪应变定义γxy=∂u/∂y+∂v/∂x，表示变形前后微元体直角的变化量，由两部分剪切变形叠加而成γxy=∂u/∂y+∂v/∂x变形协调条件3个几何方程联系3个应变分量与2个位移分量，方程数多于未知数，应变分量不能完全独立——必须满足变形协调条件3方程·2未知位移普适性特征几何方程仅反映变形的几何关系，与材料弹性常数E、μ无关，对所有连续介质均成立∀连续介质Chapter02平面问题基本理论物理方程：广义胡克定律的平面简化物理方程将应力分量与应变分量通过材料弹性常数（E、μ）联系起来，是广义胡克定律在平面条件下的简化形式。平面应力与平面应变的物理方程形式相似但系数不同，反映了材料力学行为对约束条件的敏感性。NormalStrain正应变耦合效应εx=(σx−μσy)/E，εy=(σy−μσx)/E，反映面内应力对正应变的耦合效应ShearRelation剪应变与剪应力γxy=2(1+μ)τxy/E=τxy/G，剪应变与剪应力成正比，比例系数为剪切模量G的倒数MathematicalUnity数学统一性将E替换为E/(1−μ²)、μ替换为μ/(1−μ)即可由平面应力方程获得，体现了两类问题的数学统一性StiffnessIncrease等效刚度提升εx=(1−μ²)/E·[σx−μ/(1−μ)·σy]，纵向约束使面内变形刚度增大，等效弹性模量提高平面应力物理方程平面应变物理方程Chapter2·FundamentalEquations三大基本方程体系的完整梳理弹性力学平面问题的三大基本方程分别从静力学、几何学和物理学三个角度约束8个未知函数，构成封闭的偏微分方程组，是求解平面问题的核心数学模型。平衡微分方程2个∂σx/∂x+∂τxy/∂y+fx=0∂σy/∂y+∂τxy/∂x+fy=0静力学条件几何方程3个εx=∂u/∂x，εy=∂v/∂yγxy=∂u/∂y+∂v/∂x运动学条件物理方程3个由广义胡克定律简化得到，包含材料常数E和μ，是联系应力与应变的本构纽带。本构条件方程组封闭8=88个方程对应8个未知函数（σx、σy、τxy、εx、εy、γxy、u、v），弹性体内部封闭，求解仍需边界条件。边界条件Chapter2·BasicTheoryofPlaneProblems一点的应力状态与主应力分析弹性体内一点的应力状态由σx、σy、τxy三个分量完全确定。通过应力张量的坐标变换，可以求出过该点任意方向截面上的应力分量；而主应力和主方向则由应力张量特征方程给出，反映了该点应力状态的本质特征。01斜截面应力公式已知σx、σy、τxy，通过方向余弦(l,m)可求任意斜面上的正应力σn和剪应力τn，建立应力分量与截面方位的完整映射关系。02主应力特征方程λ²−(σx+σy)λ+(σxσy−τxy²)=0，解得两个主应力σ₁和σ₂，对应剪应力为零的特殊截面方向。03应力不变量σx+σy=σ₁+σ₂（第一不变量），σxσy−τxy²=σ₁σ₂（第二不变量），不随坐标系旋转而改变，是应力状态的内在属性。04最大剪应力τmax=(σ₁−σ₂)/2，发生在与主应力方向成45°角的截面上，是强度分析和失效判据的关键参数。CHAPTER03边界条件与圣维南原理理解弹性体边界上的力学约束及其简化处理方法CHAPTER2·平面问题的基本理论三类边界条件的定义与数学表达弹性力学边界条件分为应力边界条件、位移边界条件和混合边界条件三种类型。边界条件反映弹性体的个性——不同的约束方式导致完全不同的应力和位移分布，是区分数学模型通解与特解的关键。应力边界条件TRACTIONB.C.l·σx+m·τxy=f̄xm·σy+l·τxy=f̄y在边界Sσ上给定面力(f̄x,f̄y)，(l,m)为边界外法线方向余弦，保证边界上内力与外力平衡。Sσ位移边界条件DISPLACEMENTB.C.u=ūv=v̄在边界Su上给定位移(ū,v̄)。固定端典型情形：ū=0、v̄=0；铰支座或滑动支座则部分分量为零。Su混合边界条件MIXEDB.C.Sσ∪Su=S边界一部分给定面力、另一部分给定位移。同一边界上也可部分方向给面力、部分方向给位移，如简支梁支座处。Sσ+SuBOUNDARYCONDITIONS应力边界条件的推导与典型情形应力边界条件通过在边界上取微元三角形并列力平衡方程推导得到。方向余弦(l,m)取边界外法线方向，对于平行于坐标轴的直线边界，可大幅简化。01推导方法在边界上取微元三角形体元，斜面上作用应力(σx,σy,τxy)，外法线方向余弦为(l,m)，对微元列x、y方向力平衡方程。02一般形式l·σx+m·τxy=f̄x，m·σy+l·τxy=f̄y，其中f̄x、f̄y为边界上给定的面力分量。03左右边界右边界(x=const,l=1,m=0)：σx=f̄x，τxy=f̄y；左边界(l=−1,m=0)：−σx=f̄x，−τxy=f̄y。04上下边界上边界(y=const,l=0,m=1)：τxy=f̄x，σy=f̄y；下边界(l=0,m=−1)：−τxy=f̄x，−σy=f̄y。Chapter2·BasicTheory圣维南原理：基本概念与核心思想圣维南原理指出：在弹性体的一小部分边界上，将面力替换为静力等效的力系（合力和合力矩相同），仅在力的作用区域附近应力分布发生显著改变，远离该区域处的应力分布几乎不受影响。该原理为简化复杂边界条件提供了理论依据。核心思想局部边界上面力的具体分布方式只影响附近的应力状态，对远处的应力分布影响可忽略。局部效应静力等效条件替换前后的力系必须具有相同的合力（两个方向的分力）和合力矩，共三个等效力条件。3个等效条件影响范围应力差异主要局限在力作用区域附近约一个特征尺寸范围内，超出此范围差异迅速衰减。≈1倍特征尺寸工程意义当边界上面力的精确分布未知时，可用静力等效的简单分布替代，大幅简化求解过程。简化边界条件APPLICATIONMETHODS圣维南原理的工程应用方法圣维南原理在工程中有两种典型应用方式：静力等效法（用简单分布替代复杂面力）和静力平衡法（通过截面合力条件建立边界约束）。这两种方法在梁端受力、杆件连接等大量工程问题中被广泛采用，是弹性力学分析中最重要的简化手段之一。简支梁弯曲实验·材料力学实验室01静力等效法在局部边界上用均匀分布或线性分布的简单应力假设替代实际复杂面力，仅需保证合力与合力矩等效02静力平衡法直接在边界截面上列写合力平衡方程∫σdA=N、∫σy·dA=M等，不关心应力的详细分布03梁端部应用精确端部面力分布难以确定，用圣维南原理将端部力系简化为轴力N、剪力Q和弯矩M三个合力分量04应用注意事项仅适用于"一小部分边界"，不能用于大面积边界或整体边界的替换；替代区域附近的解不可信BOUNDARYCONDITIONS边界条件的本质与求解意义基本方程描述弹性体的共性——同一类问题方程相同；边界条件体现弹性体的个性——每个具体问题约束不同。只有二者结合才能获得唯一解。圣维南原理通过允许局部边界的静力等效替换，在保证远场精度的前提下大幅降低了求解难度。必要条件：基本方程是弹性力学问题的必要条件——所有满足基本方程的解构成通解族，但通解包含待定函数，不能直接用于工程充分条件：边界条件是充分条件的关键——通过边界条件确定通解中的待定函数或待定常数，从而获得满足具体约束的特解桥梁作用：圣维南原理的桥梁作用——当精确边界条件难以建立或求解过于复杂时，用等效边界条件替代可获得足够精度的远场解完整定义：区域Ω内的基本方程+边界∂Ω上的约束条件+材料的弹性常数(E,μ)，三者缺一不可Chapter04平面问题的基本解法位移解法与应力解法——两条经典的弹性力学求解路径DISPLACEMENTMETHOD位移解法：以位移为基本未知函数位移解法将8个未知函数缩减为2个位移分量(u,v)，通过几何方程代入物理方程再代入平衡微分方程，得到Lamé-Navier方程求解。消元路径先将位移(u,v)代入几何方程得应变，再将应变代入物理方程得用位移表示的应力表达式8→2Lamé-Navier方程将应力表达式代入平衡微分方程，得到(λ+G)∂θ/∂x+G∇²u+fx=0等两个方程，其中θ=∂u/∂x+∂v/∂y为体积应变2个PDE边界条件处理位移边界条件下直接给定u和v，天然适配；应力边界条件需转化为用位移表示的形式，处理较复杂天然适配求解步骤解Navier方程得位移→代入几何方程得应变→代入物理方程得应力，三步顺序完成全部求解3步求解STRESSAPPROACH应力解法：以应力为基本未知函数应力解法以3个应力分量为基本未知函数，但平衡微分方程仅有2个，方程不足。需引入应变协调方程（变形连续条件）作为补充——将物理方程代入几何方程消去位移，得到应变之间的相容性约束，再转化为用应力表示的协调方程，从而使方程组封闭可解。013:2方程不足问题：3个应力分量需要3个方程，但平衡微分方程仅2个，必须从几何方程中再提取1个补充方程02ε→γ应变协调方程：∂²εx/∂y²+∂²εy/∂x²=∂²γxy/∂x∂y——这是3个几何方程消去2个位移后的唯一独立相容条件03∇²σ应力协调方程：将物理方程代入应变协调方程，得到用应力表示的∇²(σx+σy)=-(1+μ)(∂fx/∂x+∂fy/∂y)（平面应力情形）04∇²=0常体力简化：当体力fx、fy为常数时，协调方程简化为∇²(σx+σy)=0，即应力第一不变量满足调和方程COMPARISON位移解法与应力解法的系统对比位移解法方程少但边界条件处理复杂，应力解法方程多但边界条件直接适配，工程中应力解法应用更广泛。位移解法22个方程基本未知量少（仅u、v两个），方程为Lamé-Navier方程（2个二阶偏微分方程）。由于方程数量少，理论推导相对简洁，但求解过程中需要处理复杂的边界条件转换。天然适配位移边界条件，但应力边界条件需转化为位移形式，数学处理较复杂。这种转换涉及微分运算和本构关系的联立，增加了求解难度。适用场景：位移边界条件明确的结构分析、有限元法中的位移元formulation。应力解法33个方程基本未知量为σx、σy、τxy（3个），需2个平衡方程加1个协调方程联立求解。虽然方程数量增加，但物理意义更加直观，便于工程理解。天然适配应力边界条件，工程中大多数问题的边界条件以面力给定为主，故应力解法更常用。可直接将已知的面力边界条件代入求解，无需复杂转换。适用场景：应力边界条件主导的弹性力学问题、平面应力/应变问题、应力集中分析。INSIGHT应力边界条件在工程中更为常见，应力解法因此成为主流选择弹性力学·平面问题基本理论从应力解法到应力函数解法的过渡应力解法需要同时求解3个偏微分方程，Airy应力函数法通过引入标量函数φ(x,y)，将问题缩减为求解双调和方程∇⁴φ=0，大幅降低求解难度。01应力解法的难点3个应力分量需满足2个平衡方程和1个协调方程，方程耦合，直接求解困难02Airy应力函数的引入构造函数φ(x,y)，使由φ导出的应力分量自动满足平衡微分方程03简化效果将3个未知函数缩减为1个应力函数，将3个方程缩减为1个双调和方程∇⁴φ=004应用前提常体力情况（如重力），工程中大量问题体力为常量；非常体力情况需另行处理Chapter05常体力下的应力函数解法Airy应力函数与双调和方程——平面问题最核心的解析求解工具第二章·平面问题基本理论Airy应力函数的定义与应力分量关系Airy应力函数φ(x,y)通过其二阶偏导数定义应力分量：σx=∂²φ/∂y²、σy=∂²φ/∂x²、τxy=-∂²φ/∂x∂y。这种定义方式使平衡微分方程在无体力或常体力时自动满足，从而将求解问题转化为在边界条件下求解φ的协调方程。01应力分量定义σx=∂²φ/∂y²，σy=∂²φ/∂x²，τxy=-∂²φ/∂x∂y，三个应力分量由φ的二阶偏导数确定。∂²φ/∂y²02平衡方程自满足将应力分量代入∂σx/∂x+∂τxy/∂y=0，得∂³φ/∂x∂y²-∂³φ/∂x∂y²=0，恒成立；y方向同理。恒成立03常体力修正体力为常量时，令σx'=σx-fx·x等，将体力吸收到应力函数的修正中，平衡方程仍可自动满足。σx'修正04物理意义诠释φ本身无直接物理意义，但其等值线可反映应力分布特征，是一种数学构造工具。等值线弹性力学·平面问题基本理论双调和方程：应力函数的协调条件将Airy应力函数代入协调方程，得到∇⁴φ=0，满足此方程的任何双调和函数均可产生满足全部基本方程的应力场。01推导过程σx+σy=∂²φ/∂y²+∂²φ/∂x²=∇²φ，代入协调方程∇²(σx+σy)=0得∇²(∇²φ)=0，即∇⁴φ=0。这一推导建立了应力函数与协调条件之间的数学联系。02展开形式∂⁴φ/∂x⁴+2·∂⁴φ/(∂x²·∂y²)+∂⁴φ/∂y⁴=0，是四阶线性齐次偏微分方程。该方程在直角坐标系下具有对称性，反映了平面问题的内在几何约束。03方程性质双调和方程的通解不能写成有限项的显式形式，一般需要采用逆解法或半逆解法求解。根据边界条件假设应力函数形式，再验证是否满足方程。04物理解读∇⁴φ=0本质上保证了由φ导出的应变场满足变形协调条件，即变形后的弹性体仍然连续无裂缝。这是弹性力学解存在的必要条件。SOLUTIONMETHODS应力函数解法的完整求解流程应力函数解法将弹性力学平面问题转化为寻找满足双调和方程∇⁴φ=0的函数φ，并在边界条件下确定待定参数。由于双调和方程通解无法写成有限项形式，工程实践中主要采用逆解法（假设φ→验证方程→匹配边界）和半逆解法（由力学直觉假设应力→反推φ）两种策略。01逆解法流程假设φ的形式→验证∇⁴φ=0→由φ求σx、σy、τxy→检查边界条件→确定待定常数→获得精确解φ→∇⁴φ=0→精确解02半逆解法流程根据力学直觉假设部分应力分量→反推φ的函数形式→代入双调和方程确定完整φ→验证全部条件直觉假设→反推φ03关键技巧多项式函数天然是双调和函数（≤3次的多项式自动满足∇⁴φ=0），是最常用的应力函数试探形式≤3次多项式≡双调和04求解难点φ的假设需要经验和物理直觉，不同问题需要不同的函数形式（多项式、三角函数、极坐标形式等）经验与物理直觉POLYNOMIALSTRESSFUN