弹性力学作业总结

1、 一、综述这学期我们有幸跟着邱老师学习了弹性力学这门课程，虽然我本科是学习机械专业的，但经过这学期的系统学习，使我对弹性力学的认识也越发的清晰，我对平面问题、空间问题等基本知识有了较为清晰的了解与掌握，会用逆解法、半逆解法、差分法、变分法和有限元法解决一些基础的弹性力学问题。弹性力学是固体力学的一个分支，研究弹性体由于外力作用或温度改变等原因而发生的应力、形变和位移。它是学习塑性力学、断裂力学、有限元方法的基础，广泛应用于建筑、机械、化工、航天等工程领域。本课程较为完整的表现了力学问题的数学建模过程，建立了弹性力学的基本方程和边值条件，并对一些问题进行了求解。弹性力学基本方程的建立为进一步的数

2、值方法奠定了基础。2、 绪论 弹性力学所依据的基本规律有三个：变形连续规律、应力-应变关系和运动(或平衡)规律，它们有时被称为弹性力学三大基本规律。弹性力学中许多定理、公式和结论等，都可以从三大基本规律推导出来。通过对弹性力学的学习，我感觉整本书就讲了十五个控制方程解十五个未知数。而剩下的问题就是如何求解这些方程的问题，这也是数学和力学结合最紧密的地方。而求解的方法无外乎有:基于位移的求解（位移法）和基于应力的求解（应力函数法），差分法、变分法。而前人的研究大部分都是如何使这些方程求解起来更方便。弹性力学思路清晰，但是方程和公式复杂。1.工程力学问题建立力学模型的过程，一般要对三方面进行简化：

3、结构简化、材料简化及受力简化。建模过程如右图：结构简化：如空间问题向平面问题的简化，向轴对称问题的简化，实体结构向板、壳结构的简化。受力简化：根据圣维南原理，复杂力系简化为等效力系。 材料简化：根据各向同性、连续、均匀等假设进行简化。 在建立数学模型的过程中，通常要注意分清问题的性质进行简化：线性化和实验验证。2弹性力学的基本内容就是研究研究弹性体由于外力作用或温度改变等原因而发生的应力、形变和位移。应用在工程中的实例有比赛斜塔，水轮机以及各种齿轮等等。但是我们需要对工程实例中的问题进行假设，具体如下：（1）连续性假设：这样物体内的一些物理量，例如应力、应变和位移等可用连续函数表示。（2）线弹

4、性假设：假定物体服从胡克定律。（3）均匀性假设：假定物体由同一材料组成，这样材料常数不随位置坐标变化。（4）各向同性假设：物体内一点的弹性性质在各个方向上相同。（5）小变形假设：假定位移和应变是微小的。这样，可以用变形前的尺寸代替变形后的尺寸，在考察物体的应变和位移时，可以略去高阶小量，这对于方程的线性化十分重要。在上述简化的基础上，我们现在的问题是要得出普遍的描述上述问题的力学和数学模型。3.几个基本概念：外力：分为体力和表面力。 体力：是分布在物体体积内的力。物体内个点受体力的情况一般不 相同。 物体在一点的体力的集度：。 面力：是分布在物体表面上的力，物体在其表面上各点受面力的情况一般也

5、不相同。物体在一点的面力的集度：。 方向：沿坐标轴正向则为正。应力：包括正应力和切应力。 方向：正面正向为正，负面负向为正。（正面：截面外法线方向为正方向的截面） 注意下标字母的含义：正应力是作用在垂直于x轴的面上，同时也是沿着x轴的方向作用的。切应力有两个下标字母，前一个表示作用面垂直于哪个坐标轴，后一个字母表示作用方向沿哪个轴。 切应力互等性：作用在两个互相垂直的面上并且垂直于该两面交线的切应力是大小相等，正负相同的。（区别于材料力学中的正负号相反）形变：包括线应变和切应变。 线应变：各线段的每单位长度的伸缩。（，）正负判定：伸长为正。 切应变：各线段之间的直角的改变。（，）正负判定：直角

