唐紀芳,張子衛(wèi)
(四川職業(yè)技術(shù)學(xué)院,四川 遂寧 629000)
關(guān)于古塔變形研究的數(shù)學(xué)模型
唐紀芳,張子衛(wèi)
(四川職業(yè)技術(shù)學(xué)院,四川 遂寧 629000)
我國是一個文明古國,各地的古塔建筑非常的多,由于長時間承受自重、氣溫、風(fēng)力等各種作用,偶然還要受地震、颶風(fēng)的影響,古塔會產(chǎn)生各種變形,諸如傾斜、彎曲、扭曲等.為保護古塔,文物部門需適時對古塔進行觀測,了解各種變形量,以制定必要的保護措施.本文通過對古塔測量數(shù)據(jù)進行分析研究,對古塔各層的中心、變形的塔身傾斜度、扭曲度、古塔彎度等給出了明確的定義,并運用數(shù)學(xué)軟件進行了相應(yīng)計算,對古塔的各個測試時間的變形程度給出了具體量化結(jié)論,并建立了古塔變形的數(shù)學(xué)模型預(yù)測其變形趨勢,給文物管理部門提出了有關(guān)古塔保護的建議.
古塔變形;各層中心點;傾斜度;彎曲度;扭曲度;變形預(yù)測
由于長時間承受自重、氣溫、風(fēng)力等各種作用,偶然還要受地震、颶風(fēng)的影響,古塔會產(chǎn)生各種變形,諸如傾斜、彎曲、扭曲等。為保護古塔,文物部門需適時對古塔進行觀測,了解各種變形量,以制定必要的保護措施.
某古塔已有上千年歷史,是我國重點保護文物.管理部門委托測繪公司先后于1986年7月、1996年8月、2009年3月和2011年3月對該塔進行了4次觀測.
請你們根據(jù)附件1提供的4次觀測數(shù)據(jù),討論以下問題:
(1)給出確定古塔各層中心位置的通用方法,并列表給出各次測量的古塔各層中心坐標(biāo).
(2)分析該塔傾斜、彎曲、扭曲等變形情況.
(3)分析該塔的變形趨勢.
(1)假設(shè)附件1中的數(shù)據(jù)是以同樣的方式對同樣的8個點進行了4次觀測.
(2)假設(shè)古塔最初建筑時是正八邊形,而它的中心就是重心.
(3)假設(shè)塔尖就是一個點.
(4)假設(shè)底層不考慮彎曲和扭曲程度.
(5)假設(shè)各層的傾斜、彎曲和扭曲均是以第1層為基礎(chǔ)作比較.
(6)假設(shè)古塔是公元1000年修建.
(7)假設(shè)古塔的變形不受強烈地質(zhì)災(zāi)害影響.
3.1 問題1的求解
(1)分析.根據(jù)附件提供的數(shù)據(jù)信息,發(fā)現(xiàn)所測數(shù)據(jù)為古塔各層底部角點處的空間坐標(biāo)數(shù)據(jù),經(jīng)過試算,1至13層為邊長相等的空間八邊形,因此我們認為,古塔修建時,1至13層為水平的平面正八邊形,正八邊形的中心就是各頂點的重心,由于時間久遠,長時間承受自重、氣溫、風(fēng)力等各種作用,偶然還要受地震、颶風(fēng)的影響,從而產(chǎn)生變形,各層底面已經(jīng)不再是正多邊形,頂點已經(jīng)變化位置,但中心應(yīng)該在各層底面各頂點的重心附近,因此我們先根據(jù)《圖論》中的知識得出空間有限個點的中心的定義,再運用matlab軟件編程在各層底部各頂點的重心附近搜尋中心位置.
(2)空間個n點Ai(i=1,2,…,n)的中心O的定義.若空間一點O滿足min﹛ max﹛ OAi﹜﹜(i=1,2,…,n),則O叫做空間n個點Ai(i=1,2,…,n)的中心.即對空間上的任一點到這n個點的距離中有一個最大值,則所有可能的最大值中那個最小的值所對應(yīng)的那個點我們就把它叫做這n個點Ai(i=1,2,…,n)的中心.
(3)各層底面中心的確定.在研究各層的中心之前,我們首先對缺失的數(shù)據(jù)進行了補充.在1986年和1996年的第十三層的第5個觀測點的數(shù)據(jù)都缺失了.我們將2009年和2011年該層的同一個觀測點的數(shù)據(jù)差值除以年份差作為每年的變化差值,然后計算出1986年和1996年的補充數(shù)據(jù).然后,各層底面中心應(yīng)該在各層8個觀測點的重心附近,所以我們先求出其重心坐標(biāo),以其為中心構(gòu)造長方體型區(qū)域,在這個區(qū)域內(nèi)編程搜索滿足min﹛ max﹛ OAi﹜﹜(i=1,2,…,n)的點O,其程序見附錄1,計算結(jié)果如下(表1):
表1 各年所測古塔各層中心坐標(biāo)
我們運用matlab軟件的三維作圖工具作出古塔在各年所測的頂點與中心的圖形,見圖一.
圖一
3.2 問題2的求解
(1)古塔傾斜角的計算.我們認為,古塔的傾斜分為塔基的傾斜與第i層的傾斜,而傾斜度的刻畫分為與鉛直方向的偏離角度和偏向水平面的具體方位,這個方位我們用塔身向量在平面上的投影與軸的正方向的夾角來確定.底層傾斜角的計算我們采用將底層各頂點擬合的平面的法向量與鉛直方向的偏離角度和偏向水平面的具體方位來確定,第i層的傾斜程度我們用第1層的中心為起點,第i層的中心為終點的向量與鉛直方向的偏離角度和偏向水平面的具體方位來刻畫.(見下圖)
1)第一層底面(塔基)的傾斜度
我們先求出各頂點的擬合平面,設(shè)擬合平面方程為Axil+Byil+Czil=1(C>0),代入八個頂點坐標(biāo)得到八個關(guān)于A,B,C的一次方程,然后用matlab軟件解超定方程的矩陣左除法得到擬合平面的法向量(程序見附錄3),m=﹛A,B,﹜C得到它與鉛直方向向量{0,0,1}的夾角度數(shù)
在水平面上偏離ox軸的正方向的夾角為
2)古塔第1層至第j層的傾斜度的計算
表2 各年豎直傾斜度α
表3 各年水平偏轉(zhuǎn)度θ
我們?nèi)「鲗拥钠骄鶅A斜度為古塔的傾斜度,結(jié)果見上表.
(2)古塔扭曲度的計算.我們在每層塔的測量數(shù)據(jù)中選定第一個觀測點與該層的中心點連線成一個空間向量,然后將古塔第j-1層與第j層兩層的這兩個向量的夾角作為古塔第j-1層的局部扭曲度,運用表2數(shù)據(jù)我們得到如下結(jié)果,見表4.
表4 扭曲度
1986年最大扭曲度在第十層,扭曲度為1.467943度,1996年最大扭曲度在第十層,扭曲度為1.445559度,2009年最大扭曲度在第五層,扭曲度為0.917536度,2011年最大扭曲在第五層,扭曲度為0.93176度.
(3)古塔彎曲度的計算.古塔在第j層的彎曲度,我們認為是由Oj-1,Oj,Oj+1這三個點確定。這三個點可確定一個圓,把這個圓的半徑的倒數(shù)作為古塔在第j層的彎曲度.
接下來我們計算彎曲度,我們以第1,2,3層為例:
表5 各年彎曲度對比表
1986年、1996年、2009年、2011年最大扭曲分別發(fā)生在第11層、第10層、第9層、第8層,彎曲度為分別為:
3.3 問題3的求解
古塔變形趨勢的預(yù)測我們假設(shè)古塔建于公元1000年,其傾斜度,扭曲度、彎曲度均為0,我們以時間t為變量,建立第t年古塔最大傾斜度,最大扭曲度、最大彎曲度的函數(shù)f(t),g(t),h(t).
首先,取古塔每年最大豎直及水平偏轉(zhuǎn)度,建立了表6,再用這5個點進行回歸,,通過回歸方法求出傾斜函數(shù)f(t),然后根據(jù)此函數(shù)計算,得出2011——2026年預(yù)測數(shù)據(jù)見表7,作出相應(yīng)的圖形,見圖二、圖三.
表6
傾斜函數(shù):
表7
圖二
圖三
取古塔每年最大扭曲度建立了表8,利用這5個點進行回歸分析,通過回歸方法求出扭曲函數(shù),即最大扭曲度函數(shù):,然后根據(jù)函數(shù)計算,得出2012——2026年預(yù)測數(shù)據(jù)見表9,作出相應(yīng)的圖形,見圖四.
表8
表9
圖四
表10
表11
圖五
本文在實際問題的基礎(chǔ)上,通過合理的假設(shè),預(yù)測出古塔的最大可能的傾倒方向,可以協(xié)助有關(guān)部門對古跡的維護.
充分利用excle來對大量的數(shù)據(jù)進行計算,完成表格的輸入、統(tǒng)計、分析等多項工作,可生成精美直觀的表格、圖表.
Matlab是一種科學(xué)計算軟件,專門以矩陣的形式處理數(shù)據(jù),利用matlab計算可以使我們所計算的數(shù)據(jù)更加精確.
[1]陳國華,韋程東.數(shù)學(xué)建模與數(shù)學(xué)建模方法[M].天津:南開大學(xué)出版社,2012.
[2]陳光亭,裘哲勇.數(shù)學(xué)建模[M].北京:高等教育出版社,2012.
[3]梁???古塔變形測量方法探討[J].城市勘測,2011,(3).
[4]吳望名.圖論及其應(yīng)用[M].北京:科學(xué)出版社,1984.
責(zé)任編輯:張隆輝
O 29
A
1672-2094(2013)06-0155-05
2013-10-28
唐紀芳(1982-),女,四川遂寧人,四川職業(yè)技術(shù)學(xué)院應(yīng)用數(shù)學(xué)與經(jīng)濟系講師,碩士。
張子衛(wèi)(1963-),男,四川遂寧人,四川職業(yè)技術(shù)學(xué)院應(yīng)用數(shù)學(xué)與經(jīng)濟系副教授。