对弹性梁形状的研究(数学建模)[].docVIP

对弹性梁形状的研究(数学建模)[1] 对弹性梁形状的研究 摘 要 本文主要针对弹性梁在静态平衡时的形状进行研究。以一悬臂梁为例，建立以固定点为原点，墙所在的直线为Y轴，且向下为正的平面直角坐标系，以利于分析描述中性线受力时的形状变化，并且可转化为数学上的曲线来分析。 对问题一，在建立无量纲化的数学模型之前，首先根据材料力学相关知识建立确定弹性梁形状的一般模型，得到度量梁变形后横截面位移的基本挠曲线方程和转角方程，然后根据材料了力学上的叠加原理对梁是否与地面接触分为两种情况讨论，引入各个以梁长为特征长度的无量纲的参数，并提出无量纲变换的控制方程，使模型最终得以量纲化而与梁长无关。 对问题二，由于不与地面接触时的挠曲线函数是与地面接触时挠曲线函数的特殊情况，故仅需要对更加普遍的情况进行讨论研究。首先对中性线的形状做定性分析，由于挠曲线与参数a,p有关，会随参数的变化而变化，从而使中性线的形状发生变化。分为两种情况分析讨论：（1）令a?24，通过MATLAB绘制三维曲面来对不同p值下的中性线形状进行分析，发现p值对中性线形变时的横向位移和纵向位移都会产生影响；（2）1，同样从三维图上定性分析对不同a值下的中性线形状的变化情况，发现a值只2 影响中性线形变时的纵向位移而不影响横向位移。然后，通过绘制三维图分别对斜率、曲率、极值点以及拐点等决定中性线形状的相关量与参数的关系进行研究，间接而又深入地分析了中性线的形状。 令p? 对问题三，由于梁是否与地面接触与参数p有关，提出判据：若p在(0,1)范围内也即0?F/gDL?1时则与地面接触，若p?1即F/gDL?0时，梁不与地面接触。根据已2121gDL。确定了参数的值，，此时F?3232 则挠曲线函数也唯一确定，绘制出中性线图像。此外，对曲率的最大、最小值进行求解， 21得到最小值为kmin(p)??1.418，负号表示方向，对应的坐标为(,?1.418)，最大值为32给出的条件，结合已有的挠曲线方程求得p?kmax(0)?3.75，对应的坐标为(0,3.75)。最后求解k?0时中性线上点的集合，为 5,0)}。 16 综上所诉，本模型充分运用数学曲线的性质，并且结合材料力学??的叠加原理来研究弹性梁受力时的发生的形变，具有一定的可靠性和科学性。 关键词：叠加原理 无量纲变换控制方程 挠曲线方程 {(1,0),( 1 1、 问题重述与分析 弹性梁等杆件受到外力的作用会发生轴向和纵向的变形，均匀薄板状的弹性梁在理想状况下（不受任何作用力包括重力）时状态为以平板。现有一端固定的弹性梁，假设其墙外部分的中性线受力时只考虑纵向弯曲，不考虑横向变形，即其长度L不变。以固定点为原点，对墙和梁建立平面直角坐标系，已知自由端的纵向位移为Y?H?0和梁在静态平衡时中性线满足的条件。 对问题一，为避免物理单位的选择对数值结果的影响，以L为特征长度，引入无量纲量，建立确定梁静态平衡时中性线的无量纲化的数学模型。度量梁变形后横截面位移的两个基本量是横截面形心（即轴线上的点）在垂直于X轴方向的线位移，称为该截面的挠度，以及横截面对其原来位置的角位移，称为该截面的转角[1]。由于梁变形后的轴线是一条光滑的连续曲线，且横截面仍与曲线保持垂直，因此，横截面的转角也就是曲线在该点处的切线与X轴之间的夹角。变形后的中性线即为挠曲线，考虑挠度和转角可以确定挠曲线方程，即为中性线的曲线方程，再引入可以取代最大挠度，最大转角等的无量纲量，便可以将模型量纲化。 另外，又考虑到梁与地面接触的情况，则除了受到重力的作用还受到地面对梁的支持力，挠曲线由于受力的改变而发生变化，可以建立新的挠曲线模型，同样引入无量纲量，得到无量纲化的数学模型。 对问题二，要求在斜率、曲率、拐点等与参数的关系等方面对中性线的几何形状做定性分析。由于在问题一中已确定中性线的挠曲线方程，并可以通过求导，来得到斜率、曲率等量，在将其与问题一中引入的参数联系起来，找出他们之间的关系。此外，可以通过对挠曲线方程求一阶导数和二阶导数来判断曲线上的极大值，极小值以及拐点，从而对中性线的形状进行定性分析。 对问题三，首先要求判断梁是否与地面接触，可以通过参数p来判断，若0?p?1，即存在一个地面对梁的支持力且小于梁的重力，则判断梁与地面接触，若p?1，即不存在一个地面对梁的支持力，F?0，则判断梁与地面不接触，p?0amp;p?1的情况是不存在的，因为梁的一端是固定在墙上的，地面对梁的支持力不可能大于也不可能等于梁的重力，出现负数的情况也是不符合实际的。然后要求画出中性线图像，便求k的最大值、最小值对应的点的坐标，可以