6、变小为正。位移：,（三个坐标的投影表示）。3 平面问题的基本理论3.1平衡微分方程 微元体尺寸：dx, dy, 1。 推导利用：泰勒公式： ；。 求得平面应力表示的平衡微分方程： ； 。说明：1.平衡方程仅反映物体内部的平衡，当应力分量满足平衡方程，则物体内部是平衡的。 2.平衡方程也反映了应力分量与体力的关系。3.2几何方程P点的位移为（u,v）。 PA的线应变为。 PB的线应变为。PA与PB之间的直角改变为。推得平面问题的几何方程： ；。说明：1.几何方程反映了位移和应变之间的关系。 2.当位移完全确定时，应变也确定；反之，当应变确定时，位移不一定确定。刚体位移：与形变无关的位移。由，得到

7、，。说明：当物体发生一定的形变时，由于约束条件的不同，它可能有不同的刚体位移，因而它的位移并不完全确定。必须有适当的刚体约束条件来确定其中不确定的常数。3.3物理方程 由广义胡克定律公式推导：在平面应力问题中，。代入胡克定律，得到物理方程：，。在平面应变问题中，。其物理方程为平面应力情况下，将方程中的，进行转换即可。3.4平面问题的边界条件应力边界条件位移边界条件混合边界条件在应力边界S上在位移边界Su上 在S上 在Su上3.5圣维南原理如果将物体的一小部分边界上的面力，变换为分布不同但静力等效的面力（主矢、主矩相同），那么，近处的应力分布将有显著的改变，但是远处所受的影响可以不计。(1) 利

8、用静力等效力系(2)局部边界用近似边界条件3.6按位移求解平面问题位移法：取位移分量为基本未知函数，从基本方程和应力边界条件中消去应力和应变分量，导出用位移表示的平衡微分方程和用位移表示的应力边界条件。并由此求出位移分量，通过几何方程求出应变分量，再由物理方程求出应力分量。按位移求解平面为题时，位移分量（u，v）必须满足以下全部条件：（1）用位移表示的平衡微分方程（平面应力问题）（2）用位移表示的应力边界条件（平面应力问题）（3）位移边界条件 对于平面应变问题，只需将（1）（2）的式子中的、分别换成和。3.7按应力求解平面问题 应力法：它是以应力分量为基本未知函数，从方程和边界条件中消去位移分

9、量和形变分量，导出只含应力分量的方程和边界条件，由此解出应力分量，然后再求形变分量和位移分量。关键问题：如何得到方程组。解决思路：平衡方程中有两个方程，三个未知数，需引入一个补充方程： 相容方程：，说明了连续体的形变分量，不 是相互独立的，而是相关的，它们之间必须满足相容方程才能保证对应的位移分量的存在。 将物理方程代入相容方程中，再利用平衡方程化简，就得到补充方程： 注意：分清所用的物理方程是关于平面应力问题还是平面应变问题的，其他方程及条件都是相同的，平面应力与应变问题的转换：，。3.8常体力情况下的简化及应力函数在体力为常量的情况下，平衡微分方程、相容方程和应力边界条件都不包含弹性常数，

10、从而对于两种平面问题都是相同的。因此，当体力为常量时，在单连体的应力边界问题中，如果两个弹性体具有相同的边界形状，并受到同样分布的外力，那么不管这两个弹性体材料是否相同和处于那种平面问题状态，应力分量， ，的分布是相同的（以及形变和位移却不一定相同）。求解体力为常数时的平衡微分方程： ， 。注意求通解的时候要利用性质：。可得到通解形如：，。其中称为应力函数。将上述通解代入相容方程，得到：，即，这就是用应力函数表示的相容方程。 （注：代表）。四、平面问题的直角坐标解答4.1逆解法与半逆解法 逆解法的基本思路：（1）设定各种形式的应力函数 ，要求：满足相容方程（2）由下式求得应力分量（3）由应力边

