悬臂梁MATLAB有限元算例注释
用有限元法对悬臂梁分析的算例
算例:
如下图所示的悬臂梁,受均布载荷q=1N/mm2作用。E=2.1×105N/mm2, μ=0.3厚度h=10mm。现用有限元法分析其位移及应力。
梁可视为平面应力状态,先按图示尺寸划分为均匀的三角形网格,共有8×10=80个单元,5×ll=55个节点,坐标轴以及单元与节点的编号如图。将均布载荷分配到各相应节点上,把有约束的节点5l、52、53、54、55视作固定铰链,建立如图所示的离散化计算模型。
程序计算框图:
(续左)
(接右)
程序中的函数功能介绍及源代码
1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym)――该函数用于计算平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)、第二个节点坐标为(xj,yj)、第三个节点坐标为(xm,ym)时的线性三角形元的单元刚度矩阵.该函数返回6×6的单位刚度矩阵k.
2. LinearTriangleAssemble(K,k,i,j,m)――该函数将连接节点i,j,m的线性三角形元的单元刚度矩阵k集成到整体刚度矩阵K。每集成一个单元,该函数都将返回2N×2N的整体刚度矩阵K.
3. LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u)-- 该函数计算在平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)第二个节点坐标为(xj,yj)、第三个节点坐标为(xm,ym)以及单元位移矢量为u时的单元应力。该函数返回单元应力矢