拱轴系数程序确定.docVIP

下载本文档

45
0
约6.65千字
约 7页
2019-10-18 发布于江西
举报
版权申诉

拱轴系数程序确定.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

在拱桥设计中，拱桥的拱轴系数m必须由试算确定，试算的次数及计算量相当繁重，只能由电算确定：现将确定空腹式拱桥的拱轴系数m，并计算其跨径l,矢高f的ORTRAN程序清单写于下面（仅供参考）：L0，F0，S1——拱的跨径，净矢高，实腹段长一半；A1,B1——截面面积，拱圈宽度；R，R1，R2——主拱圈材料，拱顶填料，腹孔墩材料的单位体积重；D，HD，Y1——主拱圈的高度，拱顶填料厚度，腹孔墩墩顶至主拱圈顶的竖向距离；NN——半拱内腹空墩个数；S（J），P（J），AI（J）——腹孔墩J离拱顶的水平距离、腹孔墩上腹孔自重、腹孔墩截面面积；? ?? ?? ? ? ? PROGRAM GQSL? ? ? ? REAL M0,L,L0,M? ? ? ? CHARACTER*12 INDAT? ? ? ? DIMENSION Y(150),Y1(150),S(150),A1(150),P(150)? ? ? ? WRITE(*,*) Please input primary data file name!? ? ? ? READ(*,(A12)) INDAT? ? ? ? OPEN(1,FILE=INDAT,STATUS=OLD)? ? ? ? READ(1,*)M0? ? ? ? READ(1,*)F0,L0,S1? ? ? ? READ(1,*)D,HD,YJ? ? ? ? READ(1,*)R,R1,R2? ? ? ? READ(1,*)NN? ? ? ? DO 5 J=1,NN,1? ? ? ? READ(1,*)S(J),P(J),A1(J)5? ? ? ? CONTINUE? ? ? ? A1=0.9? ? ? ? B1=1? ? ? ? AM1=M0? ? ? ? BL=L0/210? ? ? ? AK1=LOG(AM1+SQRT(AM1**2-1))? ? ? ? QJ=020? ?? ?D1=(F0+D*(1-COS(QJ))/2)/(BL+D*SIN(QJ)/2)? ? ? ? D1=D1*AK1*(EXP(AK1)-EXP(-AK1))*0.5/(AM1-1)? ? ? ? QI=ATAN(D1)? ? ? ? ? ? ? ? QQ=ABS(QI-QJ)? ? ? ? IF(QQ.LT.0.0001)THEN? ? ? ?? ?L=L0+D*SIN(QJ)? ? ? ?? ?F=F0+D*(1-COS(QJ))/2? ? ? ?? ?A2=L/4? ? ? ?? ?AL=L/2? ? ? ?? ?TG=AK1*F/(AL*(AM1-1))? ? ? ?? ?A=0? ? ? ?? ?B=S1? ? ? ?? ?CALL JIFEN1(B1,HD,D,TG,AK1,AL,F,AM1,R1,A2,B,A,X1)? ? ? ?? ?B=A2? ? ? ?? ?A=0? ? ? ?? ?CALL JIFEN2(A1,TG,AK1,AL,R,A2,B,A,X2)? ? ? ?? ?A=0? ? ? ?? ?B=S1? ? ? ?? ?CALL JIFEN3(B1,HD,D,TG,AK1,AL,F,AM1,R1,B,A,SJ1)? ? ? ?? ?B=AL? ? ? ?? ?A=0? ? ? ?? ?CALL JIFEN4(A1,TG,AK1,AL,R,B,A,SJ2)? ? ? ?? ?X3=0? ? ? ?? ?DO 30 J=1,NN,1? ? ? ?? ???Y(J)=F/(AM1-1)*((EXP(AK1*S(J)/AL)+EXP(-AK1*S(J)/AL))/2-1)? ? ? ?? ???TG=AK1*F/(AL*(AM1-1))? ? ? ?? ???TG1=TG*(EXP(AK1*S(J)/AL)+EXP(-AK1*S(J)/AL))/2? ? ? ?? ???SC=1/SQRT(1+TG1**2)? ? ? ?? ???Y1(J)=Y(J)+D/2+YJ-D/2/SC? ? ? ?? ???IF(AL/2-S(J).GE.0)THEN? ? ? ?? ?? ? X3=X3+(P(J)+Y1(J)*A1(J)*R2)*(AL/2-S(J))? ? ? ?? ???ELSE ? ?? ?? ?? ???X3=X3+0? ? ? ?? ???END IF30? ?? ? ? ?? ?CONTINUE? ? ? ?? ?SJ3=0? ? ? ?? ?DO 50 J=1,NN,1? ? ? ?? ???PP=P(J)+Y1(J)*A1(J)*R2? ? ? ?? ???SJ3=SJ3+PP*