11、界条件式和弹性体的边界形状找到应力分量对应的面力， 从而得知所选取的应力函数 可以解决的问题。半逆解法的基本思路：（1）针对所要求解的问题，根据边界形状和受力情况，假设部分或全部应力分量的函数形式；（2）推出应力函数的形式；（3）代入相容方程，求出应力函数的具体表达形式；（4）再按相容方程式由应力函数求得应力分量；（5）考查应力分量是否满足全部边界条件（多连体还要满足位移单值）；（6）满足是问题的解，不满足重新假设求解。注（1）线性应力函数对应于无体力、无面力、无应力的状态；（2）把平面问题的应力函数加上一个线性函数，并不影响应力。4.2矩形梁的纯弯曲如右图所示，建立应力函数（1）上下两个主要

12、边界：（2）左右两个次要边界切应力：（3）左右两个次要边界正应力（引入圣维南原理）：因此应力分量为：小结：（1）弹性力学的解和材料力学的解完全一致；（2）只有在纯弯曲的情况下，解才是完全精确的；（3）两端的面力按其它分布时，根据圣维南原理，在梁两端附近时 有显著误差，在离开梁较远处，误差可以忽略不计。4.3位移分量的求出（1）把应力分量带入物理方程求得应变分量；（2）把应变分量带入几何方程得到用应变分量表示的位移方程组；其中： 是表示刚体位移的常量，可由位移约束条件求出。小结：（1）弹性力学验证了材料力学的平截面假定是正确的；铅直线段的转角（2）梁纵向纤维的曲率 和材料力学完全相同；（3）相对

13、于材料力学而言，弹性力学的位移分量更全面的表示了梁的变形情况；（4）对于平面应变情况，只要。4.4简支梁受均布载荷弹性力学结果材料力学结果 4.5楔形体受重力和液体压力五、用差分法和变分法解平面问题5.1差分公式的推导差分法：是微分方程的近似解法，具体的讲，差分法就是把微分用差分来代替，把导数用差分商来代替，从而把基本方程和边界条件（微分方程）近似用差分方程来表示，把求解微分方程的问题变成求解代数方程问题。差分法的数学基础：泰勒公式在节点0处的：x方向的上的差分公式y方向的上的差分公式混合二阶导数的差分公式四阶导数的差分公式5.2应力函数的差分解用差分法解弹性平面问题时，可按下列步骤进行： （

14、1）在边界上任意选定一个结点作为基点A，取，然后由面力的矩及面力之和算出边界上所有各结点处 的值，以及所必需的一些 及 值，即垂直于边界方向的导数值。（2）应用以下公式，将边界外一行虚结点处的值用边界内的相应结点处的值来表示。（3）对边界内的各结点建立差分方程（如下），联立求解这些结点处的 值。（4）按照下面的公式，算出边界外一行的各虚结点处的值。（5）按照下面公式计算应力的分量。说明： 如果一部分边界是曲线的，或是不与坐标轴正交，则边界附近将出现不规则的内结点。对于这样的结点，对（3）中的差分方程必须加以修正。 5.3弹性体的变形势能和外力势能变形势能：注意：变形势能是变形分量或位移分量的二

15、次泛函，叠加原理不再适用； 变形或位移发生时，变形势能总是正的。 外力势能：5.4位移变分方程（1）位移变分方程（拉格朗日变分方程）假想，位移分量发生了位移边界条件所容许的微小改变，即虚位移，或位移变分。同时给出弹性体的限制条件，即没有温度改变（热能没变）和没有速度改变（动能没变）。二维：三维：（2）虚功方程二维：三维：虚功方程表示：如果在虚位移发生前，弹性体是处于平衡状态，那么，在虚位移过程中，外力在虚位移上所做的虚功，等于应力在虚应变上所做的虚功。（3）极小势能原理 在给定外力作用下，在满足位移边界条件的所有各组位移中间，实际存在的一组位移应使总势能成为极值，对于稳定平衡状态，这个值是极小

