弹性棒及弯曲变形

1、弹性棒的弯曲变形研究案例背景 对于某一均匀圆柱形细长弹性棒，在棒的两端施加方向相反大小相等（等于F ）的轴向压力. 实验表明，当且仅当 F 大于某力 F1 时，棒才会发生弯曲. 这个 力 F1 称为临界力. 而当 F = F2 2F1 时，棒的平衡状态如示意图 1 那样：棒的两个端点重合（假设棒仍处在弹性限度之内）. 试建立数学模型来验证实验的结果.并且完成以下任务：一、计算力学特征：1. 临界力 F1 的表达式；2 F2 / F1 的精确的理论比值.二、计算 F = F2 时棒的平衡状态曲线的几何特征：1宽度与棒长之比;2高度与棒长之比;3端点重合处的夹角(单位：度).示意图 1一问题的分析

2、仔细思考问题的条件和要求，思考如下：1. 工作前奏概念理解F2 / F1 的精确的理论比值，通过对材料力学和数学的理解，并查阅了相关资料，认为它是指相对于给定的,这里的精确不是指误差为0，而是指在允许的范围内较2更为精确的比值，比值越精确，模型越是合理。2. 中心任务建立在作用下弹性棒的的曲线方程并求解（1） 问题转化。问题要求确定一个精确比值，以及计算出 F = F2 时棒的平衡状态曲线的几何特征；即使要求建立一个合理的曲线方程模型来表示弹性棒的变形，并通过解此方程球的相应的几何特征。（2） 具体方法。具体方法是否得当，不仅仅是该模型建立本身的合理性、逻辑严密性，而且要求该模型具有可解性，即

3、该模型能够通过现有的数学手段求解。从F1的求解入手，层层递进，使模型逐步优化。1） 按照纯弯曲理论，对细长杆小变形压杆失稳用梁挠曲线近似微分方程式按照材料力学公式求解即欧拉公式。2） 在1）的基础上，对大变形理想弹性棒进行分析，据结构力学对称性分析，将F2作用下的弹性棒简化为一端固定，一端自由的弹性棒，类比欧拉公式，得出其挠曲线微分方程。3） 对该方程进行进一步分析和简化，发现其是不可积方程，故对该方程利用进行优化，通过在固定端的近似求解，将弹性棒的变形曲线分解成两段。通过边界条件，近似解出该方程。3. 目标工作得到该方程后，代入相应条件，得出曲线的几何特征。4. 模型反思在支座附近曲线求解上

4、采取了近似处理，将整个棒上弯矩取为常数自由端弯矩，事实上，并非常数。截断点离固端越近，求解越是精确。 二 模型假设1. 棒是均匀的，圆柱形的，细长弹性的。2. 棒的截面尺寸相对于其长度为无穷小。3. 在模型、中，弹性棒发生的是纯弯曲变形，轴力和剪力对弹性棒的变形不起作用。4. 在模型中，在作用下，棒处于临界状态，可认为杆发生了极其微小的侧向位移。5. 在模型中,不同杆长所形成的曲线是相似的,即不同杆长所形成的曲线的相应几何特征之比应等于杆长之比. 三 符号约定：弹性棒的截面抗弯刚度（其中为弹性模量，I为截面惯性矩）。M： 截面弯矩,规定以逆时针为正。：曲率半径。F：压力。:曲线上任一点的斜率。

5、 四 模型的设计和求解根据问题分析，问题先转化为平衡状态曲线方程的确立，再以此确定相应的F值和几何特征。1. 模型（1） 模型的建立和求解根据问题分析，可将该弹性细长棒等效为两端用球形铰支座支承的细长杆，在它的两端作用轴向压力，当F的的大小到达时杆因水平方向的微小干扰而产生微小弯曲，并在此微弯状态下保持平衡。因属于小变形，可以梁挠曲线的近似微分方程式描述，即 图一 临界状态条件下 (1-1)由图：在任意截面处杆的玩具M(X)有 代入(1-1)式得： ，令=,则 (1-2)其通解为，代入边界条件X=0,Y=0 X=L,Y=0解得 B=0, KL=N（N为整数）， 代入=,得 （1-3）取N=1，

6、即得临界荷载 （1-4）即为欧拉公式，此公式由欧拉导出。2，模型模型只能解决小变形问题，在作用下棒的弯曲属于大弯曲问题，须从梁在纯弯曲时绕曲线的曲率公式从头推起。(1) 模型的建立根据结构力学对称性分析可将细长弹性棒简化为如图所示，取上半段进行研究。梁在纯弯曲时绕曲线的曲率公式为 （1-5）图二 状态下在这里，M,已不是常量，而是X的函数，可将上式改写为 （1-6）而 （1-7）代入得 考虑到M（X）0,<0 （从曲线大致形状得出）取 （1-8） 将M(X)=代入上式，并令K=,并积分得（C为常数）代入边界条件 ,得到C=,代入上式得 （1-9）将进行整理得 从而得曲线表达式： （1-1

7、0）利用棒长为,得（1-11） 上述方程组仅含两个未知数，理论上可解。（2） 解模上述积分方程均没有原函数，故不能用常规方法求解，考虑到条件F2 2F1，则可以通过计算机搜索法求解.具体如下:由F2 2F1, K=,再通过曲线大致图像得=0，及在 恒成立得出，有假设5取L=1，理论上则可通过MATLAB编程得到K,D的解。具体算法如下：开 始将k，D代入方程中，得出曲线方程不成立，则返回结束判断等式（1-11）是否成立成立，则输出K，D的值给定K，D的值，代入方程组（1-11（3） 通过编程运行发现，得不出想要的结果。通过对方程组（1-11）进行分析发现，由于被积函数在X=D的奇异性，使得函数

