应用有限元分析工程实例-1.ppt

三角形常应变单元，单元刚度矩阵是6×6的对称矩阵，其显式形式如下： 例2.1 平面应力直角三角形单元ijm，直角边长分别为a, b，厚度为 t，弹性模量为 E，泊松比为 ，求单元的刚度阵。 三角形面积为 刚度阵中的 单元刚度矩阵 当单元为等腰直角三角形单元时，其单元刚度矩阵与单元的边长无关。 对于平面问题，如果材料的物理特性相同，则几何形状相似的单元都具有相同的单元刚度矩阵。 在程序中利用这一特点，可以减少计算量。 当a=b时，为等腰直角三角形单元，其刚度矩阵为 （a) （b) （c) （d) 讨论 (a) (b) (c) (d) 三角形ABC 和 三角形 abc的单元刚度阵是否相同？ 四边形 ABCD和四边形abcd 的单元刚度阵是否相同？ 在图（a）中，若边AB、ab固定、边BC和边bc 仅限制其竖直方向位移，并在点 C 和 点c 施加同样的力 F。 假定两个三角形材料相同，并按线弹性考虑，则： 1、点 C 和 点c的水平位移是否相等？若不相等，哪个点的位移大？ 为什么？ 2、三角形ABC 和 三角形 abc的位移场和应力场是否相同？ 2.2.4 单元等效节点载荷 作用在结构上的外载荷:集中力、表面力、体积力等。 集中力作用点划分分节点。 表面力和体积力采用虚位移原理移植到单元的节点上。 移植后的节点载荷和移植前的载荷在约束允许的任意虚位移上所做的虚功相等。 设单元的虚位移： （1）体积力:自重、离心力等。特点：体积力分布在整个单元上。 令单位体积的力 令与体积力等效的节点载荷： 等效节点载荷的虚功 = 移植前载荷的虚功 体积力可以是常数，也可以是坐标的函数。 如果体积力是重力，且重力方向为负 y方向，则： 单元体积力等效到单元节点上的计算公式为： （2）表面力 表面力作用在单元边界上，令 表面力矢量为： 令 与表面力等效的单元节点载荷为： 等效节点载荷的虚功 = 移植前表面力的虚功 单元表面力的等效节点载荷计算公式 积分沿作用力边界。 表面力可以是沿单元边界的均布力、非均布力 均布力：以右图所示为例 是边 的长度。 相当于将表面力均分到两个节点上。 主要参考文献： [1] 王勖成,邵 敏 编著. 有限元法基本原理和数值方法 [M]. 北京： 清华大学出版社,1999 [2] [美]K.J. Bathe著(付子智译). 工程分析中的有限元法[M]. 北京： 机械工业出版社, 1991 [3][英]O.C.监凯维奇著. 有限元法 [M]. 北京：科学出版社,1985 [4] 赵经文,王宏钰编著. 结构有限元分析 [M] . 北京：科学出版社,2001 [5]梁醒培,王辉编著. 应用有限元分析[M]. 北京：清华大学出版社,2010 [6]梁醒培,王辉著. 基于有限元法的结构优化设计—原理与工程应用[M]. 北京：清华大学出版社,2010 [7]博弈创作室 编著. ANSYS9.0 经典产品高级分析技术与实例详解[M]. 北京：中国水利水出版社,2005 [8]ANSYS程序电子文档 本讲义系为学习班专门编写，未经许可，不得翻印。 梁醒培 2012.8.25 第一讲 1 绪论 1 有限元法简介 1.1 有限元法发展概况 （1）起源： 飞机的结构分析，用三角形单元求得了平面应力问题的解答 中国科学院冯康教授 1960年克拉夫（Clough）第一次提出“有限单元法”这一名称 数学家和工程师共同奠定了有限元法的理论基础和框架 （2）总体发展概况: 国外20世纪70年代就有了商业有限元程序系统 国内： 1980年以前 高校、科研院所、大型厂矿企业开始研究， 特点是功能单一、重复、不能解决较为复杂的工程问题 1980-1990大型有限元程序的