古塔变形的数学模型大学生数学建模竞赛C题全国二等奖.docx
《古塔变形的数学模型大学生数学建模竞赛C题全国二等奖.docx》由会员分享，可在线阅读，更多相关《古塔变形的数学模型大学生数学建模竞赛C题全国二等奖.docx（80页珍藏版）》请在冰豆网上搜索。
古塔变形的数学模型大学生数学建模竞赛C题全国二等奖
2013高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）：
C
我们的参赛报名号为（如果赛区设置报名号的话）：
13289002
所属学校（请填写完整的全名）：
湖北职业技术学院
参赛队员（打印并签名）：
1.罗周希
2.金超
3.李婉丽
指导教师或指导教师组负责人（打印并签名）：
数模组
日期：
2013年9月13日
赛区评阅编号（由赛区组委会评阅前进行编号）：
2013高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号（由赛区组委会评阅前进行编号）：
赛区评阅记录（可供赛区评阅时使用）：
评
阅
人
评
分
备
注
全国统一编号（由赛区组委会送交全国前编号）：
全国评阅编号（由全国组委会评阅前进行编号）：
古塔变形的数学模型
摘要
古塔被誉为中国古代杰出的高层建筑物，历史悠久，值得并需要我们的保护。
本文研究了关于古塔变形的问题，古塔的变形与塔身的中心紧密相关，具体分析了古塔倾斜、弯曲、扭曲的变形情况及趋势。
对于问题1，建立中心位置模型，采用多边形组合形心的算法，求的结果是表6，7，8，9中的数据；
在问题2-1中，研究塔身的倾斜建立了古塔自身倾斜角的模型和古塔相对倾斜角模型两个数学模型，模型2-1.1利用三角函数相关知识确定倾斜角，结果是1986年的塔身倾斜了1.5308°，1996年的塔身倾斜了1.5558°，2009年的塔身倾斜了1.5564°，2011年的塔身倾斜了1.5339°；模型2-1.2对各年份各层的中心点数据进行空间直线拟合，采用空间向量法计算两直线的夹角，结果是古塔1996年相对于1986年倾斜了0.049°，2009年相对于1996年未发生倾斜，2011年相对于2009年倾斜了0.0245°。
在问题2-2中，研究塔身的弯曲建立古塔弯曲的数学模型，采用三次多项式拟合的算法，得到了古塔每年每层的弯曲率。
在问题2-3中，研究塔身的扭曲建立古塔扭曲的数学模型，采用两点之间的距离公式和扭曲加权法，得出1996年总扭曲距离11.9399,2009年总扭曲距离8.4095，2011年总扭曲距离为3.8508。
在问题3中，采用层次分析法建立古塔变形趋势的数学模型，对古塔的变形程度和趋势做出综合量化评价，分析得到塔楼变形的趋势为从1986年到1996年变形程度加深，所以塔楼变形的趋势为从1996年到2009年变形程度减轻，从2009年到2011年变形程度稍有加深。
关键词：
组合图形的形心拟合中心变形曲线层次分析法扭曲加权
一．问题重述
古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响会产生各种变形，诸如倾斜、弯曲、扭曲等。
为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。
现有一古塔已上千年历史，是我国重点保护文物。
管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。
请根据题目提供的4次观测数据，讨论以下问题：
1.给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。
2.分析该塔倾斜、弯曲、扭曲等变形情况。
3.分析该塔的变形趋势。
二、问题分析
塔的变形趋势与塔身的中心有关，因此本题的主要目标是计算出塔中心的变化规律。
研究题中数据发现x轴坐标值逐渐偏大，y轴坐标值逐渐偏小，z轴坐标值逐渐偏小，可以由此联想到塔正向着x轴偏大的方向、y轴偏小的方向倾斜，且下沉。
说明x轴与y轴决定倾斜方向；z轴和y轴数值变化越大，倾斜程度越大，反之则越小。
因此在求解层面中心点的z坐标点取八个点z值的平均。
问题1属于计算多边形组合形心的数学问题。
因为1986年受自然、社会影响度较小，最接近古塔的原始图样，所以图形较为有参考价值，用MATLAB画出1986年塔身的大致图样为八边形（图1-1），
图1-1
通过假设得出正多边形的形心就是中心，因此可运用组合图形形心公式计算出古塔各层中心坐标。
问题2-0塔的倾斜、弯曲、扭曲都与塔的中心有关。
从数据点上看，四年13组中心点的数据大致分布在空间的一个平面上（1986年图2-0-1，1996年图2-0-2，2009年图2-0-3，2011年图2-0-4），所以用空间直线的最小二乘法来拟合中心线。
图2-0-1图2-0-2
图2-0-3图2-0-4
问题2-1塔的倾斜程度用倾斜角W表示，塔身拟合后的中心线与地面的夹角为倾斜角（图2-1-1）。
图2-1-1
问题2-1.1塔自身倾斜情况。
观察图2-1-1得到Rt△OPQ，利用三角形正切公式确定自身倾斜角W。
1986、1996年题目分别给有4组塔尖坐标，分析数据结合实际生活我们将4组数据求平均值，得到1组数据用于计算塔的倾斜角W。
问题2-1.2年份的相对倾斜角△W可表现塔倾斜幅度。
两直线的方向向量的夹角（小于90O）叫做直线的夹角，现古塔的倾斜度一定小于90O，所以各年份中心线夹角可转化为两直线的方向向量的夹角，将每年各层中心和塔尖的数据用空间直线的最小二乘法进行拟合成中心线。
用法向量法得到各中心直线的方向向量，利用两向量的夹角的余弦公式确定相对倾斜角。
问题2-2塔的弯曲与空间直角坐标系的z轴有关。
以塔第一层楼面的中心为原点建立空间直角坐标系，拟合各年各层中心点竖坐标的直线，多次拟合后发现三次多项式的拟合效果较好（附录图2-2-1，图2-2-2，图2-2-3），
问题2-3塔的扭曲在空间直角坐标系中的x,y轴有关。
扭曲是相邻年份相对的物理量，以每层中心点偏移的平面距离作为扭曲距离，并用层次分析法对各层的扭曲距离赋予权重，得到总扭曲距离来衡量不同年份古塔的扭曲程度
问题3运用层次分析法对古塔的变形程度做综合量化评价，并根据组合权向量对古塔变形趋势做出大致分析。
三、模型假设
1，参照中国古塔的塔身样式假设塔身水平面是正八角形；
2，假设古塔是均质物体，形心就是中心；
3，以地平面为X轴，建立空间直角坐标系X，Y，Z；
4，假设塔尖面积极小，近视为一个点；
5，假设模型4中古塔所在地不存在地震等这种偶然事件；
6，假设古塔所在的地平面是水平面。
四．符号说明
Ci--第i个三角形的形心（i=1,2,……6）；
Si--第i个三角形的面积（单位：
m2）（i=1,2,……6）；
Ai--第i个三角形的形心坐标（单位：
m）；
Wt--第t年塔身的倾斜角（t=1986,1996,2009,2011）；
Dt--第t年塔身在地面上的水平投影（t=1986,1996,2009,2011）；
Ht--第t年塔上部与下部的相对位移（t=1986,1996,2009,2011）；
Lt--第t年塔的拟合直线（t=1986,1996,2009,2011）；
R--塔的倾斜率；
mt--第t年中心直线x轴的法向量（t=1986,1996,2009,2011）；
nt--第t年中心直线y轴的法向量（t=1986,1996,2009,2011）；
Jt--第t年的中心直线方向向量（t=1986,1996,2009,2011）；
K--塔的曲率；
五．模型的建立与求解
5.1、模型的准备
（1）给数据用图像形式表现出来题目所给的4次观测数据1986年、1996年塔13层的第五个点数据残缺。
参照其10层、11层、12层的数据知道第五个点的数值在第四个点和第六个点的数值之间，另外分析后发现这两组数据对题目结果的影响不大。
因此结合实际和理论，采取取相邻年份该层的第四个点和第六个点的平均值将两组数据补齐。
（2）用MATLAB软件将题中所给数据及补充的数据画图表现出来。
5.2问题1中心位置模型
问题一求古塔各层中心坐标，通过假设知道实质就是计算正多边形的中心，运用组合图形求形心的方法。
将每层正八边形的塔面都分为六个三角形（图1-2）S1、S2、S3、S4、S5、S6，
图5-1-2
每年各层六个三角形的形心：
C1（
，
）、
C2（
）、
C3（
）、
C4（
）、
C5（
）、
C6（
），
代入测量数据，运行MATLAB软件得出：
1986年各层分面积的形心
塔层
坐标（m）
A1
A2
A3
A4
A5
A6
1
565.6707，526.9707
564.3160，524.7823
565.1133，523.6693
567.8320，525.6177
570.8977，523.7220
570.3780，520.4057
2
565.7107，526.7563
564.4383，524.6417
565.2177，523.5770
567.8177，525.4710
570.7933，523.6700
570.3203，520.4773
3
565.7503，526.5460
564.5587，524.5040
565.3203，523.4870
567.8037，525.3273
567.6070，524.0603
567.1797，520.9890
4
565.7813，526.3787
564.6540，524.3940
565.4017，523.4147
567.7920，525.2127
570.6087，523.5783
570.2180，520.6037
5
565.8150，526.1987
564.7567，524.2760
565.4893，523.3370
567.7800，525.0893
570.5210，523.5347
570.1693，520.6640
6
565.8483，526.0230
564.8583，524.1610
565.5757，523.2620
567.2693，524.9697
570.4320，523.4920
570.1233，520.7217
7
565.9323，525.8013
565.0030，524.0360
565.6087，522.1230
568.3457，522.8950
568.9073，520.9863
567.9333，519.2110
8
566.0143，525.5857
565.1437，523.9140
565.7230，522.1233
568.2910，522.8597
568.8363，521.0700
567.9370，519.4030
9
566.0963，525.3700
565.2862，523.7907
565.8408，522.1207
568.2398，522.8213
568.7700，521.1480
567.9455，519.7718
10
566.1781，525.1621
565.4221，523.6731
565.9518，522.1205
568.1901，522.7871
568.7034，521.2268
567.9501，519.7718
11
566.2173，525.1621
565.3411，522.7775
565.9963，522.1686
568.1775，522.8075
568.6384，521.1750
567.9245，519.9055
12
566.4234，524.7179
565.7663，523.4303
566.2073，522.0879
568.1537，522.6799
568.5723，521.3238
567.8977，520.0395
13
566.5415，524.5037
565.9325，523.3134
566.3304，522.0724
568.1358，522.6281
568.5088，521.3708
568.5314，521.3482
表1
1996年各层分面积的形心
塔层
坐标（m）
A1
A2
A3
A4
A5
A6
1
565.6709，526.9704
564.3162，524.7820
564.8909，522.3034
568.4069，523.1387
568.9919，520.6390
567.6595，518.4570
2
565.7116，526.7554
564.4392，524.6499
564.9537，522.2315
568.3784，523.0668
567.0461，520.8848
567.7212，518.5824
3
565.7519，526.5444
564.5603，524.5024
565.1606，522.2207
568.4056，523.0441
568.9876，520.7551
567.7816，518.2057
4
565.7835，526.3765
564.6562，524.3918
565.2669，522.1878
568.4048，523.0065
568.9855，520.8005
567.8299，518.8035
5
565.8178，526.1958
564.7594，524.2731
567.2741，521.0509
566.4548，524.1346
568.5384，524.6946
570.4308，523.5249
6
563.6743，521.6294
565.4214，520.6144
567.3460，521.1294
566.5451，523.9446
568.4659，524.4789
570.2377，523.4745
7
563.8114，521.6232
565.5112，520.6525
565.6125，522.1191
568.3495，522.8911
568.9111，520.9825
567.9371，519.2071
8
566.0185，525.5814
565.1478，523.9097
565.7272，522.1190
568.2953，522.8554
568.8405，521.0657
567.9412，519.3987
9
564.2483，521.6049
565.7973，520.7773
567.4596，521.2459
566.6890，523.6426
568.2713，522.7973
567.9455，519.5893
10
564.4428，521.5880
565.9130，520.7977
567.4925，521.2496
566.7358，523.5298
568.327，524.0022
569.2540，523.2009
11
564.6655，521.5956
566.0338，520.8531
567.5067，521.2830
566.6366，523.4134
568.3226，523.8619
569.7207，523.1170
12
564.8876，521.6031
566.1544，520.9094
567.5207，521.3166
566.9068，523.2978
568.2923，523.7225
569.5873，523.0337
13
565.1091，521.6032
566.2779，520.9556
567.5416，521.3413
566.9817，523.1791
568.2704，523.5812
569.4667，522.9465
表2
2009年各层分面积的形心
塔层
坐标（m）
A1
A2
A3
A4
A5
A6
1
562.4995，521.6767
564.6728，520.3384
567.1635，520.9321
566.2995，524.4607
568.7379，525.0660
570.9705，523.7300
2
562.7235，521.6673
564.8138，520.3847
567.1945，520.9648
569.5556，521.5574
569.5059，523.8796
570.8206，523.6727
3
562.9424，521.6587
564.9515，520.4308
567.2241，520.9976
566.3847，524.2815
568.6415，524.8619
570.6734，523.6173
4
563.1175，521.6513
565.0616，520.4670
567.2484，521.0232
566.4188，524.2100
568.5920，524.7807
570.5565，523.5727
5
563.3050，521.6439
565.1797，520.5062
567.2741，521.0509
566.4548，524.1346
568.5384，524.6946
570.4308，523.5249
6
563.6743，521.6294
656.4214，520.6144
567.3460，521.1294
566.5451，523.9446
568.4659，524.4789
570.2377，523.4345
7
563.8114，521.6232
565.5112，520.6525
567.3735，521.1561
566.5799，523.8915
568.4390，524.3948
570.1656，523.3980
8
564.0303，521.6142
565.6546，520.7151
567.4165，521.2012
566.6344，523.7572
568.3961，524.2642
570.0504，523.3433
9
564.2483，521.6049
565.7973，520.7773
567.4596，521.2459
566.6890，523.6426
568.3533，524.1329
569.9357，523.2876
10
564.4428，521.5880
565.9130，520.7977
567.4925，521.2496
566.7658，523.5298
568.3527，524.0022
569.8540，523.2007
11
564.6655，521.5955
566.0338，520.8535
567.5067，521.2830
566.8366，523.4134
568.3226，523.8619
569.7207，523.1170
12
564.8876，521.6031
566.1544，520.9094
567.5207，520.3166
566.9068，523.2978
568.2923，523.7225
569.5873，523.0337
13
565.1091，521.6032
566.2779，520.9556
567.5416，521.3413
566.9817，523.1791
568.2704，523.5812
569.4667，522.9465
表3
2011年各层分面积的形心
塔层
坐标（m）
A1
A2
A3
A4
A5
A6
1
562.49957，521.6766
564.6856，5204421
564.6730，520.3383
567.1637，520.9320
567.1413，520.9436
568.4800，523.1363
2
562.7237，521.6671
564.8140，520.3844
567.1947，520.9645
566.3432，524.3695
568.7042，524.9621
570.8208，523.6725
3
562.9428，521.6591
564.9518，520.4312
567.2245，520.9979
566.3851，524.2819
568.6419，524.8622
570.6737，523.6176
4
563.1179，521.6509
565.0620，520.4666
567.2488，521.0227
566.4192，524.2096
568.5924，524.7802
570.5569，523.5722
5
563.3056，521.6433
565.1803，520.5056
567.2747，521.0503
566.4554，524.1339
568.5391，524.6939
570.4314，523.5243
6
563.6750，521.6286
565.4221，520.6136
567.3467，521.1287
566.5458，523.9439
568.4666，524.4782
570.2384，523.4338
7
563.8122，521.6224
565.5120，520.6517
567.3743，521.1553
566.5807，523.8707
568.4398，524.3940
570.1664，523.3972
8
564.0311，521.6133
565.6554，520.7142
567.4173，521.2003
566.6352，523.7563
568.3969，524.2633
570.0512，523.3424
9
564.2493，521.6040
564.7983，520.7764
567.4607，521.3884
566.6899，523.1730
568.3543，523.6633
569.9366，522.8180
10
564.4439，521.5
升级会员 