《古塔变形的数学模型大学生数学建模竞赛C题全国二等奖8669874.doc》由会员分享，可在线阅读，更多相关《古塔变形的数学模型大学生数学建模竞赛C题全国二等奖8669874.doc（56页珍藏版）》请在三一文库上搜索。
1、高教社杯全国大学生数学建模竞赛古塔变形的数学模型承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛
2、章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C 我们的参赛报名号为（如果赛区设置报名号的话）： 13289002 所属学校（请填写完整的全名）： 湖北职业技术学院 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 数模组 日期： 2013 年 9 月13 日赛区评阅编号（由赛区组委会评阅前进行编号）：高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编
3、号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：古塔变形的数学模型摘要 古塔被誉为中国古代杰出的高层建筑物，历史悠久，值得并需要我们的保护。本文研究了关于古塔变形的问题，古塔的变形与塔身的中心紧密相关，具体分析了古塔倾斜、弯曲、扭曲的变形情况及趋势。对于问题1，建立中心位置模型，采用多边形组合形心的算法，求的结果是表6，7，8，9中的数据；在问题2-1中，研究塔身的倾斜建立了古塔自身倾斜角的模型和古塔相对倾斜角模型两个数学模型，模型2-1.1利用三角函数相关知识确定倾斜角
4、，结果是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年总扭曲距离
5、11.9399,2009年总扭曲距离8.4095 ，2011年总扭曲距离为3.8508。在问题3中，采用层次分析法建立古塔变形趋势的数学模型，对古塔的变形程度和趋势做出综合量化评价，分析得到塔楼变形的趋势为从1986年到1996年变形程度加深，所以塔楼变形的趋势为从1996年到2009年变形程度减轻，从2009年到2011年变形程度稍有加深。 关键词： 组合图形的形心 拟合 中心变形曲线 层次分析法 扭曲加权一 问题重述古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定
6、必要的保护措施。现有一古塔已上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。请根据题目提供的4次观测数据，讨论以下问题：1. 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。2. 分析该塔倾斜、弯曲、扭曲等变形情况。3. 分析该塔的变形趋势。二、问题分析塔的变形趋势与塔身的中心有关，因此本题的主要目标是计算出塔中心的变化规律。研究题中数据发现x轴坐标值逐渐偏大，y轴坐标值逐渐偏小，z轴坐标值逐渐偏小，可以由此联想到塔正向着x轴偏大的方向、y轴偏小的方向倾斜，且下沉。说明x轴与
7、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-
8、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得到RtOPQ，利用三角形正切公式确定自身倾斜角W。1986、1996年题目分别给有4组塔尖坐标，分析数据结合实际生活我们将4组数据求平均值，得到1组数据用于计算塔的倾斜角W。问题2-1.2年份的相对倾斜角W可表现塔倾斜幅度。两直线的方向向量的夹角（小于90O）叫做直线的夹角，现古塔的倾斜度一定小于90O，所以各年份中心线夹角可转化为两
9、直线的方向向量的夹角，将每年各层中心和塔尖的数据用空间直线的最小二乘法进行拟合成中心线。用法向量法得到各中心直线的方向向量，利用两向量的夹角的余弦公式确定相对倾斜角。 问题2-2塔的弯曲与空间直角坐标系的z轴有关。以塔第一层楼面的中心为原点建立空间直角坐标系，拟合各年各层中心点竖坐标的直线，多次拟合后发现三次多项式的拟合效果较好（附录图2-2-1，图2-2-2，图2-2-3），问题2-3塔的扭曲在空间直角坐标系中的x,y轴有关。扭曲是相邻年份相对的物理量，以每层中心点偏移的平面距离作为扭曲距离，并用层次分析法对各层的扭曲距离赋予权重，得到总扭曲距离来衡量不同年份古塔的扭曲程度问题3运用层次分析
10、法对古塔的变形程度做综合量化评价，并根据组合权向量对古塔变形趋势做出大致分析。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年塔身
11、在地面上的水平投影(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
12、年塔13层的第五个点数据残缺。参照其10层、11层、12层的数据知道第五个点的数值在第四个点和第六个点的数值之间，另外分析后发现这两组数据对题目结果的影响不大。因此结合实际和理论，采取取相邻年份该层的第四个点和第六个点的平均值将两组数据补齐。（2） 用MATLAB软件将题中所给数据及补充的数据画图表现出来。5.2 问题1 中心位置模型问题一求古塔各层中心坐标，通过假设知道实质就是计算正多边形的中心，运用组合图形求形心的方法。将每层正八边形的塔面都分为六个三角形（图1-2）S1、S2、S3、S4、S5、S6，图5-1-2每年各层六个三角形的形心：C1(，)、C2(,)、C3(,)、C4(,)、C
13、5(,)、C6(,)， 代入测量数据，运行MATLAB软件得出：1986年各层分面积的形心 塔层 坐标(m)A1A2A3A4A5A61565.6707，526.9707 564.3160，524.7823 565.1133，523.6693567.8320，525.6177 570.8977，523.7220570.3780，520.40572565.7107，526.7563564.4383，524.6417565.2177，523.5770567.8177，525.4710570.7933，523.6700570.3203，520.47733565.7503，526.5460564.558
14、7，524.5040565.3203，523.4870567.8037，525.3273567.6070，524.0603567.1797，520.98904565.7813，526.3787564.6540，524.3940565.4017，523.4147567.7920，525.2127570.6087，523.5783570.2180，520.60375565.8150，526.1987564.7567，524.2760565.4893，523.3370567.7800，525.0893570.5210，523.5347570.1693，520.66406565.8483，526.02
15、30564.8583，524.1610565.5757，523.2620567.2693，524.9697570.4320，523.4920570.1233，520.72177565.9323，525.8013565.0030，524.0360565.6087，522.1230568.3457，522.8950568.9073，520.9863567.9333，519.21108566.0143，525.5857565.1437，523.9140565.7230，522.1233568.2910，522.8597568.8363，521.0700567.9370，519.40309566.09
16、63，525.3700565.2862，523.7907565.8408，522.1207568.2398，522.8213568.7700，521.1480567.9455，519.771810566.1781，525.1621565.4221，523.6731565.9518，522.1205568.1901，522.7871568.7034，521.2268567.9501，519.771811566.2173，525.1621565.3411，522.7775565.9963，522.1686568.1775，522.8075568.6384，521.1750567.9245，519.
17、905512566.4234，524.7179565.7663，523.4303566.2073，522.0879568.1537，522.6799568.5723，521.3238567.8977，520.039513566.5415，524.5037565.9325，523.3134566.3304，522.0724568.1358，522.6281568.5088，521.3708568.5314，521.3482表11996年各层分面积的形心 塔层 坐标(m)A1A2A3A4A5A61565.6709，526.9704564.3162，524.7820564.8909，522.3034
18、568.4069，523.1387568.9919，520.6390567.6595，518.45702565.7116，526.7554564.4392，524.6499564.9537，522.2315568.3784，523.0668567.0461，520.8848567.7212，518.58243565.7519，526.5444564.5603，524.5024565.1606，522.2207568.4056，523.0441 568.9876，520.7551567.7816，518.20574565.7835，526.3765564.6562，524.3918565.266
19、9，522.1878568.4048，523.0065568.9855，520.8005567.8299，518.80355565.8178，526.1958564.7594，524.2731567.2741，521.0509566.4548，524.1346568.5384，524.6946570.4308，523.52496563.6743，521.6294565.4214，520.6144567.3460，521.1294566.5451，523.9446568.4659，524.4789570.2377，523.47457563.8114，521.6232565.5112，520.65
20、25565.6125，522.1191568.3495，522.8911568.9111，520.9825567.9371，519.20718566.0185，525.5814565.1478，523.9097565.7272，522.1190568.2953，522.8554568.8405，521.0657567.9412，519.39879564.2483，521.6049565.7973，520.7773567.4596，521.2459566.6890，523.6426568.2713，522.7973567.9455，519.589310564.4428，521.5880565.9
21、130，520.7977567.4925，521.2496566.7358，523.5298568.327，524.0022569.2540，523.200911564.6655，521.5956566.0338，520.8531567.5067，521.2830566.6366，523.4134568.3226，523.8619569.7207，523.117012564.8876，521.6031566.1544，520.9094567.5207，521.3166566.9068，523.2978568.2923，523.7225569.5873，523.033713565.1091，52
22、1.6032566.2779，520.9556567.5416，521.3413566.9817，523.1791568.2704，523.5812569.4667，522.9465表22009年各层分面积的形心 塔层 坐标(m)A1A2A3A4A5A61562.4995，521.6767564.6728，520.3384567.1635，520.9321566.2995，524.4607568.7379，525.0660570.9705，523.73002562.7235，521.6673564.8138，520.3847567.1945，520.9648569.5556，521.55745
23、69.5059，523.8796570.8206，523.67273562.9424，521.6587564.9515，520.4308567.2241，520.9976566.3847，524.2815568.6415，524.8619570.6734，523.61734563.1175，521.6513565.0616，520.4670567.2484，521.0232566.4188，524.2100568.5920，524.7807570.5565，523.57275563.3050，521.6439565.1797，520.5062567.2741，521.0509566.4548，
24、524.1346568.5384，524.6946570.4308，523.52496563.6743，521.6294656.4214，520.6144567.3460，521.1294566.5451，523.9446568.4659，524.4789570.2377，523.43457563.8114，521.6232565.5112，520.6525567.3735，521.1561566.5799，523.8915568.4390，524.3948570.1656，523.39808564.0303，521.6142565.6546，520.7151567.4165，521.2012
25、566.6344，523.7572568.3961，524.2642570.0504，523.34339564.2483，521.6049565.7973，520.7773567.4596，521.2459566.6890，523.6426568.3533，524.1329569.9357，523.287610564.4428，521.5880565.9130，520.7977567.4925，521.2496566.7658，523.5298568.3527，524.0022569.8540，523.200711564.6655，521.5955566.0338，520.8535567.50
26、67，521.2830566.8366，523.4134568.3226，523.8619569.7207，523.117012564.8876，521.6031566.1544，520.9094567.5207，520.3166566.9068，523.2978568.2923，523.7225569.5873，523.033713565.1091，521.6032566.2779，520.9556567.5416，521.3413566.9817，523.1791568.2704，523.5812569.4667，522.9465表32011年各层分面积的形心 塔层 坐标(m)A1A2A3
27、A4A5A61562.49957，521.6766564.6856，5204421564.6730，520.3383567.1637，520.9320567.1413，520.9436568.4800，523.13632562.7237，521.6671564.8140，520.3844567.1947，520.9645566.3432，524.3695568.7042，524.9621570.8208，523.67253562.9428，521.6591564.9518，520.4312567.2245，520.9979566.3851，524.2819568.6419，524.862257
28、0.6737，523.61764563.1179，521.6509565.0620，520.4666567.2488，521.0227566.4192，524.2096568.5924，524.7802570.5569，523.57225563.3056，521.6433565.1803，520.5056567.2747，521.0503566.4554，524.1339568.5391，524.6939570.4314，523.52436563.6750，521.6286565.4221，520.6136567.3467，521.1287566.5458，523.9439568.4666，5
29、24.4782570.2384，523.43387563.8122，521.6224565.5120，520.6517567.3743，521.1553566.5807，523.8707568.4398，524.3940570.1664，523.39728564.0311，521.6133565.6554，520.7142567.4173，521.2003566.6352，523.7563568.3969，524.2633570.0512，523.34249564.2493，521.6040564.7983，520.7764567.4607，521.3884566.6899，523.17305
30、68.3543，523.6633569.9366，522.818010564.4439，521.5869565.9141，520.7966567.4936，521.2485566.7669，523.5287568.3538，524.0011569.8551，523.199611564.666，521.5944566.0349，520.8524567.5078，521.2819566.8377，523.4123568.3237，523.8607569.7218，523.115912564.8889，521.6018566.1557，520.9081567.5220，521.3153566.908
31、1，523.2965568.2936，523.7212569.5886，523.032413565.1104，521.6018566.2792，520.9542567.5429，521.3398566.9830，523.1777568.2717，523.5797569.4680，522.9450表4每层六个三角形的面积S1、S2、S3、S4、S5、S6，代入测量数据，运行MATLAB软件得出：13层塔中各层划分后的三角形面积 塔层 面积(m2)S1S2S3S4S5S616.106714.639720.606821.002514.70746.256325.730713.640219.066219
32、.390213.58645.839735.370612.696217.612117.869412.535.446145.099511.970216.501616.705911.72115.140554.813711.214115.350915.505810.87774.822864.538710.502914.256414.389210.14.523974.059.267212.639112.73328.92594.002883.6.268.13911.158211.21717.85133.528293.17897.08179.77769.80146.84873.0811102.80146.1
33、6058.55718.54015.96732.6792112.47635.34497.53127.4125.23612.3588122.16844.58736.56476.36134.558372.0604131.89313.90795.69685.42383.93922.6792表5多边形的形心（xc,yc,zc)公式： 将表1，2，3，4，5中的数据代入形心公式中，MATLAB软件运行的结果如下： 1986年塔的中心坐标 塔层 坐标(m)XcYcZc1567.1157524.37331.78742567.1355524.25837.32033566.38524.267212.7553456
34、7.6101523.857317.07835567.1858523.987921.72056567.2786523.784626.23517567.0062522.346129.83698566.9845522.506633.35099567.0212522.486836.854910567.0562522.470740.172111567.0363522.34644.440912567.1611522.389348.711913567.312522.402452.8293表61996年塔的中心坐标 塔层 坐标(m)XcYcZc1566.6531522.76061.7832566.355852
35、2.69047.31463566.7792522.631712.75084566.6888523.089217.07515567.0364523.466621.7166567.04522.557926.22957566.9986522.132829.83238566.9059522.101133.34549566.8944521.920636.848310567.071522.304540.167611567.1643522.340144.435412567.2123522.299748.707413567.5044522.404552.8253表72009年塔的中心坐标 塔层 坐标(m)Xc
36、YcZc 1566.7412522.70051.76452567.7195521.77797.3093566.8095522.643712.73234566.8368522.623417.06985566.8656522.59721.70946566.9541522.547226.2117966.9871522.530529.82468567.0394522.493233.42949567.0915522.460436.853810567.1463522.589540.636811567.1897522.366844.432612567.2313522.069348.699813567.279
37、9522.282252.8184表82011年塔的中心坐标 塔层 坐标(m)XcYcZc 1565.8577520.91661.76332566.7785522.67067.29053566.8098522.644112.72694566.8374522.619717.0525566.8662522.596421.70396566.9548522.546826.20457566.9879522.524629.8178567.0402522.492333.33669567.0925522.260236.822310567.1474522.405740.269111567.1901522.3656
38、44.42512567.2326522.32648.683913567.241522.32652.8131表9找出古塔各层中心位置的通用方法：将每层平面多边形划分为6个相互独立、互不相交的三角形，并求出每层的6个三角形的形心和面积，运用组合图形的形心坐标公式 得出古塔各层中心坐标（xc，yc,zc），即中心位置；表6为1986年古塔各层的中心坐标，表7为1996年古塔各层的中心坐标，表8为2009年古塔各层的中心坐标，表9为2011年古塔各层的中心坐标。5.3-1 问题2-1.1自身倾斜模型类别塔身第一层的中心塔尖 坐标(m) 年份 xoyoZOXYZ1986567.1157524.37331
39、.7874567.2473522.24455.12331996566.6531522.76061.783567.2544522.236755.11982009566.7412522.70051.7645567.336522.214855.0912011565.8577520.91661.7633567.3375522.213555.087表10各年塔尖平面坐标(x,y）,各年塔身第一层的中心平面坐标(xo,yo)；塔身在地面上的水平投影D: 塔上部与下部的相对位移H: 在直角三角形OPQ中存在： tanWt即为倾斜度；倾斜角W: 倾斜率： r=arctanWt 代入表10的数据运算得出，W19
40、86=1.5308,W1996=1.5558,W2009=1.5564,W2011=1.5339；r1986=2.1334,r1996=0.7975,r2009=0.7679,r2011=1.9677；即1986年的塔身倾斜了1.5308，1996年的塔身倾斜了1.5558，2009年的塔身倾斜了1.5564，2011年的塔身倾斜了1.5339。5.3-2 问题2-1.2 相对倾斜模型对各年份各层的中心点数据（表6、7、8、9）进行空间直线拟合，利用MATLAB软件拟合各年份塔中心线，结果如下： L1986 x=0.0028z+567.0229 y=-0.04739z+524.57 L1996 x=0.015z+566.50