《数学模型》课件——第八章 古塔的变形.pptxVIP

古塔的变形;历史文物古塔由于建筑年代久远，长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。

某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。;请根据提供的4次观测数据，讨论以下问题：

1.给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。

2.分析该塔倾斜、弯曲、扭曲等变形情况。

3.分析该塔的变形趋势。

;1.问题的重要性分析

古建筑是人类重要的文化遗产，由于建筑年代久远，经受长期的自然损坏和人为损坏，所以对其进行维修和加固是一个必要而艰巨的任务。对现存历史文物古塔的形变情况进行研究，是进行维修和加固的前提，也是对国家不可再生的文化资源进行保护的基础。针对研究出的形变情况采取合理的结构加固措施，可以延长古建筑结构的寿命，保护人类共同的精神文化遗产。;2.问题思路的分析

对问题一的分析：

由已知数据中给出的每层的八个点坐标，可以画出任意一层的二维和三维图，确定古塔的横截面图形，找出重心即为中心点坐标；也可以根据中心点坐标的意义，由中心点到八个测量点任一点的距离相等，建立关系，使中心点到各个点的距离平方和最小找出中心点坐标；还可以根据同一层的塔角点在圆曲线上的关系，通过拟合圆曲线的方法找出中心点的坐标。;2.问题思路的分析

对问题二的分析：

若分析该塔的倾斜、弯曲、扭曲等变形情况，则需要考虑塔身各层内部八个点的变化，塔身各段的变化，和塔身整体的变化（也要特别关注突出点的变化）。由此可以建立每一个点关于坐标轴的方向角关系，进行横向分析，纵向分析，局部分析，整体分析。即可知塔的变形情况。;2.问题思路的分析

对问题三的分析：

为了求得塔的变形趋势，我们可以在问题二的基础上把每年的变形度和时间（年份）进行曲线拟合，得出函数关系式，来预测未来。;3.1模型的假设

1.假设测绘者测量的数据真实，可靠，没有误差；

2.假设2009年、2010年中塔顶层的数据为该塔层的中心点坐标；

3.假设四年的测量数据来源于同一测绘基点；;3.1模型的假设

4.假设该塔每层成圆柱状，每层数据的八点在外边缘面上。

5.假设误差范围在时，为可接受误差区间；

6.假设最初修建该塔时，塔身垂直于水平面。

;3.2符号说明

：相对于轴的方向角；

：相对于轴的方向角；

：相对于轴的方向角；;4.1.问题（1）的模型建立与求解

4.1.1求解中心的方法一

首先任意选择附录表中某个层面上测量的点的坐标，画出二维图和三维图，确定塔身横截面近似于一个正八边形，;四模型的建立与求解;四模型的建立与求解;我们假设测绘者在每一个层面上测量的点的坐标都是在塔的边缘，以此进行描点画图，发现每个层面的测量点的都可以围成圆形曲面，则整体可以判定出每一层的测量点将塔围成圆柱。则同一圆曲面上，所以通过拟合圆形曲线的方法[3]来获得中心点坐标的。;设（A，B）中点坐标，（）为测量点。

通常的圆曲线可表示为：

展开得：

令：

;可得圆曲线方程的另一种表达形式为：

当有个点（）时，所拟合的圆方程可以表示为：

即所拟合圆曲线的平面坐标为：

再以此编录程序，并加以回归分析。

;在回归分析中，发现在1986年和1996年有异常点出现，而且都在第四层数据中的第一个点和第十二层中第四个点为异常点（对每层中的八个点进行绘图，重新整理坐标点顺序，确定两两相邻点的坐标）。我们可以考率这几个点在变化中进行了单一的变化，并没有随同塔的整体趋势进行变化，因此，在求曲面的中心点坐标中为保证数据的误差较小，先将异常点进行剔除，再从新求得对应层次的新的中心点点坐标，确定中心点坐标。

;确立的中心点坐标中缺少值，因此可以根据同一圆柱形曲面的的坐标相同的原理，计算方法与方法一相同。

根据上述算法，即可求出各次测量的古塔各层中心坐标。;;;4.2．问题（2）的模型建立与求解

由于塔的变形情况分为很多种，且每层、每个时期的变形情况都不相同，我们可以逐步分析每一层之间的变化，上一层对下一层的变化，以及每个时期的整体变化。分别从倾斜，扭曲，和旋转来分析。;4.2.1塔体倾斜度

本文通过两种指标来度量塔体倾斜度。

由4.1.1得到的每一层的中心点坐标，根据倾斜后的两中心点和自身的高度确定偏转角度，定义（1）：

倾斜角=;两中心点的距离=

高度==;利用matlab软件求解（具体程序见附录（3））

得出每一年的中心坐标的倾斜角度：;中心点的倾斜角度;定义倾斜度方法二，建立每一层中心坐标点关于对Z坐标轴的方向余弦，进而得出坐标点关于坐标