刚架有限元程序设计

《刚架有限元程序设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1N1m1m1m计算各节点位移计算各约束反力计算各节点最大应力刚架结构有限元程序设计1Nm2N1m单元应力程序设计流程一、前期准备工作节点：编号、坐标单元：编号、节点组成、材料属性载荷：作用节点、方向、大小约束：作用节点、方向二、中间实战环节单元刚度矩阵：直接选用整体刚度矩阵：总自由度编号、对号入座节点力列向量：对号入座约束处理：去处奇异性三、求解及后处理解线性方程组得位移代入约束方程得反力代入单元本构方程得应力编程要点一、输入信息及其计算机语言描述节点：编号、坐标单元：编号、节点组成、材料属性载荷：作用节点、方向、大小约束：作用节点、方向

2、1、节点定义一维数组：编程要点一、输入信息及其计算机语言描述2、单元定义两维数组：i表示单元编号，j包含单元节点数、节点号、材料属性号、截面特性号等例：有限元计算编程要点一、输入信息及其计算机语言描述两维数组：i表示载荷编号，j包含载荷作用节点号、自由度号、大小例：3、载荷定义有限元计算编程要点一、输入信息及其计算机语言描述两维数组：i表示约束编号，j包含约束作用节点号、自由度号、值例：4、约束定义编程要点1、总自由度自动编号程序二、核心程序总自由度编号整体刚度矩阵组装节点力列向量形成约束处理i表示节点号，j位置表示节点局部自由度编号，初始

3、置1,运行后数值为总自由度编号。节点码——局部自由度码——总自由度码对应关系数组节点1节点2节点3节点4节点5节点6目标函数(输出信息):确定每个节点的每个自由度对应整体自由度的编号已知条件(输入信息):节点总数nTotalNode每个节点自由度数总数nNodeFreedomi表示节点号，j位置表示节点自由度顺序，数值表示总自由度号。整体编码数组充“1”整体编码数组生成2、单元刚度矩阵生成程序单元信息单元刚度矩阵3、总刚度矩阵生成程序总刚置0单刚装配单元信息单元节点总自由度编号4、节点力向量生成程序总节点力向量置0填充非0节点力5、约束(边

4、界条件)处理程序充大数法:6、解线性方程组——高斯消去法7、求约束反力