16、值。5.5位移变分法位移变分法：（1）设定一组包含若干待定系数的位移分量表达式；（2）使它们满足位移边界条件；（3）令其满足位移变分方程（代替平衡微分方程核应力边界条件）并求出待定系数，就同样地能得出实际位移解答。 弹性力学的计算实例及应用ANSYS分析对比结果实例1（差分法与ANSYS有限元法）如图正方形板，单位厚度，边长为L，网格步长h=l/4,q=p/l p=1 h=1（差分法）解：（1）取A为几点，=0 =0 =0 （2）边界上 延顺时针方向积分由式 B:=0 =0 =0 C:=0 =0 =p/2D:=0 =0 =p/20E-G与C、D一样H:=p/2h/2= =p/4 =0I:=ph

17、/2 =p/4 =0应为对称性右边的点不用算（3） 虚节点值 （4） 列差分方程 =借助计算机MATLAB解出- ANSYS有限元法求解 X轴方向应力示意图Y轴方向应力示意图XY平面上的切应力示意图有限元求解结果误差分析位移变分法求解ANSYS有限元求解位移变分法求解ANSYS有限元求解位移变分法求解ANSYS有限元求解A（0,0）-0.104绝对误差0.0070.0610.0141相对误差0.2100.4420.157实例2（位移变分法与ANSYS有限元法）如图铅直平面内的正方形薄板，边长为2a，四边固定如图所示，只受重力作用。设=0，试取位移分量的表达式为板厚0.1，长a=1; E=211

18、e9; 密度8500求点（0，0），（1，0），（0.4，0.4）的用瑞里-里兹法求解： 设=0，试取位移分量的表达式为 用瑞里-里兹法求解： （1）在位移表达式u,v中各取第一项 (a) 将=0代入 得 （b) 其中 ， ，代入（b)得 （c)确定系数。板四周边界上的位移为零，板的体力分量所以可将 化简为 （d) 将式（c)代入（d)得 对称性 积分的 可求得 位移分量的解为 应力分量表达式为 点（0,0） 点（1,0） 点（0.4,0.4） 用ANSYS有限元方法求解 建立有限元模型模型受力变形X方向应力示意图Y方向受力示意图XY面切应力示意图ANSYS后处理结果点（0,0） 点（1,0）

19、 点（0.4,0.4） 位移变分法求解ANSYS有限元求解位移变分法求解ANSYS有限元求解位移变分法求解ANSYS有限元求解A（0,0）000000B（1,0）000358674.5353922C（0.4,0.4）2389.2882255.923630.3124832108329.399245.7绝对误差A：0B：0C:133.388A：0B：C:1201.69A：0B：4752.5C:9083.6相对误差A：B：C：0.0558A：C：0.05A：0.013C：0.08误差分析分析误差产生原因：有限元法分析的概念是用较简单的问题代替复杂问题，得到的是一个近似解，非准确解，存在误差是不可避免

20、的。基于有限元法的理论，用ANSYS软件分析结构，对结果的主要影响因素有：模拟的单元类型选择，节点选择，建模方式，网格划分精度和约束边界条件简化以及荷载施加条件，另外，网格质量，材料参数等都对结果有一定的影响。在静力分析当中，网格划分越密，应力分析结果准确度越高。对于简单的结构，经过ANSYS简化分析，得出的结果一般比较准确，与弹性力学结果吻合得较好。而对于较复杂的结构，弹性力学计算与实际值相差就较大，ANSYS分析结果的准确性很大程度上就取决于建模和边界条件等是否合理了。ANSYS分析的结果值误差通常在10左右，如果不选择其他单元，而只选择在节点上施加荷载或应力集中处的单元，误差值有可能达到50或以上。六、结束语经过一个学期的弹性力学学习，虽然学起来还真的比较的抽象不过经过弹性力学的学习，还是了解到了一些相关的基本理论和一些解题思想。弹性力学，是固体力学的一个分支，研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、形变和位移。弹性力学的研究对象是完全弹性体，弹性体是变形体的一种，在外力作用下物体变形，当外力不超过某一限度时，出去外力