有限元基础和CAE.ppt

有限元与CAE技术 江苏科技大学 机械设计教研室 CH1绪论 一、有限元的概念 2.单元分析 对于弹性力学问题，单元分析，就是建立各个单元的节点位移和节点力之间的关系式。 由于将单元的节点位移作为基本变量，进行单元分析首先要为单元内部的位移确定一个近似表达式，然后计算单元的应变、应力，再建立单元中节点力与节点位移的关系式。 单元分析的步骤可表示如下： 第二章 弹性力学基础 一、弹性力学与材料力学 第二章弹性力学基础1、研究的内容：基本上没有什么区别。 弹性力学也是研究弹性体在外力作用下的平衡和运动，以及由此产生的应力和变形。 2、研究的对象：有相同也有区别。 材料力学基本上只研究杆、梁、柱、轴等杆状构件，即长度远大于宽度和厚度的构件。弹性力学虽然也研究杆状构件，但还研究材料力学无法研究的板与壳及其它实体结构，即两个尺寸远大于第三个尺寸，或三个尺寸相当的构件。 3、研究的方法：有较大的区别。 虽然都从静力学、几何学与物理学三方面进行研究，但是在建立这三方面条件时，采用了不同的分析方法。材料力学是对构件的整个截面来建立这些条件的，因而要常常引用一些截面的变形状况或应力情况的假设。这样虽然大大简化了数学推演，但是得出的结果往往是近似的，而不是精确的。而弹性力学是对构件的无限小单元体来建立这些条件的，因而无须引用那些假设，分析的方法比较严密，得出的结论也比较精确。所以，我们可以用弹性力学的解答来估计材料力学解答的精确程度，并确定它们的适用范围。 五、应力与外力之间的平衡方程 根据微元体受合力为零的条件可以得出三维平衡方程： 六、两种平面问题 弹性力学可分为空间问题和平面问题，严格地说，任何一个弹性体都是空间物体，一般的外力都是空间力系，因而任何实际问题都是空间问题，都必须考虑所有的位移分量、应变分量和应力分量。但是，如果所考虑的弹性体具有特殊的形状，并且承受的是特殊外力，就有可能把空间问题简化为近似的平面问题，只考虑部分的位移分量、应变分量和应力分量即可。 1、平面应力问题 1）应力分量 只需要考虑三个应变分量即可，于是应变矩阵简化为： 第三章 杆系结构的有限元法 §3.1 连续梁问题的有限元法 1.离散化 在杆系结构中各杆具有自然划分，成为杆单元，记为e。图1所示的连续梁可划分为两个单元，分别称为单元①、②，各个铰接点称为节点，统一编号为1、2、3称为节点总码，而称为节点力矩载荷，如图2所示： 2.单元分析 单元分析就是对已划分的单元进行力学分析，没有必要对每一个单元都进行分析，而只需要对一个典型单元进行分析即可。为此取任意单元e进行分析，首先对该单元的两个端点重新编码为 i、j，称为局部码。在该单元的两个端点分别作用有 两个力矩，称为杆端力矩。在这两个力矩的作用下杆发生变形，如图3所示，在节点处分别产生了转角 ，称为节点转角。 3.整体分析 整体分析的主要任务就是在单元分析的基础上得到整体刚度矩阵。 厚度为t的很薄的均匀木板。只在边缘上受到平行于板面且不沿厚度变化的面力，同时，体力也平行于板面且不沿厚度变化。 以薄板的中面为xy面，以垂直于中面的任一直线为Z轴。由于薄板两表面上没有垂直和平行于板面的外力，所以板面上各点均有： 另外由于平板很薄，外力又不沿厚度变化，可认为在整个薄板内各点均有： 于是，在六个应力分量中，只需要研究剩下的平行于XOY平面的三个应力分量， 所以称为平面应力问题。 2）物理方程 物理方程中后两式可见，这时的剪应变： 由物理方程的第三式可见： 如图1所示为一个连续梁，在铰点处分别作用有力矩载荷M1、M2、M3，求连续梁的内力。 两种情况的叠加 将上述两式写成矩阵形式为： 1）暂时不引入支承条件和载荷情况，先建立整体刚度矩阵K，写出M 与θ之间的转换关系 2) 在节点1和2引入载荷值，在固定端引入支承条件θ3=0，将上式修改为 为了求解未知转角θ1、θ2，将上式展开： 为了便于编程，我们希望修改后的矩阵仍然保留原矩阵的阶数和排列顺序。为此扩大如下形式： 矩阵形式 刚度方程修改方法： 1）把节点力向量M换成节点载荷向量。 P 2）在矩阵K与零转角对应的行列中，主对角线元素改为1，其它元素改为0。 3）在载荷向量中，与零转角对应的元素改为0。 5.非节点载荷的引入 1）求等效节点载荷 图10 所示的连续梁为例： Ⅰ、 求单元1与2产生的固端力矩：（加角标0表示固端力矩，杆端采用单元局部编码） Ⅱ、各节点的约束力矩分别为该节点的相关单元固端力矩之和。 III、各节点的约束力矩的反作用力矩为原非节点载荷的等效节点载荷，即 2）求各杆端弯矩 连续梁在非节点载荷作用下的杆端弯矩由两部分组成：一部