浅谈数值分析在机械工程领域的应用

1、浅浅谈谈数数值值分分析析在在机机械械工工程程领领域域的的应应用用摘摘要要： MATLAB 是目前国际上最流行的科学与工程计算的软件工具 , 它具有强大的数值分析、矩阵运算、信号处理、图形显示、模拟仿真和最优化设计等功能。本文浅谈 MATLAB 在机械设计优化问题的几点应用。关键词关键词： MATLABMATLAB 约束条件约束条件 机械设计优化机械设计优化 数值分析数值分析引言引言：在线性规划和非线性规划等领域经常遇到求函数极值等最优化问题，当函数或约束条件复杂到一定程度时就无法求解，而只能求助于极值分析算法，如果借助计算器进行手工计算的话，计算量会很大，如果要求遇到求解极值问题的每个人都去用

2、 BASIC,C 和 FORTRAN 之类的高级语言编写一套程序的话，那是非一朝一日可以解决的，但如用 MATLAB 语言实现极值问题的数值解算，就可以避免计算量过大和编程难的两大难题，可以轻松高效地得到极值问题的数值解，而且可以达到足够的精度。 数值分析是一门研究如何在计算机上求解数学问题算法的学科,主要内容有:误差分析,插值法,数值微积分,数值代数, 矩阵计算和微分方程数值解法等, 是工科各专业大学本科及研究生中开设的一门计算量大,算法多,实践性比较强的专业课。在长期的教学实践中,数值分析课程常采用C语言进行教学和实验, 要求学生既要对算法有充分了解,又要熟练掌握C语言的语法和编程技巧,

3、导致学生和教师将大量的时间和精力都花在繁琐的数值计算以及对各种结果绘图上面，学习效果往往令人不满意。M a t l a b 是M a t h W o r k s 公司开发的一款以数值计算为主要特色的数学工具软件, 在数值计算领域独领风骚。其所带强大的符号运算功能, 几乎包括高等数学所涉及的运算, 如求极限、导数、微分、积分、函数的级数展开、解常微分方程等等, 并且样条工具箱中的命令调用格式极为简单方便, 对工科学生来说, 掌握起来无需费多大力气, 而对机械系等理工科系的同学,通过初步了解M a t l a b还可以进一步挖掘其强大的功能, 对学习其他课程也有帮助。本文讨论基于matlab在机械

4、方面的数值分析。一数值分析方法的研究数值分析方法的研究1、数值分析方法意义、数值分析方法意义数学是一种工具，用于解决日常生活、工业工程上的相关问题。针对于数值分析中的数学方法，我们小组将主要内容概括分解，将使用到的方法进行对比分析。2、数值分析主要部分。、数值分析主要部分。1 各类插值方法我们讲过拉格朗日插值、牛顿插值、分段插值、样条插值。2 函数逼近及拟合。3 数值积分、欧拉法解常微分方程、龙格-库塔法解常微分方程、方程组。【1】插值对于牛顿插值相对于拉氏插值增加一个节点，所有的插值基本多项式要重新取、重新算.2 而牛顿插值，节点增加，次数增加，即高次插值函数计算量大，有剧烈震荡，数值稳定性

5、较差（例如龙格现象） ；分段插值在分段点上仅连续（即函数值相等） ，但是有尖点，不光滑（尖点导数不连续） ；样条函数可以解决以上问题：使插值函数既是低次阶分段函数，又是光滑的函数。【2】理解逼近问题与拟合问题：1）逼近问题：函数 f(x)在区间a,b具有一阶光滑度，求多项式 p(x)是f(x)-p(x)在某衡量标准下最小的问题。 2）拟合问题：从理论上讲 y=f(x)是客观存在的，但在实际中，仅仅从一些离散的数据（xi,yi）(i=1,2)是不可能求出 f(x)的准确表达式，只能求出其近似表达式 （x） 。【3】 Romberg（龙贝格）求积法和 Gauss 求积法的基本思想：(主要研究方法)

6、1）复化求积公式精度较高，但需要事先确定步长，欠灵活性，在计算过程中将步长逐次减半得到一个新的序列，用此新序列逼近 I 的算法为 Romberg 求积法。2）对插值型求积公式，若能选取适当的 xk.Ak 使其具有 2n+1 阶代数精度，则称此类求积公式为 Gauss 型。【4】Runge-Kutta 方法的基本思想：借助于 Taylor 级数法的思想，将 yn+1=yn+hy()中的 y()（平均斜率）表示为 f 在若干点处值的线性组合，通过选择适当的系数使公式达到一定的阶。1. 常微分方程初值问题数值求解的优缺点分析Euler 方法计算简单但精度差；向后 Euler 方法与 Euler 方法

