既有曲線整正

1、表表6-5 6-5 漸伸線長度的實用計算公式漸伸線長度的實用計算公式3.3.撥距的大小和方向撥距的大小和方向b1o測量始點(diǎn)b1b2b2既有曲線設(shè)計曲線a1a2始切線圖圖6-17 6-17 既有曲線整正既有曲線整正 jsEE 為正：曲線內(nèi)壓，向圓心方向撥動為正：曲線內(nèi)壓，向圓心方向撥動 為負(fù)：曲線外挑，向切線方向撥動為負(fù)：曲線外挑，向切線方向撥動 嚴(yán)格來說，撥動方向是沿著漸伸線方向撥動的。但是因為撥嚴(yán)格來說，撥動方向是沿著漸伸線方向撥動的。但是因為撥動量很小，在既有曲線附近很小的范圍內(nèi)，完全可以認(rèn)為漸伸線動量很小，在既有曲線附近很小的范圍內(nèi)，完全可以認(rèn)為漸伸線方向與既有曲線的法線方向重合。所以，

2、實際工作中，既有曲線方向與既有曲線的法線方向重合。所以，實際工作中，既有曲線撥正系沿法線方向撥動。撥正系沿法線方向撥動。 4.4.撥距計算的步驟撥距計算的步驟1 1）計算既有曲線各測點(diǎn)的漸伸線長度）計算既有曲線各測點(diǎn)的漸伸線長度 2 2）選定設(shè)計曲線半徑）選定設(shè)計曲線半徑 3) 3) 計算設(shè)計曲線各測點(diǎn)的漸伸線長度計算設(shè)計曲線各測點(diǎn)的漸伸線長度4) 4) 求出撥動距離求出撥動距離 jESRsEjsEE ExcelExcel撥距計算撥距計算一、既有曲線資料輸入一、既有曲線資料輸入（一）在（一）在A6A6、A23A23分別輸入分別輸入“不不0 0、不、不1313、不、不21”21”；在；在B6-B

3、27B6-B27單元格中，分別輸入單元格中，分別輸入“K169+980K169+980、K170+000320”K170+000320”；在；在A28A28單元格中，輸入單元格中，輸入“前視點(diǎn)前視點(diǎn)”；在在C7-C28C7-C28單元格中，輸入單元格中，輸入“180180、1803”1803”；在；在D7-D28D7-D28單元格中，輸入單元格中，輸入“0 0、 5716”5716”；在；在E7-E28E7-E28單元格中，單元格中，輸入輸入“3737、5649”5649”；在；在L7-L19L7-L19單元格中，分別輸入單元格中，分別輸入“2020、40260”40260”；在；在L7-L1

4、9L7-L19單元格中，分別輸入單元格中，分別輸入“2020、40160”40160”。如圖。如圖2-12-1所示：所示：（二）在第（二）在第1111列的列的K7K7單元格中，輸入單元格中，輸入“=I7=I7* *J7”J7”，回車，然后復(fù)制此公式到，回車，然后復(fù)制此公式到K28K28單元格，如圖所示：單元格，如圖所示：（三）既有曲線漸伸線長度（三）既有曲線漸伸線長度 計算計算jE在第在第1212列列L7L7單元格中，輸入單元格中，輸入“=K7”=K7”，回車，然后復(fù)制，回車，然后復(fù)制此公式到此公式到K19K19單元格；在單元格；在L20L20單元格中，輸入單元格中，輸入“=K20+57.26

5、4”=K20+57.264”，回車，然后復(fù)制此公式到，回車，然后復(fù)制此公式到K27K27單元格；單元格；如圖所示：如圖所示：Fortran程序計算撥距程序計算撥距一、已知數(shù)據(jù)一、已知數(shù)據(jù) 已知數(shù)據(jù)有、曲線偏角資料已知數(shù)據(jù)有、曲線偏角資料、置鏡點(diǎn)里程、置鏡點(diǎn)里程zjlczjlc、里程、里程lclc、前、前視角度視角度qsjqsj。 讀入已知數(shù)據(jù)，利用讀入已知數(shù)據(jù)，利用open語句來實現(xiàn)數(shù)據(jù)的讀入。語句來實現(xiàn)數(shù)據(jù)的讀入。 open(1,file=shuru.dat,status=old)二、曲線偏角資料二、曲線偏角資料度轉(zhuǎn)換成弧度度轉(zhuǎn)換成弧度 a=bd(i)b=(bd(i)-a*1.0)*100c

6、=int(bd(i)-a)*100-b)*100)bhd(i)=(a+b/60.0+c/3600.0)*3.1415926/180 因為我們輸入的角度數(shù)據(jù)是度分秒的形式，所以必因為我們輸入的角度數(shù)據(jù)是度分秒的形式，所以必須把它轉(zhuǎn)換為弧度，在程序中默認(rèn)的角度是弧度。須把它轉(zhuǎn)換為弧度，在程序中默認(rèn)的角度是弧度。A A為度，為度，b b為分，為分，c c為秒為秒三、里程數(shù)組的計算三、里程數(shù)組的計算 這個和第一個置鏡點(diǎn)里程有關(guān)系，利用這個關(guān)系，用一個循環(huán)這個和第一個置鏡點(diǎn)里程有關(guān)系，利用這個關(guān)系，用一個循環(huán)語句來實現(xiàn)。語句來實現(xiàn)。 lc(1)=zjdlc(1) do 15 i=2,23 lc(i)=l

