平面桁架的计算原理及编制解剖.pptVIP

计算杆件内力： 其中 杆件的伸长量为 表示为矩阵形式 杆单元内力为 计算杆件内力的流程表示如下： {已知总位移向量DP} {开设局部变量} BEGIN {调出STIF(NT,L,C,B,RD);} {调出I0J0(NT,II);} {计算杆件内力} END; 计算第NT号杆件的内力过程是： SUBROUTINE INTFORCE(NDISP,NT,DP,T) USE PDATA INTEGER NT,I,J,II(4) REAL DP(120),C(4,4),B(4),T,RD,U ! GIVE NT – THE ELEMENT_NO, RETURN TENSION. DP IS GLOBAL DISP CALL STIF(NT,L,C,B,RD) ! GET RD AND B CALL I0J0(NT,II) ! GET ELEMENT DISP NUMBERING VECTOR II T=0.0 ! SET 0.0 BEFORE ACCUMULATE DO 10 I=1,4 J=II(I) IF ((J .GT. 0).AND.(J .LT. NDISP)) THEN U=DP(J) ! THE RESPECTIVE DISP DEPENDS ON IS IT FIXED OF NOT ELSE U=0.0 END IF T=T+B(I)*U ! ACCUMULATE THE ENLONGATION OF BAR 10 CONTINUE T=T*RD RETURN END 计算和输出杆件内力： WRITE(2,*) BAR HL HR EF L TENSION‘ DO 90 NT=1,IESG L=BAR_TENSION(NT) CALL INTFORCE(NDISP,NT,DP,BAR_TENSION(NT)) WRITE(2,*) NT,BAR_HL(NT),BAR_HR(NT),BAR_EF(NT),L, BAR_TENSION(NT) 90 CONTINUE 本章结束！ 矩阵位移法求解一般步骤： （1）结构计算简图。（节点、单元编号，建立一个统一的坐标系等） （2）分析节点位移的力学特性，确定位移未知数。 （3）建立每根杆件两端位移和内力的关系。（单根杆件的刚度矩阵） （4）根据每根杆件上的上述关系建立结构可动节点的平衡方程。（结构总刚度矩阵） （5）求解平衡方程，得到节点位移。 （6）根据求得的位移，利用每根杆件位移与内力关系计算各杆的内力。 程序流程 平面桁架计算程序的总功能表示 {依据平面桁架及计算模型的有关数据，进行结构静力计算，并输出计算结果} 根据平面桁架的特性进一步展开为： {总体数据结构的设计} {输入描述平面桁架本体结构和性质的数据} {输出结构图形} {计算各杆件的单元刚度矩阵} {形成结构总刚度矩阵} {总刚度矩阵三角化} {输入荷载数据，形成右端项总外力向量} {回代求解，求出总位移向量} {打印位移} {计算各杆内力} {输出各杆内力} 根据平面桁架的特性进一步展开为： {总体数据结构的设计} {输入描述平面桁架本体结构和性质的数据} {输出结构图形} {计算各杆件的单元刚度矩阵} {形成结构总刚度矩阵} {总刚度矩阵三角化} {输入荷载数据，形成右端项总外力向量} {回代求解，求出总位移向量} {打印位移} {计算各杆内力} {输出各杆内力} {总体数据结构的设计} 平面桁架程序的总体数据结构设计，包括全局量（标识符）的存储安排与说明。 一个平面桁架结构的数据有三个方面： 结构本体描述数据（NW, IESG, NU, X, Y, HL, HR） 性质数据（EF） 荷载数据（PX, PY） NW为节点总数 IESG杆件总数 NU可动节点总数 X, Y 节点坐标 HL, HR 每根杆件两端节点编号 EF 性质数据 PX, PY 外载荷数据 MODULE PDATA REAL, PUBLIC :: BAR_HL(100),BAR_HR(100), BAR_TENSION(100),BAR_EF(100) REAL, PUBLIC :: NODP_X(80),NODP_Y(80) END MODULE REAL B(4),DP(120),R(1200),C(4,4) INTEGER II(4),V(120),NW,IESG,NU {总体数据结构的设计} 4 ‘ NW 4 ‘ IESG 2 ‘ NU 0 1