刚架有限元程序设计.pptVIP

* 1N 1m 1m 1m 计算各节点位移 计算各约束反力 计算各节点最大应力 刚架结构有限元程序设计 1Nm 2N 1m 单元应力 程序设计流程 一、前期准备工作 节点：编号、坐标 单元：编号、节点组成、材料属性 载荷：作用节点、方向、大小 约束：作用节点、方向 二、中间实战环节 单元刚度矩阵：直接选用 整体刚度矩阵：总自由度编号、对号入座 节点力列向量：对号入座 约束处理：去处奇异性 三、求解及后处理 解线性方程组得位移 代入约束方程得反力 代入单元本构方程得应力 编程要点 一、输入信息及其计算机语言描述 节点：编号、坐标 单元：编号、节点组成、材料属性 载荷：作用节点、方向、大小 约束：作用节点、方向 1、节点定义 一维数组： 编程要点 一、输入信息及其计算机语言描述 2、单元定义 两维数组： i 表示单元编号，j 包含单元节点数、节点号、材料属性号、 截面特性号等 例： 有限元计算编程要点 一、输入信息及其计算机语言描述 两维数组： i 表示载荷编号，j 包含载荷作用节点号、自由度号、大小 例： 3、载荷定义 有限元计算编程要点 一、输入信息及其计算机语言描述 两维数组： i 表示约束编号，j 包含约束作用节点号、自由度号、值 例： 4、约束定义 编程要点 1、总自由度自动编号程序 二、核心程序 总自由度编号 整体刚度矩阵组装 节点力列向量形成 约束处理 i 表示节点号，