8、无论通过符号还是数值都无法求解。故模型失败。3 模型模型无法求解，是因为被积函数存在奇异性，且在在理论上是不可积函数。故而须用近似的模型代替模型进行求解，但是近似模型必须具有足够的精度。故而有了模型.（1） 模型的建立 奇异点的处理图三 D-D段的近似计算 在B点附近（如D点），弯矩M已经足够大，从D到B点M增量相对于总量可忽略不计，即将M近似为常数。设=L,（可通过模型求出），在D-D段，由纯弯曲理论不难得 （1-12）其中，同理可解得 （1-13）则D-D的长度 为 （1-14），可积理论上，只要D点越靠近B点，就越精确。在A-D区域，仍采用模型的函数。（2） 解模算法A-D的曲线长度为，

9、D点纵坐标=，弹性棒总长度为L=2+ （1-15）方程（1-15）虽含有两个未知数K、D，但还是可以通过计算机算法求解。因为这里在建模时就已经认为A，B点在同一水平面，这是隐含的一个条件。同模型，可用MATLAB编程进行循环搜索。算法流程图如下：开 始将k，D代入方程中，得出曲线方程不成立，则返回结束判断等式（1-15）是否成立成立，则输出K，D的值给定K，D的值，代入方程组（1-11初步计算结果经过编制程序1计算，初步得出结果如下：组数KDL119.478327481851911.00926702107444 219.578327481851910.40805716292769 1.0066

10、8043730346319.678327481851910.40701902364968419.778327481851910.40598876766673519.878327481851910.404966295709280.99903857352555619.978327481851910.403951510249060.99652975593795720.078327481851910.402944315460010.994039771649438可以发现K值得精确解应该在19.420.2之间，D的精确解应该在0.40左右，的精确解应该在0，110.12之间.计算结果优化编制程序2，在K

11、在9.420.2之间和D 在0.4附近对结果进行进一步优化。优化结果如下（部分数据）组数KDL119.223 0.3970.00945530068819 0.99991892959063219.2240.3970.009483004072000.99999655550612319.2250.3970.009510709471311.00007420031572419.2260.3970.00953841688707519.2270.3970.009566126320241.00022954664339上面五组数据是所有数据最合理的解。选择最优解，故最后结果可取K=19.224 D=0.397将

12、其带入原函数解得将其带入原函数解得=0.2315，= = 从而得到原问题的解1) 2) F2 / F1= 1.948,3) 宽度与棒长之比,即D，为0.397，4）高度与棒长之比即2=0.253436，取0.2535)端点重合处的夹角(单位：度) 2=进一步优化由假设可知，到D点越靠近B点时，计算结果越是精确，故在程序二的基础上将D点向B点靠近，取L=0.385，得到程序三。部分优化结果如下：组数KDL1 20.981 0.3960.0472652366350 1.00001343189113220.8520.3970.048692716554931.00003448116483320.471

13、0.4000.052452069591460.99997896825555420.5970.3991.000015610564725 20.8510.3970.048659791347410.99994593507208 上面五组数据是所有数据最精确的解，最后结果可取上面五组数据的平均值，即得K=20.7504， D=0.3978将其带入原函数解得将其带入原函数解得=0.2487，= = 从而得到原问题的解1) 2) F2 / F1= 2.1024,取2.103) 宽度与棒长之比,即D，为0.3978，4）高度与棒长之比即2=0.252218，取0.2525)端点重合处的夹角(单位：度) 2=

14、即79.86度。模型的最后优化上面得到D=0.3978，故可取=0.390，同上，将程序三做微小改变，得程序四，得如下数据：组数KDL120.7645 0.3975 0.03688617058714 0.99999176154476220.76500.39750.036902766125491.00003587991195320.76550.39750.036919362743201.00008000291076420.82850.39700.035865992906080.99992892320689520.82900.39700.03588256886970.99997298214455故

15、得最终结果K=20.7645， D=0.3975将其带入原函数解得将其带入原函数解得=0.2484，=从而得到原问题的解1) 2) F2 / F1= 2.1039,取2.103) 宽度与棒长之比,即D，为0.39754）高度与棒长之比即2=0.252584，取0.25265)端点重合处的夹角(单位：度) 2=79.65度。此时精度已足够高，故可以此为最优解。 五 .问题的进一步分析上述模型解决了在纯弯曲条件下弹性棒在作用下的变形曲线。但是实际情况下，弹性棒变形时受多方面影响的，从以下几个方面说明：1. 轴力和剪力对变形曲线的影响。弯矩能引起结构的弯曲变形，在轴力和剪力作用下，构件会发生相应的拉

16、压和剪切变形，对曲线将会有一定的影响，但是如果考虑其影响，模型的建立无疑是非常复杂的，甚至可能没法求解。但是在理想细长弹性棒的条件下，模型是合理的。2. 构件截面的不均匀性实际的构件，由于加工、材料等各方面特性，各截面抗力会在一定范围内波动。如果假设弹性棒沿长度方向是非均匀的，求解将更加复杂。3. 任意F作用下弹性棒的变形在任意F作用下，弹性棒会发生相应的变形，同的求解类似，通过建立合理的坐标系，选择正确的边界条件，总是可以求出相应的曲线方程的。4. 材料抗弯能力对曲线的影响。实际构件中，材料都有一定的强度，当变形过大时，局部应力超过了材料抗力，材料可能会脆断或发生塑性变形，不在保持弹性状态。 六 模型的评价1. 模型的优点（1） 原理简单。整个模型建立过程中都是用的最基本的材料力学