分数阶导数型粘弹性结构动力学方程数值算法研究.pdfVIP

第34卷第2期林业大学学报(自然科学版)Vol.34,No.2

2010年3月JournalofNanjingForestryUniversity(NaturalScienceEdition)Mar.,2010

http://www.nldxb.com

分数阶导数型粘弹性结构动力学方程的数值算法

12311

银花,,,

(1.林业大学土木,江苏210037;2.林业大学机械学院,江苏210037)

:基于精细积分方法,提出了具有分数阶导数型本构关系的粘弹性结构动力响应的一种新的

数值计算方法。该方法首先将系统的动力学微积分方程转化为含分数阶导数项的一阶常微积分

方程组,然后采用精细积分法对方程进行积分计算得到系统响应。数值计算结果与解析法及

Zhang2Shimizu算法的结果相吻合,并显示随计算步长减小其计算的收敛性更好。

:分数阶导数;粘弹性;动力学;状态方程;精细积分

号:O322文献标志码:A文章编号:1000-2006(2010)02-0115-04

Anumericalalgorithmforthedynamicsequationofthefractional

derivativeviscoelasticitystructure

YINHua1,CHENNing23,ZHAOChen1,WANGDa2ming1

(1.CollegeofCivilEngineering,NanjingForestryUniversity,Nanjing210037,;2.CollegeofElectronicand

MechanicalEngineering,NanjingForestryUniversity,Nanjing210037,)

:Accordingtothepreciseintegralmethod,anewnumericalalgorithmisproposedfordynamicalresponsesofvis2

coelasticitystructurewiththefractionalderivativeconstitutiverelation.Inthisapproach,thefractionaldifferentialequation

forthedynamicofasystemistransformedintoasetoffirstorderordinarydifferentialequationswhichcontainfractional

derivativeterms.Thepreciseintegralmethodisusedtointegratethesetermsandobtaintheresponseofthesystem.Nu2

mericalresultsobtainedusingthisschemeagreewellwiththoseobtai