有限元第四章整体分析.pptVIP

哈尔滨工业大学 土木学院 王焕定 第四章 杆系结构的整体分析 4.1坐标转换 4.2结构整体刚度方程 结构整体刚度方程性质 整体分析的物理实质 4.3边界条件处理 4.4单元内力的计算 4.5程序调试中关键量的速算方法 4.6杆系结构静力分析程序简介 ★ 3) 关于斜边界的处理 如图示意的斜支座情况，有多种处理方案。 3-1) 通过单元的坐标转换来处理 x y 图示有斜支座单元，r 结点处 以倾角? - ? 来进行坐标转换，也 即在r 结点处整体坐标为图示 xy 。 r 3-2) 通过增加一个单元来处理 图示有斜支座单元，r 结点处沿 y 方向增加一个 刚结的单元，此单元有“无穷大”的抗拉刚度、但 没有抗弯刚度。单元长度可任意。 3-3) 对整体刚度矩阵进行处理（参见教材P.83） 4.3 边界条件的处理 4.4 单元内力的计算 4.4.1 杆端内力计算公式 解方程的结果可以得到结点位移，有了位移就 可以进一步求各单元的内力。 解算步骤： 从[?]根据定位向量取出单元杆端位移。 由单元倾角确定是算 还是 。 减去等效结点荷载或加上固端力矩阵。 4.4 单元内力的计算 4.4.2 杆中任意截面内力计算公式 注意事项： 由图示隔离体图，列杆段的平衡方程即可得到 任意截面的内力计算公式。请自行写出。 按上述隔离体图所求得的内力是精确的。 求得单元结点位移后，代入单元位移场、求应 变、求内力，这样所得的内力一般不满足平衡条 件。只是近似解。 4.5 程序调试中关键量的速算方法 4.5.1 总刚度矩阵元素的确定 1) 总刚度矩阵元素的物理意义 整体刚度方程为 如果 ，则可见 总刚度矩阵元素 的物理意义为：当且仅当 时，在 处所需施加对应于 的广义力。或理解 为：当且仅当 时，在限制 位移的约束上 所产生的约束反力。 2) 指定总刚度矩阵元素的速算方法 注意：（1）要牢记总刚元素的物理意义。 （2）仅仅产生 。 （3）实质上这里纯粹是用位移法来求解。 4.5 程序调试中关键量的速算方法 根据总刚元素 的物理意义，令仅仅产生 , 利用位移法中的形常数作弯矩图，象位移法一样即 可求得指定总刚元素值，为校核总刚集成的正确性 提供测试数据。 举例 4.5.2 综合等效荷载元素的确定 1) 综合等效荷载元素的组成 综合等效荷载为 也即，它由直接结点荷载和等效结点荷载组成。 2) 综合等效荷载元素的确定 直接结点荷载只需将外荷载坐标方向投影即可， 因此关键是确定等效结点荷载。 由单元分析可知，单元刚度方程和象位移法用叠 加所得力-位移关系（转角位移方程）一样，因此 4.5 程序调试中关键量的速算方法 只要熟记载常数，即可将单元荷载转化为单元等效 结点荷载，再经过往坐标方向的投影，即可获得作 用在结点的 元素。由此和直接结点荷载相加 即得到需求的综合等效荷载元素。 注意：（1）坐标正向的“荷载”为正。 （2）建议先按载常数确定固端力的实际方向和数值, 然后反方向得到等效荷载实际方向(局部坐标方向)。 （3）将所有单元荷载的等效荷载作用到结点，同时 考虑直接结点荷载（斜杆需投影）即可得需求值。 4.5 程序调试中关键量的速算方法 举例 4.5.3 单元杆端内力元素的确定 1) 单元杆端位移的确定 整体刚度方程求解结果，所得到的是整体坐标下 的结点位移，为求单元杆端内力，需作两件事： 从整体位移矩阵中根据定位向量驱除单元位移； 将整体坐标的位移往单元局部坐标方向投影。 这样即可获得单元局部坐标下的单元杆端位移。 2) 由单元刚度方程来求 “固端力” 4.5 程序调试中关键量的速算方法 如果只需求某指定内力，实际并不需要作整个矩 阵乘。 按杆端力方程求要作矩阵运算，为避免它，可在 获得局部坐标位移后，利用形、载常数通过叠加来 得到某指定内力值。 3) 由形、载常数叠加来求 注意：（1）如果要求整体坐标下的内力该怎麽办？ （2）“固端力”符号规定和位移法有区别。 （3）建议用叠加法求。 4.5 程序调试中关键量的速算方法 最大 半带宽 ** 杆系结构整体分析其他内容 1） 总刚度矩阵的存储与对应解法 因为总刚度矩阵对称、带状稀疏，利用这一特 点可减少存