资源描述：
《有限单元法及程序设计实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、——连续梁及平面桁架实验报告学院：_____机械学院_______姓名：_______朱琳_________学号：_____201107812______班级：______78班__________实验一：连续梁静力分析一、连续梁结构分析程序总框图开始输入原始数据形成结点载荷列阵形成整体刚度矩阵引入支撑条件解方程并打印位移计算并打印杆端力结束二、连续梁静力分析源程序PROGRAMCBSAPCSTRUCTURALANALYSISPROGRAMFORCONTINUOUSBEAMDIMENSIONGC(20),GX(20),AMF(20,2),AM(20,2),#
2、P(21),AK(21,21),XG(20)CHARACTER*80TITLEOPEN(2,FILE='CBSAP.IN')OPEN(6,FILE='CBSAP.OUT')READ(2,*)TITLEWRITE(6,*)TITLEREAD(2,20)NJ,IL,IR20FORMAT(3I5)WRITE(6,30)NJ,IL,IY30FORMAT(/3X,'NUMEROFNODE='I2,5X,'BIND(L)='I1,2X,'BIND(R)='#I1,3X,'1=FIXED0=SIMPLYSUPPORTED')NNE=NJ-1READ(2,40)(GC(I)
3、,GX(I),(AMF(I,J),J=1,2),I=1,NNE)40FORMAT(4F10.3)WRITE(6,50)(I,GC(I),GX(I),(AMF(I,J),J=1,2),I=1,NNE)50FORMAT(/3X,'NO.E',7X,'LENGHT',6X,'STIFFNESS',6X,'MOMENT(L)',#6X,'MOMENT(R)'/(I5,4F15.4))DOI=1,NNEXG(I)=GX(I)/GC(I)ENDDOP(1)=-AMF(1,1)P(NJ)=-AMF(NNE,2)DOJ=2,NNEP(J)=-AMF(J-1,2)-AMF(J
4、,1)ENDDODOI=1,NJDOJ=1,NJAK(I,J)=0ENDDOENDDOAK(1,1)=4.0*XG(1)AK(NJ,NJ)=4.0*XG(NNE)DOJ=2,NNEAK(J,J)=4.0*(XG(J-1)+XG(J))ENDDODOJ=2,NJAK(J,J-1)=2.0*XG(J-1)AK(J-1,J)=2.0*XG(J-1)ENDDOIF(IL.NE.0)THENAK(1,1)=1.0AK(1,2)=0.0AK(2,1)=0.0P(1)=0.0ELSEENDIFIF(IR.NE.0)THENAK(NJ,NJ)=1.0AK(NNE,NJ)=0
5、.0AK(NJ,NNE)=0.0P(NJ)=0.0ELSEENDIFDOK=1,NNEDOI=K+1,NJC=AK(I,K)/AK(K,K)P(I)=P(I)-C*P(K)DOJ=1,NJAK(I,J)=AK(I,J)-C*AK(K,J)ENDDOENDDOENDDOP(NJ)=P(NJ)/AK(NJ,NJ)DOI=NNE,1,-1DOJ=I+1,NJP(I)=P(I)-AK(I,J)*P(J)ENDDOP(I)=P(I)/AK(I,I)ENDDOWRITE(6,60)(I,P(I),I=1,NJ)60FORMAT(/10X,'ANGLE(RADIAN)OF
6、FOTATION'//(5(I3,F10.4)))DOI=1,NNEAM(I,1)=4.0*XG(I)*P(I)+2.0*XG(I)*P(I+1)+AMF(I,1)AM(I,2)=2.0*XG(I)*P(I)+4.0*XG(I)*P(I+1)+AMF(I,2)ENDDOWRITE(6,70)(I,(AM(I,J),J=1,2),I=1,NNE)70FORMAT(/17X,'MOMENTSOFBEAMENDS'//3X,#'NO.BEAM',6X,'1(LEFT)',5X,'2(RIGHT)'/(I8,2F16.4))ENDPROGRAMCBSAP三、连续梁实
7、验题目试用连续梁静力分析程序CBSAP计算图示连续梁的杆端弯矩并给弯矩图四、连续梁计算输入数据标点EXAMPLE---------(3-1)结点数4左支承0右支承1单元号长度抗弯刚度固端弯矩左端右端①612-6060②1010-250250③88-187.5112.5五、连续梁输入及输出文件六、连续梁弯矩图实验二：平面桁架静力分析一、平面桁架静力分析程序总框图开始(1)标题及数组说明结束NJ是否为零？(8)子程序STIFFN(2)输入并打印原始数据(9)子程序TRANS(3)打算并打印半宽带(10)子程序MULVA(4)组集整体刚度矩阵(11)子程序CALM
8、(5)形成结点载荷列阵(12)子程序CSL(6)解方
