Nghiên cứu ứng dụng công nghệ gps trong trắc địa công trình ở việt nam

-2- Lời cam đoan Tôi xin cam đoan công trình nghiên cứu riêng Các số liệu, kết nêu luận án trung thực cha đợc công bố công trình khác Tác giả luận án -3- Mục lục Trang Trang phụ bìa Lời cam đoan.2 Mục lục Danh mục ký hiệu, chữ viết tắt Danh mục bảng .6 Danh mục hình vẽ, đồ thị .8 Mở đầu ..11 Chơng 1- tỉng quan vỊ øng dơng cđa c«ng nghƯ GPS 1.1 ứng dụng công nghệ GPS lĩnh vực trắc địa giới Việt nam 15 1.2 Các kết đ nghiên cứu ứng dụng GPS trắc địa công trình19 1.3.Định hớng nghiên cứu Luận án25 Chơng - Nghiên cứu khả ứng dụng công nghệ GPS trắc địa công trình Việt nam 2.1.Khả ứng dụng công nghệ GPS giai đoạn quy hoạch - thiết kế công trình 26 2.2.Khả ứng dụng công nghệ GPS giai đoạn thi công xây dựng công trình .32 2.3.Khả ứng dụng công nghệ GPS quan trắc chuyển dịch ngang công trình40 Chơng - Một số biện pháp nâng cao chất lợng lới GPS trắc địa công trình 3.1.Nâng cao chất lợng lới GPS biện pháp tối u hoá thiết kế lới GPS .59 -4- 3.2.Nâng cao chất lợng lới GPS biện pháp đo ngoại nghiệp 82 3.3 Nâng cao chất lợng lới GPS biện pháp xử lý số liệu đo GPS 91 Chơng - Phơng pháp xử lý số liệu lới GPS trắc địa công trình 4.1 Tính chuyển toạ độ điểm đo GPS hệ toạ độ thi công công trình 94 4.2.Tính chuyển độ cao điểm đo GPS hệ độ cao thi công công trình 119 4.3 Thành lập phần mềm thiết kế tính chuyển tọa độ lới GPS trắc địa công trình 129 Kết luận kiến nghị Danh mục công trình tác giả Tài liệu tham khảo Phụ lục Phụ lôc Phô lôc Phô lôc -5- Danh mục ký hiệu chữ viết tắt GPS : Global Positioning system DGPS : Differential GPS GPS RTK : Real Time Kinematic GPS GPS PPK : Post Processing Kinematic GPS DOP : Dilution of Precision PDOP : Position Dilution of Precision GDOP : Geometric Position Dilution of Precision HDOP : Horizontal Position Dilution of Precision VDOP : Vertical Position Dilution of Precision 10 C/A code : Coarse / Acquisition code 11 P- Code : Precision code 12 NCKH : Nghiên cứu khoa học 13.KTXD : Kỹ thuật xây dựng -6- Danh mục bảng biểu B¶ng 2.1 Sai số vị trí điểm lới khống chế đo vẽ đồ địa hình 28 Bảng 2.2 Độ xác cho phép trắc chuyển dịch ngang công trình 41 Bảng 2.3 Chiều dài cạnh lới đo thực nghiƯm…………………… 48 B¶ng 2.4 Thêi gian thu tÝn hiƯu vƯ tinh ngày 17/11/2007.51 Bảng 2.5 Toạ độ độ cao điểm A 52 Bảng 2.6 Kết đo thực nghiệm chuyển dịch ngang công nghệ GPS.52 Bảng 2.7 Kết tính chuyển dịch ngang theo phơng án đo GPS.55 Bảng 2.8 So sánh phơng án quan trắc chuyển dịch ngang công nghệ GPS56 Bảng 2.9 So sánh biến đổi toạ độ độ cao điểm cố định ca đo.57 Bảng 3.1 Kết ớc tính độ xác đánh giá độ xác sau bình sai lới GPS cầu B i Cháy 67 Bảng 3.2 Kết ớc tính độ xác đánh giá độ xác sau bình sai lới GPS khu du lịch sinh thái tỉnh Hoà Bình 68 Bảng 3.3 Kết ớc tính độ xác đánh giá độ xác sau bình sai lới GPS cầu Thanh Trì 69 Bảng 3.4 Giá trị sai số trung phơng cho phép đào thông hầm76 Bảng 3.5 Kết thiết kế lới GPS hầm thuỷ điện Quảng Trị 82 Bảng 3.6 Độ xác đo chiều cao anten máy thu GPS 86 Bảng 4.1 So sánh chiều dài cạnh lới khống chế thi công cầu B i Cháy96 Bảng 4.2 So sánh kết đo cạnh lới GPS khu công nghiệp Yên Phong.98 Bảng 4.3 Toạ độ điểm GPS lới khống chế thi công cầu B i Cháy 115 Bảng 4.4 Toạ ®é c¸c ®iĨm ®o GPS mói 108 ( sau tính chuyển ) 116 Bảng 4.5 So sánh chiều dài cạnh lới GPS theo phơng án đ tính chuyển.116 Bảng 4.6 So sánh kết bình sai lới khống chế thi công cầu B i Cháy 117 -7- Bảng 4.7 So sánh kết đo cạnh lới GPS khu công nghiệp Yên Phong (sau tính chuyển ) 118 Bảng 4.8 Bảng toạ độ độ cao điểm đo GPS cầu Thanh Trì126 Bảng 4.9 Kết tính toán thực nghiệm tính chuyển độ cao lới GPS127 Bảng 4.10 So sánh độ xác tính chuyển độ cao GPS hệ độ cao thi công.128 -8- Danh mục hình vẽ, đồ thị -Hình 2.1.ứng dụng công nghệ GPS PPK để đo vẽ đồ địa hình tỉnh Thái Nguyên 30 Hình 2.2 Sơ đồ lới GPS thi công xây dựng đờng hầm 36 Hình 2.3 Sơ đồ lới khống chế thi công cầu vợt 38 Hình 2.4 Sơ đồ lới GPS thi công xây dựng cầu vợt 39 Hình 2.5 Sử dụng máy thu GPS R3 để lập lới khống chế thi công.39 Hình 2.6 Toàn cảnh đập thuỷ điện Sông Đà - Hoà Bình.42 Hình 2.7 Chuyển dịch ngang qi điểm quan trắc i.43 Hình 2.8 Lới đo chuyển dịch ngang công nghệ GPS 45 Hình 2.9 Mốc quan trắc chuyển dịch ngang đập thủy điện Hoà Bình 46 Hình 2.10 Sơ đồ đo thực nghiệm quan trắc chuyển dịch ngang công nghệ GPS47 Hình 2.11 Cấu tạo mốc định tâm bắt buộc.48 Hình 2.12 Cấu tạo chuyển dịch máy thu GPS49 Hình 2.13 Sơ đồ lới GPS đo thực nghiệm chuyển dịch ngang 50 Hình 2.14 Sơ đồ tính chuyển dịch ngang từ hai điểm cố định máy thu GPS54 Hình 2.15 Sơ đồ tính chuyển dịch ngang từ ba điểm cố định- máy thu GPS 55 Hình 3.1 Sơ đồ lới GPS cầu B i Cháy 67 Hình 3.2 Lới GPS khu du lịch sinh thái tỉnh Hoà Bình68 Hình 3.3 Sơ đồ lới GPS cầu Thanh Trì.69 Hình 3.4 Quy trình thiết kế lới khống chế đợc thành lập công nghệ GPS theo h−íng thiÕt kÕ tèi −u…………………………………… 71 H×nh 3.5 ThiÕt kÕ lới GPS xây dựng cầu vợt 75 -9- Hình 3.6 Sơ đồ hệ thống sở trắc địa thi công xây dựng đờng hầm .78 Hình 3.7 Sơ đồ hầm thuỷ điện Quảng Trị 79 Hình 3.8 Sơ đồ thiết kế lới khống chế thi công hầm thuỷ điện Quảng Trị 80 Hình 3.9 Sơ đồ thiết kế lới GPS (Phơng án ban đầu) 81 Hình 3.10 Quan hệ sai số đo chiều cao Antenna đại lợng đo GPS 86 Hình 3.11 Cửa sổ đặt điều kiện quan sát máy thu GPS 89 Hình 3.12 Cửa số xác định thời gian quan sát ngày 89 Hình 4.1 Lới GPS thi công cầu B i Cháy.96 Hình 4.2 Sơ đồ lới GPS khu công nghiệp Yên Phong Bắc Ninh97 Hình 4.3 Hệ toạ độ vuông góc không gian địa tâm 100 Hình 4.4 Hệ toạ độ trắc địa 101 Hình 4.5 Hệ tọa độ địa diện chân trời 103 Hình 4.6 Quan hệ hệ toạ độ địa diện xích đạo hệ toạ độ địa diện chân trời 104 Hình 4.7 Hệ quy chiếu toàn cầu WGS-84106 Hình 4.8 Chọn mặt chiếu trắc địa công trình 107 Hình 4.9 Nguyên lý tính chuyển độ cao GPS theo mặt Ellipsoid thi công 122 Hình 4.10 Sơ đồ lới GPS cầu Thanh Trì.126 Hình 4.11 Sơ đồ lới GPS cầu B i Cháy 127 Hình 4.12 Giao diện phần mềm GP.AGE129 Hình 4.13 Modul ớc tính độ xác lới GPS 132 Hình 4.14 Giao diện hình ớc tính lới GPS133 Hình 4.15 Kết ớc tính độ xác lới GPS thiết kế 133 Hình 4.16 Giao diện hình tính chuyển toạ độ lới GPS 134 Hình 4.17 Modul tính chuyển toạ độ GPS 134 Hình 4.18 Kết tính chuyển toạ độ GPS 135 - 10 - Hình 4.19 Giao diện hình hiệu chØnh d÷ liƯu thiÕt kÕ l−íi GPS………135 - 11 - Mở đầu Tính cấp thiết đề tài Hiện công nghệ GPS với đặc tính u việt ® ®−ỵc øng dơng rÊt réng r i ë nhiỊu lĩnh vực giới.Trong trắc địa, công nghệ GPS đ đợc ứng dụng để thành lập lới toạ độ liên lục địa, lới toạ độ quốc gia đo chi tiết Công nghệ GPS đ đợc ứng dụng trắc địa công trình ba giai đoạn : quy hoạch thiết kế, thi công xây dựng quan trắc an toàn công trình Việt nam, trắc địa công trình, công nghệ GPS đ đợc ứng dụng để thành lập lới trắc địa thi công cho số công trình xây dựng lớn nh công trình cầu B i cháy, cầu Thanh trì, hầm đờng Hải vân, thuỷ điện Sơn La.v.v Tại công trình này, công nghệ GPS đ đợc ứng dụng để thành lập lới khống chế thi công, nhng cha đề cập đến vấn đề thiết kế lới GPS nhằm bảo đảm yêu cầu công tác trắc địa công trình.Việc tính chuyển toạ độ độ cao điểm đo GPS hệ tọa độ độ cao thiết kế công trình cha đợc nghiên cứu cách đầy đủ nên thực tế sản xuất gặp nhiều lúng túng, khó khăn Về lý thuyết : Hiện nớc ta có số đề tài nghiên cứu khoa học cấp lĩnh vực với nghiên cứu ứng dụng công nghệ GPS trắc địa công trình nói chung [16], ứng dụng công nghệ GPS xây dựng nhà cao tầng nói riêng [28] Trong trắc địa công trình, lới khống chế thi công công trình quan trắc chuyển dịch biến dạng công trình có yêu cầu độ xác rÊt cao.Trong ®ã ë n−íc ta, th−êng míi sư dụng lịch vệ tinh quảng bá mà sử dụng lịch vệ tinh xác để xử lý số liệu đo GPS Việt Nam phần mềm xử lý xác số liệu đo GPS nh phần mềm GAMIT đợc sử dụng hạn - 190 - Bảng 13 Chênh lệch toạ độ độ cao điểm D ca quan trắc TT X(mm) Y(mm) ∆H(mm) +1.4 -2.8 -4.6 +0.1 +2.0 +5.3 +0.4 -0.3 - 1.1 -1.8 +2.6 +11.3 0.0 0.0 + 0.3 0.0 0.0 - 0.3 +0.1 +0.1 +1.2 Bảng 14 Toạ độ độ cao ®iĨm MDC thêi gian quan tr¾c TT X(m) Y(m) H(m) mx my mP (mm) (mm) (mm) 2330307.6389 580364.9121 -3.6617 0.141 0.132 0.193 2330307.6335 580364.9082 -3.6299 0.090 0.130 0.158 2330307.6368 580364.9042 -3.6102 0.077 0.135 0.155 2330307.6343 580364.9064 -3.6304 0.545 0.761 0.936 2330307.6505 580364.9102 -3.6478 0.142 0.144 0.202 2330307.6527 580364.9092 -3.6526 0.063 0.065 0.090 2330307.6538 580364.9085 -3.6462 0.389 0.352 0.524 2330307.6536 580364.9087 -3.6409 0.042 0.038 0.056 - 191 - Bảng15 Chênh lệch toạ độ độ cao điểm MDC ca quan tr¾c TT ∆X(mm) ∆Y(mm) ∆H(mm) -5.4 -3.9 -31.8 +3.3 - 4.0 -19.7 - 2.5 +4.0 -20.2 +16.2 +3.8 -17.4 + 2.2 -1.0 +4.8 + 1.1 - 0.7 - 6.4 - 0.2 +0.2 -5.3 - 192 - Phụ lục Chơng trình tính øng dông - 193 - $DEBUG $LARGE C PROGRAME THIET KE LUOI GPS Implicit real *8(A-H,O-Z) Character FF*30,FF1*10,FF2*10,FF3*10,FF4*10,FF5*10,sl*50 Character * 10 TEN(30) Character TCT*30 Common/H1/X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/V(100),QF(100),STP(100),Fa(100) Dimension mauso(100),St(100),STPa(100),Va(100) PI=3.141592653589793D0 RO=180*3600/PI PI2=PI*2 open(11,file='c:\dl',status='old') read(11,'(a)') sl Open(1,FILE=sl,STATUS='OLD') Open(4,FILE='c:\tkchung.kq',STATUS='NEW') Read(1,1)TCT Format(A) Write(4,2)TCT Format(//10x'KET QUA UOC TINH DO CHINH XAC:',A30) Write(4,3) Format(18x'==================&==================') Write(4,*) Read(1,*) nc,npv,ndm,ndg Read(1,*) as,bs,apv,bpv Do j=1,(ndm+ndg) read(1,*) j,x(j),y(j) continue Do i = 1,(ndm+ndg) Read(1,6)TEN(i) Format(A) Continue l=2*ndm mt= ndm+ndg Open(2,File = 'Dcx.dat',Status = 'New') Do i=1,nc read(1,*) i,id,ic,dn Call dshcc(l,id,ic,si) Write(2,'(3i3,f12.3,(50f10.4)') i,id,ic,si,(a(j),j=1,l) dms=(as**2+(bs*(si/1000))**2)/dn dms1 = dsqrt(dms) p =1/dms Call dptc(l,p) continue Write(4,*) Write(2,*) Do 10 i=1,npv Read(1,*) i,ia,ib,dn Do 43 j=1,l IF(j.EQ.ia) Then Call dshcgpv(Ro,l,ia,ib,si1) Fa(j*2-1) = a(j*2-1) Fa(j*2) = a(j*2) Else Fa(j*2-1)= 0.d0 Fa(j*2) = 0.d0 End IF - 194 - 43 Continue cds1 = si1/1000 dma=((apv*apv)+((bpv/cds1)*(bpv/cds1)))/dn dma1 =dsqrt(dma) p=1/dma Write(2,'(3i3,50f10.3)')i,ia,ib,(Fa(j),j=1,l) Call dptc(l,p) 10 continue Call DNDMT(l) Write(2,*) Write(4,*) Write(4,15) 15 format(20x'I/CHI TIEU KY THUAT CUA LUOI' */20x, ' & -') Write(4,16)mt 16 format(18x'1-Tong so diem luoi :',i3) Write(4,17)ndg 17 format(18x'2-So diem goc :',i3) Write(4,18)ndm :',i3) 18 format(18x'3-So diem moi Write(4,19)nc 19 format(18x'5-So canh :',i3) Write(4,20)npv 20 Format(18x'6-So phuong vi :',i3) Write(4,*) Write(4,*) Write(4,21) 21 format(24x 'II/DIEM GOC LUOI GPS') Write(4,22) 22 format(15x'+ -+ -+' */15x, '| | TOA DO |' */15x, '| DIEM | -+ -|' */15x, '| | X(m) | Y(m) |' */15x, '+ -+ -+ -|') Do 23 i=(ndm+1),(ndm+ndg) 23 Write(4,24)TEN(i),X(i),Y(i) 24 format(15x'|'A8,1x,'|'F13.3'|'F13.3'|') Write(4,25) 25 format(15x'+ -+ -+ -+') Write(4,*) Write(4,*) Write(4,26) 26 format(27x'III/TOA DO THIET KE') Write(4,27) 27 format(18x '+ + +' */18x, '| | TOA DO |' */18x, '| DIEM | -+ |' */18x, '| | X(m) | Y(m) |' */18x, '| + -+ |') 28 j=1,ndm Write(4,29)TEN(j),X(j),Y(j) 29 format(18x'|'A6'|'F13.3'|'F14.3'|') 28 continue Write(4,30) 30 format(18x '+ + -+ +') write(4,*) Write(4,31) 31 format(20x'IV/DO CHINH XAC DO DAC DU KIEN', - 195 - */20x, ' ') Write(4,32)as,bs 32 format(6x'SSTP canh luoi GPS : Ms = ' *,f3.1,2x'(mm) +',f3.1,2x'ppm') Write(4,*) Write(4,34) 34 format(16x'V/SAI SO TRUNG PHUONG XAC DINH VI TRI DIEM' */16x, ' *** ') Write(4,35) 35 format( * 8x,'+ + -+ -+' */8x,'| | TOA DO THIET KE | SAI SO VI TRI DIEM |' */8x,'| DIEM| -+ -+ + + -|' */8x,'| | X(m) | Y(m) | Mx | My | Mp(cm)|' */8x,'| + -+ -+ + + -|') smmax = 0.d0 36 j=1,ndm j1=(2*j)-1 j2=(2*j) smx= dsqrt(Q(j1,j1)) smy= dsqrt(Q(j2,j2)) sm=sqrt(SMX**2+SMY**2) write(4,37)TEN(j),X(j),Y(j),smx,smy,sm 37 format(8X'|'A6'|'F13.3'|'F13.3'|'F6.2'|'F6.2'|'F6.2' |') IF (sm.GT.smmax) Then smmax = sm FF1 = TEN(j) End IF 36 Continue Write(4,38) 38 format( *8x'+ + + + -+ -+ -+') Write(4,*) Write(4,39) 39 Format(25x'VI/BANG TUONG HO VI TRI DIEM') Write(4,40) 40 Format(25x' -*** ') Open(2,File ='Dcx.dat',Status='Old') msomin = 100000000 Do 45 i= 1,nc Read(2,50) i,id,ic,si,(a(j),j=1,l) St(i) = si 50 Format(i3,i3,i3,f12.3,50f10.4) Do 46 k = 1,l S = 0.d0 Do 47 m = 1,l 47 S = S +Q(k,m)* a(m) 46 V(k) = S S1 = 0.d0 Do 150 i1 = 1,l 150 S1 = S1 + a(i1)* V(i1) QF(i)= S1 STP(i)= dsqrt(QF(i))/100 mauso(i) = INT(St(i)/STP(i)) IF(mauso(i).LE.msomin) Then msomin = mauso(i) FF2 = TEN(id) FF3 = TEN(ic) - 196 - End IF Continue stpamax = 0.d0 Read(2,*) Do 61 i= 1,npv Read(2,62) i,ia,ib,(Fa(j),j=1,l) 62 Format(i3,i3,i3,50f10.4) Do 63 k = 1,l Sa = 0.d0 Do 64 m = 1,l 64 Sa = Sa +Q(k,m)* Fa(m) 63 Va(k) = Sa Call dnhan(i,Fa,Va,l) Rewind(2) Write(4,66) 66 format( *8x, '+ + -+ -+ -+ -+ + +' */8x,'| TT | Diem | Diem | S | Ms | Ms/s | Ma |' */8x,'| | dau | cuoi | (m) | (m) | | (") |' */8x,'| + -+ -+ -+ -+ + |') Do 67 i= 1,nc Read(2,69)i,id,ic 69 Format(i3,i3,i3) Write(4,68)i,TEN(id),TEN(ic),St(i),STP(i),mauso(i),STPa(i) 68 format(8x'|'i3,1x'|' A7'|'A7'|'f7.0'|'F6.4,1x'|','1:',i8'|' *f5.2,1x'|') 67 Continue Write(4,70) 70 Format( *8x'+ + -+ -+ -+ -+ + +') Write(4,*) Write(4,80) 80 Format(1x,'VII/KET QUA UOC TINH DCX LUOI GPS' */1x, ' ') Write(4,81) 81 Format( *1x,'1/ Sai so trung phuong so don vi : M =1.0') Write(4,*) Write(4,82) FF1,smmax 82 Format( *1x,'2/SSTP vi tri diem yeu nhat:' */10x,'Diem -',2x,A7,5x,': Mp =',f7.4,'(cm)') Write(4,83) FF2,FF3,msomin 83 Format( *1x,'3/SSTP tuong doi canh yeu nhat:', */10x,'Canh:',A7,'-',A7,': 1/',i8) Write(4,84) FF4,FF5,stpmax 84 Format( *1x,'4/SSTP phuong vi yeu nhat:', */10x,'Canh:',A7,'-',A7,': Ma =',f8.2,'"') Write(4,*) Write(4,*) Write(4,85) 85 Format(10x,'Tinh theo chuong trinh "GP.AGE"', */10x, ' -') Write(*,*) Write(*,*) Write(*,71) 45 - 197 - 71 c c c c Format(20x,'DA CHAY XONG CHUONG TRINH') Write(*,*) Write(*,72) 72 Format( *' Copyright(C) 2006 by:'/ *' Tran Viet Tuan'/ *' Khoa Trac Dia - Dai hoc Mo -Dia chat'/ *' -*** ') STOP END -subroutine dshcc(l,id,ic,si) implicit real *8(a-h,o-z) common/H1/ X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/V(100),QF(100),STP(100),QFa(100) i=1,l A(i)=0.d0 continue px=X(ic)-X(id) py=Y(ic)-Y(id) si=dsqrt(px**2+py**2) A(id*2-1)= px/si A(id*2)= py/si A(ic*2-1)= -px/si A(ic*2)= -py/si return end -subroutine dshcgpv(Ro,l,ia,ib,si1) implicit real *8(a-h,o-z) common/H1/X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/V(100),QF(100),STP(100),QFa(100) i=1,l A(i)=0.d0 continue px=X(ib)-X(ia) py=Y(ib)-Y(ia) si1=dsqrt((px*px)+(py*py)) A(ia*2-1)= Ro*py/((si1**2)*100) A(ia*2)= -Ro*px/((si1**2)*100) A(ib*2-1)= -Ro*py/((si1**2)*100) A(ib*2)= Ro*px/((si1**2)*100) return end -subroutine dptc(l,p) implicit real *8(a-h,o-z) common/H1/X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/V(100),QF(100),STP(100),QFa(100) i=1,l j=i,l R(i,j)= R(i,j)+A(i)*A(j)*p continue continue return end Subroutine DNDMT(l) - 198 - 100 11 10 12 13 15 18 17 Implicit real *8(a-h,o-z) Common/H1/X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/V(100),QF(100),STP(100),QFa(100) logical W(100) Dimension B(100),C(100,100),Y1(100) Do i = 1,l W(i) = True S = 0.0 i1 = i-1 i2 = i+1 IF(i.EQ.1) goto Do k=1,i1 IF (W(i)) goto S = (S+ C(k,i)*C(k,i))*1 Goto S = (S - C(k,i)* c(k,i))*1 Continue S= S+R(i,i) IF (S.GE.0) Goto S = - S W(i)= False IF (S.LT.1E-25) Goto Goto Write(4,100) Format(/10x'MA TRAN SUY BIEN'/) Goto 99 C(i,i) = dsqrt(S) IF (i.EQ.l) Goto Do k=i2,l S=0.d0 IF(i.EQ.1) Goto 10 Do j=1,i1 IF (W(i)) Goto 11 S = (S+C(j,i)*C(j,k))*1 Goto S = (S-C(j,i)*C(j,k))*1 Continue S = S + R(i,k) IF (.NOT.W(i)) Goto 12 Goto 13 S = -S IF (C(i,i).EQ.0.d0) Goto C(i,k)= S/C(i,i) Continue Continue Do 14 L1= 1,l Do 15 i = 1,l X1(i) = 0.d0 Y1(i) = 0.d0 B(i) = 0.d0 B(L1) = Do 16 i = 1,l S = 0.d0 IF(i.EQ.1) Goto 17 i1 = i- Do 18 j = 1,i1 S = ( S - C(j,i) * Y1(j)) * IF (C(i,i).EQ.0) Goto 16 - 199 - Y1(i) = (B(i) + S)/ C(i,i) Continue Do 19 i2 = 1,l i = l +1 -i2 i1 = i +1 S = 0.d0 IF (i.EQ.n) Goto 20 Do 21 j = i1,l 21 S = ( S -C(i,j)* X1(j))*1 20 IF (C(i,i).EQ.0) Goto 14 X1(i) = (Y1(i) + S)/C(i,i) 19 Q(i,L1)= X1(i) 14 Continue Goto 999 99 Stop 999 Return End c Subroutine dnhan(i,Fs,Vs,l) implicit real *8(a-h,o-z) common/H1/ X(100),Y(100),A(200),R(100,100),Q(100,100),X1(100) Common/H2/ V(100),QF(100),STP(100),Fa(100) Dimension Fs(100),Vs(100) S=0.d0 Do k = 1,l S = S + Fs(k)*Vs(k) QF(i) = S Return End 16 - 200 - DEBUG $LARGE C PROGRAME TINH CHUYEN TOA DO LUOI GPS Implicit Real*8(a-h,o-z) Character FF*30,FF1(50)*10 COMMON/H1/ X(100),Y(100),Z(100),B(100),XL(100),XH(100) DIMENSION b1(50),b2(50),b3(50),a1(50),a2(50),a3(50) DIMENSION x2(50),y2(50),z2(50) Common/v1/ HV(8000),ptc(30000) bb=1 pi=4*datan(bb) ro=180*3600/pi pi2=2*pi Write(*,1) Format( *18x,'TINH CHUYEN TOA DO GPS VE HE TOA DO CONG TRINH'/ *18x,'=================***********=================='/) aa=6378137 bb=aa*(1-1/298.257223563) e2=(aa*aa-bb*bb)/(aa*aa) eph2=(aa*aa-bb*bb)/(bb*bb) Write(*,'(a\)')' * Nhap ten file so lieu:' Read(*,'(a)')FF Open(1,file=ff,status='old') Write(*,'(a\)')' * Nhap ten file ket qua:' Read(*,'(a)')FF Open(4,file = FF,status='new') Read(1,41) FF 41 Format(A) Write(4,5) FF Format(/20x'TEN CONG TRINH:',A) Write(4,6) Format(20x'============***============') Read(1,*) ndt Do i = 1,ndt Read(1,*) i,X(i),Y(i),Z(i) xt=xt+x(i) yt=yt+y(i) zt=zt+z(i) Continue Write(4,30)xt,yt,zt 30 Format(3f20.6) Write(4,*) Write(4,9) Format(19x,'KET QUA TINH CHUYEN TOA DO GPS ') Write(4,13) 13 Format(19x' -& -') Write(4,*) Write(4,12) ndt 12 Format(10x,'I/So luong diem tinh :',i5,' diem') Write(4,*) Write(4,14) 14 Format(10x'II/TOA DO VUONG GOC KHONG GIAN DIA TAM' */2x,'+ -+' */2x,'| No |Ten diem| X (m) | Y (m) | Z (m) |' */2x,'+ + + -+ -+ -+') Do 15 i= 1,ndt Read(1,'(a)') FF1(i) - 201 - Write(4,16)i,FF1(i),X(i),Y(i),Z(i) Format(2x'|'i3,1x'|',2x,A6'|'f15.4'|',f13.4,'|'f13.4'|') Continue Write(4,17) 17 Format( *2x '+ -+') Xo=xt/ndt Yo=yt/ndt Zo=zt/ndt Write(4,*) Write(4,31)Xo,Yo,Zo 31 Format(3f20.6) p=dsqrt(Xo*Xo +Yo*Yo) bt=datan(Zo/p/(1-e2)) 27 w=dsqrt(1-e2*(dsin(bt)**2)) xhh=(p/cos(bt))- xn bs=datan((Zo)/p*(xn+xhh)/(xn+xhh-e2*xn)) if(ep.le.0.00001/ro) Goto 28 bt=bs Goto 27 28 Bo=bs XLo=datan(Yo/Xo) IF(XLo.lt.0) XLo = XLo + 180 *3600/ro XHo =(p/dcos(Bo))-(aa/dsqrt(1-e2*(dsin(Bo)**2))) Write(4,*) Write(4,32)Bo,XLo,XHo 32 Format(2x'Bo=',f12.6,'Xlo=',f12.6,'XHo=',f12.6) Do 43 i=1,ndt p=dsqrt(x(i)*x(i)+y(i)*y(i)) bt=datan(z(i)/p/(1-e2)) 20 w=dsqrt(1-e2*(dsin(bt)**2)) xhh=(p/cos(bt))-xn bs=datan(z(i)/p*(xn+xhh)/(xn+xhh-e2*xn)) ep=bs-bt if(ep.le.0.00001/ro) Goto 19 bt=bs Goto 20 19 B(i)=bs Write(4,*) Write(4,34)bs 34 Format(2x'Bs='f12.8) XL(i)=datan(Y(i)/X(i)) if(XL(i).lt.0) XL(i)=XL(i)+180*3600/ro w1 = 0.d0 xn1 = 0.d0 w1=dsqrt(1-e2*(dsin(Bo)**2)) XH(i)=(p/cos(bs))- xn1 a1(i) = -sin(bs)* cos(XL(i) a2(i) = -sin(bs)* sin(XL(i) a3(i) = cos(bs) b1(i) = X(i) - (xn1 + XHo)* cos(Bo)* cos(XLo) b2(i) = Y(i) - (xn1 + XHo)* cos(Bo)* sin(XLo) b3(i) = Z(i) - (xn1*(1-e2)+ x2(i)=-sin(bs)*cos(xl(i))*(x(i)-(xn1+XHo))*cos(Bo)*cos(XLo) *-sin(bs)*sin(xl(i))*(y(i)-((xn1+XHo)*cos(Bo)*sin(XLo)))+cos(bs) **(z(i)-((xn1*(1-e2)+XHo)*sin(Bo))) y2(i)=-sin(XL(i))*(X(i)-(xn1+XHo)*cos(Bo)*cos(XLo))+cos(XL( *i))*(Y(i)-(xn1+XHo)*cos(Bo)*sin(XLo)) 16 15 - 202 - z2(i)=cos(bs)*cos(XL(i))*(X(i)-(xn1+XHo)*cos(Bo)*cos(XLo)) *+cos(bs)*sin(XL(i))*(Y(i)-(xn1+XHo)*cos(Bo)*sin(XLo)) *+cos(bs)*sin(XL(i))*(Z(i)-((xn1*(1-e2)+ XHo)*sin(Bo))) 43 Continue Write(4,*) 33 Format(5x'w1=',f14.8,'xn1=',f14.8) Write(4,21) 21 Format(16x' TOA DO VUONG GOC KHONG GIAN DIA DIEN' */10x, '+ -+' */10x, '| No |Ten diem| x (m) | y (m) | z (m) |' */10x, '+ -+ + + + +') Do 22 i =1,ndt Write(4,23)i,FF1(i),x2(i),y2(i),z2(i) 22 Continue 23 Format(10x'|'i3,1x'|',2xA6'|'f12.4'|'f12.4'|'f12.4'|') Write(4,24) 24 Format(10x'+ +') Write(4,25) 25 Format( 20x 'TOA DO THI CONG ' */10x, '+ +' */10x, '| No |Ten diem| x (m) | y (m) |' */10x, '+ -+ + + +') Write(4,26) 26 Format(10x'+ -+') STOP END - 203 - - 204 - ... nghiên cứu Luận án25 Chơng - Nghiên cứu khả ứng dụng công nghệ GPS trắc địa công trình Việt nam 2.1.Khả ứng dụng công nghệ GPS giai đoạn quy hoạch - thiết kế công trình 26 2.2.Khả ứng dụng công nghệ. .. bao gồm: trắc địa công trình thành phố; trắc địa công trình công nghiệp; trắc địa công trình đờng sắt, đờng bộ; trắc địa công trình cầu vợt; trắc địa công trình ngầm; trắc địa công trình thuỷ... ngang công nghệ GPS, đo đạc thực nghiệm để kiểm định khả ứng dụng GPS quan trắc chuyển dịch ngang công trình - 26 - Chơng nghiên cứu khả ứng dụng công nghệ GPS trắc địa công trình việt nam 2.1.Khả