能量原理与变分法.ppt

对上式中每一项进行分部积分，有：因为在边界上，有：在域内，所以，有，第95页,共123页，星期六，2024年，5月§11-8应力变分法应于扭转问题1.扭转应力变分方程等截面直杆的扭转问题中，存在应力函数?，横截面剪应力可表示为：形变余能及其变分式中：函数?为Prandtl应力函数。将其代入形变余能计算式：∵应力函数?仅为x、y的函数，可将上述积分变为：其中：L为杆的长度；G为剪切弹性模量。第96页,共123页，星期六，2024年，5月（a）外力的功及其变分扭转杆件侧面上无外力，因而不存在面力的功。在两端作用有方向相反的两扭矩M，两端的相对转角为：KL，则面力在位移上的功为：W=MKL由上一章的结果：得外力的功为：外力的功的变分为：（b）扭转变分方程将式（a）（b）代入变分方程，有：第97页,共123页，星期六，2024年，5月扭转变分方程将式（a）（b）代入变分方程，有：或：（c）以上两式即为扭转问题的变分方程或最小余能原理。（1）式（c）中：——扭转问题的总余能说明：（2）式（c）中的应力函数?已满足了两端的边界条件。第98页,共123页，星期六，2024年，5月2.扭转问题的变分方法由于扭转应力函数?要求在边界上的值等于零，其形式可设为：其中：Am为互不相关的m个系数。为使应力函数?在边界上的值等于零，必须要求函数?m都在横截面的边界上的值为零。将?代入扭转变分方程，注意到其变分是由系数Am的变分来实现的，所以有：（11-26）得到一m阶的线性方程组，恰好可用来m个系数Am。第99页,共123页，星期六，2024年，5月例：图示矩形扭转杆，材料的剪切弹模为G。试求其单位长度扭转角，剪应力等。yxOAa/2a/2解：设定扭转应力函数矩形四根边界线的方程：为满足?在边界上的值为零，可取：（d）由截面的对称性，或薄膜比拟，应力函数?应为x、y的偶数，所以，式中m、n都只需取为偶数。对正方形截面杆（b=a），若在（d）中只取一项（m=n=0）,则有第100页,共123页，星期六，2024年，5月对正方形截面杆（b=a），若在（d）中只取一项（m=n=0）,则有yxOAa/2a/2代入式（11-26）有：（11-26）代入扭转变分方程确定待定系数第101页,共123页，星期六，2024年，5月经积分运算，得：从而，有：（e）由公式（10-5）：有：由此求得：对照式（10-21）：有：?的精确值：0.141，相差：0.14%。两者仅将求得的K代入式（e），有：第102页,共123页，星期六，2024年，5月（10-2）对照式（10-22）：由公式（10-2），可求得应力分量：精确值为：0.208,相差：6.8%。两者如果要得更精确的解，需在式（d）中较多的系数项。如：第103页,共123页，星期六，2024年，5月进行与上相同的运算，得到：由此算出的单位扭转角K比精确值只小0.14%；最大剪应力?max比精确值只大出4%。第104页,共123页，星期六，2024年，5月§11-9解答的唯一性1.问题的提出弹性力学问题的求解方法与途径：（1）解析解：就所取未知量，有：按应力求解；按位移求解；就所用坐标系，有：直角坐标求解；极坐标求解；就解的函数形式，有：多项解；级数解；其它函数解；复变函数解；（2）数值解：有限差分解（FDM）；有限单元法（FEM）；边界单元法（BEM）；离散单元法（UDEC）；不同方法、不同途径得到的不同形式的解，其数值是否唯一？第105页,共123页，星期六，2024年，5月2.解答的唯一性及其证明相应于一定的体力和边界条件，某一弹性力学问题的解是唯一的。——也称解的唯一性定理。解的唯一性定理证明：（反证法）假设：在一定的体力、面力、边界条件下，某个弹性力学问题存在两组解：（1）（2）考察这两组解是否相同？∵它们都为同一问题的解，∴应满足相的平衡方程和边界条件。对于第一组解，有：第106页,共123页，星期六，2024年，5月对于第一组解，类似有：将上述两组不同的解方程两边分相减，有：第107页,共123页，星期六，2024年，5月可见：两组不同的解的差，对应的状态为