劉靜

泰山職業(yè)技術(shù)學(xué)院,山東泰安271000

基于模型構(gòu)建的古塔變形研究

劉靜

泰山職業(yè)技術(shù)學(xué)院,山東泰安271000

根據(jù)在1986年、1996年、2009年、2011年對古塔的實(shí)際觀測數(shù)據(jù)，通過Matlab作圖、編程等功能，以及Excel軟件計(jì)算分析相關(guān)數(shù)據(jù)，給出數(shù)據(jù)分析的結(jié)果并作出了相應(yīng)的圖形。首先利用Matlab軟件繪制出1986年與2011年古塔形狀的圖形，然后利用模型構(gòu)建并求解，累加模型確定古塔各層次中心位置，夾角模型描繪古塔的傾斜、彎曲、扭曲等一系列變形情況，最后給出古塔變形趨勢。本文采用模型構(gòu)建方法研究古塔變形，計(jì)算分析合理，提供了較精確的結(jié)果，對于維護(hù)古建筑物的安全性有著重要的意義。

古塔變形;Matlab和Excel軟件;模型構(gòu)建;夾角模型

文物是國家的不可再生的人類創(chuàng)造的文化資源，它們不僅具有非常高的科學(xué)藝術(shù)文化價(jià)值，而且是一種珍貴的實(shí)物資料。為了有效的保護(hù)文物古塔，要對古塔的平面、高度和傾斜度等進(jìn)行定期監(jiān)測，使管理部門及時(shí)掌握古塔變形情況，從而制定出針對性的措施，對古塔進(jìn)行更好的高效率的保護(hù)。古塔高13層，每層有8個(gè)觀測點(diǎn)，4次觀測數(shù)據(jù)，為了對古塔科學(xué)系統(tǒng)性的研究與保護(hù)，不同學(xué)者從曲線擬合、最小二乘原理、橢圓等擬合模型、灰色預(yù)測模型、曲率撓率等不同角度進(jìn)行了分析預(yù)測。本文主要采用模型構(gòu)建的思想，建立了數(shù)學(xué)模型并借助于數(shù)學(xué)軟件Matlab與Excel對數(shù)據(jù)計(jì)算、編程、作圖、分析等比較全面地研究問題。首先做出了古塔的三維立體圖形及塔底的平面圖形，直觀地掌握古塔的大體形狀；其次比較分析了確定古塔各層次中心位置的方法，應(yīng)用計(jì)算精度較高的平面薄板的重心坐標(biāo)公式構(gòu)建累加模型，結(jié)合Matlab編程，確定平面中心坐標(biāo)；接著構(gòu)建夾角模型，求出古塔傾斜、彎曲與扭曲的角度；最后進(jìn)行垂直位移檢測，給出變形趨勢。模型的構(gòu)建具有較高的擬合精度和適用性。

1　模型準(zhǔn)備與假設(shè)

1.1缺失數(shù)據(jù)的處理

處理方法：（1）刪除法：直接刪除有缺失數(shù)據(jù)的記錄；（2）填補(bǔ)法：包括平均值填充法，人工填寫法，就近補(bǔ)齊法等。（3）不處理：對數(shù)據(jù)直接分析或構(gòu)建模型。對于第1次和第2次測量中第13層第5點(diǎn)坐標(biāo)的缺失，采用了人工填寫和就近補(bǔ)齊相結(jié)合的方法，用第4個(gè)監(jiān)測點(diǎn)的發(fā)展趨勢，結(jié)合第3次第13層第5點(diǎn)的坐標(biāo)填補(bǔ)第1次和第2次中第13層第5個(gè)監(jiān)測點(diǎn)的坐標(biāo)[1]。

1.2模型假設(shè)

模型的假設(shè)有（1）不考慮測量數(shù)據(jù)過程中所產(chǎn)生的誤差；（2）忽略分析計(jì)算中出現(xiàn)的誤差;(3)以東西方向?yàn)閄軸，其正向?yàn)闁|向；南北方向?yàn)閅軸，其中北向?yàn)檎颉?/p>

2　古塔及其底部的形狀

利用給出的測量數(shù)據(jù)，使用數(shù)學(xué)軟件Matlab繪制圖像，做出了1986年和2011年古塔的形狀如圖1、圖2所示；圖3給出了1986年塔底的形狀為凸多邊形[2]。

圖1　1986年古塔的形狀Fig.1 Ancient pagoda’s shape in 1986

圖2　2011年古塔的形狀Fig.2 Ancient pagoda’s shape in 2011

圖3　1986年古塔底部的形狀Fig.3 The shape of ancient pagoda’s bottom in 1986