7、c(i-1)+2015 continue這里這里23是個跟測是個跟測點(diǎn)偏角個數(shù)有關(guān)點(diǎn)偏角個數(shù)有關(guān)的變量，測得的的變量，測得的偏角越多，這個偏角越多，這個值越大！值越大！四、四、l長度計算（長度計算（20、40。）。） n=1i=1m=(zjdlc(i+1)-zjdlc(i)/20 do 30 j=n,ml(j)=20*j30continue 35i=i+1if(i.le.3) then n=n+mm=(zjdlc(i+1)-zjdlc(i)/20 do 40 j=n,m+n-1k=j-n+1l(j)=20*k40continueend ifif(i.le.2) goto 35do 50 i=1

8、,22ej1(i)=bhd(i)*l(i)50 continue這塊，我們需要搞清這塊，我們需要搞清楚，置鏡點(diǎn)里程和測楚，置鏡點(diǎn)里程和測點(diǎn)里程的關(guān)系，由測點(diǎn)里程的關(guān)系，由測點(diǎn)里程減去置鏡點(diǎn)里點(diǎn)里程減去置鏡點(diǎn)里程來獲得。程來獲得。3和置鏡點(diǎn)的個數(shù)和置鏡點(diǎn)的個數(shù)有關(guān)系有關(guān)系 五、計算既有漸伸線長度五、計算既有漸伸線長度 i=1n=1k=0m=(zjdlc(i+1)-zjdlc(i)/20 do 60 j=n,mej(j)=ej1(j)60continue65 k=m+n-1p=ej(k)n=n+m i=i+1m=(zjdlc(i+1)-zjdlc(i)/20 do 70 j=n,n+m-1ej(j

9、)=ej1(j)+p70continue if(i.le.2) goto 65這個和上一個有相這個和上一個有相似之處。需要在置似之處。需要在置鏡點(diǎn)以后的數(shù)據(jù)都鏡點(diǎn)以后的數(shù)據(jù)都要進(jìn)行累加處理。要進(jìn)行累加處理。 六、計算曲線轉(zhuǎn)角六、計算曲線轉(zhuǎn)角a 在公式中，我們知道曲線轉(zhuǎn)角在公式中，我們知道曲線轉(zhuǎn)角a a等于最后一個置鏡點(diǎn)處的累計偏角加上等于最后一個置鏡點(diǎn)處的累計偏角加上前視角度，這個計算可以通過賦值語句來完成。前視角度，這個計算可以通過賦值語句來完成。alfa=qsjd+bd(nd)a=alfab=(alfa-a*1.0)*100c=int(alfa-a)*100-b)*100)if(c.gt.

10、60) thenb=b+1c=0endifalfa1=(a+b/60.0+c/3600.0)*3.1415926/180 計算測量終點(diǎn)到計算測量終點(diǎn)到QZ點(diǎn)的距離點(diǎn)的距離X X=EJ(nd)/ALFA1C 計算計算QZ里程里程 QZ=ZJDLC(nz)-XC 計算計算QZ對應(yīng)的點(diǎn)號對應(yīng)的點(diǎn)號 七、計算七、計算qzqz里程以及對應(yīng)的樁號里程以及對應(yīng)的樁號QZ=QZ=終點(diǎn)里程終點(diǎn)里程- -測量終點(diǎn)的漸伸線長度測量終點(diǎn)的漸伸線長度/ /曲線轉(zhuǎn)角曲線轉(zhuǎn)角設(shè)設(shè)X=X=各測點(diǎn)里程各測點(diǎn)里程- -曲中里程曲中里程, ,取取X X最小時的測點(diǎn)號為最小時的測點(diǎn)號為曲中的點(diǎn)號。曲中的點(diǎn)號。Nd為偏角總個數(shù)，為偏角

11、總個數(shù)，nz為置鏡點(diǎn)的總個數(shù)為置鏡點(diǎn)的總個數(shù) 八、計算曲線半徑八、計算曲線半徑PM 1M2M1PM 2Pi1Pi2PinMi1Mi2MinA1A2A3AnAn -1計算三級差計算三級差DO 75 J=1,nd-1DO 75 J=1,nd-1 YJC(J)=EJ(J+1)-EJ(J) YJC(J)=EJ(J+1)-EJ(J)75 CONTINUE 75 CONTINUE DO 80 J=2,nd-1DO 80 J=2,nd-1 RJC(J)=YJC(J+1)-YJC(J) RJC(J)=YJC(J+1)-YJC(J)80 CONTINUE80 CONTINUEDO 85 J=3,nd-1DO 8

12、5 J=3,nd-1 SJC(J)=RJC(J+1)-RJC(J) SJC(J)=RJC(J+1)-RJC(J)85 CONTINUE85 CONTINUE程序語句程序語句計算計算QZQZ對應(yīng)的點(diǎn)號對應(yīng)的點(diǎn)號 CHA(1)=ABS(QZ-LC(1)CHA(1)=ABS(QZ-LC(1)ZXZ=CHA(1) ZXZ=CHA(1) DO 95 J=2,nd+1DO 95 J=2,nd+1CHA(J)=ABS(QZ-LC(J)CHA(J)=ABS(QZ-LC(J)IF (CHA(J).LT.ZXZ) THEN IF (CHA(J).LT.ZXZ) THEN ZXZ=CHA(J)ZXZ=CHA(J)M

13、=JM=JENDIFENDIF95 CONTINUE95 CONTINUE因為我們要從因為我們要從qzqz向兩邊搜索向兩邊搜索m1m1和和m2m2，所以要確定，所以要確定qzqz的位置。的位置。計算計算M1,M2 M1,M2 w=0.05w=0.05DO 98 I=M,3,-1DO 98 I=M,3,-1IF (ABS(SJC(I).LT.w) M1=I IF (ABS(SJC(I).LT.w) M1=I 98 CONTINUE98 CONTINUEDO 105 I=M,nd-2DO 105 I=M,nd-2IF (ABS(SJC(I).LT.w) M2=I IF (ABS(SJC(I).LT

14、.w) M2=I 105 CONTINUE105 CONTINUE計算曲線半徑Rs S=0DO 110 I=M1,M2 S=RJC(I)+S110 CONTINUERS=400/S*(M2-M1+1)RS=INT(RS/10)*10 計算未加設(shè)緩和曲線圓曲線長度計算未加設(shè)緩和曲線圓曲線長度LSLSLS=RSLS=RS* *ALFA1ALFA1ZY=QZ-LS/2.0ZY=QZ-LS/2.0YZ=QZ+LS/2.0 YZ=QZ+LS/2.0 曲線半徑優(yōu)化曲線半徑優(yōu)化* 計算梯形面積TXMJ=(EJ(M2)-(LC(M2+1)-ZY)*2/(2*RS)+ $ (EJ(M1)-(LC(M1+1)-Z

15、Y)*2/(2*RS)*(LC(M2)-LC(M1)/2.0C 計算小梯形面積之和STXMJ=0.0DO 112 I=M1,M2-1STXMJ=STXMJ+(EJ(I)-(LC(I+1)-ZY)*2/2/RS)+ $ (EJ(I+1)-(LC(I+2)-ZY)*2/2/RS)*20/2.0 112 CONTINUEf=txmj-stxmj if (f.gt.0) thencall bjyhd(f,txmj,stxmj,alfa1,rs,m1,m2,ej,lc,zy) elsecall bjyhx(f,txmj,stxmj,alfa1,rs,m1,m2,ej,lc,zy) end if IF (

16、ALFA1.GT.50*3.1415926/180) THENRS=anint(RS*10)/10ELSE IF (ALFA1.GT.30*3.1415926/180) THENRS=anint(RS)ELSE IF (ALFA1.GT.20*3.1415926/180) THENRS=anint(RS/5)*5ELSE IF (ALFA1.GT.10*3.1415926/180) THENRS=anint(RS/10)*10ELSE IF (ALFA1.GT.0) THENRS=anint(RS/25)*25ENDIF 半徑調(diào)整PY ZPZYPM 2PM 1Y ZM2M1ZY 計算未加設(shè)緩和

17、曲線圓曲線長度計算未加設(shè)緩和曲線圓曲線長度LSLS LS=RS LS=RS* *ALFA1ALFA1ZY=QZ-LS/2.0ZY=QZ-LS/2.0YZ=QZ+LS/2.0YZ=QZ+LS/2.0 計算緩和曲線長度計算緩和曲線長度l0l0S1=0S1=0DO 115 I=M1+1,M2-1DO 115 I=M1+1,M2-1S1=EJ(I-1)-(LC(I)-ZY)S1=EJ(I-1)-(LC(I)-ZY)* * *2/(22/(2* *RS)+S1RS)+S1115 CONTINUE115 CONTINUEP=S1/(M2-M1-1) P=S1/(M2-M1-1) L0=SQRT(24L0=SQRT(24* *RSRS* *abs(P)abs(P)L0=INT(L0/10)L0=INT(L0/10)* *10 10 這里只有等長緩和曲線的計算，還需考慮不等長的情況。計算五大樁里程計算五大樁里程zh=zy-l0/2zh=zy-l0/2hy=zy+l0/2hy=zy+l0/2yh=yz-l0/2yh=yz-l0/2hz=yz+l0/2 hz=yz+l0/2 計算設(shè)計漸伸線長度do 120 i=2,nd+1if (lc(i).gt.hz) thenX