关于古塔变形的数学模型.docVIP

关于古塔变形的数学模型 摘 要 本文主要研究古塔在自重、气温、风力等因素的影响下产生变形的问题。采用中垂线求解外切圆圆心的模型以及多次平均除误差的方法，找到了确定古塔中心的通用方法，并用多元线性回归模型及插值拟合等方法对倾斜、弯曲、扭曲等变形情况进行分析，从而通过残差拟合得出预测数据对古塔变形趋势进行描述。 针对问题一：论文采用古塔八个角点中任意三个角点构成的两两连线，取其中垂线的交点得到外接圆圆心，已知正八边形的中心与外接圆圆心一致，但古塔八角点构成的八边形存在轻微不规则，所以我们采用多次取点求外接圆圆心，并用其平均值消除误差，最后对不同取点方式进行了精度分析（答案详见表一）。 针对问题二：首先是古塔倾斜分析，根据测量学本文取塔尖和塔底的中心连线作为倾斜角计算的倾斜方程，算出塔顶在水平面投影与塔底中心的间距，引入实测高程数据，可以得到古塔四次测量的倾斜角（），对其倾斜情况经行描述；然后是弯曲情况分析，根据问题一中古塔各层中点坐标，本文对其进行多元回归分析及多项式拟合，得出函数曲线，并将其和倾斜方程进行比较得到最大差值即挠度（材料力学中对弯曲的描述量）；最后是扭曲分析，本文分垂直和水平两个方向进行讨论，垂直方向上涉及高程Z，即对各层中心点多元线性回归得到的拟合值与实测值进行残差分析，得到扭曲描述量（）。水平方向，本文参考材料力学中扭转角的计算，对古塔各层间的轴向扭转进行分析，得到扭转角对古塔扭曲情况进行描述。 针对问题三：在分析了四次观测值中倾斜、弯曲，扭曲的情况下，本文采用加权平均的方法各产生影响数据进行处理后，进行残差拟合，得到下一次观测的模拟数据，对古塔的变形进行变形趋势描述 关键词 多边形中心确定 多元回归分析 多项式拟合 残差分析 一、问题重述 由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。 请你们根据附件1提供的4次观测数据，讨论以下问题： 1.给出确定古塔各层中心位置的通用方法，列表给出各次测量的古塔各层中心坐标。 2. 分析该塔倾斜、弯曲、扭曲等变形情况。 3. 分析该塔的变形趋势。 二、模型假设 根据测量学对倾斜的描述，倾斜角的计算为塔顶与塔底中心连线的倾角，假设各层中心的偏移不影响塔尖和底层的倾斜计算。 2、假设古塔的各层面是一个稳定整体，平面转动是同时进行的。 3、假设附件给出的数据是八角古塔各角点的坐标。 三、符号说明 1、：圆上任意三点构成的弦的斜率 2、：塔顶在底面的投影与塔底中心的距离 3、：塔顶的高程、塔层间距离 4、：位移变化量、残差值 5、：单位扭转角 6、：多元线性回归系数 7、：i次项变量 8、：参与计算的古塔层数 9、：扭曲描述量 四、问题分析 对于问题一，求解该层的中心，已知数据是八边形古塔的各角点，根据正多边形中心的定义，各对角线的交点以及外接圆的圆心是正多边形的中心，但是对于不规则的多边形，其对角线的交点并不重合，同时外切圆的圆心也会偏离。于是在不同选点情况下，圆心位置可能存在的偏差，我们采取多次取点计算平均值的方式来消除。由于八边形存在的种选点情况，于是我们对不同角点选取对圆心精度变化的情况进行了验证。 对于问题二，倾斜分析，在第一问得出的古塔各层中心的基础上，首先按照测量学对倾斜观测的方法，根据塔顶三维坐标与底层中心，可以得出各层高差以及在平面投影上的中心差，由此可得倾斜角：； 弯曲分析，在倾斜分析中我们已经得出塔尖和底层中心的连线，可以用各层中心点与该连线的偏离程度来描述古塔的弯曲程度，根据问题一得出的古塔各层中点坐标，通过多项式插值和多元回归分析可以得到中心点连线的回归方程，将此方程分别投影到，平面后，和塔的倾斜方程进行比较得到最大差值即挠度（材料力学中对弯曲的描述量）和差值总量可以得到古塔的弯曲描述； 扭曲分析，按照结构力学中扭曲是空间概念的说明，我们分垂直和水平两个方向对古塔的扭曲状况进行讨论，垂直方向上涉及高程Z，即对各层中心点多元线性回归得到的拟合值与实际值进行残差分析，得到扭曲描述量（）。水平方向，参考材料力学中扭转角的计算，对古塔各层间的轴向扭转进行分析，计算得到平面内位移量后通过计算古塔各层高差得到扭转角，对古塔平面扭曲情况进行描述。 对于问题三，古塔的变形趋势，首先对前四次测得数据中对倾斜，弯曲，扭曲有影响的变量根据观测年份间距进行加权平均，然后通过残差拟合得出下一次测量的预测值，再对预测数据进行倾斜、弯曲、扭曲的分析，比对前四次测