计算结构力学程序..doc

计算结构力学编程大作业 时间： 2007年6月 !!!**************************************************************************** !!! 关于程序的说明 !!!**************************************************************************** !一、功能： ! 1、可计算包括节点力，一般非节点力，支座沉降、温度荷载作用、制造误差的平 ! 面桁架、梁、刚架及其组合结构的节点位移与杆端力； ! 2、可同时计算多种工况下的节点位移与杆端力。 !***************************************************************************** !****************************************************************************** ! ! 二、变量说明： ! NE——单元数； ! N——结构中自由度数； ! NJ——节点数； ! NS——特殊节点数，包括支座节点、主从节点（1节点不做主节点）、连接桁架的铰节点（没有转角）； ! NAI——结构的单元截面类型数； ! MT——单元截面类型号； ! NL——荷载工况数； ! H——截面高度； ! E——弹性模量； ! JC——单元定位向量数组； ! X(NJ),Y(NJ)——节点的X，Y坐标值； ! JE(NE,2)——单元两端节点码数组； ! AI(NAI,2)——按单元类型顺序存放A与I，AI（I，1）—第I类单元的截面积，AI（I，2）—第I类单元的 ! 惯性矩； ! MT(NE)——单元所属单元类型号； ! JS(NS,4)——特殊节点信息，JS(I,1)—结点码；JS(I,2)，JS(I,3)，JS(I,4)—U，V，CETA约束信息， ! 有约束为1，没有约束为0；从节点某位移同主节点时位移时，该位移约束信息填主节点码； ! ! PJ(NP,3)——节点荷载信息数组；PJ（I，1）—节点力所在节点号；PJ（I，2）—节点力作用坐标方向： ! 坐标方向U，V，M分别为1，2，3； PJ（I，3）—节点力的大小(含正负号)；U，V方向集中力时， ! 与坐标轴正向同向为正，M按右手法则为正；本程序推导过程取y轴向下为正。 ! ! PF(NF,4)——非节点荷载数组，并给出以下类型说明： ! 前6类型数据输法(梯形等可以用叠加法计算)： ! PF(I,1)-单元码；PF(I,2)-类型；PF(I,3)-荷载大小；PF(I,4)-c值； ! 1——垂直于单元的均布力,大小为q，以坐标轴正向为正,c为荷载末端距i节点距离; ! 2——非节点集中力P，c为荷载距i节点距离; ! 3——非节点集中力距M，c为荷载距i节点距离，右手法则判正负; ! 4——三角形荷载，c为荷载距i节点距离，i端为0，距离i端c时力为q; ! j端为0的三角形，可按叠加法处理。 ! 5——沿杆轴向均布力,大小为q,c为荷载末端距i节点距离; ! 6——沿杆轴向集中力,大小为q,c为荷载末端距i节点距离; ! ! 从第7到第9类型（支座沉降）数据输法：PF(I,1)-单元码；PF(I,2)-类型；PF(I,3)-位移大小（含正负），坐标轴正向位为正，转角按右手法则；PF(I,4)-沉降所在的单元位移分量，i端为1-3，j端为4-6； ! ! 7——沿轴向支座沉降； ! 8——垂直于轴向支座沉降； ! 9——支座转动； ! 10——制造误差,PF(I,1)—制造误差所在单元，PF(I,2)-类型；PF(I,3)-误差大小(含正负),正负取决于消除 ! 误差时端点的运动方向，PF(I,4)—误差所在坐标号； ! 11——