本科毕业论文古塔的变形—数学建模论文.doc

2013高教社杯全国大学生数学建模竞赛 承 诺 书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛 赛区评阅编号（由赛区组委会评阅前进行编号）： 2013高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）:  C题 古塔的变形 摘要 文物保护是一件非常重要的事情，如何根据对古塔的测量数据了解其倾斜、弯曲、扭曲等情况并能够确定变形趋势是本文要解决的问题。针对问题一，本文将各层八个观测点看作同一平面上一个八边形的角点，利用公式，并编写C++程序计算其中心坐标，有较高的精确度。问题二要求确定古塔的倾斜、弯曲、扭曲等变形情况，结合测量专业知识，利用公式计算得到四年各层中心相对于下一层中心的偏心距（合位移）与倾斜角度以及四年古塔整体倾斜的合位移与倾斜角度： 1986年 1996年 2009年 2011年 合位移（m） 0.3779 0.3739 0.3778 0.3285 倾斜角（弧度） 0.0129 0.0121 0.0113 0.0103 利用公式计算反映各段的弯曲，利用；及 最终确定出各年古塔整体扭曲情况： 1986年 1996年 2009年 2011年 扭曲方向（度） -34.3 -31.66 -36.52 -34.34 扭曲位移（m） 0.3752 0.3159 0.3035 0.2713 数据分析表明四年中底部1-3层及顶部11-13层扭曲现象严重，而中段4-10扭曲较平缓。应加强对底部及顶部的观测与维护。弯曲方面，四年相对变化不大，主要体现在各年不同层次的弯曲有一些变异，整体倾斜度不大且减小趋势。 关键词：重心（中心），倾斜、弯曲、扭曲 C++ MATLAB EXCEL 一 问题重述 由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。 请根据附件1提供的4次观测数据，讨论以下问题： 1.给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。 2. 分析该塔倾斜、弯曲、扭曲等变形情况。 3. 分析该塔的变形趋势。 二 模型假设 1.古塔变形测量的监测点选取符合变形观测的原则，测量数据有效。 2．楼层平面质地均匀。 3. 在古塔各层选定的8（或7个）个监测点在同一平面内。 三 符号约定 1．：第个观测点的坐标。 2. △x：各层中心相对于下一层中心在x轴方向的位移量。 2. △y：各层中心相对于下一层中心在y轴方向的位移量。 3． d：各层中心相对于底层中心的合位移。 4．：各层中心相对于底层中心的扭曲方向。 5．：各层中心扭曲方向的平均值，反映建筑物主体的扭曲方向。 6．h：各层中心相对于第一层中心的高度值。 四 模型分析建立与求解 问题一：确定古塔各层中心 分析:附件给出了测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对古塔进行4次观测结果，即古塔十三层及塔尖各监测测点的坐标值。 在建筑物的变形测量中每次都选取固定监测点观测。按附件提供的测测点的坐标值作图（图一），显现出正八边形形状，所以可以假设古塔为八角形，各层所设8个（7个）监测点为各层八角形的塔角。各层所设8个（7个）监测点坐标值的z坐标值虽有所不同，但差值最大仅为0.194m，可以忽略高度变化，故各层中心计算时假设8个（7个）塔角在同一平面内。 查阅测量相关资料知道：对于外轮廓为正八边形的塔形垂直的建筑物，需对8条棱进行观测，每层观测8个角点（在同一平面内），从8个角点构成八边形得到的几何中心就是该层中心。 在附件给出的4次观测数据中，还发现1986年和1996年数据中十三层却监测点5。根据古塔已上千年历史，可推断古塔1986年和1996年进行测量时，13层塔角5处于损毁状态，到2009年第三次测量前已经修复。 13层塔角5在两次测量无数据，故对1986年和1996年13层数据计算得到中心坐标误差较大，这对后面计算会产生影响。 模型建立：在高等数学中已知物理学中求平面薄片重心的公式为： ， 当薄片质量分布均匀，即密度为一常数时，分式上下可约分，此时分母转化为平面图形的面积，公式简化