发明名称 一种基于机械加工过程的动力吸振器参数调谐优化方法
摘要 本发明公开了一种基于机械加工过程的动力吸振器参数的调谐优化方法,该方法包括以下步骤:首先建立振动系统两自由度的数学模型,然后求出动力学微分方程稳态解的幅值;将主结构和动力吸振器的相对位移带入粘性阻尼消耗能量公式;绘制能量曲线,然后统计单一因素变化对能量峰值的关系;根据统计数据拟合出对应各个质量比的最优频率比和最优阻尼比;根据外界干扰信号的频率范围对动力吸振器能够吸收的能量进行积分;然后绘制能量柱状图找到最大值,进而得到最优质量比;最后得到动力吸振器的各个最优参数。本发明对外界干扰引起的振动进行抑制可以使响应速度提高大约50%。
申请公布号 CN103455728B 申请公布日期 2017.01.18
申请号 CN201310420434.2 申请日期 2013.09.13
申请人 天津大学 发明人 陈东;林滨;张彦斌;李新月
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 李素兰
主权项 一种基于机械加工过程的动力吸振器参数的调谐优化方法,其特征在于,该方法包括以下步骤:步骤一、建立振动系统两自由度的数学模型,求出两自由度振动系统的运动微分方程,即<maths num="0001"><math><![CDATA[<mrow><msub><mi>M</mi><mi>m</mi></msub><msub><mover><mi>x</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mi>m</mi></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>.</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>.</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>F</mi><mn>0</mn></msub><mi>sin</mi><mi>&omega;t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000011.GIF" wi="1635" he="75" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>M</mi><mi>a</mi></msub><msub><mover><mi>x</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mn>2</mn></msub><mo>+</mo><msub><mi>K</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>.</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>.</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000012.GIF" wi="1626" he="75" /></maths>式中:M、K和C分别表示质量、刚度和阻尼,脚标m和a分别代表主结构和吸振器;x<sub>1</sub>和x<sub>2</sub>分别表示主结构和吸振器的振动位移,F<sub>0</sub>和ω分别表示激振力幅值和频率;求出动力学微分方程稳态解的幅值,包括主结构和吸振器的位移幅值,并进行无量纲化,得到<maths num="0003"><math><![CDATA[<mrow><mfrac><msub><mi>X</mi><mn>1</mn></msub><msub><mi>&delta;</mi><mi>st</mi></msub></mfrac><mo>=</mo><mfrac><mrow><msup><mi>f</mi><mn>2</mn></msup><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&zeta;g</mi></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>f</mi><mn>2</mn></msup><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><msup><mi>&mu;f</mi><mn>2</mn></msup><msup><mi>g</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&zeta;g</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&mu;g</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000013.GIF" wi="1651" he="154" /></maths><maths num="0004"><math><![CDATA[<mrow><mfrac><msub><mi>X</mi><mn>2</mn></msub><msub><mi>&delta;</mi><mi>st</mi></msub></mfrac><mo>=</mo><mfrac><mrow><msup><mi>f</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mo>&bull;</mo><mn>2</mn><mi>&zeta;g</mi></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>f</mi><mn>2</mn></msup><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><msup><mi>&mu;f</mi><mn>2</mn></msup><msup><mi>g</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mo>&bull;</mo><mn>2</mn><mi>&zeta;g</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&mu;g</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000014.GIF" wi="1634" he="147" /></maths>式中:X<sub>1</sub>和X<sub>2</sub>分别表示主结构和吸振器的位移幅值,δ<sub>st</sub>是主结构的静位移,μ、f和ζ分别是吸振器与主结构的质量比、频率比和阻尼比,g是激振力频率与主结构频率的比值;步骤二、将主结构和动力吸振器的相对位移代入粘性阻尼消耗能量公式,求出动力吸振器粘性阻尼的能量吸收方程,具体处理包括:假设简谐运动为x(t)=Xsinω<sub>d</sub>t,其中X为该运动的振幅,ω<sub>d</sub>为阻尼振动频率,根据粘性阻尼吸收能量原理,系统在一个振动周期内消耗的能量ΔW为:<maths num="0005"><math><![CDATA[<mrow><mi>&Delta;W</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>&pi;</mi><mo>/</mo><msub><mi>&omega;</mi><mi>d</mi></msub></mrow></msubsup><mi>c</mi><msup><mrow><mo>(</mo><mfrac><mi>dx</mi><mi>dt</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mi>dt</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><mi>c</mi><msup><mi>X</mi><mn>2</mn></msup><msub><mi>&omega;</mi><mi>d</mi></msub><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&omega;</mi><mi>d</mi></msub><mi>td</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>d</mi></msub><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>c&pi;&omega;</mi><mi>d</mi></msub><msup><mi>X</mi><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000015.GIF" wi="1619" he="154" /></maths>主结构在受到简谐激振力的作用下发生振动,而附着在上面的动力吸振器可以起到抑制作用;根据公式(3)和公式(4)得到动力吸振器的相对振幅如下:<maths num="0006"><math><![CDATA[<mrow><mi>&Delta;X</mi><mo>=</mo><msub><mi>X</mi><mn>2</mn></msub><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mo>=</mo><mfrac><msup><mi>g</mi><mn>2</mn></msup><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>f</mi><mn>2</mn></msup><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><msup><mi>&mu;f</mi><mn>2</mn></msup><msup><mi>g</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&zeta;g</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&mu;g</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msub><mi>&delta;</mi><mi>st</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000016.GIF" wi="1635" he="155" /></maths>将公式(6)代入公式(5),整理得一个振动周期内吸收的能量:<maths num="0007"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfrac><mi>&Delta;W</mi><mrow><mn>2</mn><mi>&pi;</mi><mo>&bull;</mo><msub><mi>K</mi><mi>m</mi></msub><msubsup><mi>&delta;</mi><mi>st</mi><mn>2</mn></msubsup></mrow></mfrac><mo>=</mo><mfrac><msup><mi>&mu;&zeta;g</mi><mn>5</mn></msup><mrow><msup><mrow><mo>[</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>f</mi><mn>2</mn></msup><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><msup><mi>&mu;f</mi><mn>2</mn></msup><msup><mi>g</mi><mn>2</mn></msup><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mn>4</mn><mi>&zeta;</mi></mrow><mn>2</mn></msup><msup><mi>g</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>g</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&mu;g</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000021.GIF" wi="1635" he="167" /></maths>步骤三、绘制能量与激振频率比的关系曲线,然后统计单一因素变化对能量峰值的关系;根据统计数据拟合出对应各个质量比的最优频率比和最优阻尼比,具体包括以下处理:根据不同动力吸振器参数情况下能量曲线的峰值数据进行曲线拟合得出峰值最大且相等时的频率比,然后以最优频率比情况下不同阻尼比对应的吸收能量峰值为基础经过拟合得出最优的阻尼比;步骤四、根据外界干扰信号的频率范围对动力吸振器能够吸收的能量进行积分;然后绘制能量柱状图找到最大值,进而得到最优质量比;最优质量比根据激振力信号的频率的范围进行选择,选择依据为吸振器在一定频率范围内吸收能量的多少,即<maths num="0008"><math><![CDATA[<mrow><mi>W</mi><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>g</mi><mi>low</mi></msub><msub><mi>g</mi><mi>up</mi></msub></msubsup><mi>&Delta;Wdg</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000381940780000026.GIF" wi="1036" he="108" /></maths>根据上述公式计算吸振器吸收激振力能量的大小,吸收能量最多的即为最优质量比;步骤五、通过最优的质量比、频率比和阻尼比得到吸振器的最优参数,即质量、刚度和阻尼。<img file="FDA0000381940780000022.GIF" wi="166" he="137" />,进而有,M<sub>a</sub>=μ·M<sub>m</sub><maths num="0009"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><mfrac><msub><mi>&omega;</mi><mi>a</mi></msub><msub><mi>&omega;</mi><mi>m</mi></msub></mfrac><mo>=</mo><msqrt><mfrac><mfrac><msub><mi>K</mi><mi>a</mi></msub><msub><mi>M</mi><mi>a</mi></msub></mfrac><mfrac><msub><mi>K</mi><mi>m</mi></msub><msub><mi>M</mi><mi>m</mi></msub></mfrac></mfrac><mo>=</mo><msqrt><mfrac><mrow><msub><mi>K</mi><mi>a</mi></msub><mo>&bull;</mo><msub><mi>M</mi><mi>m</mi></msub></mrow><mrow><msub><mi>M</mi><mi>a</mi></msub><mo>&bull;</mo><msub><mi>K</mi><mi>m</mi></msub></mrow></mfrac></msqrt></msqrt></mrow>]]></math><img file="FDA0000381940780000023.GIF" wi="600" he="289" /></maths>,进而有,<maths num="0010"><math><![CDATA[<mrow><msub><mi>K</mi><mi>a</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>f</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msub><mi>M</mi><mi>a</mi></msub><mo></mo><msub><mrow><mo>&CenterDot;</mo><mi>K</mi></mrow><mi>m</mi></msub></mrow><msub><mi>M</mi><mi>m</mi></msub></mfrac></mrow>]]></math><img file="FDA0000381940780000024.GIF" wi="374" he="147" /></maths><maths num="0011"><math><![CDATA[<mrow><mi>&zeta;</mi><mo>=</mo><mfrac><msub><mi>C</mi><mi>a</mi></msub><msub><mi>C</mi><mi>c</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>C</mi><mi>a</mi></msub><mrow><mn>2</mn><mo>&CenterDot;</mo><msub><mi>M</mi><mi>a</mi></msub><mo>&CenterDot;</mo><mo></mo><msub><mi>&omega;</mi><mi>m</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000381940780000025.GIF" wi="408" he="137" /></maths>,进而有,C<sub>a</sub>=2·ζ·M<sub>a</sub>·ω<sub>m</sub>。
地址 300072 天津市南开区卫津路92号