有关拱桥拱轴系数m的确定

拱轴系数m：拱顶处横载强度gd与拱脚处横载强度gj的比值。即m（拱轴系数）=gd/gj。

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/2

10 AK1=LOG(AM1+SQRT(AM1**2-1))

QJ=0

20 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 IF

30 CONTINUE

SJ3=0

DO 50 J=1,NN,1

PP=P(J)+Y1(J)*A1(J)*R2

SJ3=SJ3+PP*(AL-S(J))

50 CONTINUE

AM4=X1+X2+X3

AMJ=SJ1+SJ2+SJ3

YF=AM4/AMJ

AM2=(1/YF-2)**2/2-1

FF=ABS(AM2-AM1)

IF(FF.GT.0.005)THEN

AM1=AM2

GOTO 10

ELSE

M=AM2

WRITE(*,*)’ GONG ZHOU XI SHU:M,KUA JING:L,SHI GAO:F’

WRITE(*,*)’ M=’,M

WRITE(*,*)’ L=’,L

WRITE(*,*)’ F=’,F

END IF

ELSE

QJ=QI

GOTO 20

END IF

CLOSE(1)

END

C

SUBROUTINE JIFEN1(B1,HD,D,TG,AK1,AL,F,AM1,R1,A2,B,A,X1) REAL K1,I0,I1

FNA(X)=B1*(HD+D/2*(1-SQRT(1+(TG*(EXP(AK1*X/AL)-

+ EXP(-AK1*X/AL))*0.5)**2))+F*(EXP(AK1*X/AL)+

+ EXP(-AK1*X/AL))*0.5/AL*(AM1-1))*R1*(A2-X)

MM=6000

H=(B-A)/2

T=0

N11=1

E1=FNA(A)+FNA(B)

EE=FNA(A+H)

I1=H*(E1+4*EE)/3

10 H=H/2

N11=2*N11

T=T+EE

EE=0

K1=A+H

DO 20 J9=1,N11,1

EE=EE+FNA(K1)

20 K1=K1+H+H

I0=H*(E1+4*EE+2*T)/3

IF(N11.GT.MM) THEN

WRITE(*,*)’ ***********************’

ELSE

IF(ABS(I1-I0).LT.0.005)THEN

C WRITE(*,*)’I01=’,I0

X1=I0

ELSE

I1=I0

GOTO 10

END IF

END IF

END

C

SUBROUTINE JIFEN2(A1,TG,AK1,AL,R,A2,B,A,X2) REAL K1,I0,I1

FNA(X)=A1*SQRT(1+(TG*(EXP(AK1*X/AL)-

+ EXP(-AK1*X/AL))*0.5)**2)*R*(A2-X) MM=6000

H=(B-A)/2

T=0

N11=1

E1=FNA(A)+FNA(B)

EE=FNA(A+H)

I1=H*(E1+4*EE)/3