研究生数学建模竞赛优秀论文(最终版)C.doc

全国第三届研究生数学建模竞赛 题 目 维修线性流量阀时的内筒设计问题（C题） 针对问题1，首先考察了内孔为四种特殊形状的情况下，“过流面积”随曲线下降距离的变化情况，得到凸凹圆曲线与严格线性面积特性曲线偏差的平方和最小，线性关系保持得比较良好。此后利用微元法证明了“过流面积”呈严格线性变化时曲线和外孔圆交点横坐标的差为定值这一性质，得出了在此种情况下曲线在两交点处的斜率应为无穷大。基于以上分析，利用最小二乘原理建立了无约束泛函极值模型，采用了变分法将其转化为微分方程，再转化为等效的变分原理，采用Ritz算法近似求解。最后通过对内筒孔曲线的合理假设，得到了满足线性关系较好的内孔曲线形状（见图11），其样本点的偏差平方和为0.064412。 针对问题2，利用最小二乘原理建立了有约束泛函极值模型。根据文中第四节中的引理，给出理想状态下的内孔形状。之后对其进行了微调，通过牺牲严格的线性关系来使其逐渐满足两个约束和，并最终找到了合适的内孔设计方案（见图13（b））。最后针对外孔磨损情况提出了基于自动控制理论和逆向工程技术等的解决办法。 本文提出的模型是从考察内孔的特殊形状中得到启发的，从而具有实际应用价值和准确性。 关键词：线性阀体 最小二乘法 泛函极值模型 变分原理 非线性规划 一、问题的提出 阀体是我们日常工作和生活中一种十分常见的工具。它种类繁多，其中线性阀体可使阀体的旋转角度和流量成正比。因而它可使人们方便地对流量进行控制。而如何设计线性阀体成为当今控制领域中研究的热点问题之一。 现在我们需要设计出一种阀体，它由两个同心圆柱筒组成。外筒固定，其侧面上有一个孔，形状为两个直径不等的圆柱体的交线。内筒和外筒轴向之间没有相对运动，内筒可以自由转动。内筒的侧面上也有一个孔，但它原来的形状未知。 要求设计出内筒孔的形状，使得“过流面积”与内筒旋转角成近似线性关系；在线性区间至少达“最大范围”区间长度的75%以上，而且主要工作区的最大“过流面积”至少要达到外筒孔面积的85%以上，并且使“过流面积”和内筒的旋转角度之间的“线性关系”尽量好的约束限制下，重新设计内筒孔的形状。并且还要考虑当外筒孔发生磨损时要采取的应对措施。 二、模型假设 1、阀体的旋转角度与内圆筒相对移动距离成正比，圆筒移动距离与“过流面积”成正比。 2．线性阀体内外筒为薄壁筒，不考虑其壁厚给设计带来的影响。 3、外圆筒直径与外圆孔直径相差很大，展开后外圆孔面积变化足够小，可近似视为圆形。 4、内筒在转动过程中，只存在周向水平运动，不存在垂直方向的运动。 5、假设内圆孔设计曲线与外圆孔曲线最多只有两个交点，可以有一段相切，且曲线连续。 6、为简化计算，假设外圆孔半径为一个单位长度。 三、变量设定 ：圆的半径，在本文中为一个单位长度1； ：待求内孔的曲线方程； ：内孔下边沿曲线方程； ：外圆孔上半圆方程，即圆的方程； ：曲线下降的距离微元； ：曲线下降到某一位置时其与初始位置的距离； ：曲线从初始位置下降至“过流面积”达到最大值时的距离； 、、、：分别表示曲线F(x)在移动过程中与曲线G(x)的交点； ，，，：分别表示点、、、的坐标值； ：曲线下降的距离与“过流面积”之间的线性比例； ：曲线下降时“过流面积”的增加量； ：“过流面积”的理想值，。 四、问题的分析 本文将内外两个圆柱筒展开为平面，得到两个长方形，于是将三维空间中物体的转动问题化简为二维平面上内孔与外孔相对移动的问题来求解，此外根据问题假设可将外筒孔近似视为圆孔。 建立如图1所示直角坐标系，用以坐标原点为圆心的单位圆来表示外圆孔，X轴与内、外筒的轴心平行。用任意曲线表示内圆孔曲线初始位置时的一部分，另一部分与其组成封闭图形，但是未画出的部分与圆不相交，如图1（a）所示。 （a） （b） 引理：若要使内孔旋转角度（称为开度）与“过流面积”满足线性关系（这种关系称为面积特性曲线），则内孔曲线必满足其与外孔圆的交点横坐标之差衡为常数，即，其中分别为内孔曲线与外孔圆的交点横坐标。或者可以说即为面积特性曲线保持线性的必要条件。 证明：假设某一时刻内孔曲线向下移动与圆相交，其方程为，当曲线向下移动微元时，“过流面积”的增加量由三部分组成，两边近似三角形面积和中间矩形面积（如图1（b）所示），并可用以下积分表示： （1） 若要使内孔旋转角度（称为开度）与“过流面积”满足线性关系（这种关系式称为面积特性曲线），则只须使曲线的向下移动距离与“过流面积”满足线性关系即可，即微元面积也与有线性关系： （2） 曲线与圆的交点坐标x由方程（表示下降时的曲线）求得： ，