学士学位论文古塔变形的数学模型数学建模竞赛c题全国二等奖论文Word文档格式.docx
- 文档编号:18386949
- 上传时间:2022-12-15
- 格式:DOCX
- 页数:69
- 大小:479.65KB
学士学位论文古塔变形的数学模型数学建模竞赛c题全国二等奖论文Word文档格式.docx
《学士学位论文古塔变形的数学模型数学建模竞赛c题全国二等奖论文Word文档格式.docx》由会员分享,可在线阅读,更多相关《学士学位论文古塔变形的数学模型数学建模竞赛c题全国二等奖论文Word文档格式.docx(69页珍藏版)》请在冰豆网上搜索。
,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运用层次分析法对古塔的变形程度做综合量化评价,并根据组合权向量对古塔变形趋势做出大致分析。
3、模型假设
1,参照中国古塔的塔身样式假设塔身水平面是正八角形;
2,假设古塔是均质物体,形心就是中心;
3,以地平面为X轴,建立空间直角坐标系X,Y,Z;
4,假设塔尖面积极小,近视为一个点;
5,假设模型4中古塔所在地不存在地震等这种偶然事件;
6,假设古塔所在的地平面是水平面。
4.符号说明
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.模型的建立与求解
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年各层分面积的形心
565.6709,526.9704
564.3162,524.7820
564.8909,522.3034
568.4069,523.1387
568.9919,520.6390
567.6595,518.4570
565.7116,526.7554
564.4392,524.6499
564.9537,522.2315
568.3784,523.0668
567.0461,520.8848
567.7212,518.5824
565.7519,526.5444
564.5603,524.5024
565.1606,522.2207
568.4056,523.0441
568.9876,520.7551
567.7816,518.2057
565.7835,526.3765
564.6562,524.3918
565.2669,522.1878
568.4048,523.0065
568.9855,520.8005
567.8299,518.8035
565.8178,526.1958
564.7594,524.2731
567.2741,521.0509
566.4548,524.1346
568.5384,524.6946
570.4308,523.5249
563.6743,521.6294
565.4214,520.6144
567.3460,521.1294
566.5451,523.9446
568.4659,524.4789
570.2377,523.4745
563.8114,521.6232
565.5112,520.6525
565.6125,522.1191
568.3495,522.8911
568.9111,520.9825
567.9371,519.2071
566.0185,525.5814
565.1478,523.9097
565.7272,522.1190
568.2953,522.8554
568.8405,521.0657
567.9412,519.3987
564.2483,521.6049
565.7973,520.7773
567.4596,521.2459
566.6890,523.6426
568.2713,522.7973
567.9455,519.5893
564.4428,521.5880
565.9130,520.7977
567.4925,521.2496
566.7358,523.5298
568.327,524.0022
569.2540,523.2009
564.6655,521.5956
566.0338,520.8531
567.5067,521.2830
566.6366,523.4134
568.3226,523.8619
569.7207,523.1170
564.8876,521.6031
566.1544,520.9094
567.5207,521.3166
566.9068,523.2978
568.2923,523.7225
569.5873,523.0337
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年各层分面积的形心
562.4995,521.6767
564.6728,520.3384
567.1635,520.9321
566.2995,524.4607
568.7379,525.0660
570.9705,523.7300
562.7235,521.6673
564.8138,520.3847
567.1945,520.9648
569.5556,521.5574
569.5059,523.8796
570.8206,523.6727
562.9424,521.6587
564.9515,520.4308
567.2241,520.9976
566.3847,524.2815
568.6415,524.8619
570.6734,523.6173
563.1175,521.6513
565.0616,520.4670
567.2484,521.0232
566.4188,524.2100
568.5920,524.7807
570.5565,523.5727
563.3050,521.6439
565.1797,520.5062
656.4214,520.6144
570.2377,523.4345
567.3735,521.1561
566.5799,523.8915
568.4390,524.3948
570.1656,523.3980
564.0303,521.6142
565.6546,520.7151
567.4165,521.2012
566.6344,523.7572
568.3961,524.2642
570.0504,523.3433
568.3533,524.1329
569.9357,523.2876
566.7658,523.5298
568.3527,524.0022
569.8540,523.2007
564.6655,521.5955
566.0338,520.8535
566.8366,523.4134
567.5207,520.3166
表3
2011年各层分面积的形心
562.49957,521.6766
564.6856,5204421
564.6730,520.3383
567.1637,520.9320
567.1413,520.9436
568.4800,523.1363
562.7237,521.6671
564.8140,520.3844
567.1947,520.9645
566.3432,524.3695
568.7042,524.9621
570.8208,523.6725
562.9428,521.6591
564.9518,520.4312
567.2245,520.9979
566.3851,524.2819
568.6419,524.8622
570.6737,523.6176
563.1179,521.6509
565.0620,520.4666
567.2488,521.0227
566.4192,524.2096
568.5924,524.7802
570.5569,523.5722
563.3056,521.6433
565.1803,520.5056
567.2747,521.0503
566.4554,524.1339
568.5391,524.6939
570.4314,523.5243
563.6750,521.6286
565.4221,520.6136
567.3467,521.1287
566.5458,523.9439
568.4666,524.4782
570.2384,523.4338
563.8122,521.6224
565.5120,520.6517
567.3743,521.1553
566.5807,523.8707
568.4398,524.3940
570.1664,523.3972
564.0311,521.6133
565.6554,520.7142
567.4173,521.2003
566.6352,523.7563
568.3969,524.2633
570.0512,523.3424
564.2493,521.6040
564.7983,520.7764
567.4607,521.3884
566.6899,523.1730
568.3543,523.6633
569.9366,522.8180
564.4439
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 学士学位 论文 变形 数学模型 数学 建模 竞赛 全国 二等奖