古塔变形的数学模型.docx

古塔变形的数学模型 摘要 古塔被誉为中国古代杰出的高层建筑物， 历史悠久，值得并需要我们的保护。本文研究了关于古塔变形的问题， 古塔的变形与塔身的中心紧密相关， 具体 分析了古塔倾斜、弯曲、扭曲的变形情况及趋势。 对于问题一， 建立中心位置模型， 采用多边形组合形心的算法， 求的结果是表 6，7，8，9 中的数据； 对于问题二， 我们将古塔的倾斜、 弯曲和扭曲等变形情况， 分别给予合理的 ，即塔尖与底层中心的水平距离与塔高的比值；对于弯曲变形，我们定义了弯曲率 K ，即用中心点所拟合出的空间曲线的曲率来描述古塔各处弯曲率； 对于扭曲变形， 我们定义了相对扭曲度，利用坐标的旋转变换角度描述古塔的扭曲变形情况。 利用空间曲线拟合、 坐标变换等方法以及 MATLAB程序，分别求出了三个变形刻画量的量化指标。 对于问题三， 我们考虑通过古塔的倾斜、 弯曲及扭曲程度来分析古塔的变形趋势，由于数据较少，我们利用灰色预测法来分析这三种变形，利用 matlab 得到的数据来分析古塔的变形趋势 关键词：关键词： 组合图形的形心 拟合 中心变形曲线 灰色预测 一． 问题重述 古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 现有一古塔已上千年历史， 是我国重点保护文物。 管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了 4次观测。 请根据题目提供的 4次观测数据，讨论以下问题： 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。 分析该塔倾斜、弯曲、扭曲等变形情况。 分析该塔的变形趋势。 二、 模型假设 1，参照中国古塔的塔身样式假设塔身水平面是正八角形； 2，假设古塔是均质物体，形心就是中心； 3，以地平面为 X 轴，建立空间直角坐标系 X，Y， Z； 4，假设塔尖面积极小，近视为一个点； 5，假设模型 4中古塔所在地不存在地震等这种偶然事件； 6，假设古塔所在的地平面是水平面。 三、符号说明 Ci -- 第 i 个三角形的形心（ i=1,2, 6）； S 2 i -- 第 i 个三角形的面积（单位： m) （ 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-- 塔的曲率； 四、模型准备 ）给数据用图像形式表现出来题目所给的 4 次观测数据 1986 年、1996 年塔 13 层的第五个点数据残缺。参照其 10 层、11 层、12 层的数据知道第五个点的数值在第四个点和第六个点的数值之间， 另外分析后发现这两组数据对题目结果的影响不大。 因此结合实际和理论， 采取取相邻年份该层的第四个点和第六个点的平均值将两组数据补齐。 ）用 MATLAB软件将题中所给数据及补充的数据画图表现出来。 五、模型的建立与求解 问题一 塔的变形趋势与塔身的中心有关， 因此本题的主要目标是计算出塔中心的变化规律。 问题 1 属于计算多边形组合形心的数学问题。因为 1986 年受自然、社会影 响度较小，最接近古塔的原始图样，所以图形较为有参考价值，用 MATLAB画出 1986 年塔身的大致图样为八边形，通过假设得出正多边形的形心就是中心，因 此可运用组合图形形心公式计算出古塔各层中心坐标。 问题一求古塔各层中心坐标，通过假设知道实质就是计算正多边形的中心， 运用组合图形求形心的方法。 将每层正八边形的塔面都分为六个三角形（图 1-2）S1、S2、S3、S4、S5、S6， 每年各层六个三角形的形心： x1 x 2 x8  y1 y2 y8 C1( ) 、 3 ， 3 x2 x3 x8 y2 y3 y8 C2( 3 , 3 ) 、 x3 x4 x8 y3 y 4 y8 C3( 3 , 3 )