資源描述：
《矩陣位移法教學(xué)資料》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、矩陣位移法教學(xué)資料一單元?jiǎng)偠染仃囆纬沙绦?單元?jiǎng)偠染仃嚫髟氐乃闶綄?duì)稱將局部坐標(biāo)系下的單元?jiǎng)偠染仃?、坐?biāo)變換矩陣代入變換式，得結(jié)構(gòu)坐標(biāo)系下單元?jiǎng)偠染仃嚫髟氐乃闶剑涸谟?jì)算機(jī)程序中，單元?jiǎng)偠染仃嚫髟匾话阒苯影雌渌闶酱_定，然后用一個(gè)二維數(shù)組KE(6,6)存放，程序中：EL─單元桿件的長(zhǎng)度EF(N,3),EF(N,4)─單元的抗拉剛度EA、抗彎剛度EI2形成單元?jiǎng)偠染仃嚨淖映绦騍UBROUTINESETKE(N,EF,EL,SI,CO,KE)REALKEDIMENSIONEF(100,4),KE(6,6)A=EF(N,3)*CO*CO/EL+12*EF(N,4)*SI
2、*SI/EL**3KE(1,1)=AKE(1,4)=-AKE(4,4)=AB=(EF(N,3)/EL-12*EF(N,4)/EL**3)*CO*SIKE(1,2)=BKE(1,5)=-BKE(2,4)=-BKE(4,5)=BC=6*EF(N,4)*SI/EL**210KE(1,3)=-CKE(1,6)=-CKE(3,4)=CKE(4,6)=CD=EF(N,3)*SI*SI/EL+12*EF(N,4)*CO*CO/EL**3KE(2,2)=DKE(2,5)=-DKE(5,5)=DE=6*EF(N,4)*CO/EL**2KE(2,3)=EKE(2,6)=EKE(3,5
3、)=-EKE(5,6)=-EF=4*EL(N,4)/ELKE(3,3)=FKE(3,6)=F/2KE(6,6)=FDO10I=1,5DO10J=I+1,6KE(J,I)=KE(I,J)10CONTINUERETURNEND二結(jié)構(gòu)剛度矩陣的形成程序1程序說(shuō)明在結(jié)構(gòu)靜力分析程序中，通過(guò)調(diào)用以下子程序可以建立起結(jié)構(gòu)剛度矩陣，程序中有關(guān)變量及數(shù)組：NE─單元總數(shù)NN─結(jié)點(diǎn)位移分量總數(shù)XY(100,2)─結(jié)點(diǎn)坐標(biāo)數(shù)組ID(100,3)─結(jié)點(diǎn)位移編碼數(shù)組EF(100,4)─單元始端﹑末端結(jié)點(diǎn)編碼、EA、EI數(shù)組KE(6,6)─結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃嚁?shù)組IJ(6)─單元定位
4、向量數(shù)組KS(NN,NN)─結(jié)構(gòu)剛度矩陣數(shù)組102形成結(jié)構(gòu)剛度矩陣的子程序SUBROUTINESETKS(NE,NN,XY,ID,EF,KS)REALKS,KEDIMENSIONXY(100,2),ID(100,3),EF(100,4),IJ(6),KE(6,6),KS(NN,NN)DO40N=1,NEME=NCALLLENG(ME,XY,EF,EL,SI,CO)CALLSETKE(N,EF,EL,SI,CO,KE)IE=EF(ME,1)JE=EF(ME,2)DO10I=1,3IJ(I)=ID(IE,I)IJ(I+3)=ID(JE,I)10CONTINUEDO30
5、I=1,6IF(IJ(I)﹒EQ﹒0)GOTO30DO20J=1,6IF(IJ(J)﹒EQ﹒0)GOTO20KS(IJ(I),IJ(J))=KS(IJ(I),IJ(J))+KE(I,J)20CONTINUE30CONTINUE40CONTINUERETURNEND三建立結(jié)點(diǎn)荷載向量的程序1程序說(shuō)明在結(jié)構(gòu)靜力分析程序中調(diào)用如下子程序來(lái)建立結(jié)點(diǎn)荷載向量，程序中有關(guān)變量如下：NP—結(jié)點(diǎn)荷載總數(shù)NF—非結(jié)點(diǎn)荷載總數(shù)xy(100,2)—結(jié)點(diǎn)坐標(biāo)數(shù)組ID(100,3)—結(jié)點(diǎn)位移編碼數(shù)組EF(100,4)—單元始端、末端結(jié)點(diǎn)編碼、EA、EI值數(shù)組IJ(6)—單元定位向量數(shù)組PJ
6、(50,3)—結(jié)點(diǎn)荷載作用結(jié)點(diǎn)號(hào)、作用方向代碼、結(jié)點(diǎn)荷載值數(shù)組PE(50,4)—非結(jié)點(diǎn)荷載作用單元號(hào)、荷載類型代碼、荷載值、荷載作用點(diǎn)距始端距離數(shù)組10P(6)—單元桿端固端力數(shù)組PS(100)—結(jié)點(diǎn)荷載數(shù)組２建立結(jié)點(diǎn)荷載向量的子程序SUBROUTINESETPS(NP,NF,XY,ID,EF,PJ,PE,PS)DIMENSIONXY(100,2),ID(100,3),EF(100,4),PJ(50,3),PE(50,4),PS(100),IJ(6),P(6)IF(NP·EQ·0)GOTO20DO10N=1,NPJJ=PJ(N,1)JD=PJ(N,2)JQ=ID(
7、JJ,JD)PS(JQ)=PS(JQ)+PJ(N,3)10CONTINUE20IF(NF·EQ·0)GOTO50DO40N=1,NFME=PE(N,1)JT=PE(N,2)EP=PE(N,3)A=PE(N,4)CALLLENG(ME,XY,EF,EL,SJ,CO)CALLSETEP(JT,EP,A,EL,P)CALLTRANST(SI,CO,P)IE=EF(ME,1)JE=EF(ME,2)DO30I=1,3IJ(I)=ID(IE,I)IJ(I+3)=ID(JE,I)30CONTINUEDO40I=1,6IF(IJ(I)·EQ·0)GOTO40PS(IJ(I))=P
8、S(IJ(
