有限单元法和程序设计

《有限单元法和程序设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、有限单元法及程序设计有限单元法及程序设计绪论1.力学分析方法:解析法,数值法有限元法——实际结构形状和所受载荷比较复杂,大多用解析法很困难,因而数值法得到不断发展,随着电子计算机的进步,而发展起来的一种新兴的数值分析方法.2基本步骤：（1）结构离散化：将结构从集合上用线或面划分为有限个单元。（2）单元分析：导出单元的节点位移和结点力之间的关系（单元刚度矩阵）。（3）整体分析：将各单元组成的结构整体进行分析，导出征个结构点位移与结点力之间的关系。3程序设计的步骤：提出问题，拟定解决方案构造数学模型画出

2、程序流程图编写程序编译调试程序试算验证程序4.根据国家标准（GB-1526-89）规定的程序流程图标准化符号及规定：EMBEDPBrush图表示程序流程图的起点和终点；图表示数据信息的输入和输出；图表示数据进行系列运算之前要完成的数据预置；图表示判断条件；图表示各种处理功能，如数学运算方式等；图表示流程的路径和指向。杆件结构的有限单元法及程序设计平面杆件单元的有限单元法有限单元法的基本概念基本思路：先分后合（先单元分析，再整体分析）基本概念：整体号：节点端点号按自然数1，2，3，……（在整体坐标系x

3、Oy下）局部号：每一个单元始末用i，j标记（在单元的局部坐标EMBEDEquation.3系下，方向与整体坐标系一致）。3．Fe=keδe其中：ke=EMBEDEquation.3单元刚度矩阵，各元素为刚度系数δe=EMBEDEquation.3单元杆段位移列阵Fe=EMBEDEquation.3单元杆端力列阵KEMBEDEquation.3=P（1-7）K=EMBEDEquation.3整体刚度矩阵EMBEDEquation.3=EMBEDEquation.3EMBEDEquation.3位移列阵

4、P=EMBEDEquation.3节点载荷列阵3.有限元位移法分析连续梁需要考虑的问题刚度集成法:①将(1-3)K扩阶,扩大的元素为0,得到单元贡献矩阵单元①:K①=EMBEDEquation.3单元②:K②=EMBEDEquation.3②将单元贡献矩阵想叠加,形成整体刚度矩阵K=K①+K②=EMBEDEquation.3两端支承条件的引入先不考虑约束条件,得到整体刚度矩阵后,将其主对角线元素kEMBEDEquation.3改为1,第i行,第j列其余元素改为0,对应的载荷元素也改为0.非结点荷载的

5、处理利用等效结点荷载进行分析:各结点(包括两端结点)加约束,阻止结点转动,其约束力矩分别为交于该结点的各相关单元的固端力矩之和,顺时针为正.去掉附加约束(相当在各结点施加外力荷载PEMBEDEquation.3,其大小与约束力矩相同,方向相反)将两部分杆端弯矩叠加起来.第二节局部坐标系中的单元刚度矩阵一般单元设单元eqoac(○,e)的弹性模量、截面惯性矩、截面积分别为E、I、A,杆长为l。单元的i、j端各有三个杆端力、EMBEDEquation.3(即轴力、剪力和弯矩)和与其相应的三个杆端位移

6、EMBEDEquation.3，如图1-7所示。图中EMBEDEquation.3为单元局部坐标系，取i点位于坐标原点，EMBEDEquation.3轴与杆轴重合，规定由i到j为EMBEDEquation.3轴的正方向，由EMBEDEquation.3轴顺时针旋转90?为EMBEDEquation.3轴正方向。力和位移的正方向如图1-7所示。EMBEDPBrush在此单元中，单元杆端力列阵和杆端位移列阵分别为EMBEDEquation.3单元杆端力列阵EMBEDEquation.3杆端位移列阵为了导

7、出一般单元杆端力与杆端位移之间的关系，我们分别考虑以下两种情况。首先分析两个杆端轴力EMBEDEquation.3与轴向位移EMBEDEquation.3的关系。根据胡克定律，有EMBEDEquation.3（a）其次考虑杆端弯矩EMBEDEquation.3与杆端剪力EMBEDEquation.3与杆端转角EMBEDEquation.3和横向位移EMBEDEquation.3的关系。根据结构力学位移法的转角位移方程，并按照本节规定的符号和正负号，可得EMBEDEquation.3(b)将(a)、(

8、b)两式合在一起，并写成矩阵形式如下EMBEDEquation.3EMBEDEquation.3(1-17)上式可简写成EMBEDEquation.3(1-18)其中单元刚度矩阵为EMBEDEquation.3（1-19）单元刚度矩阵的性质每个元素代表单位杆端位移引起的杆端力,任一元素EMBEDEquation.3(r、s取1至6)的物理意义是第s个杆端位移分量等于1时,所引起的第r各杆端力分量值.是对称矩阵,其元素EMBEDEquation.3.是奇异矩阵,它的元素