质量_弹簧振子系统的ANSYS有限元分析
第24卷第2期2011年4月
常州工学院学报
JournalofChangzhouInstituteofTechnology
Vol．24No．2
Apr．2011
质量—弹簧振子系统的ANSYS有限元分析
黄兆梁
（常州工学院光电工程学院，江苏常州213002）
摘要：对质量—弹簧振子系统进行的ANSYS有限元模态分析，得出了与理论计算特征频率完
全一致的结果。通过这种相互一致的印证，可以确信关于弹簧有效质量的理论计算公式是准确、可靠的。
关键词：模态频率；有效质量；质量—弹簧振子；振动
中图分类号：O326；O241.82文献标志码：A文章编号：1671－0436（2011）02－0038－04
TheANSYSFiniteElementAnalysisofMass-springOscillatorSystem
HUANGZhao-liang
（SchoolofPhotoelectricEngineering，ChangzhouInstituteofTechnology，Changzhou213002）
Abstract：ThispaperappliesANSYSmodalanalysistothemass-springoscillatorsystemtoobtainthemodalfrequency，whoseresultsareinexactaccordancewiththoseeigenfrequencycalculatedbythetheory．Theconsen-susofsupporteachothercanmakesurethattheeffectivequalityoftheformulaisaccurateandreliable．
Keywords：modalfrequency；effectivemass；mass-springoscillator；vibration
条件所限，实验精度都不足以给出与理论值之间
高度吻合的有力证据，因此，对弹簧的有效质量问题并没有定论。有些文献对于弹簧质量1/3修正问题存有异议，或认为1/3仅仅是0.346的近，也有从各自的实验数据的分析中得出了该
［8－10］
，质量修正值的不同结论正所谓众说纷纭，似
莫衷一是。本文拟采用ANSYS有限元分析工具
来分析质量—弹簧振子的振动模态，获得了与文6］献［高度吻合的结果，从而对解决弹簧的有效质量问题提供了具有相当说服力的佐证。
［7］
0引言
对于弹簧质量不可忽略的质量—弹簧振子系
［1－5］
，统振动问题的研究国内外有许多报道多数
都认为只要将弹簧质量的1/3加到振子质量上去，就可以将该振动系统等效看作无弹簧质量的理想振子系统，其振动周期为：T=2πM为振子的质量，m为弹簧质量，k为其中，
弹簧的弹簧系数。但进一步的理论研究揭示，当弹簧质量不可忽略时，质量—弹簧振子系统的振
也不属于严格的周期运动，而是动并非简谐振动，
其振动的模式有无穷多，各模式的振准周期运动，
［6］
动频率之间也都不成有理数的倍数关系。大量文献报道了各自从理论或实验的不同角度分析该振动系统而得出的不同结论。由于受实验技术等
收稿日期：2011-04-02作者简介：黄兆梁（1956—
）男，副教授。
1ANSYS有限元建模分析
有限元建模中质量单元采用了mass21，并定
义为二维质量单元，弹簧采用combin14，并定义为二维纵向弹簧单元。作为对照与校验，同时也采用了beam3梁单元来模拟线性弹簧，不过梁
第2期黄兆梁：质量—弹簧振子系统的ANSYS有限元分析
39
beam3单元并不能直接输入弹性系数k，而需要通过输入弹性模量E，梁截面面积A，以及梁的长度L来实现梁单元的纵向（沿梁的长度方向）的弹性系数为k的目的。由于这里的梁并不用来作承重分析，且不需要用到梁的截面惯量矩，因此截面形状可以随意选取，方便起见，截面形状取正方
2
边长都设为1m，则面积为1m，梁的长度也形，
设为1m，在输入梁的弹性系数k时，在数值上已
都只给出了一阶振动模态，这意味着在这种情况
下确实只存在唯一的一阶振动模态，而不存在其这清楚地表明了理想弹簧振子的他的振动模态，
运动是属于简谐振动。ANSYS有限元分析所给出的5位有效数字与理论计算值是完全一致的。
对于一质量为m、弹簧系数为k的连续弹簧每在作ANSYS分析时需将其分成均等的n段，一小段连续弹簧都被看作一质量为m/n的集中质量与刚度系数为n·k的无质量小弹簧构成，由n个这样的离散小质量—弹簧振子组成的系统来代替原来的连续弹簧，经这样离散化处理后的情况与一个实际的连续的有质量的弹簧在力学特性上是基本等效的。当然，这里所分的段数n越多就越能精确模拟实际的连续弹簧，稍有影响的是这个小的集中质量是该放置在小弹簧的什么位置，为此在分析中作了2种极端情况的处理，即将小弹簧的集中质量分别放置在小弹簧离振子较远
并分别称之为远振子端和近振一端和较近一端，
子端。分析结果显示，在这2种情况下相同阶的
模态频率存在着少许差异，这种差异可以作为ANSYS的分析精度的一种衡量尺度，可以反映出分析的精度，而正确的模态频率应该介于2种极端情况之间。为了能够清楚地展示ANSYS的分
本文将弹簧析精度随划分的单元数的变化情况，
1000和10000个单元的离散质划分成了100、
量—弹簧振子，每种情况下又可分为“远振子端”“近振子端”、与两种极端情况，并分别用“+”“－”附在所划分的单元数之后作为标记。
经与弹性模量E相同，从而可以将弹性系数k的
数值直接输入到梁单元的弹性模量EX中，而泊松比直接以0填入到泊松比输入提示符PRXY相应的输入框中去。至于梁单元被设置为正方体及梁单元不作承重而作弹簧似乎不合常理的担忧是没有必要的，这样做只是为了ANSYS的材料特性参数设置的便利并确保分析结果的正确，况且由于有combin14弹簧单元作为参考印证，可以保可靠。弹簧质量m可以通过证分析结果的准确、
对材料密度来进行设置，在这里可以直接将弹簧质量m的数值输入到密度输入提示框DENS中。分析中约束了除纵向振动之外的所有自由度，而固定端则需作全约束。为了检验ANSYS建模分析的有效性与可靠性，首先分析无质量弹簧的理想振子的情况，并与理论计算的结果相比较以检
假定振子质验两者相符的程度。为了计算方便，
弹簧的弹性系数k=1N/m，由理论量M=1kg，
公式可计算得出其振动频率为：
11ω
f===≈2π2πM2π10．159154946（Hz）
用梁单元beam3模拟无质量（m=0）及弹性系数k=1N/m，振子质量M=1kg的理想弹簧振子，无论将弹簧划分为多少个网格单元，其第一阶振动模态的频率都为0.15915Hz，这与理论值完全吻合。
用ANSYS的弹簧单元combin14来模拟弹簧就要简单得多，可以直接设置弹簧质量m和弹簧系数k，用其进行模态分析得出的理想弹簧振子的第一阶振动模态的频率也同样是0.15915Hz，两者没有任何差异。
虽然在模态分析中要求展开至30阶，而实际结果无论是采用combin14单元还是beam4单元
2模态分析结果
采用ANSYS的APDL语言进行简单的编程，将振子质量M、弹簧质量m、弹簧弹性系数k、弹簧划分单元数n均作为输入参数，通过在输入界面上更改这些参数就可以很方便地对不同参数的情况进行反复计算。经过上述的ANSYS仿真分析后，即可求解得出各模态的振动频率。这里针对combin14和beam3单元分别编写了两段程1000和10000等序，在将弹簧分别划分为100、
份的情况下，对振动模态进行有限元仿真分析，其中对弹簧单元combin14还需区分远振子端与近
其前30阶振动模态的振动频率振子端2种情况，
40
以及理论计算值见表1。
常州工学院学报2011年
表1ANSYS计算值与理论值的比较
（M=1kg，m=1kg，k=1N/m，combin14与beam3单元）
阶次[***********][***********][**************]
combin14单元
100+0．137180．545371．02451．51612．01132．50753．00403．50013．99574．49064．98445．47725．96866．45876．94727．43407．91918．40228．88329．36219．838710．31310．78411．25311．72012．18312．64313．10113．55514．005
100－0．136680．544981．02421．51602．01112．50743．00393．50013．99574．49054．98445．47715．96866．45876．94727．43407．91908．40218．88329．36219．838610．31310．78411．25311．72012．18312．64313．10113．55414．005
1000+0．136950．545221．02451．51662．01262．51013．00843．50714．00624．50555．00495．50436．00396．50347．00317．50278．00238．50209．00169．501310．00110．50111．00011．50011．99912．49912．99813．49813．99714．497
1000－0．136900．545181．02451．51662．01252．51003．00843．50714．00624．50555．00495．50436．00396．50347．00307．50278．00238．50209．00169．501310．00110．50111．00011．50011．99912．49912．99813．49813．99714．497
10000+0．136930．545211．02451．51662．01262．51013．00843．50724．00634．50565．00515．50466．00426．50397．00367．50348．00328．50309．00289．502710．00310．50211．00211．50212．00212．50213．00213．50214．00214．502
10000－0．136920．545201．02451．51662．01262．51013．00843．50724．00634．50565．00515．50466．00426．50397．00367．50348．00328．50309．00289．502710．00310．50211．00211．50212．00212．50213．00213．50214．00214．502
1000．136930．545231．02471．51722．01392．51273．01293．51434．01694．52065．02575．53206．03986．54927．06027．57308．08778．60449．12329．644310．16810．69411．22211．75412．28812．82513．36613．90914．45615．007
beam3单元
10000．136930．54521．02451．51662．01262．51013．00853．50734．00644．50585．00535．50496．00466．50437．00427．50418．00408．50409．00409．504110．00410．50411．00411．50512．00512．50513．00613．50614．00614．507
100000．136930．54521．02451．51662．01262．51013．00843．50724．00634．50565．00515．50466．00426．50397．00367．50348．00328．50309．00289．502710．00310．50211．00211．50212．00212．50213．00213．50214．00214．502
理论值0．1369260．5452041．0245281．5166412．0125602．5100783．0084123．5072174．0063194．5056205．0050595．5046006．0042186．5038947．0036167．5033758．0031658．5029799．0028139．50266510．0025310．5024111．0023011．5022012．0021112．5020313．0019513．5018814．0018114．50175
由表1可见：
1）ANSYS计算值与理论值非常吻合，尤其是ANSYS所给出的5位有效当划分单元数较多时，
数字与理论计算值无差异。
2）当弹簧质量偏置于远振子端时频率偏高；而当弹簧质量偏置于近振子端时频率偏低；其平均值更为接近理论值。
3）相比而言，远振子端情况的数据要比近振
第2期黄兆梁：质量—弹簧振子系统的ANSYS有限元分析
41
子端情况的数据略为准确一些。这是由于弹簧质量的一级修正值为弹簧质量的1/3。
4）采用梁单元beam3来计算振动模态比用combin14单元更为准确，其收敛于理论值的速度更快、效果更好，也不用区分远振子端还是近振子端。但在无理论值的情况下不便于精度估计。
6］中的式（8）表1中的理论值是根据文献［
表2
列表项目一阶频率/Hz相对误差/%
无修正0．1591516
ωωm
tan=，，其中ωi=2πfi=用数学软件
mωiωiM
Mathematica7．0通过编程计算得到的前30阶振
i=1，2，…，30，动模态频率fi，单位Hz。
对于一阶振动模态，表2给出了采用Math-ematica7．0根据文献［6］各级修正公式计算得出的近似频率、理论值及其相应的相对误差。
一阶振动模态各级修正公式的精度
一级修正0．137830．7
二级修正0．13670－0．2
三级修正0．136970．03
四级修正0．136926－0．00025
理论值0．13692634
———
随着计算公式的修正级次的增加，通常计算
，精度也越高，但修正级次过高会出现“振荡”适当选取修正级次可获得较好效果。表2中的理论
6］值是根据文献［中的式（8）计算得到的，其五级
（原文献［6］式（11）中ω1（4）为ω1（5）之误，即所给公式为一阶振动模态ω1的五级修正而非四级修
）振动频率修正是根据文献［6］正。中的式（11）的角频率计算得出，即：
ω1（5）=
而四级修正为：
ω1（4）=
其余可依次类推。
1m1mm
1+－M+
15M63M3[[
()()
()()
2
11m
+
4725M
()
3
m1
+
7425M
()]
4
1m1mm
－1+M+
15M63M3
2
11m
+
4725M
()]
3
值得注意的是，当采用beam3单元来对上述
如果将弹质量—弹簧振子系统进行模态分析时，簧单元仅仅划分为一个单元时，这时的质量—弹
簧振子系统就仅仅只有一个自由度，因此其得出的振动模态就仅仅只有一阶，其频率为0.13783Hz，正好与表2中的一级修正公式也即1/3弹簧质量修正公式得到的振动频率完全相同，这并不是偶然的巧合。事实上，如果将弹簧单元仅仅划分为一个单元，这就意味着弹簧被限制为只能作
而不能作疏密相间的弹性均匀的线性伸缩变化，
1/3的弹簧质量修正的波动的变化。众所周知，
公式的推导正是建立在弹簧仅仅作均匀的线性弹
性伸缩变化的假定之下推导而得到的结论，因此，这里两者结果完全相同也是在意料之中。一般情况下弹簧的伸缩并不是线性变化的，而是呈现出
比较复杂的波动，这正是存在无穷多振动模态的
原因所在。
3结论
由ANSYS仿真得出的各阶振动模态的频率
［6］
与由本征值方程ytany=β求解得出的对应特征频率完全一致，这就充分说明将连续弹簧视作均匀的连续介质，用波动方程来描述弹簧的振动是行之有效的，并且具有相当高的精度，能够相当精确地描述弹簧振子在振动过程中由弹簧引起的波动现象。另外，在认可ANSYS分析的结果是6］准确、可靠的基础上，那么也就能验证文献［中的式（11）所求出的振动模态的正确性，而作为其中一阶振动模态的弹簧有效质量的计算公式自然也就是准确无误的。
（下转第73页）
第2期2.4
朱建国：基于SWOT模型的江苏省高校体育产业现状分析与对策
73
广施“双赢”战略，在校内外建立携手共进
的伙伴关系
相对而言，社会体育事业产业化之路起步早，发展机制灵活，除拥有较强资金实力外，或拥有政府支持的背景，或拥有垄断体育竞赛资源的优势，然而它们与高校体育具有相似性和互补性，高校体育产业的最佳市场营销战略应定位于与这些社会体育事业建立“双赢”的伙伴关系。同时还必在体育产业须充分利用校内多学科的人才优势，发展理念、市场建设、营销策略以及体育产品开发等方面全面加强合作，携手共进，这样才能实现体育产业健康快速发展。
培育校内外体育消费市场与需求有机地结合起来，走创新发展之路；在策略上，以品牌建设和经努力打造出国内第三产业的优势济效益为中心，
品牌。江苏省高校体育产业的发展一定能跃上一个新台阶，成为一个新的经济增长点。
［参考文献］
［1］鞠成军，廉俊颖，张孔军，等．美国高校体育产业发展的启
J］．体育学刊，2003（2）：140－142．示［
［2］史松涛，韩建．江苏省大学生体育消费状况的调查研究［J］．
2007（10）：136－138皖西学院学报，
［3］江苏人民政府．江苏人民政府关于加快发展体育产业的实施
Z］．苏政发〔2010〕110号，2010－09－10．意见［
［4］江苏省统计局．江苏城镇居民文教娱乐消费简析［EB/OL］．
（2007－12－18）．http：//WWW．jssb．gov．cn/jssb/tjfxzl/[1**********]36．htm．
［5］“大拇指经济”［EB/OL］．（2000－叶岱夫．“注意力经济”与
06－22）．http：//www．people．com．cn/GB/channel7/498/20000622/113912．html．
3结语
体育产业是新兴产业、绿色产业、朝阳产业，
是现代服务业的重要组成部分。江苏省高校在体育产业发展应当充分发挥自身的优势，优化劣势，抓住机遇，取长补短。在认识上，理清高校体育产业发展的思路，以培育经济新的增长点为目标，把
责任编辑：张秀兰
櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚櫚（上接第41页）
通过理论分析与有限元分析结果吻合一致的相互印证，可以有充分理由来终结有关弹簧有效质量到底该如何准确计算的长期悬而未决的争议，也能充分理解为何通过实验来进行弹簧有效质量检验的种种努力始终都未能给出令人信服的结果。
［参考文献］
［1］．大学物理，丁履成，司明扬．有质量的弹簧的振动问题［J］
1985，4（9）：15－18．
［2］罗蔚茵．关于弹簧振子固有频率的进一步讨论［J］．大学物
1985（11）：9－11．理，
［3］LeclercM．振动弹簧的有效弹性常数和有效质量［J］．李国
1990，9（6）：10－11．强，译．大学物理，
［4］刘大鹏，关荣华．弹簧质量对弹簧谐振子圆频率的影响［J］．
1995，4（10）：22－24．大学物理，
［5］NunedaSilvaJM．负荷弹簧的重整化振动研究［J］．孙厚谦，
1995，14（10）：28－30．洪林，译．大学物理，
［6］．大学物理，1998，3（17）：黄兆梁．弹簧质量对振动的影响［J］
12－16．
［7］陈美华，蒋毅．弹簧振子计及弹簧质量的误差分析［J］．嘉应
2006，24（6）：27－30．学院学报，
［8］刘伟．用振动法研究弹簧振子的有效质量和倔强系数［J］．实
2006，10（23）：38－39．验技术与管理，
［9］．汉中师范孙彦清．小振动系统弹簧质量影响因子的测量［J］
2003，6（21）：50－52．学院，
［10］肖波齐．基于Matlab的弹簧振子简谐振动研究［J］．陕西科
2009，12（27）：116－120．技大学学报，
责任编辑：张秀兰