7、误差相似；梯形方法比 Euler 方法精度高但算法复杂、计算量很大；改进 Euler 方法结合了 Euler 方法和梯形法的优点；龙格-库塔方法二机械运动的相关简单概念机械运动的相关简单概念1机械运动 在物理学中，把一个物体相对于另一个物体位置的变化称作为机械运动，简称运动。 2参照物 要判断一个物体是否在运动，必须选择另一个物体作为标准，这个作为标准的物体叫做参照物。对于同一个物体的运动，选择的参照物不同，得出的结论也有可能是不同的。 3运动和静止的相对性 自然界中一切物体都在运动，因为地球本身在自转，所以绝对静止的物体是不存在的。通常所描述的物体的运动或静止都是相对于某一个参照物而言的。三

8、机械运动的数值仿真三机械运动的数值仿真现举一例：螺栓组联结的优化设计螺栓组联结的优化设计 如图 1 所示的压力容器螺栓组联接中，已知 D1= 400mm,D2 = 250mm，缸内工作压力为 p=1.5 MPa，螺栓材料为 35 号钢，s=320Mpa,安全系数 S=3,取残余预紧力 Qp=1.6F,采用铜皮石棉密封垫片。现从安全、可靠、经济的角度来选择螺栓的个数 n 和螺栓的直径 d 。3 3 1 1 设设计计问问题题分分析析若从经济性考虑，螺栓数量尽量少些、尺寸小些，但这会使降低联结的强度和密封性，不能保证安全可靠的工作；若从安全、可靠度考虑，螺栓数量应多一些、尺寸大一些为好，显然经济性差

9、，甚至造成安装扳手空间过小，操作困难。为此，该问题的设计思想是：在追求螺栓组联结经济成本最小化的同时，还要保证联结工作安全、可靠。 3 3 2 2 设设计计变变量量 目目标标函函数数 约约束束条条件件 3 3 2 2 . .1 1 设设计计变变量量 选取螺栓的个数 n 和直径 d(mm)为设计变量:T21Txxdn X 3 3 2 2 . .2 2 目标函数目标函数 追求螺栓组联结经济成本 Cn最小为目标。而当螺栓的长度、材料和加工条件一定时，螺栓的总成本与 nd 值成正比，所以本问题优化设计的目标函数为min F(X) = Cn = nd = x1x2 强度约束条件强度约束条件 为了保证安全

10、可靠地工作，螺栓组联结必须满足强度条件 32 . 521dQca ； 其中MpaSs106. 3320； nnpnDFFFFQQp6093742505 . 16 . 246 . 26 . 26 . 1222 N；对于粗牙普通螺纹：由文献3推荐，小径 d1=0.85d 所以，强度约束条件为：0106146192106146192106105624)(2212211xxndndXg 密封密封约束条件约束条件 为了保证密封安全，螺栓间距应小于 10d，所以，密封约束条件为：01040010)(2112xxdnDXg 安装扳手空间安装扳手空间约束条件约束条件 为了保证足够的扳手空间，螺栓间距应大于5d

11、，所以，安装安装约束条件为：040055)(1213xxnDdXg 边界约束条件边界约束条件 0)(14xXg ；0)(25xXg3 3 3 3 3 3 建建立立数数学学模模型型综上所述，本问题的数学模型可表达为：设计变量：T21xxX 目标函数：min F(X) = x1x2约束条件： s.t. 0)(Xgi ( i = 1, 2, 3, 4, 5,)现运用MATLAB的优化函数进行求解 : 先编写M文件 function c,ceq=mynas(x)c(1)=146192/(x(1)*x(2)2)-106; % 非线性不等式约束 c(2)=400*pi/x(1)-10*x(2);c(3)=

12、-400*pi/x(1)+5*x(2);ceq=; % 非线性等式约束在MATLAB命令窗口输入:fun=x(1)*x(2); % 目标函数x0=4,6; % 设计变量初始值A=-1,0;0,-1; % 线性不等式约束矩阵b=0;0; Aeq=; % 线性等式约束矩阵beq=;lb=; % 边界约束矩阵ub=;x,fval=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,mynlsub) % 调用有约束优化函数运行结果如下:x 所以，该问题优化结果为： n =11.4499 ,d =，目标函数最小值：F(X)=。根据实际问题的意义取整、标准化： n =12 ,d = 12。由

13、此例可以看出，与其它编程语言相比， MATLAB语言可以简化编程。图2是调用MATLAB绘图函数自动对上例绘制的数学模型要素图（标注数字的曲线为目标函数的等值线） ，为此在MATLAB命令窗口输入: x1=0.1:20;y1=146192./(106.*x1.2);y2=400.*pi./(10.*x1);y3=400.*pi./(5.*x1);plot(y1,x1,y2,x1,y3,x1,x(1),x(2),o)y4=0.1:0.1:20;y4,x1=meshgrid(y4,x1);Q=y4.*x1;hold on;c,h=contour(y4,x1,Q); hold on;clabel(c,