侯學慧

新疆警察學院信息安全工程系， 新疆 烏魯木齊 830011

自東漢初年建成第一座塔約兩千年以來，中華大地上出現(xiàn)的塔數(shù)逾萬，雖然歷經(jīng)自然淘汰和人為損毀，據(jù)說尚余三千，僅此數(shù)也號稱世界之最.古塔由于年代久遠，長時間承受自重、氣溫、風力等各種作用，偶然還要受地震、颶風和洪澇的影響，古塔會產(chǎn)生各種變形，諸如傾斜、彎曲和扭曲等[1～4].本文探究古塔受外界環(huán)境影響而產(chǎn)生變形、彎曲、扭曲等情況并將其量化得到具體的數(shù)據(jù)以便預測未來許多年內(nèi)，古塔的變形情況或者考慮古塔新建立未變形時的具體位置和高度等.

開封鐵塔已有上千年歷史，是我國重點保護文物.管理部門委托測繪公司先后于1986年7月、1996年8月、2009年3月和2011年3月對該塔進行了4次觀測(圖1).

根據(jù)圖1提供的4次觀測數(shù)據(jù)，本文主要討論以下3個問題，并通過數(shù)據(jù)理論分析該古塔的變形情況：

(1)給出確定古塔各層中心位置的通用方法，并列表給出各次測量的古塔各層中心坐標；

(2)分析該塔的傾斜、彎曲和扭曲等變形情況；

(3)分析該塔的變形趨勢.

1 問題分析

中國古塔成近千年，歷盡滄桑，且不具備再生性.近千年來遭受到地震、冰雹、風災和水患等自然災害，尤其是在戰(zhàn)爭時期曾用飛機、大炮進行轟炸等一系列人為因素，歷經(jīng)滄桑，不少古塔身患重病、年久失修，導致古塔發(fā)生了各種形變，例如傾斜、彎曲、扭曲、下沉、構(gòu)件斷裂等，但這里我們目前只考慮傾斜、彎曲、扭曲這3種主要的形變.

1.1 確定古塔各層中心位置的通用方法

確定三種解決方案：

圖1 用Mathmatic 計算每一層中心(圓心)坐標
Fig.1 Use Mathmatic to calculate the center coordinate of each layer

(1)我們利用Mathmatica作圖可知，給出的每層8個點是出于頂點的位置，因此，我們可以用8個點的平均值點來代替中心位置(采用平均值法).

(2)中國古代建筑是對稱結(jié)構(gòu)，塔面多數(shù)采用偶數(shù)邊的正多邊形.因此采用求外接圓的圓心代替中心.

(3)用軟件計算出圓心.

1.2 古塔的變形程度

對于變形問題我們只考慮3個主要因素：傾斜，扭曲和彎曲.解決時量化這3個量，求其變化角度.注意：給出的數(shù)據(jù)在2009年到2011年的坐標對應發(fā)生了變化，需要調(diào)整.

傾斜程度：構(gòu)造每層中心點與古塔初建造時地面中心點的直角三角形，在直角三角形當中求解傾斜角.

扭曲程度：我們利用軟件量化扭曲程度，將古塔同一標號的13個點在地面上的投影做函數(shù)擬合成直線(或者用一元線性回歸方程)，從而分析擬合后的投影直線與變化后的投影直線的夾角，即為扭轉(zhuǎn)角度.

1.3 分析古塔的變化趨勢

我們采用數(shù)據(jù)直線擬合的方法，根據(jù)古塔的變化角度，分析其變化量，建立擬合函數(shù)方程，根據(jù)古塔的變形程度得到的數(shù)據(jù)，確定函數(shù)，從而預測n年以后，古塔的變化趨勢.

2 古塔變形問題的解決

首先假設：

(1)古塔的傾斜程度，即傾斜角記為θ；

(2)古塔的扭曲程度，即傾斜角記為α；

(3)古塔同一標號位置13個點在地面上投影點集，擬合成的直線記為L.

2.1 確定古塔各層中心

2.1.1 確定古塔各層中心位置的通用方法

(1)方法一：求平均值法(用8個點的平衡點來代替中心點)

