采用极坐标求解弹性力学平面问题基本问题.doc

采用极坐标弹性力学平面问题基本问题一内容介绍 在弹性力学问题的处理时，坐标系的选择从本质上讲并不影响问题的求解，但是坐标的选取直接影响边界条件的描述形式，从而关系到问题求解的难易程度。对于圆形，楔形，扇形等工程构件，采用极坐标系统求解将比直角坐标系统要方便的多。本章的任务就是推导极坐标表示的弹性力学平面问题基本方程，并且求解一些典型问题。 二重点 1基本未知量和基本方程的极坐标形式；2双调和方程的极坐标形式；3轴对称应力与厚壁圆筒应力；4曲梁纯弯曲、楔形体和圆孔等典型问题学习思路: 选取极坐标系处理弹性力学平面问题，首先必须将弹性力学的基本方程以及边界条件通过极坐标形式描述和表达。 本节的主要工作是介绍基本物理量，包括位移、应力和应变的极坐标形式；并且将基本方程，包括平衡微分方程、几何方程和本构关系转化为极坐标形式。 由于仍然采用应力解法，因此应力函数的极坐标表达是必要的。 应该注意的是坐标系的选取与问题求解性质无关，因此弹性力学直角坐标解的基本概念仍然适用于极坐标。 学习要点： 1；2；3；4；5；6Laplace算符；7。 为了表明极坐标系统中的应力分量，从考察的平面物体中分割出微分单元体ABCD，其由两个相距d? 的圆柱面和互成 d? 的两个径向面构成， 在极坐标系中，用???表示径向正应力，用???表示环向正应力,????和????分别表示圆柱面和径向面的切应力，根据切应力互等定理，????＝????。 首先推导平衡微分方程的极坐标形式。考虑到应力分量是随位置的变化，如果假设AB面上的应力分量为???和???,则CD面上的应力分量为 如果AD面上的应力分量为???和????，则BC面上的应力分量为 。 同时，体力分量在极坐标径向??和环向?? 方向的分量分别为Fb???和Fb??。 设单元体的厚度为1， 考察其平衡首先讨论径向的平衡，注意到 ，可以得到 简化上式，并且略去三阶微量，则 同理，考虑微分单元体切向平衡，可得 简化上式，可以得到极坐标系下的平衡微分方程，即 以下推导极坐标系统的几何方程。 在极坐标系中，位移分量为u?，u??，分别为径向位移和环向位移。 极坐标对应的应变分量为：径向线应变??，即径向微分线段的正应变；环向线应变???为环向微分线段的正应变；切应变???为径向和环向微分线段之间的直角改变量。 首先讨论线应变与位移分量的关系，分别考虑径向位移环向位移u?，u?? 所引起的应变。如果只有径向位移u?，借助于与直角坐标同样的推导，可以得到径向微分线段AD的线应变 为；环向微分线段AB=? d? 的相对伸长为 ；如果只有环向位移u?? 时，径向微分线段线没有变形， 环向微分线段的相对伸长为 ；将上述结果相加，可以得到正应变分量 , 下面考察切应变与位移之间的关系。设微分单元体ABCD在变形后变为ABCD， 因此切应变为 ??? =??+ (??- ?)? 上式中??表示环向微分线段AB向???方向转过的角度，即 ；??表示径向微分线段AD向??方向转过的角度，因此 ；而 ??角应等于A点的环向位移除以该点的径向坐标?，即 。 将上述结果回代，则一点的切应变为。综上所述，可以得到极坐标系的几何方程为 由于讨论的物体是各向同性材料的，因此极坐标系的本构方程与直角坐标的表达形式是相同的，只要将其中的坐标 x 和 y 换成???和??就可以了。对于平面应力问题，有 对于平面应变问题，只要将上述公式中的弹性常数E，??分别换为就可以。Laplace算符 平面问题以应力分量形式表达的变形协调方程在直角坐标系中为 。由于??x+??y= ???+????为应力不变量，因此对于极坐标问题，仅需要将直角坐标中的Laplace转换为极坐标的形式。 因为，x=? cos?， y=? sin? ，即。将? 和? 和分别对x和y求偏导数，可得 根据上述关系式，可得以下运算符号 则 将以上两式相加，简化可以得到极坐标系的Laplace算符。 另外，注意到 ，因此在极坐标系下，平面问题的由应力表达的变形协调方程变换为 如果弹性体体力为零，则可以采用应力函数解法求解。不难证明下列应力表达式是满足平衡微分方程的 这里(?，?)是极坐标形式的应力函数，假设其具有连续到四阶的偏导数。将上述应力分量表达式代入变形协调方程，可得 显然这是极坐标形式的双调和方程。 总而言之，用极坐标解弹性力学的平面问题，与直角坐标求解一样，都归结为在给定的边界条件下求解双调和方程。 在应力函数解出后，可以应用 求解应力，然后通过和求解应变分量和位移分量学习思路: 如果弹性体的结构几何形状、材料性质和边界条件等均对称于某一个轴时，称为轴对称结构。轴对称结构的应力