結(jié)構(gòu)力學(xué)電算實(shí)習(xí)報(bào)告及論大學(xué)生寫作能力

結(jié)構(gòu)力學(xué)電算實(shí)習(xí)報(bào)告摘要:現(xiàn)實(shí)生活中往往存在一些很復(fù)雜的連續(xù)梁模型，它們的支撐性質(zhì)和受力特點(diǎn)都比較復(fù)雜，結(jié)構(gòu)力學(xué)電算則是將它們有規(guī)律的進(jìn)行簡化，最后將其所有已知的性質(zhì)表達(dá)為一組數(shù)據(jù)，將這些數(shù)據(jù)輸入到我們所編好的程序中，既可以得到我們所要求一些量如：桿端力，各截面內(nèi)力，某界面的內(nèi)力影響線等。關(guān)鍵詞：連續(xù)梁、電算程序、截面內(nèi)力程序功能：連續(xù)梁程序可對連續(xù)梁結(jié)構(gòu)進(jìn)行相關(guān)的計(jì)算。能夠計(jì)算在任意形式均布荷載、跨中集中力荷載、跨中力偶荷載、節(jié)點(diǎn)力偶荷載組合下連續(xù)梁單元的桿端力，也可計(jì)算連續(xù)梁各截面內(nèi)力，從而繪制連續(xù)梁的彎矩圖。影響線程序可以計(jì)算單位1集中力荷載作用下的指定截面中內(nèi)力在單位1載荷按指定距離移動(dòng)情況下的變化情況，從而提取數(shù)據(jù)繪制影響線。使用說明：（1）連續(xù)梁內(nèi)力程序在連續(xù)梁計(jì)算程序（.exe文件）所在文件夾中建立一新的txt文件，并命名（1.txt）。在文件其中依次輸入單元數(shù)、支承類型、節(jié)點(diǎn)荷載個(gè)數(shù)、非節(jié)點(diǎn)荷載個(gè)數(shù)、指定截面?zhèn)€數(shù)、彈性模量、桿長&慣性矩（依次按組輸入）、節(jié)點(diǎn)荷載大小&對應(yīng)節(jié)點(diǎn)編號（依次按組輸入）、非節(jié)點(diǎn)荷載值&距該單元左端距離&單元號&荷載類型號（依次按組輸入）、所求截面所在單元號&距該單元左端距離（依次按組輸入），然后保存該文件。運(yùn)行連續(xù)梁內(nèi)力程序，按提示輸入起始文件即初始數(shù)據(jù)的文件（要求全名），再輸入結(jié)果文件名，在相應(yīng)文件夾下打開結(jié)果文件即可查看結(jié)果。（2）影響線程序在影響線程序（.exe文件）所在文件夾中建立一新的txt文件，并命名。在其中一次輸入單元數(shù)、支承類型、要求影響線的數(shù)量、彈性模量、單位載荷每次移動(dòng)的距離、桿長&慣性矩、截面所在單元號&距該單元左端距離，然后保存該文件。運(yùn)行影響線程序，輸入起始文件即初始數(shù)據(jù)的文件，再輸入結(jié)果文件名，在相應(yīng)文件夾下打開結(jié)果文件即可查看結(jié)果。提取期中的數(shù)據(jù)，利用excel或其它工具繪制影響線即可。程序注釋：支承類型：兩端鉸支為1，左固右鉸為2，左鉸右固為3，兩端固定為4。非節(jié)點(diǎn)荷載類型：均布荷載為1，集中力為2，集中力偶為3；節(jié)點(diǎn)荷載只有外力偶。非節(jié)點(diǎn)荷載作用距離：指到單元左端的距離，其中均布荷載指均布荷載結(jié)束位置到單元左端的距離。荷載正負(fù)號：力偶逆時(shí)針為正，力向上為正。算例1：求圖示結(jié)構(gòu)的彎矩圖及指定截面內(nèi)力?！咀⑨尅繄D中形如的粗線所示斷面為任意斷面，程序會(huì)求出任意斷面的剪力及彎矩根據(jù)此連續(xù)梁的計(jì)算簡圖，將其性質(zhì)及受力特征表現(xiàn)為一組數(shù)據(jù)，并輸入運(yùn)行程序進(jìn)行計(jì)算，如下：輸入數(shù)據(jù):442571365104861210344-6311-5322-43314131-1234311.52123313342.9943.01得出如下數(shù)據(jù)結(jié)果輸出結(jié)果:*************連續(xù)梁內(nèi)力計(jì)算*****************單元數(shù)=4支承類型=4節(jié)點(diǎn)荷載個(gè)數(shù)=2非節(jié)點(diǎn)荷載個(gè)數(shù)=5指定截面?zhèn)€數(shù)=7彈性模量=1.0000桿長,慣性矩GC(NE),GX(NE)3.0006.0005.00010.0004.0008.0006.00012.000節(jié)點(diǎn)荷載大小,對應(yīng)未知數(shù)序號PJ(I,1),PJ(I,2)10.0003.0004.0004.000非結(jié)點(diǎn)荷載值,距離,單元號,荷載類型號-6.0003.0001.0001.000-5.0003.0002.0002.000-4.0003.0003.0001.0004.0001.0003.0001.000-12.0003.0004.0003.000截面所在單元號，距離1.0001.5002.0001.0002.0003.0003.0001.0003.0003.0004.0002.9904.0003.010：：：：：：：：：位移：；：：：：：：：：結(jié)點(diǎn)號=1.0000結(jié)點(diǎn)號=2.0109結(jié)點(diǎn)號=3.4815結(jié)點(diǎn)號=4.5463結(jié)點(diǎn)號=5.0000.................各單元桿端內(nèi)力....................單元號=1左端彎矩=4.543右端彎矩=-4.413單元號=2左端彎矩=4.413右端彎矩=.296單元號=3左端彎矩=9.704右端彎矩=2.630單元號=4左端彎矩=1.370右端彎矩=-.815;;;;;;;;;;;;;;;;;;;;;所求截面內(nèi)力;;;;;;;;;;;;;;;;;;;;截面1彎矩=2.272剪力=-.043截面2彎矩=-1.471剪力=-2.942截面3彎矩=4.412剪力=-2.942截面4彎矩=-2.621剪力=-7.083截面5彎矩=3.546剪力=.917截面6彎矩=-7.073剪力=1.907截面7彎矩=4.888剪力=1.907======================計(jì)算結(jié)束======================根據(jù)此得出數(shù)據(jù)，繪制連續(xù)梁彎矩圖：如下：算例2：求圖示結(jié)構(gòu)1、2截面的剪力影響線。輸入數(shù)據(jù)：312114455661221輸出結(jié)果：***********連續(xù)梁彎矩剪力影響線計(jì)算*****************單元數(shù)=3支承類型=1要求影響線數(shù)量=2彈性模量=1.0000單位荷載每次移動(dòng)距離1.0000桿長,慣性矩GC(NE),GX(NE)4.0004.0005.0005.0006.0006.000截面所在單元號，距離1.0002.0002.0001.000單位力所在單元號1距離1.0000截面1彎矩=.375剪力=.313截面2彎矩=-.187剪力=-.062單位力所在單元號1距離2.0000截面1彎矩=.800剪力=-.400截面2彎矩=-.300剪力=-.100單位力所在單元號1距離3.0000截面1彎矩=.325剪力=-.163截面2彎矩=-.262剪力=-.087單位力所在單元號1距離4.0000截面1彎矩=.000剪力=.000截面2彎矩=.000剪力=.000單位力所在單元號2距離1.0000截面1彎矩=-.160剪力=.080截面2彎矩=.512剪力=-.832單位力所在單元號2距離2.0000截面1彎矩=-.200剪力=.100截面2彎矩=.216剪力=-.616單位力所在單元號2距離3.0000截面1彎矩=-.160剪力=.080截面2彎矩=.064剪力=-.384單位力所在單元號2距離4.0000截面1彎矩=-.080剪力=.040截面2彎矩=.008剪力=-.168單位力所在單元號2距離5.0000截面1彎矩=.000剪力=.000截面2彎矩=.000剪力=.000單位力所在單元號3距離1.0000截面1彎矩=.051剪力=-.025截面2彎矩=.000剪力=.102單位力所在單元號3距離2.0000截面1彎矩=.074剪力=-.037截面2彎矩=.000剪力=.148單位力所在單元號3距離3.0000截面1彎矩=.075剪力=-.038截面2彎矩=.000剪力=.150單位力所在單元號3距離4.0000截面1彎矩=.059剪力=-.030截面2彎矩=.000剪力=.119單位力所在單元號3距離5.0000截面1彎矩=.032剪力=-.016截面2彎矩=.000剪力=.065單位力所在單元號3距離6.0000截面1彎矩=.000剪力=.000截面2彎矩=.000剪力=.000======================計(jì)算結(jié)束====================== 1截面彎矩影響線1截面彎矩影響線11截面剪力影響線截面2的彎矩剪力影響線同理，不再做出。附：程序清單（程序中含有注釋）（1）連續(xù)梁內(nèi)力程序PROGRAMLXLDIMENSIONGC(20),GX(20),PJ(20,2),PF(10,4),sqjm(20,2),DK(2,2),P(45),1,F0(2),WY(2),F(2),ZK1(45),ZK2(45),gd(20,2),qdnl(20,2) character*12dat0,dat1 write(*,*)'輸入:初始數(shù)據(jù)文件名' read(*,10)dat0 write(*,*)'輸入:計(jì)算結(jié)果數(shù)據(jù)文件名' read(*,10)dat110format(A12)open(8,file=dat0,status='old') open(9,file=dat1,status='new') write(9,11)11format(8x,'*************連續(xù)梁內(nèi)力計(jì)算*****************') READ(8,*)NE,NZ,NP,NF,njm,E0 WRITE(9,12)NE,NZ,NP,NF,njm,E0 12FORMAT(1X,'單元數(shù)=',6x,I3,5x,'支承類型=',I3,5x,'節(jié)點(diǎn)荷載個(gè)數(shù)='I3,/,11x,'非節(jié)點(diǎn)荷載個(gè)數(shù)=',I3,5x,'指定截面?zhèn)€數(shù)=',i3,5x,'彈性模量=',F12.,14)NJ=NE+1CALLSRSJ(NE,NP,NF,GC,GX,PJ,PF,sqjm,njm)CALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,CALLFCQJ(NJ,ZK1,ZK2,P)CALLSCWY(NJ,P)CALLGDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,gd)calljmnl(gd,sqjm,njm,qdnl,gc,ne,nf,pf) ENDC************************************c輸入初始數(shù)據(jù)SRSJ子程序C************************************SUBROUTINESRSJ(NE,NP,NF,GC,GX,PJ,PF,sqjm,njm)DIMENSIONGC(NE),GX(NE),PJ(NP,2),PF(NF,4),sqjm(njm,2)c sqjm:所求截面njm:截面數(shù)目sqjm(njm,1)=所屬單元sqjm(njm,2)=距左端距離CGC(NE)：桿長GX(NE)：慣性矩C結(jié)點(diǎn)荷載:PJ(NP,1)=荷載大小PJ(NP,2)=對應(yīng)未知數(shù)序號C非結(jié)點(diǎn)荷載:PF(NE,1)=荷載大小PF(NE,2)=距離PF(NE,3)=單元號PF(NE,4)=荷載類型號WRITE(9,11)11FORMAT(1X,'桿長,慣性矩GC(NE),GX(NE)')READ(8,*)(GC(I),GX(I),I=1,NE) write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,2F8.3)IF(NP.GT.0)THENWRITE(9,22)22FORMAT(1X,'節(jié)點(diǎn)荷載大小,對應(yīng)未知數(shù)序號PJ(I,1),PJ(I,2)')READ(8,*)(PJ(I,1),PJ(I,2),I=1,NP) write(9,23)(PJ(I,1),PJ(I,2),I=1,NP)23 FORMAT(1X,2F8.3)ENDIFIF(NF.GT.0)THENWRITE(9,33)33FORMAT(1X,'非結(jié)點(diǎn)荷載值,距離,單元號,荷載類型號')READ(8,*)((PF(I,J),J=1,4),I=1,NF)WRITE(9,34)((PF(I,J),J=1,4),I=1,NF)34FORMAT(1X,4F8.3) ENDIF if(njm.gt.0)then write(9,35)35format(1x,'截面所在單元號，距離')read(8,*)((sqjm(i,j),j=1,2),i=1,njm) write(9,36)((sqjm(i,j),j=1,2),i=1,njm)36format(1x,2f8.3) endifENDC*******************************#***************c計(jì)算第NHF個(gè)非結(jié)點(diǎn)荷截引起的等效結(jié)點(diǎn)荷截列陣F0C***********************************************SUBROUTINEDJH(NHF,NE,NF,PF,F0,GC)DIMENSIONPF(NF,4),F0(2),GC(NE)CG—荷載值,C—作用距離,NT--單元號,ID—荷載類型號G=PF(NHF,1)C=PF(NHF,2)NT=INT(PF(NHF,3)+0.1)ID=INT(PF(NHF,4)+0.1)BL=GC(NT)D=BL-CC1=C/BLC2=C1*C1GOTO(10,20,30),IDC均布荷載作用下的等效荷載列陣F0(2)10F0(1)=G*C*C*(6.0-8.0*C1+3.0*C2)/12.0F0(2)=-G*C*C*C*(4.0-3.0*C1)/12.0/BLGOTO200C集中力作用下的等效荷截列陣P0(2)20F0(1)=G*C*D*D/BL/BLF0(2)=-G*D*C2 goto200C 跨中力偶作用下的等效荷載列陣Q0(2)30 F0(1)=+(6.0*G*c*c/bl/bl-6.0*g*c/bl+(2.0*g*c*bl-3.0*g*c*c)/bl/bl+g) F0(2)=-(2.0*g*c*bl-3.0*g*c*c)/bl/bl200RETURNENDC******************************************c計(jì)算第NE0個(gè)單元的單元?jiǎng)偠染仃嘍K(2,2)C******************************************SUBROUTINEDG(NE0,NE,E0,GC,GX,DK)DIMENSIONGC(NE),GX(NE),DK(2,2)DO15I=1,2DO10J=1,2DK(I,J)=0.010CONTINUE15CONTINUECDL—桿長,DI-慣性矩,S--線剛度DL=GC(NE0)DI=GX(NE0)S=E0*DI/DLDK(1,1)=4.0*SDK(1,2)=2.0*SDK(2,1)=2.0*SDK(2,2)=4.0*SENDc*********************************************************C集成總體剛度矩陣,ZKl,ZK2分別存主對角元素和付對角元素c*********************************************************SUBROUTINEJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)DIMENSIONZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE)DO10I=1,NJZK1(I)=0.0ZK2(I)=0.010CONTINUEDO100I=1,NECALLDG(I,NE,E0,GC,GX,DK)ZK1(I)=ZK1(I)+DK(1,1)ZK2(I)=ZK2(I)+DK(1,2)ZK1(I+1)=ZK1(I+1)+DK(2,2)100CONTINUEENDC*****************************C形成總荷裁矩陣C*****************************SUBROUTINEXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)DIMENSIONP(NJ),PJ(NP,2),PF(NF,4),F0(2),GC(NE)DO10I=1,NJP(I)=0.010CONTINUEIF(NP.GT.0)THENDO20I=1,NPJ=INT(PJ(I,2)+0.1)P(J)=PJ(I,1)20CONTINUEENDIFIF(NF.GT.0)THENDO100NF0=1,NFCALLDJH(NF0,NE,NF,PF,F0,GC)ND=INT(PF(NF0,3)+0.1)P(ND)=P(ND)+F0(1)P(ND+1)=P(ND+1)+F0(2)100CONTINUEENDIFENDc**********************************c進(jìn)行支承條件處理c**********************************SUBROUTINEZCCL(NZ,NJ,ZK1,ZK2,P)DIMENSIONZK1(NJ),ZK2(NJ),P(NJ)GOTO(10,20,30,20),NZ10GOTO10020ZK1(1)=1.0P(1)=0.0ZK2(1)=0.0IF(NZ.EQ.4)GOTO30GOTO10030ZK1(NJ)=1.0ZK2(NJ-1)=0.0P(NJ)=0.0100CONTINUEENDc******************************c解方程求節(jié)點(diǎn)位移Pc******************************SUBROUTINEFCQJ(NJ,ZK1,ZK2,P)DIMENSIONZK1(NJ),ZK2(NJ),P(NJ)DO10I=1,NJ-1ZK1(I+1)=ZK1(I+1)-ZK2(I)*ZK2(I)/ZK1(I)P(I+1)=P(I+1)-ZK2(I)*P(I)/ZK1(I)10CONTINUEP(NJ)=P(NJ)/ZK1(NJ)DO20I=1,NJ-1P(NJ-I)=(P(NJ-I)-ZK2(NJ-I)*P(NJ-I+1))/ZK1(NJ-I)20CONTINUEENDc*******************************c輸出位移c*******************************SUBROUTINESCWY(NJ,P)DIMENSIONP(NJ)WRITE(9,10)10FORMAT(1X,'：：：：：：：：：位移：；：：：：：：：：')DO100I=1,NJWRITE(9,20)I,P(I)20FORMAT(1X,'結(jié)點(diǎn)號=',I2,5X,F10.4)100CONTINUEENDc******************************c計(jì)算單元桿端內(nèi)力c*****************SUBROUTINEGDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,gd)DIMENSIONP(NJ),WY(2),F(2),F0(2),DK(2,2),PF(NF,4)DIMENSIONGC(NE),GX(NE),gd(ne,2)WRITE(9,2)2FORMAT(1X,'.................各單元桿端內(nèi)力....................')DO200NE0=1,NECALLDG(NE0,NE,E0,GC,GX,DK)WY(1)=P(NE0)WY(2)=P(NE0+1)DO35I=1,2F(I)=0.0DO30J=1,2F(I)=F(I)+DK(I,J)*WY(J)30CONTINUE35CONTINUEIF(NF.GT.0)THENDO50I=1,NFIF(INT(PF(I,3)+0.1).EQ.NE0)THENCALLDJH(I,NE,NF,PF,F0,GC)DO40J=1,2F(J)=F(J)-F0(J)40CONTINUEENDIF50CONTINUEENDIFWRITE(9,150)NE0,F(1),F(2) gd(ne0,1)=f(1) !將桿端力數(shù)據(jù)存入內(nèi)存以備調(diào)用 gd(ne0,2)=f(2)150FORMAT(1X,'單元號=',I2,5X,'左端彎矩=',F9.3,2X,'右端彎矩=',F9.3)200CONTINUE endc*********************************************c 計(jì)算任意截面的彎矩和剪力 （利用疊加法）c************************ subroutinejmnl(gd,sqjm,njm,qdnl,gc,ne,nf,pf) dimensiongd(ne,2),sqjm(njm,2),gc(ne),qdnl1(nf,njm,2),pf(nf,4),qdn， 1l(njm,2)cqdnl(njm,1)=彎矩qdnl(njm,2)=剪力c取左側(cè)梁為隔離體，彎矩逆時(shí)針為正，剪力向上為正， write(9,20)20 format(1x,';;;;;;;;;;;;;;;;;;;;;所求截面內(nèi)力;;;;;;;;;;;;;;;;;;;;') donjm0=1,njm !按要求的截面順序循環(huán) dyh=sqjm(njm0,1) qdnl(njm0,2)=-(gd(dyh,2)+gd(dyh,1))/gc(dyh) !單元上無荷載情況 qdnl(njm0,1)=-(gd(dyh,1)+qdnl(njm0,2)*sqjm(njm0,2)) donf0=1,nf c=pf(nf0,2) d=gc(dyh)-c if((abs(pf(nf0,3)-dyh)<1.0e-6).and.(pf(nf0,2)>=sqjm(njm0,2)))then if(pf(nf0,4)==1)then!單元上有均布荷載且荷載結(jié)束位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=pf(nf0,1)*c*(2.0*gc(dyh)-c)/2.0/gc(dyh)-pf(nf0,1, 1)*sqjm(njm0,2) qdnl1(nf0,njm0,1)=-(pf(nf0,1)*c*(2.0*gc(dyh)-c)/2.0/gc(dyh)*sqjm(n, 1jm0,2)-pf(nf0,1)*sqjm(njm0,2)*sqjm(njm0,2)/2.0) elseif(pf(nf0,4)==2)then!單元上有集中力荷載且位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=pf(nf0,1)*d/gc(dyh) qdnl1(nf0,njm0,1)=-qdnl1(nf0,njm0,2)*sqjm(njm0,2) elseif(pf(nf0,4)==3)then!單元上有集中力偶且位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=-pf(nf0,1)/gc(dyh) qdnl1(nf0,njm0,1)=-qdnl1(nf0,njm0,2)*sqjm(njm0,2) endifelseif((abs(pf(nf0,3)-dyh)<1.0e-6).and.(pf(nf0,2)<sqjm(njm0,2)))& &then if(pf(nf0,4)==1)then!單元上有均布荷載且荷載結(jié)束位置在所求截面右側(cè) qdnl1(nf0,njm0,2)=-pf(nf0,1)*c*c/2.0/gc(dyh) qdnl1(nf0,njm0,1)=-pf(nf0,1)*c*c/2.0/gc(dyh)*(gc(dyh)-sqjm(njm0,2), 1) elseif(pf(nf0,4)==2)then!單元上有集中力荷載且位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=-pf(nf0,1)*c/gc(dyh) qdnl1(nf0,njm0,1)=-pf(nf0,1)*d*c/gc(dyh)-qdnl1(nf0,njm0,2)*(sqjm(n, 1jm0,2)-c) elseif(pf(nf0,4)==3)then!單元上有集中力偶且位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=-pf(nf0,1)/gc(dyh) qdnl1(nf0,njm0,1)=pf(nf0,1)/gc(dyh)*(sqjm(njm0,2)-gc(dyh)) endif else qdnl1(nf0,njm0,2)=0 !其他情況下需疊加的內(nèi)力為0 qdnl1(nf0,njm0,1)=0 endif qdnl(njm0,1)=qdnl(njm0,1)+qdnl1(nf0,njm0,1) !疊加求內(nèi)力 qdnl(njm0,2)=qdnl(njm0,2)+qdnl1(nf0,njm0,2) enddowrite(9,30)njm0,qdnl(njm0,1),qdnl(njm0,2)30 format(1x,'截面',i3,2x,'彎矩=',f9.3,2x,'剪力=',f9.3) enddoWRITE(9,11)11FORMAT('======================計(jì)算結(jié)束======================') end（2）影響線程序PROGRAMLXLyxxDIMENSIONGC(20),GX(20),PF(1,2),sqjm(20,2),DK(2,2),P(45),F0(2),WY, 1(2),F(2),ZK1(45),ZK2(45),gd(20,2),qdnl(20,2),fs(20),fs0(20) character*12dat0,dat1 write(*,*)'輸入:初始數(shù)據(jù)文件名' read(*,10)dat0 write(*,*)'輸入:計(jì)算結(jié)果數(shù)據(jù)文件名' read(*,10)dat110format(A12)open(8,file=dat0,status='old') open(9,file=dat1,status='new') write(9,11)11format(8x,'***********連續(xù)梁彎矩剪力影響線計(jì)算*****************') READ(8,*)NE,NZ,njm,E0,hz WRITE(9,12)NE,NZ,njm,E0,hz 12FORMAT(1X,'單元數(shù)=',6x,I3,5x,'支承類型=',I3,5x,'要求影響線數(shù)量=',1i3,5x,'彈性模量=',F12.4,6x,'單位荷載每次移動(dòng)距離',f12.4,6x)NJ=NE+1CALLSRSJ(NE,GC,GX,sqjm,njm)done0=1,ne !讓單位載荷按指定距離移動(dòng)的循環(huán) fs(ne0)=gc(ne0)/hz fs0(ne0)=int(fs(ne0)+0.1) enddo ne0=0 done0=1,ne pf(1,2)=ne0 dofs00=1,fs0(ne0) pf(1,1)=hz*fs00 write(9,14)ne0,pf(1,1)14 format(/8x,'單位力所在單元號',i3,12x,'距離',f8.4)CALLXCP(NJ,NE,P,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLFCQJ(NJ,ZK1,ZK2,P)CALLGDL(NE,NJ,E0,P,WY,F,F0,DK,PF,GC,GX,gd)calljmnl(gd,sqjm,njm,qdnl,gc,ne,pf) enddo enddoWRITE(9,13)13FORMAT('======================計(jì)算結(jié)束======================') ENDC************************************c輸入初始數(shù)據(jù)SRSJ子程序C************************************SUBROUTINESRSJ(NE,GC,GX,sqjm,njm)DIMENSIONGC(NE),GX(NE),sqjm(njm,2)c sqjm:所求截面njm:截面數(shù)目sqjm(njm,1)=所屬單元sqjm(njm,2)=距左端距離CGC(NE)：桿長GX(NE)：慣性矩C非結(jié)點(diǎn)荷載單位1:PF(1,1)=距離PF(1,2)=單元號WRITE(9,11)11FORMAT(1X,'桿長,慣性矩GC(NE),GX(NE)')READ(8,*)(GC(I),GX(I),I=1,NE) write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,2F8.3) if(njm.gt.0)then write(9,35)35format(1x,'截面所在單元號，距離')read(8,*)((sqjm(i,j),j=1,2),i=1,njm) write(9,36)((sqjm(i,j),j=1,2),i=1,njm)36format(1x,2f8.3) endifENDC*******************************#***************c計(jì)算單位1力引起的等效結(jié)點(diǎn)荷截列陣F0C***********************************************SUBROUTINEDJH(NE,PF,F0,GC)DIMENSIONPF(1,2),F0(2),GC(NE)CC—作用距離,NT--單元號C=PF(1,1)NT=INT(PF(1,2)+0.1)BL=GC(NT)D=BL-CC1=C/BLC2=C1*C1C集中力作用下的等效荷截列陣P0(2)20F0(1)=-1.0*C*D*D/BL/BLF0(2)=1.0*D*C2ENDC******************************************c計(jì)算第NE0個(gè)單元的單元?jiǎng)偠染仃嘍K(2,2)C******************************************SUBROUTINEDG(NE0,NE,E0,GC,GX,DK)DIMENSIONGC(NE),GX(NE),DK(2,2)DO15I=1,2DO10J=1,2DK(I,J)=0.010CONTINUE15CONTINUECDL—桿長,DI-慣性矩,S--線剛度DL=GC(NE0)DI=GX(NE0)S=E0*DI/DLDK(1,1)=4.0*SDK(1,2)=2.0*SDK(2,1)=2.0*SDK(2,2)=4.0*SENDc*********************************************************C集成總體剛度矩陣,ZKl,ZK2分別存主對角元素和付對角元素c*********************************************************SUBROUTINEJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)DIMENSIONZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE)DO10I=1,NJZK1(I)=0.0ZK2(I)=0.010CONTINUEDO100I=1,NECALLDG(I,NE,E0,GC,GX,DK)ZK1(I)=ZK1(I)+DK(1,1)ZK2(I)=ZK2(I)+DK(1,2)ZK1(I+1)=ZK1(I+1)+DK(2,2)100CONTINUEENDC*****************************C形成總荷裁矩陣C*****************************SUBROUTINEXCP(NJ,NE,P,PF,F0,GC)DIMENSIONP(NJ),PF(1,2),F0(2),GC(NE)DO10I=1,NJP(I)=0.010CONTINUECALLDJH(NE,PF,F0,GC)ND=INT(PF(1,2)+0.1)P(ND)=P(ND)+F0(1)P(ND+1)=P(ND+1)+F0(2)ENDc**********************************c進(jìn)行支承條件處理c**********************************SUBROUTINEZCCL(NZ,NJ,ZK1,ZK2,P)DIMENSIONZK1(NJ),ZK2(NJ),P(NJ)GOTO(10,20,30,20),NZ10GOTO10020ZK1(1)=1.0P(1)=0.0ZK2(1)=0.0IF(NZ.EQ.4)GOTO30GOTO10030ZK1(NJ)=1.0ZK2(NJ-1)=0.0P(NJ)=0.0100CONTINUEENDc******************************c解方程求節(jié)點(diǎn)位移Pc******************************SUBROUTINEFCQJ(NJ,ZK1,ZK2,P)DIMENSIONZK1(NJ),ZK2(NJ),P(NJ)DO10I=1,NJ-1ZK1(I+1)=ZK1(I+1)-ZK2(I)*ZK2(I)/ZK1(I)P(I+1)=P(I+1)-ZK2(I)*P(I)/ZK1(I)10CONTINUEP(NJ)=P(NJ)/ZK1(NJ)DO20I=1,NJ-1P(NJ-I)=(P(NJ-I)-ZK2(NJ-I)*P(NJ-I+1))/ZK1(NJ-I)20CONTINUEENDc*******************************c輸出位移c******************************c計(jì)算單元桿端內(nèi)力c*****************SUBROUTINEGDL(NE,NJ,E0,P,WY,F,F0,DK,PF,GC,GX,gd)DIMENSIONP(NJ),WY(2),F(2),F0(2),DK(2,2),PF(1,2)DIMENSIONGC(NE),GX(NE),gd(ne,2)DO200NE0=1,NECALLDG(NE0,NE,E0,GC,GX,DK)WY(1)=P(NE0)WY(2)=P(NE0+1)DO35I=1,2F(I)=0.0DO30J=1,2F(I)=F(I)+DK(I,J)*WY(J)30CONTINUE35CONTINUEIF(INT(PF(1,2)+0.1).EQ.NE0)THENCALLDJH(NE,PF,F0,GC)DO40J=1,2F(J)=F(J)-F0(J)40CONTINUEENDIF50CONTINUE gd(ne0,1)=f(1) gd(ne0,2)=f(2)200CONTINUE endc*********************************************c 計(jì)算任意截面的彎矩和剪力c************************ subroutinejmnl(gd,sqjm,njm,qdnl,gc,ne,pf) dimensiongd(ne,2),sqjm(njm,2),gc(ne),qdnl1(1,njm,2),pf(1,2),qdnl(, 1njm,2)cqdnl(njm,1)=彎矩qdnl(njm,2)=剪力c取左側(cè)梁為隔離體，彎矩逆時(shí)針為正，剪力向上為正， donjm0=1,njm dyh=sqjm(njm0,1) qdnl(njm0,2)=-(gd(dyh,2)+gd(dyh,1))/gc(dyh) qdnl(njm0,1)=-(gd(dyh,1)+qdnl(njm0,2)*sqjm(njm0,2)) c=pf(1,1) d=gc(dyh)-c if((abs(pf(1,2)-dyh)<1.0e-6).and.(pf(1,1)>=sqjm(njm0,2)))then qdnl1(1,njm0,2)=-1.0*d/gc(dyh) qdnl1(1,njm0,1)=-qdnl1(1,njm0,2)*sqjm(njm0,2)elseif((abs(pf(1,2)-dyh)<1.0e-6).and.(pf(1,1)<sqjm(njm0,2)))then if(d<1.0e-6)thenqdnl1(1,njm0,2)=0 qdnl1(1,njm0,1)=1.0*d*c/gc(dyh)-qdnl1(1,njm0,2)*(sqjm(njm0,2)-c) elseqdnl1(1,njm0,2)=1.0*c/gc(dyh) qdnl1(1,njm0,1)=1.0*d*c/gc(dyh)-qdnl1(1,njm0,2)*(sqjm(njm0,2)-c) endifelse qdnl1(1,njm0,2)=0 qdnl1(1,njm0,1)=0 endif qdnl(njm0,1)=qdnl(njm0,1)+qdnl1(1,njm0,1) qdnl(njm0,2)=qdnl(njm0,2)+qdnl1(1,njm0,2) write(9,30)njm0,qdnl(njm0,1),qdnl(njm0,2)30 format(1x,'截面',i3,2x,'彎矩=',f9.3,10x,'剪力=',f9.3) enddo end

論大學(xué)生寫作能力寫作能力是對自己所積累的信息進(jìn)行選擇、提取、加工、改造并將之形成為書面文字的能力。積累是寫作的基礎(chǔ)，積累越厚實(shí)，寫作就越有基礎(chǔ)，文章就能根深葉茂開奇葩。沒有積累，胸?zé)o點(diǎn)墨，怎么也不會(huì)寫出作文來的。寫作能力是每個(gè)大學(xué)生必須具備的能力。從目前高校整體情況上看，大學(xué)生的寫作能力較為欠缺。一、大學(xué)生應(yīng)用文寫作能力的定義那么，大學(xué)生的寫作能力究竟是指什么呢？葉圣陶先生曾經(jīng)說過，“大學(xué)畢業(yè)生不一定能寫小說詩歌，但是一定要寫工作和生活中實(shí)用的文章，而且非寫得既通順又扎實(shí)不可?！睂τ诖髮W(xué)生的寫作能力應(yīng)包含什么，可能有多種理解，但從