杆梁结构的有限元分析原理.pptVIP

ANSYS应用实例2 第62页，共92页，编辑于2022年，星期四 4.4梁单元及其坐标变换 一般平面梁单元的描述 第63页，共92页，编辑于2022年，星期四 纯弯梁单元 第64页，共92页，编辑于2022年，星期四 由于单元有四个位移分量，可设梁单元的位移模式v(x)为包含4个待定常数的三次多项式： 由材料力学知,各截面的转角: 第65页，共92页，编辑于2022年，星期四 第30页，共92页，编辑于2022年，星期四 第31页，共92页，编辑于2022年，星期四 第32页，共92页，编辑于2022年，星期四 第33页，共92页，编辑于2022年，星期四 第34页，共92页，编辑于2022年，星期四 4.3杆单元及其坐标变换-局部坐标 由于杆单元只有两个节点位移，故可以设杆单元的位移模式为之包含两个待定常数的形式 u(x)=a1+a2x 根据有限元法的基本思路，将弹性体离散成有限个单元体的组合，以结点的位移作为未知量。弹性体内实际的位移分布可以用单元内的位移分布函数来分块近似地表示。在单元内的位移变化可以假定一个函数来表示，这个函数称为单元位移函数、或单元位移模式。 第35页，共92页，编辑于2022年，星期四 回代得 写成矩阵形式为 其中Ni，Nj是形函数。 根据位移条件有u(0)=u0， u(l)=ul，从而得 第36页，共92页，编辑于2022年，星期四 根据几何方程得 根据物理方程得 从而，根据单元分析结果，进行整体分析，求解整体方程组，进行结果分析 第37页，共92页，编辑于2022年，星期四 4.3.2杆单元的坐标变换 第38页，共92页，编辑于2022年，星期四 规定：杆端位移和杆端力取在截面形心上，符号以与单元系坐标正向相同为正，相反为负。下面讨论整体坐标系下与局部坐标系下的转换关系式。整体坐标系单元杆端位移和杆端力仍定义在截面形心上，符号以与坐标正向同向为正反之为负。 局部坐标系 整体坐标系 4.3.2杆单元的坐标变换-平面问题 第39页，共92页，编辑于2022年，星期四 其中是一个单位正交矩阵，单位正交矩阵的逆即等于其转置。 从上图可以得出，整体坐标系逆针旋转α角后与单元系相重合。 写成矩阵形式为 第40页，共92页，编辑于2022年，星期四 第41页，共92页，编辑于2022年，星期四 由于单元的势能是一个标量(能量)，不会因坐标系的不同而改变，因此，可将节点位移的坐标变换关系代入原来基于局部坐标系的势能表达式中， 整体坐标系下的刚度方程 第42页，共92页，编辑于2022年，星期四 根据 得 其中 单刚的性质： 是对称矩阵。 是奇异矩阵。 坐标变换并不改变矩阵的奇异性质。 第43页，共92页，编辑于2022年，星期四 第44页，共92页，编辑于2022年，星期四 1 结构的离散化与编号 第45页，共92页，编辑于2022年，星期四 2各个单元的矩阵描述 结构包括有斜杆，所以必须在总体坐标下对节点位移进行表达，所推导的单元刚度矩阵也要进行变换 第46页，共92页，编辑于2022年，星期四 3 建立整体刚度方程 1.将所得到的各个单元刚度矩阵按节点编号进行组装，可以形成整体刚度矩阵;2.同时将所有节点载荷也进行组装。 第47页，共92页，编辑于2022年，星期四 4 边界条件的处理及刚度方程求解 第48页，共92页，编辑于2022年，星期四 第49页，共92页，编辑于2022年，星期四 5 各单元应力的计算 第50页，共92页，编辑于2022年，星期四 6 支反力的计算 将节点位移的结果代入整体刚度方程中 第51页，共92页，编辑于2022年，星期四 基于MATLAB平台求解该 （1） 结构的离散化与编号 （2）计算各单元的刚度矩阵 1.建立一个工作目录，将所编制的用于平面桁架单元分析的四个MATLAB函数 （1.单元刚度；2.总刚矩阵的组装；3.单元应力的求解；4.支反力的求解） 2.在MATLAB环境中，输入弹性模量E、横截面积A，各点坐标、角度 3.调用四次单元刚度矩阵计算函数，得到各个单元的刚度矩阵 第52页，共92页，编辑于2022年，星期四 单元的刚度矩阵的计算 function k=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2,alpha) %该函数计算单元的刚度矩阵 %输入弹性模量E，横截面积A %输入第一个节点坐标（x1,y1），第二个节点坐标（x2,y2），角度alpha（单位是度） %输出单元刚度矩阵k(4X4)。 %--------------