普通圆柱螺旋弹簧的最优化设计6页

1、设计弹簧时，除选择材料及规定热处理要求外，主要是根据最大工作载荷、最大变形以及结构要求等来确定弹簧的钢丝直径、中径、工作圈数、节距或螺旋升角和高度等，通常取弹簧钢丝直径、中径、工作圈数为最优化设计的设计变量，即 （1）目标函数可根据弹簧的工作特点和对它的专门要求来建立。例如，由于因工作特点极易导致疲劳损坏的弹簧，则应以疲劳安全系数最大作为最优化设计的目标；对于受到高速运转机构变载作用的弹簧，则应以其一阶自振频率最大或最小作为最优化设计的目标，使自振频率值远离载荷变化频率值，以避免共振；对于安装空间很紧、要求尽量减少轮廓尺寸的弹簧，则应以其外径或高度最小，从而得到最小安装尺寸作为最优化设计的目标

2、；当价格成为主要问题时，也可以以弹簧的成本最小作为目标；还有按满应力原则建立目标函数的。对于一般弹簧，通常以质量或钢丝的体积最小作为最优化设计的目标，这时目标函数可表达为： （2）式中，为弹簧钢丝材料的密度，kg/mm3将值及式（1）代入式（2），得以弹簧工作部分（除支撑圈外）的质量为目标的函数表达式： （3）约束条件可根据对弹簧功能的要求和结构限制列出：（1）根据对弹簧刚度的要求范围：（），得约束条件 （4） （5）式中为弹簧材料的剪切弹性模量。（2）根据弹簧钢丝的产品尺寸规格，给出弹簧钢丝直径的限制范围：，从而得约束条件 （6） （7） （3）根据弹簧安装空间对其中径的限制而有 （8） （

3、9）（4）根据对工作圈数的规定范围：而有 （10） （11）（5）根据旋绕比（弹簧指数）的范围：，得 （12） （13）（6）根据弹簧在最大载荷下不碰圈的要求：式中，弹簧自由高度，当支撑圈数且弹簧两端磨平时； 节距，计算时可取； 弹簧在最大工作载荷下的变形量，弹簧并紧高度，当支撑圈数且弹簧两端磨平时，得约束条件： （14）（7）根据弹簧的强度条件： （15）式中，在最大工作载荷作用下或在压并状态下钢丝截面内侧所产生的最大扭转应力；曲度系数：许用扭转应力，视弹簧材料及受载情况而定。对于用碳素弹簧钢丝制成的1类弹簧：；对于2类弹簧：；对于3类弹簧：。进行预压处理后。碳素钢以外材料的值可查有关手册或

4、图书，亦可按上式选取。弹簧钢丝抗拉强度极限，与材料品质及钢丝直径有关，设计手册中列表给出。利用表中数据通过线性回归分析，可建立弹簧钢丝的回归方程，其形式为式中，A，B为线性回归分析所求得的已知常数。将及代入式（15），得 （16）由此得约束条件 （17）（8）根据压缩弹簧的稳定性条件：得约束条件 （18）式中，临界高径比，根据弹簧的支撑方式不同而异；当两端固定时 当一端固定另一端可转动时 当两端均可转动时 （9）根据弹簧的疲劳强度要求给出约束条件当弹簧承受交变载荷并在及之间作交替变化时，在钢丝截面内侧产生的相应扭转应力分别为：对于为常数的弹簧，其疲劳强度安全系数为式中，许用安全系数，当弹簧的设

5、计计算及材料试验数据精确性较高时，取=1.31.7；当精确性较低时，取=1.82.2。弹簧材料的脉动疲劳极限，可按表选取载荷作用次数1041051061070.450.350.330.30将=0.33及代入上式，经整理后得约束条件 （19）（10）根据承受高速交变载荷的弹簧不发生共振的要求，弹簧的自振频率应远离其受载的载荷变化频率，它们之间应满足下列关系：对阀门弹簧 对减振弹簧 （20）弹簧的一阶自振频率(Hz):式中，弹簧刚度：； 弹簧工作部分的质量； 弹簧钢丝材料的密度，kg/mm3。对于圆截面钢丝弹簧，其一阶自振频率为当两端固定时 当一端固定另一端可转动时 当支撑一质量时 式中弹簧钢丝直

6、径、中径的单位均为mm，工作圈数为整数。将有关弹簧的自振频率值代入式（20），则可得约束条件。综合上述各式，可得到以工作部分质量为目标函数的圆柱螺旋压缩弹簧最优化设计的数学模型为：min s.t. 举例：要求：设计一个普通圆柱螺旋拉伸弹簧，要求其体积最小。已知该弹簧在一般载荷条件下工作，并要求外径不超过mm，弹簧的许用应力为MPa。当弹簧拉伸变伸量=7.5mm时，拉力=180N；拉伸变形量mm时，拉力=340N（最大载荷）。1. 确定设计变量取弹簧钢丝直径、中径、工作圈数为最优化设计的设计变量，即 2. 确定目标函数以质量最轻作为追求的目标，建立优化设计目标函数： （2）式中，为弹簧钢丝材料的密度，kg/mm3； 为弹簧支撑圈数，取。根据要求，只需要优化3. 建立约束条件（1）弹簧刚度的约束由已知条件可知，弹簧刚度为N/mm；而(其中)，可得约束条件(2)根据弹簧的强度条件，取最大拉力=340N，得约束条件：（3）根据旋绕比的范围，得约束条件：（4）根据弹簧安装空间的要求，得约束条件：程序编写：（1）所求目标函数：function f=myfun(x)f=x(1)2*x(2)*x(3);（2）约束函数：function c,ceq=nonlcon(x)c(1)=1438.2*x(2)0.84/x(1)2.84-637.5;ceq=7895*x(1)4/(x(2)3*x(3)-16