板杆式薄壁结构-有限元大作业(共23页)

1、南京航空航天大学有限元大作业报告学院：航空宇航学院题目：板杆式薄壁结构组号：15组组员： 日期：2013.10一、问题重述15、板杆式薄壁结构PPPP4aPP8a已知：P=2000N，a=30cm，各杆截面积f=10mm2，壁板厚t=0.1cm材料参数均相同：E=7×106N/cm2，m=0.3求：各杆轴力和板剪流二、初步分析1、本题中所给结构为杆板式薄壁结构，并含有开口，一般用于飞行器机翼结构设计上，可以从一般的板杆结构考虑，做有限元结构分析，同时亦可保持一定准确性。2、该结构为对称结构，受力对称，由此我们可以初步确定该结构上下各相应单位受力相同，方向相反。同此我们也可据此考查所得

2、结果精确与否。3，结构含有两种单元，完成单元分析后，我们先写出杆的刚度矩阵，再累加上板的刚度矩阵，从而得到总刚度矩阵，最后带入边界条件和受力情况用改进平方根法求解编写计算模型。 4，计算模型的检验，以及求解。三、建立计算模型应用直接刚度法对结构进行分析，主要包括以下几方面的工作：1、结构离散化建立有限元的计算模型；2、单元分析建立单元刚度矩阵Ke；3、结构分析形成结构总刚度矩阵"K" ；4、采用“置大数法”对结构进行约束处理；5、利用改进平方根算法求解线性方程组，求得节点位移；6、根据求得的节点位移求出各单元的内力。四、 计算模型设计开始输入控制信息(NP,NE,NL,NR

3、,ND,NF)N=NF×NP NDF=ND×NF输入原始数据形成LD(N)数组调用FLD子程序计算局部坐标下的单元刚度矩阵坐标变化矩阵总体坐标下的单元刚度矩阵调用FKE子程序计算局部坐标下的单元刚度矩阵坐标变换矩阵计算单元IS(NDF)数组形成单元的IS(NDF)数组（调用FIS子程序）取出单元在总体坐标下的节点位移单元刚度矩阵向总刚度结构矩阵叠加（调用FASUM子程序）计算在局部坐标下单元节点位移进行约束处理调用FR子程序计算局部坐标下单元节点力和单元内力输出单元内力求解线性方程组求得结构节点位移结束输出节点位移五、模型符号数据说明1、 节点编号：共44个节点2、 杆的编

4、号：共56根杆3、 板的编号：共28块板六、模型原理1标识符说明【1】控制信息NP：结构节点总数；NE：结构离散杆单元总数；NEE：结构纯剪板单元总数；NL：载荷组数；ND：每个单元的节点数；NF：每个节点的自由度数；N=NF x NP ：结构的节点位移总数；NDF=NF x ND ：每个单元的节点自由度数；【2】输入原始数据1)节点坐标X (NP) ：节点的x方向坐标；Y (NP) ：节点的y方向坐标；2)单元的特征数据ME (NDNE) ：单元节点的总体编号； AE (NE) ：单元横截面面积；E0(NE) ：单元材料弹性模量；3)边界约束数据NRR (NR) ：约束为零的位移所对应的总体

5、位移编号；P (N，NL) ：NL组载荷向量；程序中的其他主要标识符LD (N) ：存放结构总刚所有主对角元素在一维数组A (NN)中的序号； IS (NDF) ：单元节点位移和节点力在结构节点位移列阵和节点载荷列阵中对应的序号；EK (NDF，NDF) ：总体坐标下的单元刚度矩阵；A (NN) ：结构刚度矩阵下三角变带宽一维压缩存储的数组； NN ：数组A的元素总数。【3】输出结果1)结构的节点位移；2)单元内力；2输入说明由于我们题目的单元数比较多，所以我们采用程序自动输入，没有输入文件，我们所有的杆单元的编号都是小号在前3， 源程序i.由于我们的程序是程序自动输入信息，所以没有输入文件，

