桥梁电算.pptVIP

交通土建工程软件分析 本课程的概述 1.本课程在交通土建工程中的地位 2.学习本课程需预先掌握的知识 3.本课程为后续课程及毕业设计服务 3.本课程的组成 理论 实践(上机) 4.学习本课程的注意点 第一章 平面杆系有限元法 §1.1 引 言 手算与电算方法之比较: 手算：力求简化方法，近似处理 电算：力求标准化，规一化 土建结构计算机分析的程序设计一般概念 程序设计：将求解问题的公式与计算过程按算法语句标准编写成机器可以运行的一个源程序。 程序设计的基础与关键问题，对分析方法，计算方法的充分理解是：编制程序的算法设计是程序设计的组成部分是关键。 1.引例:如下图，求解节点3位移 引例 以整体坐标计，结构的节点3受 到 的节点力，发生 的位移。 引例(单元1) 以局部坐标计，以①杆件为 轴，节点1至3方向为正，则节点3 在 方向的位移有： 杆件轴力为： 右图中位移箭头方向为正，轴力拉为正 引例（单元1） 以整体坐标计，①杆件的受力为： 单 元 ① 引例（单元2） 在 方向的位移有： 轴力： 单 元 ② 引例（综合） 整理得： 单 元 ① 单 元 ② 引例（综合） 矩阵形式： 单 元 ① 单 元 ② 引例（整体分析） 取节点3力的平衡，得： 即： 将单元分析结果代入： 引例（整体分析） 即： 则有： 引例（作业） 以整体坐标计，结 构的节点3受到 的节点力，发生 的位移，求解 。 §1.2 单 元 分 析 一、单元刚度矩阵: 单元刚度矩阵即用杆端位移表示杆端力的系数矩阵。 以下根据结构力学的转角位移方程与材料力学的轴力计算公式来推导平面杆单元在局部坐标下的刚度矩阵： 结构力学公式1: 结构力学公式2: 结构力学公式3: 矩阵表示 写成矩阵形式（按分量顺序排列）： 用矩阵表达式写成： 其中矩阵[KE]6×6 就是单元 在局部坐 标系里的单元刚度矩阵。 二、局部坐标和整体坐标 1．局部坐标：如右图杆e（任意杆，左端点号i，右端点号j），取杆左端点为原点，杆轴线为轴，从轴正方向逆时针旋转90o得轴。 2．整体坐标：平面内取一计算较方便的点作为参考原点，一般以水平线为X轴，竖直线为Y轴。 三、单元杆端力与杆端位移 任意杆件的杆端力与杆端位移分量表示如局部坐标系图与总体坐标系图。 1．元素顺序： 杆端位移分量顺序：左、右端杆 件位移分别按方向位移，方向位移与转 动角（局部坐标）或X方向位移u ，Y 方向位移v与转动角（总体坐标）进行 排列。 杆端力分量顺序：与杆端位移分 量排列顺序对应，左、右端杆端力分别 按方向分力， 方向分力与转动力矩局 部坐标）或X方向分力FX，Y方向分力 FY与转动力矩M（总体坐标）进行排列。 局部坐标下的分量列阵表达式： 杆端力 ； 杆端位移 总体坐标下的分量列阵表达式： 杆端力 ； 杆端位移 四、单元杆端力与杆端位移的坐标转换 如图坐标系（局部 ，总体xoy），杆e 与X轴夹角为α，根据力的分解与合成，有： 写成矩阵形式为： 例题 五、整体坐标单元刚度矩阵 §1.3总 体 刚 度 矩 阵 的 集 成 一、单元刚度矩阵的扩大矩阵（全元素） 例子刚度矩阵 练习 二、集成总刚的一般概念 三、全元素存放集成总刚的程序片断 程序片断 引入支撑条件： 划去与位移为零对应的行和列 “主1，副0” 法 单元荷载的转化 四、简单平面杆系分析主程序 单元刚度矩阵的性质和特点 a.对称矩阵 b.奇异矩阵 c.单元刚度矩阵的主对角元素都大于零 d.单元刚度矩阵与单元的材料，长度和截面 整体刚度矩阵的性质和特点 对称性 奇异性 矩阵中存在大量的零子块 1、总框图 输入基本数据 各单元分析 建立总刚矩阵 形成荷载列阵 求各杆单元内力 解方程求结点位移 输出内力与位移 结 束 边界条件处理 不同坐标系中的力与位移转换图 不同坐标系中的力与位移转换图 可表达为： 其中矩阵 [T]6×6就称为杆单