与封堵溃口有关重物落水后运动过程数学建模.docVIP

全国第七届研究生数学建模竞赛 题 目 与封堵溃口有关的重物落水后运动过程的数学建模 摘 要 本文通过分析物体在水中运动过程的会所受到各种影响因素如物体形状大小、水流速度，运用牛顿第二定律，对与封堵溃口有关的重物落水的运动过程建立数学模型，模拟运动轨迹，并将其推广到情况更为复杂的情况。 问题一：根据流体动力学在流体中运动的重物会受到压差阻力和粘性阻力，分析重物落水后的受力情况，利用牛顿第二定律及冲量定理建立数学模型，在MATLAB环境下用龙格-库塔法求解所列方程组并模拟重物落水后的运动轨迹，与试验轨迹相比误差较小。 问题二：在问题一的基础上，通过分析影响重物在水中运动的因素，把影响因素参数化，从而将原有的数学模型推广到可以模拟更为复杂情况的模型。 问题三：根据模型建立的过程，误差的产生来源于实验过程的误差，试验影响因素的误差，数值计算的误差，通过求均值误差来估计模型的合理性，并做出猜想和验证。 问题四：根据几何相似和重力相似，将所建立的模型推广到复杂流体上。在假定物体形状和投放方式的前提下，根据所建模型，带入参数，求出运动轨迹，从而确定理想的投放距离。 关键词：牛顿第二定律、压差阻力、冲量定理、粘性阻力、龙格-库塔法、相似准则 参赛队号 队员姓名 中山大学承办 与封堵溃口有关的重物落水后运动过程的数学建模 摘 要 本文通过分析物体在水中运动过程的主要影响因素，如水流速度，物体形状，物体下落高度，运用牛顿第二定律，对与封堵溃口有关的重物落水的运动过程建立数学模型，模拟运动轨迹，并将其推广到情况更为复杂的情况。 问题一：根据流体动力学在流体中运动的重物会受到压差阻力和粘性阻力，分析重物落水后的受力情况，利用牛顿第二定律及冲量定理建立数学模型，在MATLAB环境下用龙格-库塔法求解所列方程组并模拟重物落水后的运动轨迹，与试验轨迹相比误差较小。 问题二：在问题一的基础上，通过分析影响重物在水中运动的因素，把影响因素参数化，从而将原有的数学模型推广到可以模拟更为复杂情况的模型。 问题三：根据模型建立的过程，误差的产生来源于实验过程的误差，试验影响因素的误差，数值计算的误差，通过求均值误差来估计模型的合理性，并做出猜想和验证。 问题四：根据几何相似和重力相似，将所建立的模型推广到复杂流体上。在假定物体形状和投放方式的前提下，根据所建模型，带入参数，求出运动轨迹，从而确定理想的投放距离。 关键词：牛顿第二定律、冲量定理、粘性阻力、龙格-库塔法、相似准则 一、问题重述 1.1 背景分析 我国经常发生洪水，溃坝溃堤进而引发泥石流灾害造成国家和人民生命财产的严重损失。历年来的洪灾，尤其是最近的溃坝、泥石流险情给了我们深刻的教训：必须有效地开展封堵溃口的研究。由于溃口水流的流量和速度会比较大，在通常情况下很难在短时间之内将溃口彻底封堵，但如果通过投放重物对尚存的坝体产生一定的保护作用，就可以延缓溃坝溃堤的过程，为人民群众的撤离争取更多的时间。 1.2有关情况 溃坝事件常在猝不及防的情况下发生，很难从容地获得既完全又可靠的原型体的实测资料，因此研究溃坝波的主要途径是模型试验和数值计算，：水的运动粘度 ：雷诺数 ：压差阻力系数 :物体完全浸入水所用的时间 ：粘性系数 ：物体刚浸入水时垂直方向的速度 ：与速度方向垂直的表面积 ：物体刚浸入水时水平方向的速度 为大实心方砖的长宽高 三、模型的分析与建立 3.1问题一的分析与解决 3.1.1分析影响重物在水中运动过程的因素 (1)影响因素分析 物体在水中运动时受到重力，浮力，压差阻力，粘性阻力等力的作用，根据上面分析，可得影响物体在水中运动过程的因素有：物体的形状，物体的重量，物体的投放高度，物体的投放方式及水流速度。 (2)相关概念 根据流体静力学，流体对运动物体的阻力，有粘性阻力、压差阻力[1]。 粘性阻力：粘稠流体的分子结构，导致不同流层之间存在一种阻碍相对运动的阻力图1(水流为0.34m/s投放高度为0cm) 图2（水流为0.34m/投放高度为12cm） (1)由图1我们可以看出物体投放方式不同，相同条件下同一重物的运动轨迹不同。 (2)对比图1、图2，水流速度相同，投放高度不同，以相同方式投掷方砖的投放高度越高方砖在水平方向上的位移越小。 图3 (水流速度为0.55m/s投放高度为0cm) 图4（水流速度为0.34m/s投放高度0cm） (3)分别对比图1与图3，可以发现大实心方砖投放高度相同，水流速度越大方砖水平方向的位移越大。 图5(水流为0.34m/s投放高度为0cm) 图6(水流为0.34m/s投放高度为12cm) (4)分别对比图