6、所以我们的算例是在程序内部人为切换的；ii.我们的刚度矩阵是采用一维存储，从第一个非零元素开始存储；iii.约束处理，采用置大数法； !*!*平面薄壁结构静力分析*!*PROGRAM MAIN!*本次课设题目*!由于本结构单元比较多，不宜采用文件读取输入，故本程序带有专门的节点，单元信息输入程序!INTEGER: ME(2,56),MEP(4,28),LD(88)=0,NRR(100),NN,N,NP,NE,ND,NF,NDF,IE,NEE,NR,ISW,NL,ME1(4)INTEGER,ALLOCATABLE:IS(:),ISS(:)REAL,ALLOCATABLE:A(:)REAL X(4

7、4),Y(44),AE,EK(4,4),E0,G0,V,EKP(8,8),H0,A0,U(88),P(88),AA(8),C,Q(28),EKE(2,2),BL,R(2,4),S(56),CI(2,1),T(4)!*输入结构信息*OPEN(8,FILE='INPUTDATA.txt')READ(8,*) NP,NE,NF,NEE,NR,NL !输入结构信息READ(8,*) AE !输入杆的截面积READ(8,*) E0 !输入弹性模量READ(8,*) V !输入vREAD(8,*) H0 !输入板的厚度READ(8,*) BL !输入杆的长度N=88;G0=E0/(2*(1

8、+V)CALL IMPORTXY(X,Y); !输入节点坐标信息（详见xy.txt）CALL INPUTME(ME); !输入杆的信息（详见STICKSME.txt）CALL PLATE(MEP); !输入板的信息（详见PLATE.txt）CALL FLD(LD,NN,MEP,N,NP,NEE,4,NF)!输入LD数组（详见LD.txt）CALL INPUTNRR(NRR,NR) !输入约束信息P=0CALL INPUT(P,N) !输入外载荷信息!P(1)=1000!*本次课设题目*!*将所有信息输入到OUTDATA.txt中*OPEN(7,FILE='OUTDATA.TXT'

9、;)WRITE(7,200) NP,NE,NEE,NL,NR,NF200 FORMAT(5X,'*输入信息*'/,8X,'NP=',I4/,8X,'NE=',I4/,8X,'NEE=',I4/,8X,'NL=',I4/,8X,'NR=',I4/,8X,'NF=',I4/)WRITE(7,201) AE,H0,BL,E0,G0201 FORMAT(5X,'*结构属性信息*'/,8X,'面积: AE=',F4.2/,8X,'壁板厚: H0='

10、;,F4.2/,8X,'杆长: BL=',F4.1/,8X,'弹性模量: E0=',F9.1/,8X,'剪切模量: G0=',F9.1/)WRITE(7,202)202 FORMAT(5X,'*节点载荷*'/,5X,'序号',8X,'PX',9X,'PY'/)WRITE(7,203) (I,P(2*I-1),P(2*I),I=1,NP)203 FORMAT(5X,I3,3X,2F10.2/)WRITE(7,204)204 FORMAT(5X,'*约束信息*'/,5X,

11、'序号',8X,'NRR'/)WRITE(7,205) (I,NRR(I),I=1,NR)205 FORMAT(5X,I3,8X,I3)WRITE(7,206)206 FORMAT(5X,'*节点坐标*'/,5X,'序号',8X,'X',9X,'Y'/)WRITE(7,207) (I,X(I),Y(I),I=1,NP)207 FORMAT(5X,I3,3X,2F10.2)WRITE(7,208)208 FORMAT(5X,'*杆的编号*'/,5X,'序号',8X,&#

12、39;ME1',9X,'ME2'/)WRITE(7,209) (I,ME(1,I),ME(2,I),I=1,NE)209 FORMAT(5X,I3,3X,I8,4X,I8)WRITE(7,210)210 FORMAT(5X,'*板的编号*'/,5X,'序号',8X,'MEP1',8X,'MEP2',8X,'MEP3',8X,'MEP4'/)WRITE(7,211) (I,MEP(1,I),MEP(2,I),MEP(3,I),MEP(4,I),I=1,NEE)211 FORMA

13、T(5X,I3,3X,I8,4X,I8,4X,I8,4X,I8)WRITE(7,212) NN212 FORMAT(5X,'*LD数组*'/,5X,'NN=',I4/,5X,'LD(N)=')WRITE(7,213) (LD(I),I=1,N)213 FORMAT(5X,5I5)ND=2NDF=ND*NFALLOCATE(IS(NDF)ALLOCATE(A(NN)WRITE(7,214)214 FORMAT(5X,'*杆刚度矩阵*'/)!*计算杆的单元刚度矩阵，并存入一维总刚中*DO I=1,NE IE=I CALL FIS(ME

14、,IS,NF,ND,NE,IE,NDF) CALL FEK(IE,NP,NE,NDF,X,Y,ME,AE,E0,EK,ND,BL) WRITE(7,"(5X,'IE=',I3/)") I WRITE(7,215) (EK(J,K),K=1,4),J=1,4)215 FORMAT(5X,4F12.2/) CALL FASUM(A,EK,IS,LD,NN,NDF,N)ENDDOOPEN(50,FILE='ANN.TXT')WRITE(50,"(2X,5F12.2)") (A(I),I=1,NN)ND=4NDF=ND*NFALL

15、OCATE(ISS(NDF)WRITE(7,216)216 FORMAT(5X,'*板刚度矩阵*'/)!*计算板的单元刚度矩阵，并存入一维总刚中*DO I=1,NEE IE=I CALL FIS(MEP,ISS,NF,ND,NEE,IE,NDF) CALL FKEP(IE,NP,NEE,X,Y,MEP,H0,G0,EKP,A0,AA,C) WRITE(7,"(5X,'IE=',I3/)") I WRITE(7,217) (EKP(J,K),K=1,8),J=1,8)217 FORMAT(5X,8F12.2/) CALL FASUM(A,EKP

16、,ISS,LD,NN,NDF,N)ENDDOWRITE(7,218)218 FORMAT(5X,'*一维刚度矩阵*'/,5X,'A(N)=')WRITE(7,219) (A(I),I=1,NN)219 FORMAT(5X,6F12.2/)OPEN(50,FILE='ANNP.TXT')WRITE(50,"(2X,5F12.2)") (A(I),I=1,NN)CALL FR(N,NN,NR,LD,A,NRR) !置大数法处理约束CALL BAND3(N,1,NN,A,P,LD,ISW) !用改进平方根法求解位移WRITE(7,2

17、20)220 FORMAT(5X,'*节点位移*'/,5X,'序号',10X,'U',10X,'V')WRITE(7,221) (I,P(2*I-1),P(2*I),I=1,NP)221 FORMAT(5X,I3,5X,F12.6,5X,F12.6/)!*计算板的剪流*DO I=1,NEE Q(I)=0.0 CALL FIS(MEP,ISS,NF,ND,NEE,I,NDF) CALL FKEP(I,NP,NEE,X,Y,MEP,H0,G0,EKP,A0,AA,C) DO J=1,8 Q(I)=Q(I)+C*AA(J)*P(ISS(

18、J)*2 ENDDOENDDOWRITE(7,222)222 FORMAT(5X,'*板的剪流*'/,5X,'序号',10X,'Q')WRITE(7,223) (I,Q(I),I=1,NEE)223 FORMAT(5X,I3,5X,F10.4/)!*计算杆的轴力*C=E0*AE/BLND=2WRITE(7,224)224 FORMAT(5X,'*杆的轴力*'/,5X,'序号',12X,'SIJ',12X,'SJI')DO I=1,NE S(I)=0.0 NI=ME(1,I) NJ=M

19、E(2,I) XIJ=(X(NJ)-X(NI) YIJ=(Y(NJ)-Y(NI) T(1)=XIJ/BL T(2)=YIJ/BL T(3)=-T(1) T(4)=-T(2) CALL FIS(ME,IS,NF,ND,NE,I,NDF) DO J=1,4 U(J)=P(IS(J) S(I)=S(I)-C*T(J)*U(J) ENDDO SIJ=0.0DO II=1,NEEKK=0DO JJ=1,4ME1(JJ)=MEP(JJ,II)IF(ME1(JJ).EQ.NI.OR.ME1(JJ).EQ.NJ) KK=KK+1 ENDDOIF(KK.EQ.2) THEN LL=-1 IF(NI.EQ.ME1

20、(1).OR.NI.EQ.ME1(3) LL=1 SIJ=SIJ+LL*Q(II)*BL/2.0ENDIF ENDDO SJI=-SIJ+S(I)SIJ=SIJ+S(I)WRITE(7,225) I,SIJ,SJI225 FORMAT(5X,I3,5X,F12.4,5X,F12.4/)ENDDOPRINT*,'结构的信息保存在OUTDATA.TXT'END!*生成LD数组*SUBROUTINE FLD(LD,NN,ME,N,NP,NE,ND,NF) INTEGER LD(N),ME(ND,NE) LD(1)=1 DO K=1,NP IG=NP+1 DO I=1,NE DO J=

21、1,ND IF(ME(J,I).EQ.K) THEN DO J1=1,ND IF(ME(J1,I).LT.IG) IG=ME(J1,I) ENDDO ENDIFENDDO ENDDO DO I=1,NFJ=NF*(K-1)+INW=NF*(K-IG)+IIF(J.NE.1) LD(J)=LD(J-1)+NW ENDDO ENDDONN=LD(N)OPEN(40,FILE='LD.TXT')WRITE(40,"(5X,'NN=',I5)") NNWRITE(40,"(5X,'LD(I)=',8I5)") (L

22、D(I),I=1,N)RETURNEND!*输入板单元节点信息*SUBROUTINE PLATE(MEP)INTEGER I,MEP(4,28)OPEN(30,FILE='PLATE.TXT')DO I=1,13 MEP(1,I)=(I-1)/4+I MEP(2,I)=MEP(1,I)+5 MEP(3,I)=MEP(1,I)+6 MEP(4,I)=MEP(1,I)+1ENDDOMEP(1,14)=19;MEP(2,14)=23;MEP(3,14)=24;MEP(4,14)=20MEP(1,15)=21;MEP(2,15)=25;MEP(3,15)=26;MEP(4,15)=22

23、MEP(1,16)=23;MEP(2,16)=28;MEP(3,16)=29;MEP(4,16)=24DO I=17,28 MEP(1,I)=(I-17)/4+I+8 MEP(2,I)=MEP(1,I)+5 MEP(3,I)=MEP(1,I)+6 MEP(4,I)=MEP(1,I)+1ENDDODO I=1,28 WRITE(30,"(I3,2X,I3,2X,I3,2X,I3)") MEP(1,I),MEP(2,I),MEP(3,I),MEP(4,I) ENDDOEND!*输入节点坐标信息*SUBROUTINE IMPORTXY(X,Y)integer i,jreal x(

24、44),y(44)open(20,file='xy.txt')x=0;y=0do i=1,20 j=(i-1)/5 K=mod(i,5) IF(K=0) K=5 x(i)=240-j*30 y(i)=120-(K-1)*30enddox(21)=120;x(22)=120;x(23)=120;x(24)=120y(21)=120;y(22)=90;y(23)=30;y(24)=0do i=25,44 j=(i)/5 x(i)=240-j*30 y(i)=120-(mod(i,5)*30enddodo i=1,44 write(20,"(f5.1,2x,f5.1)&qu

25、ot;) x(i),y(i)enddoend !*输入杆单元节点编号*SUBROUTINE INPUTME(ME)INTEGER ME(2,56)DO I=1,7IF(I.LE.4) THENME(1,I)=IME(2,I)=I+1ELSE ME(1,I)=I-3ME(2,I)=I+2ENDIFENDDODO I=8,26ME(1,I)=ME(1,I-7)+5ME(2,I)=ME(2,I-7)+5ENDDO ME(1,27)=19ME(2,27)=23ME(1,28)=21ME(2,28)=22ME(1,29)=23ME(2,29)=24ME(1,30)=22ME(2,30)=26ME(1,3

26、1)=23ME(2,31)=28DO I=32,38ME(1,I)=ME(1,I-31)+24ME(2,I)=ME(2,I-31)+24ENDDODO I=39,56ME(1,I)=ME(1,I-7)+5ME(2,I)=ME(2,I-7)+5ENDDOOPEN(24,FILE="STICKSME.TXT")DO I=1,56WRITE(24,"(I3,2X,I3)") ME(1,I),ME(2,I)ENDDOEND!*形成IE单元的自由度编号IS*SUBROUTINE FIS(ME,IS,NF,ND,NE,IE,NDF)INTEGER ME(ND,NE)

27、,IS(NDF)DO I=1,ND DO J=1,NF II=NF*(I-1)+J IS(II)=NF*(ME(I,IE)-1)+J ENDDOENDDORETURNEND!*矩阵相乘*SUBROUTINE JZH(A,B,C,M1,N1,N2)REAL: A(M1,N1),B(N1,N2),C(M1,N2),I,J,KC=0DO I=1,M1 DO J=1,N2 DO K=1,N1 C(I,J)=C(I,J)+A(I,K)*B(K,J) ENDDO ENDDOENDDOEND!*总体坐标系下的单元刚度矩阵*SUBROUTINE FEK(IE,NP,NE,NDF,X,Y,ME,AE,E0,EK

28、,ND,BL)INTEGER: NP,NE,ND,NDF,ME(ND,NE),NI,NJREAL: K(2,2),C,EK(NDF,NDF),AE,L(4,4)REAL: X(NP),Y(NP),R(2,4),RT(4,2)=0,BLNI=ME(1,IE)NJ=ME(2,IE)R(1,1)=(X(NJ)-X(NI)/BLR(1,2)=(Y(NJ)-Y(NI)/BLR(2,3)=(X(NJ)-X(NI)/BLR(2,4)=(Y(NJ)-Y(NI)/BLRT(1,1)=(X(NJ)-X(NI)/BLRT(2,1)=(Y(NJ)-Y(NI)/BLRT(3,2)=(X(NJ)-X(NI)/BLRT(4

29、,2)=(Y(NJ)-Y(NI)/BLC=AE*E0/BLK(1,1)=C;K(1,2)=-C;K(2,1)=-C;K(2,2)=CCALL JZH(RT,K,L,4,2,2)CALL JZH(L,R,EK,4,2,4)END!*实现总刚*SUBROUTINE FASUM(A,EK,IS,LD,NN,NDF,N)INTEGER IS(NDF),LD(N)REAL EK(NDF,NDF),A(NN)DO I=1,NDF DO J=1,NDF NI=IS(I) NJ=IS(J) IF(NJ.LE.NI) THEN IJ=LD(NI)-(NI-NJ) A(IJ)=A(IJ)+EK(I,J) ENDI

30、FENDDOENDDORETURNEND!*计算板的剪流*SUBROUTINE FKEP(IE,NP,NEE,X,Y,MEP,H0,G0,EKP,A0,AA,C)REAL X(NP),Y(NP),AA(8),EKP(8,8),H0,CINTEGER MEP(4,NEE),M1(4)DO II=1,4 M1(II)=MEP(II,IE) ENDDOI0=M1(1)J0=M1(2)L0=M1(3)M0=M1(4)AA(1)=-(X(L0)-X(I0)AA(2)=-(Y(L0)-Y(I0)AA(3)=X(M0)-X(J0)AA(4)=Y(M0)-Y(J0)DO JJ=1,4 AA(JJ+4)=-AA

31、(JJ) ENDDOA0=ABS(AA(1)*AA(4)-AA(2)*AA(3)/2.0 !A0:梯形的平面面积C=G0*H0/(4*A0)DO I=1,8 DO J=1,8 EKP(I,J)=C*AA(I)*AA(J) ENDDOENDDORETURNEND!*约束处理*SUBROUTINE FR(N,NN,NR,LD,A,NRR)INTEGER LD(N),NRR(NR)REAL A(NN)DO I=1,NR N1=NRR(I) NI=LD(N1) A(NI)=1.0E35ENDDORETURNEND!*输入NRR*SUBROUTINE INPUTNRR(NRR,NR)INTEGER NR

32、,NRR(NR),IDO I=1,5 NRR(I-1)*2+1)=5*(I-1)*2+1 NRR(I-1)*2+2)=5*(I-1)*2+2ENDDODO I=6,9 NRR(I-1)*2+1)=(25+(I-6)*5-1)*2+1 NRR(I-1)*2+2)=(25+(I-6)*5-1)*2+2ENDDODO I=10,13 NRR(I-1)*2+1)=(5+(I-10)*5-1)*2+1 NRR(I-1)*2+2)=(5+(I-10)*5-1)*2+2ENDDODO I=14,18 NRR(I-1)*2+1)=(24+(I-14)*5-1)*2+1 NRR(I-1)*2+2)=(24+(I

33、-14)*5-1)*2+2ENDDOEND!*解线性方程组*SUBROUTINE BAND3(N,M,NN,A,B,LD,ISW) DIMENSION A(NN),B(N,M),LD(N) DO 20 I=1,N LDI=LD(I) IF(I .EQ. 1) GO TO 10 I0=LDI-I IM1=I-1 MI=1-I0+LD(IM1) IF(MI .GT. IM1) GO TO 10 DO 8 J=MI,IM1 J0=LD(J)-J IJ=I0+J IF(J.EQ. 1) GOTO 6 JM1=J-1 MJ=1-J0+LD(JM1) MIJ=MAX0(MI,MJ) IF(MIJ .GT.

34、 JM1) GO TO 6 S=0.0 DO 2 K=MIJ,JM1 IK=I0+K JK=J0+K 2 S=S+A(IK)*A(JK) A(IJ)=A(IJ)-S 6 DO 8 K=1,M 8 B(I,K)=B(I,K)-A(IJ)*B(J,K) 10 ALDI=A(LDI) IF(I .EQ. 1 .OR. MI .GT. IM1) GO TO 13 DO 12 J=MI,IM1 IJ=I0+J LDJ=LD(J) S=A(IJ) A(IJ)=S/A(LDJ) 12 ALDI=ALDI-S*A(IJ) 13 A(LDI)=ALDI IF(ALDI .NE. 0.0) GO TO 14 IS

35、W=0 RETURN 14 DO 20 K=1,M 20 B(I,K)=B(I,K)/ALDI DO 70 LDI=2,N I=N-LDI+2 I0=LD(I)-I MI= 1-I0+LD(I-1) IM1=I-1 IF(MI.GT.IM1) GO TO 70 DO 30 J=MI,IM1 IJ=I0+J DO 30 K= 1 ,M 30 B(J,K)=B(J,K)-A(IJ)*B(I,K) 70 CONTINUE ISW=1 RETURN END !*输入载荷* SUBROUTINE INPUT(P,N)INTEGER NL,NREAL: P(N)P=0 DO I=2,4 P(2*I-1)=

36、2000 ENDDO DO I=41,43 P(2*I-1)=-2000 ENDDO END七、模型计算结果输入文本：44 56 2 28 36 6 0.170000000.30.130输出文本：（部分）*输入信息* NP= 44 NE= 56 NEE= 28 NL= 6 NR= 36 NF= 2*结构属性信息* 面积: AE=0.10 壁板厚: H0=0.10 杆长: BL=30.0 弹性模量: E0=7000000.0 剪切模量: G0=2692307.8*LD数组* NN=1008 LD(N)= 1 3 6 10 13 17 20 24 27 31 42 54 67 81 94 108

37、121 135 148 162 173 185 198 212 225 239 252 266 279 293 304 316 329 343 356 370 383 397 410 424 435 447 460 474 483 493 504 516 525 535 546 558 561 565 576 588 601 615 626 638 651 665 678 692 705 719 732 746 757 769 782 796 809 823 836 850 863 877 888 900 913 927 940 954 967 981 994 1008*节点位移* 序号 U

38、V (cm) 1 0.000000 0.000000 2 0.051879 0.027576 3 0.085524 0.000000 4 0.051879 -0.027576 5 0.000000 0.000000 6 0.000000 0.000000 7 0.005511 -0.013434 8 0.015616 0.000000 9 0.005511 0.013434 10 0.000000 0.000000 11 0.000000 0.000000 12 -0.002485 -0.011726 13 -0.009225 0.000000 14 -0.002485 0.011726 15

39、 0.000000 0.000000 16 0.000000 0.000000 17 -0.004414 -0.002768 18 -0.007015 0.000000 19 -0.004414 0.002768 20 0.000000 0.000000 21 0.000000 0.000000 22 0.000000 0.001403 23 0.000000 -0.001403 24 0.000000 0.000000 25 0.000000 0.000000 26 0.004414 -0.002768 27 0.007015 0.000000 28 0.004414 0.002768 29

40、 0.000000 0.000000 30 0.000000 0.000000 31 0.002485 -0.011726 32 0.009225 0.000000 33 0.002485 0.011726 34 0.000000 0.000000 35 0.000000 0.000000 36 -0.005511 -0.013434 37 -0.015616 0.000000 38 -0.005511 0.013434 39 0.000000 0.000000 40 0.000000 0.000000 41 -0.051879 0.027576 42 -0.085524 0.000000 4

41、3 -0.051879 -0.027576 44 0.000000 0.000000 *板的剪流* 序号 Q 1 -73.4997 2 -12.2935 3 12.2935 4 73.4997 5 -21.2407 6 -22.7590 7 22.7590 8 21.2407 9 -9.2384 10 1.7191 11 -1.7191 12 9.2384 13 1.0913 14 -1.0913 15 -1.0913 16 1.0913 17 9.2384 18 -1.7191 19 1.7191 20 -9.2384 21 21.2407 22 22.7590 23 -22.7590 24 -21.2407 25 73.4997 26 12.2935 27 -12.2935 28 -73.4997 *杆的轴力* 序号 SIJ SJI 1 -1745.9437 459.0470 2 459.0465 827.8501 3 827.8503 459.