平面桁架程序计算原理及程序编制.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