第五节支座位移和温度变化问题的处理.ppt

第五节 支座位移和温度变化问题的处理 （1）将支座结点位移转换成与该结点连接的各单元在单元坐标系中的杆端位移。 将得到的每个杆端位移分量视为一个 PF(I,3)── PF(I,4)──由支座位移转换成的单元杆端位移的方位代码。 （2）将单元视为两端固定梁，按表3- （3）利用式（2-27）( ) （二）方法二 1.增设(118页) 2.增设变量和数组的输入处 4.引入支承条件的方法 二.温度变化问题的处理 (1)对于发生温度变化的单元，视为有一个非结点荷载 (2)将单元视为两端固定梁，按表3-3(b)中的计算 关于温度变化单元截面高度h 三.修改后的子程序EFF * * 一.支座位移的处理 （一）方法一 将支座位移转化为等效结点荷载处理。 上一节 返回上级菜单 该方法可用于先处理法和后处理法程序中，其具体做法如下 ： 例如： 例如：图(a)所示结构，已知支座结点4下沉2cm，将此支座位移分解成单元③始端沿 、 轴方向的位移， 如图(b)所示 即 其中： PF(I,1)──第I个非结点荷载所在的单元编号； PF(I,2)──第I个非结点荷载的类型码，见表3-3(a)； 非结点荷载（广义荷载），计入变量NF（非结点荷载的个数）中，有关信息仍存放在非结点荷载信息数组PF(NF,4)中， 由支座位移转换成的单元在单元坐标中的杆端位移值。其正负号规定：杆端线位移分量u、v与 、 轴同向为正，杆端转角θ顺时针方向为正。 单元始端i位移分量 、 、 的方位代码分别为 1 、2 、3 单元始端j位移分量 、 、 的方位代码分别为 4 、5 、6 例如： 如图(b)所示，对于 PF(I,1)=3， PF(I,2)=7， PF(I,3)=-0.016m， PF(I,4)=1 对于 PF(I,1)=3， PF(I,2)=8， PF(I,3)=-0.012m， PF(I,4)=2 3(a)中的计算公式，求得由给定杆端位移产生的单元固端力{FO}（e） 因此，在计算单元固端力的子程序EFF 中要增加三种非结点荷载类型，类型号分别为： 7（杆端轴向位移） 8（杆端横向位移） 9（杆端角位移） 参见表3-3(a) 求出单元等效结点荷载{PE}(e) ，并利用单元定位向量把它叠加到结构的综合结点荷载列阵{P}中去。 支座位移作为支承条件，用乘大数法来处理。 此法仅适用于后处理法。 为使后处理法程序PFL能直接计算支座位移问题，可作如下修改： 1）变量NC──已知的支座位移分量的个数 2）数组WY(NC,3) 其中 WY(I,1)──第I个支座位移分量所在的结点编号； WY(I,2)──第I个支座位移分量的方位代码，u、 v、θ的方位代码分别为1、2、3； WY(I,3)──第I个支座位移分量的代数值，u、v 与x、y轴正向一致时为正，θ以顺时 针方向为正。 3.其他相关修改 1)哑元表 由于增设了变量NC和数组WY，有关子程序的哑元表也应作相应修改。 2)数组WY的维数说明 在主程序和有关子程序中，对数组WY的维数作说明。 变量NC在主程序中输入； 数组WY的数据在子程序INPUT中输入。 用乘大数法引入支承条件的框图 J=JR(I,1) K=1,3 JR(I,K+1)=0 是 否 I=1,NC I=1,NR L=3*(J-1)+K TK(L,L)=TK(L.L) * 1.0E15 J=INT(WY(I,1)) K =INT(WY(I,2)) L=3 *(J-1)+K P(L)=TK(L,L) * WY(I,3) NC=0 是 否 对支座结点循环 取支座结点号 对支座结点的三 个位移分量循环 判断某位移分 量是否被约束 被约束的位移分量编号 总刚主元乘以大数1.0E15 判断是否有支 座位移问题 对已知的支座 位移分量循环 发生支座位移的结点号 已知支座位移分量方位代码 已知支座位移分量的编号 修改{P}中的相应分量 （119页图3-23） ──乘大数法 将温度变化作用转化为结构的非结点荷载（广义荷载）处理。 如表3-3(b)所示，其类型号为10。 该方法可用于先处理法和后处理法程序中，其具体做法如下： t2t1 PF(I,1) ──其数据形式为：×××.1 其中小数点前的数字是第I个非结点荷载（温度变化）作用所在的单元编号，小数点后的数字1表示温度变化作用，程序将自动处理为第10 类非结点荷载计算单元固端力； PF(I,2) ──温度变化值t1； PF(I,3) ──温度变化值t2； PF(I,4) ──温