三角形单元的悬臂梁圣维南原理的数值验证.docxVIP

SHANGHAIUNIVERSITY工程中的数值分析方法课程论文题 目:验证圣维南原理学 院 土木工程系专 业 建筑与土木工程学 号 生姓名 邓高林指导教师 朱杰江日 期 2017.3.16三角形单元的悬臂梁圣维南原理的数值验证摘要：由三角形单元组成的悬臂杆一端固定一端受集中拉伸外载作用为例, 数值验证了圣维南原理的正确性: 作用在物体一端( 次要边界或是小边界)的荷载, 如果只改变应力分布而不改变合成, 那么就只会显著改变该端附近的应力, 在距端部较远的地方相差甚微。关键词：三角形单元；集中拉伸外载； 圣维南原理Numerical verification of cantilever triangular element of the Saint Venant principleAbstract: the unit composed of triangle cantilever bar end fixed under concentrated tensile loads as an example, numerical simulation verified the correctness of the Saint Venant principle: the role of the object is (secondary edge or boundary) load, if only change the stress distribution without changing the synthesis, then it will only change the stress the near end, from the end of the distant place is little.Keywords: triangular element; concentrated tensile load; Saint Venant principle0 引言在实际问题中圣维南原理有着广泛的应用。大量的工程实践表明圣维南原理的正确性。本文以三角形单元悬臂梁为模型一端固定一端受集中拉伸外载作用为例，数值验证验证了圣维南原理的正确性。1 圣维南原理圣维南原理是这样叙述的: 如果把物体的一小部分边界上的面力, 变换为分布不同但静力等效的面力(主矢量相同, 对于同一点的主矩也相同), 那么, 近处的应力分布将有显著的改变, 但是远处所受的影响可以不计。它也可以这样来陈述: 如果物体一小部分边界上的面力是一个平衡力系(主矢量和主矩都等于零), 那么, 这个面力就只会使得近处产生显著的应力, 远处的应力可以不计。 这样的陈述和上面的陈述完全等效, 因为静力等效的两组面力, 它们的差异是一个平衡力系。 在应用圣维南原理时, 绝不能离开“ 静力等效”的条件。2 模型构建一个只受水平力的悬臂梁模型，假设其宽为h，长为16a，高为2a，弹性模量为E不变，泊松比μ为0，只受节点集中力P，如图1 所示。将本结构离散化，划分为由64个直角等腰三角形单元组成的结构，含有51个节点，其中节点1、18和35都是约束节点，没有位移。3 计算单元刚度阵由有限元知识理论可知三角形单元刚度矩阵其中A为单元面积(r , s=i , j , m)（i , j , m按顺序换可的j , m 对应的数）表1. 单元节点和值（i、j、m）单元刚度阵为由于平面问题的单元刚度矩阵不随单元（或坐标轴）的平行移动或作 角度（n为整数）的转动而改变，所以本结构中各个三角形单元的刚度矩阵是相同的。4 总刚合成以1单元和2单元为例其他单元按此方法“按号入座”形成总的单元刚度K（见附录1）5 计算各节点的位移P=hF(见附录2)通过编程计算得到各个节点位移（见附录3）6 计算应力由物理方程得S称为应力矩阵，由于S是常数矩阵，所以三角形单元内的应力是一样的。将表1的数据带入S矩阵中得求各个单元的应力只需将节点i、j、k 所对应的节点位移带入即可。由于应力矩阵是常数，所以单元内应力都相等。 在2、3节点中间的截面总的合力为0.97482Fh，总的合力接近P=Fh，但误差较大，但也符合圣维南原理。