3　模型的建立與求解

3.1古塔各層次中心位置模型

各層次中心位置的確定方法：（1）擬合法：形狀規(guī)則的塔，可通過各點(diǎn)的坐標(biāo)先計(jì)算出多邊形的角度，再進(jìn)行擬合確定出中心位置；不規(guī)則形狀的塔，可由各點(diǎn)的坐標(biāo)先畫出多邊形的各條邊，再擬合成多邊形，進(jìn)而求得中心位置。（2）圓心法：利用各點(diǎn)坐標(biāo)到中心的距離相等，找出中心位置。（3）平均值法：設(shè)，為平面上的N個(gè)離散數(shù)據(jù)點(diǎn)，可由公式計(jì)算出多邊形重心的坐標(biāo)，此法簡單直接，利用了離散數(shù)據(jù)點(diǎn)的坐標(biāo)求圖形重心，由于對數(shù)據(jù)點(diǎn)所圍圖形沒做處理，缺點(diǎn)是精度較差[3]。（4）利用求平面薄板重心公式求中心：它是將一個(gè)多邊形劃分為n個(gè)小區(qū)域，設(shè)分別為每一個(gè)小區(qū)域面積及其重心坐標(biāo)，利用重心坐標(biāo)公式，由積分得到累加模型:

圖4　多邊形分解Fig.4 Polygon decomposition

其中，若已知三角形三條邊的長度分別為a, b, c，由海倫公式

可以求出三角形的面積；結(jié)合使用Matlab編程，可求出每年各層的中心位置，所得結(jié)果如表1所示。

表1　各層次每年的中心坐標(biāo)Table 1 Ancient pagoda’s center coordinate of each floor in each year

3.2利用模型計(jì)算研究古塔的傾斜、彎曲、扭曲與位移

3.2.1古塔傾斜的計(jì)算模型常用的古塔傾斜觀測方法有：（1）塔頂定位控制點(diǎn)法；（2）底部幾何中心的定位控制點(diǎn)法；（3）傾斜觀測法：應(yīng)用水平角觀測方法，由塔底部與頂部的中心坐標(biāo)，計(jì)算得出塔頂部與底部中心的水平距離ΔD及塔的高度ΔH，可得出傾斜角的數(shù)學(xué)模型為[5]。

采用方法（3），計(jì)算古塔傾斜角度，其計(jì)算步驟如下：

表2　古塔傾斜角度(單位:弧度)Table 2Ancient pagoda’s tilt angle(Unit:radian)

由表2可知，1986-2011年間古塔幾乎沒有發(fā)生傾斜。

3.2.2古塔彎曲的計(jì)算模型在傾斜變形的基礎(chǔ)上彎曲變形是將直線變成曲線。

表3　古塔彎曲角度(單位:弧度)Table 3Ancient pagoda’s bending angle(Unit:radian)

圖5　1996年(-*)與2009年中心坐標(biāo)比較Fig.5 The comparison of pagoda’s center coordinate in 1996(-*)and 2009

圖6　螺旋扭曲示意圖Fig.6 The spirally twisted diagram

觀察表3及圖5得出，1986～1996年古塔彎曲情況很小，但2009年彎曲情況比較明顯，尤其是第2層至第6層有比較明顯的彎曲，第10層也有比較明顯的彎曲，第13層有復(fù)原趨勢，塔尖沒有繼續(xù)彎曲。因此，2009年較1986年古塔部分發(fā)生了彎曲，主要出現(xiàn)在1996年至2009年間，2009年至2011年幾乎沒有彎曲，但仍會有繼續(xù)彎曲趨勢。

3.2.3古塔扭曲的計(jì)算模型扭曲變形是在彎曲變形的基礎(chǔ)上旋轉(zhuǎn)，扭成螺旋狀態(tài)，如圖6。選取了四年古塔各層第一個(gè)監(jiān)測點(diǎn)的原始數(shù)據(jù)作為分析扭曲變形數(shù)據(jù)，設(shè)，，為古塔各層次第1個(gè)監(jiān)測點(diǎn)坐標(biāo)，i,j分別代表給出的年份和層數(shù)；求出向量；m也代表年份，計(jì)算步驟如下：

（1）先分別求出第i,m不同兩年的相鄰兩層第1個(gè)監(jiān)測點(diǎn)坐標(biāo)構(gòu)成的向量為：

