直杆受集中拉伸外载情形下SAINT-VENANT等效原理的.docVIP

/week114 直杆受集中拉伸外载情形下Saint-Venant等效原理的 数值验证1 冯美生，张红珠，王来贵 辽宁工程技术大学力学与工程科学系，辽宁阜新（123000） E-mail：daokedaoke@163.com 摘 要：在 ansys 平台上，采用有限元方法对 Saint-Venant 等效原理进行数值验证，定义截 面应力分布的不均匀度，定量分析不均匀度的变化规律。数值实验表明，集中荷载作用于梁 端很有限的距离内，导致截面应力不均匀度的急剧变化，这个作用范围仅占梁全尺度的 0.0011%，超过这个尺度范围不均匀度几乎为定值，即“所受影响可以不计”。验证了圣维南 原理的精确性和正确性。 关键词：圣维南原理，应力分布，平面应力，有限元法 1 引 言 在求解弹塑性力学问题时，不仅要使应力分量、位移分量在求解域内完全满足基本方程， 而且在边界上要满足给定的边界条件。但是，在工程实际中物体所受外荷载往往比较复杂， 一般很难完全满足边界条件。当所关心的不是荷载作用部分的局部应力时，圣维南原理可以 用来简化边界条件。圣维南原理[1]：将物体的一小部分边界上的面力，变换为分布不同但静 力等效的面力（主矢和主矩相同），那么，尽管在力的作用区域的附近，应力分布将有显著 的改变，但在离力的相当远处，所受的影响可以不计。圣维南原理提出至今已有 100 多年的 历史，虽然目前还没有确切的数学表示和严格的理论证明，但无数的实际计算和实验测量都 证实了它的正确性[2]。文章拟在ansys数值平台上[3]验证等效原理的可靠性及适用范围。 2 计算模型 如图 1 所示：一直杆左端固定，长L=100，高H=10，厚t=1，弹性模量E=2.1×107，右端 中部受拉力P=100 作用。简化为平面应力问题，采用beam3 单元分析，受力简图如图 1 所示。 Y  H  L 图 1 计算简图 Fig.1 Calculating sketch map  P  X  1  本课题得到国家自然科学基金面上项目的资助。 -1- /week114 3 各截面?? x 分布不均匀度定义及其变化规律研究  定义任一截面?? x 沿 y 方向的不均匀度 r 为： r?? max?? x (d , y) ? x  ，其中：?? x??  P A  ， ? x (d , y) 为 x=d 截面上沿 x 方向的应力分布，如图 2 所示。  ? x ( y)  图 2  P L－d 不均匀度计算简图 Fig.2 Calculating sketch map of non-uniform degree H=10 时，有限元计算模型图、受拉变形图及受拉应力云图见图 3。 图 3 有限元计算模型图、受拉变形图及受拉应力云图 Fig.3 Calculating model and deformation pattern and stress envelope as pulled H 等于 20，30 时的有限元计算模型及受拉变形图、拉应力云图略。经 FEM 计算出 不同 H 值时，对应的 d~r 曲线见图 4，其曲线方程为： ?r?? 1, d?≤ 99.9989; ? 2 其中，r0=0.48789；A=0.1387；w=0.00158；dc=100.00148。 由图 3、4 分析知：在 d≤99.9989 时，r 值几乎为定值 1，就是说在 d≤99.9989 的各 截面上应力分布很均匀，各应力值相等，并且和名义应力等值；在 d＞99.9989 的各截 面上应力分布变化剧烈，见图 4，其曲线斜率急剧增大，不同的梁高 H，表现了相同的 变化趋势，即应力不均匀度 r 值的变化在这一区段内积聚的改变，但在离开该区段不远 后就趋于定值 1。以上分析就应证了圣维南原理，“尽管在力的作用区域的附近，应力 分布将有显著的改变，但在离力的相当远处，所受的影响可以不计”。这里“应力值的显 著改变”突出表现在应力不均匀度 r 值由 1 剧变为 12.63。这里的“相当远处，所受的影 -2- /week114 响可以不计”，就可以认为在 d≤99.9989 时，所受的影响可以不计。 4 误差分析 以 H=10 为 例 ， 在 d ＞ 99.9989 时 ， 其 截 面 最 大 应 力?? x 的 相 对 误 差 为 ：  ??? ? x???? ?  ?100%??  12.63??10 10  ?100%?? 26.3% 。在 d≤99.9989 时，? 值几乎为零。 ? x 变化剧烈区占总长度的百分比?