我們采用近似值算法，因為中國古代的古塔多采用對稱結(jié)構(gòu)，多數(shù)塔層都是由偶數(shù)邊的正多邊形構(gòu)成的，這樣的正多邊形是相對對稱均勻的.因為已知點的坐標，我們通過Mathmatice作圖得到圖2，其中給出的點是相對的頂點.因為已知的8個點是對稱的，我們用這8個點坐標的平均值來代替每層塔面的中心位置.

(2)方法二： 將內(nèi)切圓的圓心視為每一層塔面的中心點

用平均值求近似值，相對的誤差較大，可改為采用內(nèi)切圓求其圓心.

因為已有的8個點是在頂點位置，古代對于一個建筑物的建造總是以圓為內(nèi)切圓或者外接圓建立的正多邊形，因此可以認為這8個點一定是在某一個圓心為(a,b)，半徑為r的圓上.這樣問題就轉(zhuǎn)化為已知圓上8個點的坐標，求出圓心即為中心.

圖2 Mathmatice軟件作出的古塔頂點在地面的投影Fig.2 To make the Graph of the vertex projection in the Ground

圖3 近似八角古塔
Fig.3 Approximate Eight-Square

(3)利用C語言中的Mathmatica軟件編寫程序

求解問題轉(zhuǎn)化為：已知圖上的多個點的坐標，計算該圓的圓心與半徑 ，從而確定正多邊形的中心.

2.1.2 列表給出各層量古塔各層的中心坐標

(1)依據(jù)上述方法，利用EXCLE、Mathmatic軟件，得到各層中心點的列表.如表 1所示.

(2)利用圓心去近似代替中心，我們給出1986年第一層的中心的計算方法，并得到中心坐標為(566.664 8,522.710 5,1.787 4)，其余各層方法相同.

2.2 分析古塔的傾斜、彎曲和扭曲等變形情況

2.2.1 古塔的傾斜程度分析

坐落在開封市鐵塔是一座歷史悠久的古塔，因此在1986年時，由于各種環(huán)境等因素，使古塔本身發(fā)生變形，造成傾斜.由于我們不了解古塔的背景與建造時間，我們只重點考慮1986年以后，并且假設在1986年時古塔的變形情況不嚴重，因此我們假設，將1986年古塔第一層中心的投影點，當作古塔建造時塔地面的中心O點(存在誤差).

表1 給出確定古塔各層中心坐標Tab.1 Give Tower the center coordinates of each layer

那么我們可以近似用坐標O(566.664 8, 522.710 5,0)(第一層的中心)來代替建筑物的塔底中心點.由于傾斜后，每層的中心點與地面的中心點O近似構(gòu)成相似的直角三角形,那么我們可以在直角三角形當中，記O1(x1,y1,z1)，O0(x1,y1,z1),從而通過公式計算傾斜角公式

已知4年的塔頂坐標，從而近似確定這4年的塔頂中心坐標，因1986年和1996年的塔頂4個坐標都存在，因此直接求出平均中心點作為塔頂中心，而2009年和2011年的塔頂4個坐標不全，因此近似用第13層的中心點來作為塔頂?shù)闹行妮S位置，我們根據(jù)算出的四年的塔頂中心點與塔底中心點的連線所在的直線與豎直方向存在傾斜角θ，利用EXCEL計算每一年的塔頂中心點坐標(表2)，根據(jù)傾斜角公式，我們近似得出1986年～2011年相應的傾斜角.表3給出傾斜角的角度值.

表2 每一年的塔頂中心點坐標Tab.2 Vertex coordinate of the top of tower in each year

表3 古塔傾斜角的度數(shù)Tab.3 Degree of inclination of tower

2.2.2 古塔的扭曲程度分析

因為我們對扭曲程度的變化，注重其點在地面上的投影，因此建立直角坐標系XOY，扭曲程度可分為2個方面：為了說明扭曲過程，我們都以1986年為基準.看1996年～2011年發(fā)生了多大程度的扭轉(zhuǎn).

(2)另一方面是1986年～2011年間古塔發(fā)生的扭轉(zhuǎn)程度.我們選擇計算出同一年當中同一個標號位置的不同層面的13個點在地面上投影所成的近似直線，考慮變化后直線的傾斜角.

