膜片弹簧的优化设计-MATLAB
基于MATLAB算法的机械优化设计
课程设计报告班级：机0801姓名：徐勤秀学号：081101225指导老师：边义祥日期：2012.2.17基于MATLAB算法的机械优化设计摘要：将MA TLAB算法应用于机械优化设计，提出了MATLAB算法的优化原理及其数学模型的建立，给出求解方法，最后结合实例，求解机械优化设计的最优化问题。
关键词：MATLAB；优化设计；非线性约束最小化1.概论自MathWorks公司1984年推出MA TLAB以来，历经20多年的发展和竞争，MATLAB 语言就成为最具吸引力、应用最为广泛的数值科学计算语言。
随着其功能的不断完善，可以说，MATLAB已成为集数值计算功能、符号计算功能和计算可视化为一身的强大的科学计算语言。
本文运用MA TLAB6.5的优化工具求解机械工程设计中的最优化问题。
在国民经济各部门和科学技术的各个领域中普遍存在着最优化问题，最优化问题就是从所有可能的方案中选择出最合理的、达到最优目标的方案，即最优方案，搜索最优方案的方法就是最优化方法。
将MATLAB运用于最优化方法，使得机械优化设计更趋于科学性，同时MATLAB不用编写复杂的运算程序和各种难于掌握的优化算法，而且通俗易学，从而使优化问题更通俗化。
MATLAB的最优化技术主要包括以下两个方面的内容：（1）建立数学模型。
即用数学方法来描述最优化问题。
模型中的数学关系反映了最优化问题所要达到的目标和各种约束条件。
（2）数学求解。
数学模型建好以后，选择合理的优化方法进行求解。
2.MA TLAB优化算法的几何描述由于机械优化设计多数是非线性约束最小优化问题，通常要将问题转换为更简单的子问题，这些子问题可以求解并作为迭代过程的基础。
早期的方法通常是通过构造惩罚函数等来将有约束的最优化问题转换为无约束最优化问题进行求解。
现在，这些方法已经被更有效的基于K-T（Kuhh-Tucker）方程解的方法所取代。
K-T方程是有约束最优化问题求解的必要条件，是非线性规划算法的基础，这些算法直接计算拉格朗日乘子，通过拟牛顿法更新过程，给K-T方程积累二阶信息，可以保证有约束拟牛顿法的线性收敛。
车辆离合器膜片弹簧的设计与优化
车辆离合器膜片弹簧的设计与优化纲要 :膜片弹簧是汽车离合器的重要零件，是由弹簧钢板冲压而成，形状呈碟形。
膜片弹簧构造紧凑且拥有非线性特征，高速性能好，工作稳固，踏板操作轻便，所以获得宽泛使用。
本文经过对膜片弹簧成立数学模型，特别经过引入加权系数同时对两个目标函数进行比率调理，并用MATLAB编程来优化设计参数。
通过举例，结果证明在压紧力稳固性，分别力及构造尺寸上优化结果较为理想。
重点词 :膜片弹簧；优化设计；MATLAB1.前言1.1 离合器膜片弹簧弹性特征的数学表达式膜片弹簧是汽车离合器中重要的压紧组件，构造比较复杂，内孔圆周表面上有均布的长径向槽，槽根为较大的长圆形或矩形窗孔，这部分称为分别指；从窗孔底部至弹簧外圆周的部分像一个无底宽边碟子，其截面为呈锥形，称之为碟簧。
膜片弹簧的构造如图1-1 所示。
图 1-1膜片弹簧构造表示图图1-2膜片弹簧构造主要参数膜片弹簧主要构造参数如图 2 所示。
R 是自由状态下碟簧部分大端半径。
R 1、r 1分别是压盘加载点和支承环加载点半径，H 是自由状态下碟簧部分的内截锥高度。
膜片弹簧在自由、压紧和分别状态下的变形如图1-3 所示。
图 1-3膜片弹簧在不一样工作状态下的变形膜片弹簧大端的压紧力F1与大端变形量1之间的关系为：E h F16 11ln R / r H1R r1R rh2（ 1）22R1H2 R1r1R1 r1r1式中， r 为自由状态碟簧部分小端半径(mm)； h 为膜片弹簧钢板厚度 (mm)。
明显，膜片弹簧大端的压紧力F1与大端变形量1的函数关系为非线性关系。
由式（ 1）能够看出膜片弹簧大端的压紧力F1分别为 R、 r 、H、h、R1、 r 1等参数相关，故膜片弹簧弹性特征较一般螺旋弹簧要复杂得多。
以某国产小轿车离合器为例，离合器主要性能构造参数为：最大摩擦力矩为700N·m。
从动盘为双片干式，摩擦片外径D=300mm，内径 d=175mm，摩擦因数取0.3 ，膜片弹簧资料为60Si 2MnA，资料弹性模量 E=21000MPa，泊松比μ=0.3 。
基于MATLAB目标函数的建立优化离合器膜片弹簧的设计研究
基于MATLAB目标函数的建立优化离合器膜片弹簧的设计研究摘要：探讨汽车离合器膜片弹簧在已知工作条件下，如何用优化设计方法，选择出一组膜片弹簧的优化结构参数。
使其弹性特性满足离合器的使用性能要求，而且弹簧强度也满足设计要求，以达到最佳的综和效果。
关键词：离合器膜片弹簧、压紧力、MA TLAB、优化设计、Study on Optimization Design of Clutch Spring Based on Objective Function in MATLABAbstract: On the known condition,how to choose a group optimized structure parameters through optimization method for the diaphragm spring of clutch was discussed.For the elastic properties and the spring strength to meet the use of clutch performance requirements,in order to achieve the best effect .Keywords: The clutch diaphragm spring; The pressing force; Optimized design一离合器膜片弹簧目前，汽车广泛采用膜片弹簧作为压紧弹簧的离合器，称为膜片弹簧离合器。
实质是一个用薄弹簧钢板制成的带有一定锥度，中心部分开有许多均布径向槽的圆锥形弹簧片。
二膜片弹簧基本参数1. 比值H/h和h的选择比值H/h对膜片弹簧的弹性特性影响极大。
正确选择该比值，以得到理想的特性曲线及获得最佳使用性能。
一般多取1.5~2.0.2. R/r比值和R、r的选择研究表明，比值越大，弹簧材料利用率越低，弹簧越硬，弹性特性曲线受直径误差的影响越大，切应力越高。
汽车离合器膜片弹簧的优化设计
膜 片 弹簧 在 各种 变 形 情 况 下 ， 蝶簧 部 分 内半 其 径处 应力 为最 大 ， 图 ２所 示 ， 如 Ｂ点 的应 力 总是 大 于
其 他 各点 。设 为切 向 压 缩 应 力 ，， 弯 曲应 力 ， 为
两者 的 最 大 值 在 离 合 器 分 离 过 程 中 互 相 垂 直 。 因
维普资讯
公 路 与 汽 运 ｕ ｏ ｔｖ ｔ ｍｏ ｉ ｅ Ａｐｐ ｉａ ｉ ｎ ｌｃ ｔｏ ｓ
９
ｊ
汽 车 离合 器 膜 片弹 簧 的 优化 设计
肖文 颖 ， 海 华 许
（ 东科 学技 术 职 业 学 院 ，广 东 珠 海 广 ５９９） １ ０ ０
‘ （ 拦 －） ６豢（ Ｐ × １２ — ｒ 一） ） ，Ｌ （ 一 ａ ｅ
『 罄 ） ＿ ・ ） ］ Ｈ （ 每警 Ｈ
图 ２ 膜 片 弹 簧 应 力最 大 点
维普资讯
Ｈ ｉ ｈ ｙ Ａｕ与 ｔ汽 运 ｔｏ ｓ ｇ 公 路 ｔ ｍｏ ｉ ｅ ｗａ ｓ ｏ ｖ Ａｐｐ ｉａ ｉｎ ｌｃ
际意 义。
图 １ 膜 片 弹 簧 小 端 受 载 时 的 变 形 示 震 图
式中： 为 弹性 模量 ， 于 钢 材料 ， Ｅ 对 Ｅ＝ ２ １× １ ． ０ ＭＰ ； 弹簧 厚 度 （ ｍ） 为 大 端 变 形 量 （ ｍ） ａ ｈ为 ａ ｒ ； ａ ｒ ； Ｒ 为碟簧 部 分外 半 径 （ 端 半 径 ， ｍ） ｒ为 碟 簧 部 大 ａ ｒ ； 分 内半 径 （ ｍ） 为泊 松 比 ， 于 钢 材 料 ， ０ ３ ａ ｒ ； 对 一 ． ； Ｌ为膜 片 弹簧与 压盘 接触 半 径 （ ｍ） Ｐ为 支 承环 平 ａ ｒ ； 均半 径 （ ｍ） 为分 离 轴 承 作 用半 径 （ ｍ） Ｈ 为 ａ ｒ ； ａ ｒ ； 碟 簧部 分 内截锥 高度 （ ｍ） ａ ｒ 。
基于MATLAB螺旋拉伸弹簧优化设计2800字
基于MATLAB螺旋拉伸弹簧优化设计2800字摘要：通过MATLAB优化工具箱可以对螺旋拉伸弹簧的相关尺寸得到优化。
结果证明，在满足变形要求和最大剪应力不超过允许值的情况下，使最终目标为弹簧丝体积达到最小。
毕业关键词：MATLAB 螺旋拉伸弹簧优化设计一、前言弹簧是一种通用机械零件，它可以在载荷作用下产生较大的弹性变形。
弹簧性能的好坏对一些机械如内燃机气缸的阀门弹簧及各种缓冲器用的弹簧等是否能正常工作有很大的影响，衡量弹簧优劣的重要指标有体积或质量、速度及刚度，而且每一个目标之间的约束都具有其复杂性，所以很难在同述性能目标和质量目标并求得总体意义上的最优解，早已经成为人们研究与探讨的问题。
文献[1]是根据弹簧的最大载荷、最大变形及结构要求等来决定弹簧直径、弹簧中径、工作圈数、弹簧的螺旋升角和长度等通过运用遗传算法对以弹簧丝体积最小为目标函数的圆柱螺旋拉伸弹簧进行优化设计。
文献[2-4]分别用不同的方法来求解螺旋弹簧的优化模型。
本文是以联合收割机割台仿形机构的拉伸平衡弹簧为例，在满足变形要求和最大剪应力不超过允许值的情况下，以弹簧丝体积最小为目标函数运用MATLAB进行优化设计。
二、螺旋拉伸弹簧数学模型的建立在设计螺旋拉伸弹簧时，通常是根据弹簧的最大载荷、最大变形以及结构要求等来决定弹簧丝直径、弹簧中径、工作圈数、弹簧的螺旋升角和长度等。
即要求弹簧刚度尽可能大，弹簧所用金属材料尽可能少[5]。
设计联合收割机割台仿形机构的拉伸平衡弹簧，弹簧最大的拉力为F=28760N，最大拉力时弹簧的变形f=221.3mm, =700MPa，对于碳钢G=83000N/mm2 。
（一）设计变量的确定除了拉力F已知以外,影响弹簧变形和应力的参数有弹簧平均直径D、弹簧钢丝直径d、弹簧的有效圈数i、弹簧的根数n，因此取设计变量为：（二）目标函数的确定在满足变形要求和最大剪应力不超过允许值情况下,使弹簧的重量(或用料体积最小),本文用弹簧丝的体积作为目标函数,即:引入设计变量x1,x2,x3,x4,整理后可得体积最小的目标函数为：（三）约束条件的确定考虑到变形要求和剪应力及设计变量等的界限，得约束条件如下：1.弹簧的剪应力公式为：（MPa）式中k―与旋绕比有关的系数，可按以下公式计算：因此2.弹簧丝的直径d不大于15mm，即：3.弹簧圈的平均直径D不大于100mm，即：4.弹簧的根数n必须大于或等于1，即：5.螺旋弹簧的变形公式为：式中G―材料的剪弹性系数，F―弹簧承受的载荷（N）；所以-221.3=0三、MATLAB优化工具箱MATLAB是由美国Math Works公司开发的以矩阵运算为基础，集通用数学运算、图形交互、程序设计和系统建模为一体，功能强、使用简单、容易扩展的科技应用软件，分总包和若干工具箱，其中的优化工具箱含有一系列的优化算法函数，机械优化设计把数学规划理论与数值方法应用于设计中,用计算机从大量可行方案中找出最优化设计方案, 从而大大提高设计质量和设计效率。
基于MATLAB的先导溢流阀调节弹簧的优化设计_谌铎文
基于独特型网络理论的车载多传感器信息融合算法
阳明， 唐新蓬， 周建刚 （ 华中科技大学 机械学院， 湖北 武汉 ４３００７４）
摘 要： 在独特型网络的基础上， 设计了一种应用于高速公路上， 装载 ＣＣＤ 和雷达的车辆多传感器信息融合模型。实
时更新的抗体库能提高识别速度， 而抗体库多样性则通过抗体间正交亲和度阈值来保证。独特型网络中抗体浓度随着
弹簧中径、有效圈数） ， 同时对它们进行优化设计， 即平行
设计， 从而获得一组优化的结构参数。与以往的从一个变
量的假设值开始计算其它变量， 然后代入检验、修正计
算 、再 检 验 的 设 计 过 程相 比 ， 平 行 设 计 具 有 设 计 速 度 快 、
设计精度高、节省人力、物力的优点。ＭＡＴＬＡＢ 具有强大
７２ 机械工程师 ２００６ 年第 ３ 期
Ｍｅ ｃ ｈ ａ ｎ ｉｃ ａ ｌ Ｄｅ ｓ ｉｇ ｎ
机械设计
参数。故在此选择 ｄ、Ｄ２、ｎ 为设计变量。 ｘ＝（ ｘ１ｘ２ｘ３） Ｔ＝（ ｄＤ２ｎ） Ｔ
于是， 目标函数可表示为
! Ｊ＝ １ ４π
Ｇ ２ρ
ｘ１ｘ－２２
ｘ－ １ ３
（ ６）
２．３ 约束条件 ①根据弹簧钢丝的产品尺寸规格， 给出弹簧钢丝直
ｘ－ １ ３
%
&&ｘ＝（ ｘ１ｘ２ｘ３） Ｔ＝（ ｄＤ２ｎ） Ｔ
&
’&&ｓ．ｔ． ｇ（ａ ｘ） ≤０ ａ＝１， ２…， １０
３ 优化函数选择
在 ＭＡＴＬＡＢ 优化工具箱中， 最小优化函数有多 个。
因为先导式溢流阀调节弹簧优化问题是有约束的非线性
函数的最小值问题， 所以选用 ｆｍｉｎｃｏｎ 函数。ｆｍｉｎｃｏｎ 函数
离合器膜片弹簧优化编程
基于MATLAB的膜片弹簧优化程序mopian.m文件：clear allclc%膜片弹簧结构参数值El=210000;%材料弹性模量miu=0.3;%泊松比ds=1.6;%磨损极限dt=2;%推力行程D=225;d=150;r0=37;%结构参数膜片弹簧小端内半径rf=38;%分离作用半径x0=[4.4 2.4 109 85.5 104.6 90 2.4];%分别是H h R r R1 r1 lamdax7=0:0.1:6;aa=pi*El.*x0(2).*x7/(6*(1-miu^2));bb=log(x0(3)./x0(4))./((x0(5)-x0(6)).^2);cc=x0(1)-x7.*(x0(3)-x0(4))./(x0(5)-x0(6));dd=x0(1)-0.5*x7.*(x0(3)-x0(4))./(x0(5)-x0(6));ee=x0(2).^2;F=aa.*bb.*(cc.*dd+ee);plot(x7,F,'b')%绘制原始膜片弹簧弹性特性图hold on%设计变量的上下界Lb=[2.5 1 101 82 100 82 2.5];%设计变量下限Ub=[6 4 109 97 110 97 4.5];%设计变量上限%线性不等式约束系数矩阵和常数% H h R r R1 r1 lamdaA=[1 -2.2 0 0 0 0 0;-1 1.7 0 0 0 0 0;1 0 -pi/15 pi/15 0 0 0;-1 0 pi/20 -pi/20 0 0 0;0 0 1 -1.35 0 0 0;0 0 -1 1.2 0 0 0;0 -50 1 0 0 0 0;0 35 -1 0 0 0 00 0 0 0 1 0 0;0 0 0 0 -1 0 0;0 0 1 0 -1 0 0;0 0 -1 0 1 0 0;0 0 0 -1 0 1 0;0 0 0 1 0 -1 0;];b=[0 0 0 0 0 0 0 0 D./2 -(D+d)./4 7 -1 6 0];%线性等式约束Aeq=[];beq=[];options=optimset('largescale','off','display','iter')[x,fval,exitflag,out]=fmincon(@objfun,x0,A,b,Aeq,beq,Lb,Ub,@confun,options);[c]=confun(x);x7=0:0.1:6;aa=pi*El.*x(2).*x7/(6*(1-miu.^2));bb=log(x(3)./x(4))./((x(5)-x(6)).^2);cc=x(1)-x7.*(x(3)-x(4))./(x(5)-x(6));dd=x(1)-0.5.*x7.*(x(3)-x(4))./(x(5)-x(6));ee=x(2).^2;F=aa.*bb.*(cc.*dd+ee);plot(x7,F,'r--')%绘制优化后膜片弹簧弹性特性图xconfun.m文件：%%%建立非线性约束条件function[c,ceq]=confun(x)miu=0.3;aa=pi*210000.*x(2).*x(7)/(6*(1-miu^2));bb=log(x(3)./x(4))./((x(5)-x(6)).^2);cc=x(1)-x(7).*(x(3)-x(4))./(x(5)-x(6));dd=x(1)-0.5.*x(7).*(x(3)-x(4))./(x(5)-x(6));ee=x(2).^2;kk=210000/((1-miu.^2).*x(4));e=(x(3)-x(4))./(log(x(3)./x(4)));%中性点半径tt=0.5.*(e-x(4));alfa=atan(x(1)./(x(3)-x(4)));%膜片弹簧锥形底角的计算fa=alfa+0.5.*x(2)./(e-x(4));%切向压应力达到最大值时的膜片转角thegatb=abs(kk.*(tt.*fa.^2-(2.*tt.*alfa+x(2)./2).*fa));%膜片弹簧危险部位切向压应力计算%%%%%%%%%%%%%%%rf=40;ff=(x(5)-x(6))./(x(6)-rf);%ff=(R1-r1)/(r1-rf)F2/F1力的比值fff=aa.*bb.*(cc.*dd+ee)*ff;%F2%%%%%%%%%%%%%%%n=18;%分离指的数目b=10;%分离把根部的宽度thegarb=abs(6*(x(4)-rf)*fff./(x(2).^2.*n.*b));%膜片弹簧危险部位弯曲应力的计算%%%%%%%%%%%%%%%T=246*1000;%离合器所要传递的最大转矩，单位化成N·mma=225/2;b=150/2;z=2;fz=0.28;c(1)=(thegarb-thegatb)-1500;%膜片弹簧危险点最大当量应力约束，非线性不等式1c(2)=T/(z.*fz.*(2./3).*(a.^3-b.^3)/(a.^2-b.^2))-(aa.*bb.*(cc.*dd+ee));%膜片弹簧产生压紧力的约束，非线性不等式2c(3)=(aa.*bb.*(cc.*dd+ee))-1.3*T/(z.*fz.*(2./3).*(a.^3-b.^3)/(a.^2-b.^2));%beta<1.75c(4)=1.4*T/(z.*fz.*(2./3).*(a.^3-b.^3)/(a.^2-b.^2))-(aa.*bb.*(cc.*dd+ee));%beta>1.2ceq=[];beta=(aa.*bb.*(cc.*dd+ee))/(T/(z.*fz.*(2./3).*(a.^3-b.^3)/(a.^2-b.^2)))objfun.m文件:function f=objfun(x)El=210000;%材料弹性模量miu=0.3;%泊松比rf=40;%分离轴承推力作用半径%%将弹性特性公式分成aa、bb、cc、dd、ee五部分表示aa=pi*El.*x(2).*x(7)/(6*(1-miu.^2));bb=log(x(3)./x(4))./((x(5)-x(6)).^2);cc=x(1)-x(7).*(x(3)-x(4))./(x(5)-x(6));dd=x(1)-0.5.*x(7).*(x(3)-x(4))./(x(5)-x(6));ee=x(2).^2;ff=(x(5)-x(6))./(x(6)-rf);%ff=(R1-r1)/(r1-rf)%%%%磨损后的公式参数变化ds=1.6;%磨损极限在1.6-2.2之间，取1.6mmaa1=pi*El.*x(2).*(x(7)-ds)/(6*(1-miu^2));bb=log(x(3)./x(4))./((x(5)-x(6)).^2);cc1=x(1)-(x(7)-ds).*(x(3)-x(4))./(x(5)-x(6));dd1=x(1)-0.5.*(x(7)-ds).*(x(3)-x(4))./(x(5)-x(6));ee=x(2).^2;%%%%分离过程公式参数变化dt=2;%分离行程取值为2mmaa2=pi*El.*x(2).*(x(7)+dt)/(6*(1-miu^2));bb2=log(x(3)./x(4))./((x(5)-x(6)).*(x(6)-rf));cc2=x(1)-(x(7)+dt).*(x(3)-x(4))./(x(5)-x(6));dd2=x(1)-0.5.*(x(7)+dt).*(x(3)-x(4))./(x(5)-x(6));ee=x(2).^2;%%%%双目标函数表达式f1=abs(aa.*bb.*(cc.*dd+ee)-aa1.*bb.*(cc1.*dd1+ee));%第一个目标函数：磨损极限内正压力的变化值f2=aa2.*bb2.*(cc2.*dd2+ee);%第二个目标函数：膜片弹簧在分离位置时的弹力fac=0.7;%加权因子f=fac.*f1+(1-fac)*f2;%总体目标函数优化结果曲线。
基于MATLAB的膜片弹簧参数对性能影响及优化
２０２０ 年 ７ 月
吉 林 化
工
学
院
学
Ｖｏｌ.３７ Ｎｏ.７
Ｊｕｌ. ２０２０
报
ＪＯＵＲＮＡＬ ＯＦ ＪＩＬＩＮ ＩＮＳＴＩＴＵＴＥ ＯＦ ＣＨＥＭＩＣＡＬ ＴＥＣＨＮＯＬＯＧＹ
文章编号:１００７￣２８５３(２０２０)０７￣００３２￣０４
基于 ＭＡＴＬＡＢ 的膜片弹簧参数对性能影响及优化
图 ２ 膜片弹簧弹性特性曲线
值时 λ １Ｂ 已经超过了其它两条曲线中的彻底分离
２ 参数对膜片弹簧性能的影响分析
自由状态膜片弹簧的尺寸参数如图 １ 中( ａ)
所示ꎬ其中 Ｒ 和 ｒ 分别为此状态下膜片弹簧碟簧
处大、小端的半径ꎬｈ 代表膜片弹簧的钢板厚度.
当膜片弹簧受到压紧力时ꎬ根据其受力大小的变
化ꎬ在子午断面绕中性点 Ｏ 转动 [３] .膜片弹簧受到
压盘和支撑环的共同作用ꎬ其加载点至圆心半径
分别用 Ｒ １ 和 ｒ １ 表示ꎬ如图 １( ｂ) 所示ꎬＨ 代表自由
Байду номын сангаас
３３
所受压紧力 Ｆ １ 的增大ꎬ离合器主从部件接触更为
紧密、工作更加稳定.同时 Ｈ 取极大值时膜片弹簧
磨损极限变形量 Δλ 较极小值时增大近一倍ꎬ膜
片弹簧的使用寿命延长. 因此可适当增加自由状
态膜片弹簧的内截锥高度 Ｈ 数值ꎬ实际生产中也
常增大碟簧圆锥底角 α [８] .
膜片弹簧钢板厚度 ｈ 为极小值时ꎬ膜片弹簧
１.０) λ １Ｈ
[７]
.λ １Ａ 为膜片弹簧的最大磨损形变的极
限值ꎬλ １Ｃ 为膜片弹簧工作分离形变. 膜片弹簧在
受到支撑环荷载 Ｆ １ 工作时其形变量由 λ １Ｂ 向 λ １Ｃ
弹簧摆matlab演示及源代码
弹簧摆写出弹簧摆的运动微分方程，并用matlab 编写程序，演示弹簧摆的运动弹簧摆的运动微分方程可以通过拉格朗日函数法求得，对于弹簧单摆系统，在极坐标系中其拉格朗日函数为222211 ()cos ()22L T Vmg m r r mgr k r l kθθ=-=++--+ 其中,r θ为自由度，,,,m k g l 分别为摆球质量，弹簧刚度，重力加速度和摆初始长。
将拉格朗日函数带入到拉格朗日方程00d L L dL L dt θθ∂∂⎧-=⎪⎪∂∂⎨∂∂⎪-=⎪∂∂⎩得到系统的运动微分方程：sin ()02sin 0k mg r r g r l g k r rg θθθθθ⎧--+-+=⎪⎨⎪++=⎩ 对于这个非线性方程组，利用matlab 里面自带的ode45常微分方程组求解函数，对上面的方程组进行数值求解，得到了摆球的运动轨迹，如下图所示0246810121416弹簧摆内共振动画模拟及自由度随时间的变化曲线时间t 单位秒对应的自由度图中给出了弹簧摆运动过程中自由度,r 随时间的变化关系，下面将给出这个程序的源代码，演示程序时，首先停顿三秒，然后开始摆动，图形上方的弹簧摆的示意图开始摆动并描绘出轨迹，同时图形下方自由度随时间变化关系也会动态的给出，从这个过程中可以明显看出弹簧摆的内共振现象。
程序源码：function springmasstheta0=1.5*pi/10;%单摆的初始角度m=20/9.8;k=80;g=9.8;%弹簧质量的参数保证弹簧固有频率是单摆固有频率的2倍L0=1;L=L0+m*g/k;%L0为弹簧原来长度，L为弹簧静止时长度[t,u1]=ode45(@weifen,[0:0.05:15],[L00theta00],[],L,k,m,g);[y1,x1]=pol2cart(u1(:,3),u1(:,1));y1=-y1;%将极坐标换为直角坐标axis([016-1.54]);%text(0,u1(1,1),'摆长l');%text(0,u1(1,3),'摆角\theta');%axis offtitle('弹簧摆内共振动画模拟及自由度随时间的变化曲线','fontsize',15)xlabel('时间t单位秒');ylabel('对应的自由度');hold on;R=0.055;%设置弹簧半径yy=-L0:0.01:0;xx=R*sin(yy./L0*30*pi);%用正弦曲线表示弹簧[a,r]=cart2pol(xx,yy);%用坐标变换来画初始位置的弹簧a=a+theta0;[xx,yy]=pol2cart(a,r);%弹簧的数据xx=2*xx+7;yy=yy+3.5;%xx(1)=x(1)+15;line([68],[3.53.5],'color','k','linewidth',4)%弹簧单摆系统初始状态for i=1:9line([6+0.2*i6+0.2*i+0.2],[3.53.7],'color','k','linewidth',1);endline([77],[3.52.3],'color','k','linestyle','-.','linewidth',1)bigball=line(xx(1),yy(1),'color','b','marker','.',...'markersize',30,'erasemode','xor');%球ball2=line(xx(1),yy(1),'color','g','linestyle','-',...'linewidth',1.3,'erasemode','none');%轨线spring=line(xx,yy,'color','k','linewidth',2,'erasemode','xor');%弹簧linex=line(t(1),u1(1,1),'color','r','marker','.',...'markersize',10,'erasemode','none');%摆长曲线初始liney=line(t(1),u1(1,3),'color','b','marker','.',...'markersize',10,'erasemode','none');%摆角曲线初始%以下为标注text(6,-0.75,'\uparrow');text(4,-1,'摆角随时间变化曲线');text(9.0,1.9,'摆长随时间变化曲线');text(11.4,1.70,'\downarrow');%在初始状态下停顿三秒pause(2)%以下程序开始实现动画，在弹簧摆开始振动时画出相应的变化曲线for i=2:length(t)yy=-u1(i,1):0.01:0;%弹簧xx=R*sin(yy./u1(i,1)*30*pi);%xx(i)=x(i)+15;[a,r]=cart2pol(xx,yy);a=a+u1(i,3);%set(linex,'XData',t(i),'YData',u1(i,1));%%摆长初始%set(ball2,'XData',2*x1(i)+7,'YData',y1(i)+3.5);%轨线%set(liney,'XData',t(i),'YData',u1(i,3));%摆角初始plot([t(i-1),t(i)],[u1((i-1),1),u1(i,1)],'xr');plot([t(i-1),t(i)],[u1((i-1),3),u1(i,3)],'*b');%legend('摆长曲线','摆角曲线')[xx,yy]=pol2cart(a,r);xx=2*xx+7;yy=yy+3.5;set(bigball,'XData',2*x1(i)+7,'YData',y1(i)+3.5);%球%plot(2*x1(i)+7,y1(i)+3.5,'b');%球%set(ball2,'XData',2*x1(i)+7,'YData',y1(i)+3.5);%轨线plot([2*x1(i-1)+7,2*x1(i)+7],[y1(i-1)+3.5,y1(i)+3.5],'-g');%轨线set(spring,'XData',xx,'YData',yy);%弹簧drawnow;end%以下为弹簧摆的运动微分方程function value=weifen(t,u,l,k,m,g)value=[u(2);u(1.)*u(4).^2+g*cos(u(3))-k/m*(u(1)-l+m*g/k);u(4);-2.*u(2).*u(4)./u(1)-g*sin(u(3))./u(1)];。
弹簧阻尼系统matlab课程设计
弹簧阻尼系统matlab课程设计一、课程目标知识目标：1. 理解弹簧阻尼系统的基本原理，掌握其数学模型；2. 学会使用MATLAB软件进行弹簧阻尼系统的建模与仿真；3. 掌握分析弹簧阻尼系统动态特性的方法。
技能目标：1. 能够运用MATLAB软件构建弹簧阻尼系统的数学模型；2. 能够运用MATLAB进行时域和频域分析，绘制系统响应曲线；3. 能够根据系统响应曲线，分析系统稳定性和性能指标。
情感态度价值观目标：1. 培养学生对工程问题的探究兴趣，提高解决实际问题的能力；2. 培养学生的团队协作精神，提升沟通与表达能力；3. 增强学生的创新意识，使其认识到科技对社会发展的推动作用。
课程性质：本课程为实践性较强的课程设计，旨在通过MATLAB软件的应用，让学生深入理解弹簧阻尼系统的理论知识，并能将其应用于实际问题。
学生特点：学生具备一定的数学基础和物理知识，对MATLAB软件有一定了解，具有较强的学习能力和动手能力。
教学要求：结合学生特点，将课程目标分解为具体的学习成果，注重实践操作，强调理论联系实际，提高学生的实际问题解决能力。
同时，注重培养学生的团队协作能力和创新精神。
在教学过程中，关注学生的个体差异，因材施教，确保课程目标的实现。
二、教学内容1. 弹簧阻尼系统基本原理回顾：包括弹簧、阻尼器的物理特性，系统的自由体图和受力分析；2. 弹簧阻尼系统的数学模型：介绍微分方程的建立，状态空间方程的表示；3. MATLAB软件入门：复习MATLAB基本操作，介绍MATLAB/Simulink环境下进行系统建模的方法；4. 弹簧阻尼系统建模与仿真：运用MATLAB/Simulink构建系统模型，进行时域和频域分析；5. 系统性能分析：通过绘制响应曲线，分析系统的稳定性、超调量、调节时间等性能指标；6. 实践操作与团队协作：分组进行课程设计，每组完成一个弹簧阻尼系统的建模、仿真与分析；7. 成果展示与交流：每组汇报课程设计成果，分享经验，进行交流与讨论。
基于MATLAB的呼气活门弹簧优化设计
基于MATLAB的呼气活门弹簧优化设计摘要：对呼气活门弹簧采用传统方法设计，时间长，工作量大，而且设计结果不一定是最优的情况。
建立呼气活门弹簧优化设计的数学模型，将弹簧的基本结构参数作为设计变量，轻量化作为目标函数，强度、刚度、空间尺寸等要求作为设计的约束条件，运用MATLAB编制弹簧优化设计程序，方便快速计算出符合设计要求且重量最轻的弹簧。
关键词：优化设计；MATLAB；弹簧；呼气活门前言21世纪航空技术得到迅猛发展，随着空中加油机、机载产氧技术的应用，飞机的续航能力得到大大增强，飞行员的飞行时间也随之增加，使得飞行员需要更长时间佩戴供氧面罩，因此飞行员对面罩性能提出了更高的要求。
供氧面罩的呼气阻力是标志供氧面罩性能高低的重要指标之一。
随着呼吸阻力上升，可导致飞行员呼吸肌疲劳，严重的甚至导致过度呼气发生而导致飞行事故。
而呼气阻力的大小与呼气活门弹簧（见图1所示）设计的性能息息相关，呼气活门性能高低，直接关系到整个面罩设计成败与否，因此呼气活门弹簧设计是整个呼气活门的重中之重。
而随着优化设计的兴起，设计出各项性能指标最优的呼气活门弹簧成为可能，因此对呼气活门弹簧采用优化设计可以在满足功能的前提下，对弹簧参数进一步优化，减轻重量达到降低呼气阻力的目的。
在进行实际问题的优化设计时，首先需要建立其数学模型。
优化设计的数学模型需要用设计变量、设计约束和目标函数等概念描述。
图1呼气活门弹簧示意图1.呼气活门圆柱压缩弹簧设计的一般要求1.1结构和几何尺寸要求呼气活门圆柱压缩弹簧由导向轴引导并放在限定的孔径的呼气活门座内，故其弹簧必须能装人直径为Dh的孔中和套在直径为Dn的销轴上，且都应有一定间隙。
自由高度、压缩高度、装配状态和工作状态的变形量、两端磨平状况、最少工作圈数、工作状态下圈与圈间的间隙和能合理制造的弹簧指数等都应满足要求。
1.2性能要求能承受最大工作载荷和具有合适刚度。
1.3强度及寿命要求压缩高度下的应力不应大于[]；并且呼气活门弹簧在受交变载荷时在规定的工作次数下不能产生疲劳损坏等。
基于MATLAB目标函数的建立优化离合器膜片弹簧的设计研究
张 昆 晓 ，何 海 明 ，吴 东盛 （ 广州市公用事业技 师学院，广州 ５０ ０ ） １１０
摘 要 ：探讨汽车离合器膜片弹簧在已知工作条件下 ，如何用优化设 计方法 ，选择 出一 组膜片弹簧 的优 化结构参数 。选 择在摩
ｅ ｇ ０１ ｙ ｈ ｏ ｃ ｒ．ｃ ｎ ００ ＠ ａ ｏ ． ｏ ｎ ｎ。
垫！： 垒 竺 ！ 塑 Ｏ７ 堕 ３
研 究与 开发
Ｈ ｈ对 弹 簧 特 性 的 影 响 ，正 确 地 选 择 该 比 值 ，以 得 到 理 想 的 特 ／
量接近 ｒ ；Ｒ． 略 小 于 且 尽 量 接 近 于 。 汽 车 离 合 器 膜 片 弹 簧 应
性曲线及获得最佳 的使用性能。一般汽车的膜片弹簧离合器 多
ＺＨＡＮＧ ｎ ｉ ｏ． ＨＥ ａｍｉ ｇ． Ｗ Ｕ ｎ ｓ ｅ ｇ Ｋｕ ｘａ Ｈ ｉ ｎ Ｄｏ ｇ ｈ ｎ
（ ｕ ｎｚｏ ｕ ｌ ｎ ｅｈ ｉ ａ ｏ ｅｅ ｕｎ ｚｏ １ １０ ｈｎ ） Ｇ ａ ｇｈ ｕＰ ｂｉ ａｄＴ ｃ ｎｃ ｎＣ ｌ ｇ ，Ｇ ａｇｈｕ５ ０ ０ ，Ｃ ｉａ ｃ ｉ ｌ
ａ ｓ ｌ ｔ ａ ｕ ｏ ｐ ｎ ｒ ｓｕ ｅ ｉ ｈ ｔ ｂ ｓ ｆ ａ ｅ ｃ ｏ ｓ ｅ ｓ ｎ ｂ ｅ ｃ ｎ ｔ ｉｔｃ ｎ ｉ ｏ ｓ ｒ ｇ ａ ｂ ａ ｇ ａ ｅ ａ ｄ ｑ ｉ ｋｙ ｇ ｔ ｂ ｏ ｕｅ ｖ ｌ ｅ ｆｒｓ ｒ ｇ ｐ ｅ ｓ ｒ ｎ ｔ ｅ Ｍａｌ ｏｔ ｒ ， ｈ ｏ ｅ ａ ｒａ ｏ ａ ｌ ｏ ｓｒ ｎ ｏ ｄ ｔ ｎ ，ｐ ｏ ｒ ｍ ｙ Ｃ ｌ ｎ ｕ ｇ ， ｎ ｕ ｃ ｌ ｅ ｉ ａ ｗ ａ ｉ ｔ ｅ ｉｅ ｌｏ ｔ ｚ ｔ ｎ ｒ ｓ ｌ ． ｈ ｄ ａ ｐｉ ａｉ ｅ ｕ ｔ ｍｉ ｏ ｓ Ｋｅ ｗｏ ｄ ：Ｔ ｅ ｃｕ ｃ ｉｐ ｒ ｇ ｐ ｎ ； Ｃ ｍｐ ｃｉｎ ｆｒ ｅ ｃ ａ ｇ ａｕ ；Ｗ ｅ ｒｒ ｎ ｅ Ｃ ｐ ｏ ｒ ｍｍｉ ｇ ｌ ｎ ａ ｅ Ｏ ｔ ｚ ｄ ｄ ｓｇ ｙ ｒ ｓ ｈ ｌ ｔｈ ｄａ ｈ ａ ｍ ｓ ｒ ｇ ｉ ｏ ａ ｔ ｏ ｃ ｈ ｎ ｅ ｖ ｌ ｅ ｏ ａ ａ ｇ ； ｒｇ ａ ｎ ａ ｇ ｇ ； ｐ ｉ ｅ ｅ ｉｎ ｕ ｍｉ
MATLAB案例一圆柱螺旋弹簧的优化设计说明书
g(1)=1.66*(8*Fmax/pi)*(x(1)/x(2))^0.16*x(2)/x(1)^3*1e6-tao g(2)=G*x(1)^4/(8*x(2)^3*x(3))*1e-3-24*1000 %刚度约束 g(3)=x(2)-14*x(1); %旋绕比约束 g(4)=4*x(1)-x(2); %旋绕比约束 g(5)=x(1)+x(2)-42; %安装空间
弹簧设计需要满足刚度、强度、稳定性、共振 性等条件
例一调压弹簧为普通圆柱螺旋压缩弹簧。阀腔直径为 42mm，弹簧最大工作压力为Fmax=1110N，弹簧的许用应
力[]=665MPa,弹簧的最大刚度kmax=24N/mm,弹簧结构要
求：工作圈数n≥0，弹簧指数4≤ C ≤ 14，弹簧压并高度 λb=1.1h=18.25mm。试在满足弹簧的强度条件、刚度条件、 稳定性条件、旋绕比条件和结构边界等约束条件下，确定 弹簧的簧丝直径d、中径D2和工作圈数n等三个设计参数， 使它的结构重量最轻。
Spring_f,x0,[],[],[],[],[],[],@ Spring_g,options); disp '*********弹簧优化设计最优解**********' fprintf(1, '簧丝直径 d=%3.4f\n',x(1)); fprintf(1, '弹簧中径 D=%3.4f\n',x(2)); fprintf(1, '弹簧圈数 n=%3.4f\n',x(3));
%体积目标函数 function f=Spring_f(x) f=1/4*pi*pi*x(1)^2*x(2)*x(3); %x(1)直径,x(2)中径,x(3)圈数
基于Matlab的超越弹簧离合器弹簧的优化设计
将各个约束条件中的参数转换成相应的设计变量 ， 简 化后的约束条件为： 1 ） G1 （ X） = 2 － X1 / X2 ≤0 2 ） G2 （ X） = X1 / X2 － 10 ≤0 3 ） G3 （ X） = 1 － X2 / X3 ≤0 4 ） G4 （ X） = X2 / X3 － 2 ≤0 X1 5 ） jc = － X3 －1 3 jc 2 X1 / X3 + 1 ln X1 / X3 － 1
弹簧本身的设计是弹簧离合器正常运行的最核心的要 素。如果采用恒定截面的话， 那么满足要求的弹簧的长度将 变得非常长。实际上， 把弹簧线圈宽度设计成按指数规律变 即中间宽、 两端窄， 便可以同时满足强度和柔性要求。理 化， 如果壳体施加于弹簧的扭矩呈指数规律变化， 那么线 论上， 圈宽度具有按指数规律变化的弹簧应呈现等应力分布。用 纵向剖视图表示的典型弹簧离合器结构如图 2 所示。 弹簧的两端与输入输出壳体的孔为过盈配合 。这样， 每当输入壳体试图比输出壳体旋转得快时 ， 就可以使离合 器立即接合。弹簧中部线圈在外径上与输入 、 输出壳体是 分离的， 以便使超越时打滑的弹簧的长度减至最少值 ， 从
基金项目：南京航空航天大学基本科研业务费专项科研项目 （ NP2011014 ） 作者简介：常斌（ 1987 —） ， 男， 江苏盐城人， 硕士研究生， 主要研究机械 CAD 及自动化。
Machine Building
基于MATLAB实现遗传算法的膜片弹簧优化设计
·４８·
交通 工程 技 术
２０１６年第２期 总第４３卷
基于 ＭＡＴＬＡＢ实现遗传 算法 的 膜 片 弹 簧 优 化 设 计
陈淑鑫 ，一，程兴旺
基于MATLAB的先导溢流阀调节弹簧的优化设计
基于MATLAB的先导溢流阀调节弹簧的优化设计
谌铎文;胡燕平
【期刊名称】《机械工程师》
【年(卷),期】2006(000)003
【摘要】以弹簧的3个独立结构参数(钢丝直径、弹簧中径、有效圈数)为设计变量,以先导溢流阀调节弹簧的固有频率最大为目标函数建立数学模型.运用MATLAB 的优化函数对其进行优化计算,从而获得一组优化的结构参数,其优化结果表明这是有效的弹簧设计方法.
【总页数】3页(P72-74)
【作者】谌铎文;胡燕平
【作者单位】湖南科技大学,机电工程学院,湖南,湘潭,411201;湖南科技大学,机电工程学院,湖南,湘潭,411201
【正文语种】中文
【中图分类】TH122
【相关文献】
1.基于MATLAB的先导式溢流阀动态特性分析 [J], 付娟娟;张磊;刘小宁
2.基于神经网络溢流阀调压弹簧模糊概率的优化设计 [J], 谢文;林国湘;叶进宝;耿兴春
3.基于Matlab／Simulink的先导式溢流阀研究 [J], 蒲昌顺;黄星德;谭宗柒
4.基于Matlab/Simulink的先导式溢流阀研究 [J], 蒲昌顺;黄星德;谭宗柒
5.溢流阀调节弹簧遗传优化设计 [J], 傅晓锦
因版权原因，仅展示原文概要，查看原文内容请购买。