#
单击此处编辑母版标题样式单击此处编辑母版文本样式第二级第三级第四级第五级主讲:简政有限元法及程序设计FINITE ELEMENT METHOD PROGRAME DESIGN第二章 FEM原理及方法2-1 FEM力学模型1.剖分(Disassembly)——离散化(Discretization)化整为零:将连续体(求解域)分成一组离散 单元
Lagrange单元族插值函数
Gauss —choleski 分解法(LDLT)有限元的成果表达方法ABAQUS用Pascal多项式i
#
1一般集合规则力方向[K]中的元素分布是稀疏的每一个子块行存在大量零元素结点越多0元素也越多这是因为每个子块行非零元素和相关结点有关所以:排除刚体位移后 [K]是正定的
#
单击此处编辑母版标题样式单击此处编辑母版文本样式第二级第三级第四级第五级工程力学系有限元程序设计重庆大学工程力学专业主干课程工程力学系课程设计任务及要求1)完成平面3结点(4和8结点)单元有限元程序设计2)利用设计的有限元程序计算分析一力学问题并与理论解或商用软件分析结果进行比较3)书面报告4)口头报告(5分钟)警告:独立完成不得抄袭工程力学系有限元程序系统的组成及分析过程前处理有限元分析本体程序
形成总体刚度矩阵开始[K]的存储约束引入求解行 号N→列号方阵形式的刚度矩阵[K]6 变带宽存贮:按列存贮方式从左到右逐列存放对每一列先存主对角线元素然后由下而上顺序存放直到顶线下第一个元素为止为避免混淆我们把存贮[K]的一维数组称为[A] 实现变带宽存贮的关键问题是:总刚中元素Kij在一维数组A中的地址是什么为此需要知道主元Kii在A中的位置和相应列高hi主元位置:采
5解线性方程组得位移代入约束方程得反力代入单元本构方程得应力有限元计算编程要点两维数组:i 表示节点号j 位置表示节点自由度顺序数值表示总自由度号总刚置0
违法有害信息,请在下方选择原因提交举报