分析表4數(shù)據(jù)可知，1986～2011年塔尖沒有發(fā)生扭曲變化；1986～1996年古塔發(fā)生了較小的扭曲，2009年則出現(xiàn)了明顯的扭曲，尤其是第5層產(chǎn)生比較明顯的扭曲，第13層發(fā)生了比較明顯的扭曲，塔尖沒有繼續(xù)扭曲，2011年較2009年古塔發(fā)生了很小的扭曲，2011年較1986年則發(fā)生了較大的扭曲，主要在1996～2009年間，而且會有繼續(xù)扭曲的趨向。

3.2.4古塔垂直位移檢測選擇沉降檢測的基準(zhǔn)點(diǎn)，由精密水準(zhǔn)的測量法，測量出古塔的負(fù)載沉降監(jiān)測點(diǎn)，依據(jù)高程變化量計(jì)算出古塔的沉降幅度[5]。因此提取出每年塔尖高度見表5。

由表5可知，塔頂高度從1986年到2011年間下降了3.46 cm，沉降幅度不大[6]。

3.3古塔的變形趨勢

運(yùn)用Excel，可求出結(jié)果見表6[7]。

表6　偏移量Table 6 Bias amount

由表6可知，1986～1996年古塔沿X軸正向（向東）偏移0.054 m，沿Y軸負(fù)向偏移0.0531 m；1996～2009年沿X軸正向偏移0.8292 m，沿Y軸負(fù)向偏移0.7843 m；2009～2011年沿X軸正向偏移0.0112 m，沿Y軸負(fù)向偏移0.0103 m；1986～2011年沿X軸正向偏移總量0.8942 m；沿Y軸負(fù)向偏移總量0.7277 m。

4　結(jié)語

構(gòu)建的模型及其求解過程簡練易懂，合理準(zhǔn)確地進(jìn)行問題研究。模型中涉及的各種數(shù)據(jù)的處理和計(jì)算都是通過Matlab與Excel軟件實(shí)現(xiàn)的，通過Matlab作出的實(shí)物模型圖，對古塔變形分析作用較大，以構(gòu)建的累加模型計(jì)算出的古塔各層次中心坐標(biāo)數(shù)據(jù)為基礎(chǔ)，傾斜角模型計(jì)算給出了古塔的傾斜情況；幾何向量夾角模型求出了古塔彎曲與扭曲的角度，比較精確地描繪出了古塔所發(fā)生的彎曲、扭曲的變形情況。對模型進(jìn)一步的討論可得到更加符合實(shí)際的信息，模型的推廣，可以為國家文物的保護(hù)提供參考。

[1]胡玄子,陳小雪,錢葉亮,等.數(shù)據(jù)處理中缺失數(shù)據(jù)填充方法的研究[J].湖北工業(yè)大學(xué)學(xué)報(bào),2013,28(5):82-84

[2]劉保柱,蘇彥華,張宏林.MATLAB7.0從入門到精通[M].北京:人民郵電出版社,2010:115

[3]李玉冰,郝永杰,劉恩海.多邊形重心的計(jì)算方法[J].計(jì)算機(jī)應(yīng)用,2005,12(25):391-393

[4]同濟(jì)大學(xué)數(shù)學(xué)系.高等數(shù)學(xué)第六版[M].北京:高等教育出版社,2009:170

[5]黃強(qiáng).古塔變形監(jiān)測的探討[J].測繪與空間地理信息,2013,6:217-220

[6]徐工,曲國慶.基于GM(1，1)模型的堤防沉降建模分析[J].山東理工大學(xué)學(xué)報(bào):自然科學(xué)版,2010,24(3):96-99

[7]郝惠馨,辛明影.基于Excel的數(shù)據(jù)處理方法[J].東北林業(yè)大學(xué)學(xué)報(bào),2010,38(10):131-132

Study on Ancient Pagoda’s Deformation Based on Model Construction

LIU Jing
Taishan Polytechnic College,Taian 271000,China

The deformation of ancient pagoda was investigated according to the observed data in 1986,1996,2009 and 2011 by using Matlab software and Excel.Firstly the approximate shapes of the ancient pagoda in 1986 and 2011 were drawn by using the Matlab,then the model construction and solution could be applied.The accumulation model was used to determine the center position of ancient pagoda’s each floor and the intersection angle model described a series of deformations such as tilting,bending and twisting.Finally,the deformation trend of the ancient pagoda was obtained.It has revealed a good result in the maintenance of ancient architectures.

Deformation;ancient pagoda;Matlab;Excel;model construction;intersection angle model

O29文獻(xiàn)標(biāo)示碼:A

1000-2324(2015)02-0265-05

2013-07-26

2013-08-22

劉靜(1973-),女,講師,碩士研究生,主要從事數(shù)學(xué)教學(xué)及其應(yīng)用研究.E-mail:Liujingteacher@163.com