這是標號均為1位置的1986年的13個頂點，畫出在地面上的投影，用Mathmatice輸入坐標，求出擬合的一次函數(shù)記為L1:y=k1x+b1；得到k1=-3.134 11.同上，我們選取標號為1的1996年的13個頂點在地面上的投影，擬合直線為L2:y=k1x+b1；得到k2x+b2,k2=-3.122 54，那么我們計算出L1,L2的夾角.

通過計算得出為1986年～1996年的第一個位置扭轉(zhuǎn)的角度為0.061 5°.用相同的方法得到1986年～2011年對應的每個頂點的扭轉(zhuǎn)度.在第一問題當中，同一年份的扭轉(zhuǎn)角度算法相同，由于計算量較大，就不一一列出.

(3)彎曲程度

我們在考慮古塔彎曲時，中心軸會出現(xiàn)彎曲，可以采用根據(jù)已知點的坐標，擬合函數(shù)，計算中心軸的曲率，由Mathmatica作圖可知，具體計算出較為明顯彎曲的位置.也可以采用求解扭轉(zhuǎn)角度的辦法，計算擬合直線的夾角.在理想狀態(tài)下，棱上的13個頂點，應在同一直線上，但是一旦出現(xiàn)明顯彎曲，那么任意兩點確定的直線應該會出現(xiàn)夾角，通過夾角判斷其彎曲的具體位置.

2.3 近似用線性關系分析未來古塔的變形情況

傾斜程度的分析：1986年～1996年這十年，傾斜角增加了0.000 18°,平均每年增加0.000 018°.這十年的變化較為平穩(wěn).1996年～2009年這十三年的變化是比較迅速的，傾斜角增加了0.001 43°，平均每年增加0.000 143°，2009年～2011年這三年可能是由于時間較短，傾斜角僅增加了0.000 03°，預測往后的變化會越來越慢.

由此分析古塔的傾斜程度在1986年～1996年這十年間變化較慢；1996年～2009年這十三年間，變化較快，傾斜程度明顯；2009年以后，基本趨于平穩(wěn)，最后會達到一個上界：是古塔的最大傾斜角度.我們做一個短期預測，假設傾斜角隨時間的函數(shù)是線性的，我們給出傾斜角的變化圖.類似的分析扭曲角度的值，預測未來古塔的扭轉(zhuǎn)程度.

圖4 古塔發(fā)生扭轉(zhuǎn)后頂點在地面投影Fig.4 Guta after the occurrence of a twisted vertex in the ground projection

圖5 古塔發(fā)生扭轉(zhuǎn)后每一層邊在地面的投影Fig.5 The projection of each side of the ground in Guta after twisting

圖6 傾斜角的變化圖Fig.6 Variation of inclination angle

圖7 扭轉(zhuǎn)角的變化圖
Fig.7 Variation of torsion angle

3 古塔模型的改進與評價

3.1 方案改進

(1)找到古塔開始建立的初始位置，計算出1986年的傾斜角；其余3個年份加入差額，即得出相對準確的傾斜角.(2)在求解古塔中心坐標時，可以利用C語言或Matlab 編程求出正八邊形外接圓的圓心，中心的選取更為精確.(3)在扭曲程度上，可以有多種不同衡量的扭曲程度的量，可以進一步研究，準確到南北東西，方向各轉(zhuǎn)動的度數(shù)，水平位移多少毫米.(4)彎曲上的測量，可以具體到量化，由于數(shù)據(jù)量過大，沒有給出明確的數(shù)據(jù)解釋.

3.2 對模型的評價

(1)優(yōu)點：更為精確地給出了傾斜角，扭轉(zhuǎn)角度的變化.(2)缺點：因為平衡位置是建塔的建造位置；但1986年已經(jīng)開始變形，所以從1986年到2011年，1986年的傾斜角度存在偏差.

模型的評價：首先我們對數(shù)據(jù)進行了相關性分析、處理和量化，并且提出多種方法來解決并優(yōu)化問題.同時利用Mathmatica軟件進行計算，可信度高.論文中的圖形與數(shù)據(jù)相結(jié)合更具有說服性.雖然此方法不是最優(yōu)方法，但經(jīng)過對模型擬合，數(shù)字量化已將誤差減少到最少.