基于分数导数模型的粘弹性桩振动分析.pdfVIP

第17卷2期 应用基础与工程科学学报 V01．17．No．2 2009 2009年4月 OFBASICSCIENCEANDENGINEERINGApril JOURNAL 文献标识码：A 文章编号：1005-0930(2009)01-0303-06中图分类号：TU459 基于分数导数模型的粘弹性桩振动分析 刘林超1， 杨骁2 (1．上海市应用数学和力学研究所，上海200072；2．上海大学土木工程系，上海200072) 摘要：研究成果表明混凝土桩具有粘弹性性质，为了准确分析粘弹性桩的振动特 性，必须建立准确的粘弹性本构模型．在分数导数理论、粘弹性理论、应力波理论 的基础上建立了基于分数导数模型的粘弹性桩的振动方程，利用Zhang-Shimizu 分数导数数值积分法得到了基于分数导数模型的粘弹性桩的振动方程数值解． 分析结果表明，分数导数微分算子的阶数和粘弹比对粘弹性桩桩端速度衰减的 快慢和衰减周期等有很大的影响． 关键词：分数导数；粘弹性；应力波；数值积分 桩基是一种重要的工程结构，已广泛地应用于土木和环境工程中，近年来吸引了许多 研究者的注意，特别是桩基的动力学行为，很多学者在这方面开展了研究，取得了许多成 果．然而，以往的研究通常将桩视为完全弹性杆件，而材料力学的研究成果表明，混凝土桩 体十分接近于粘弹性体．因此，将其从弹性体扩展到粘弹性，运用粘弹性理论研究桩的振 动速度及应力分布等问题，应更适合工程实际【l引．然而，由于对于粘弹性材料本构关系 的选取将直接影响其力学行为．因此，本构关系选取的恰当与否是岩土工程数值模拟需要 解决的重要问题，经典粘弹性本构模型的核函数(蠕变柔量、松弛模量)常是指数函数的 组合，这是由整数阶微分算子的性质决定的，欲精确描述实验的数据，常常不得不取消高 J． 阶的微分项或者以降低本构模型的应用范围为代价p 众所周知，具有指数函数内核的遗传力学积分方程等同于常系数线性微分方程．但实 验数据显示，仅具几个指数项的内核难以精确地拟合实验数据．这个事实表明，引进新函 数以描述粘弹性材料的特性十分必要．对这个问题本学科众多的专家学者一直在做探讨． 其中，基于分数导数理论的建模方法逐渐让人认识到是一种能较好地描述一大类粘弹性 材料力学行为的有效方法，而且研究表明利用分数导数粘弹性模型来描述混凝土的力学 行为是十分合理的H】．本文就是在采用能够较精确描述粘弹性材料力学行为的分数导数 本构模型来对粘弹性桩进行振动特性分析． 1分数导数粘弹性本构模型 分数导数模型能够精确地刻划某些粘弹性材料的力学行为，且仅需少数几个试验参 数即能在较宽阔的频带内拟合试验数据．早期的材料试验就表明，很大一部分的粘弹性材 收稿日期：2007-04-10：修订日期：2008-07-01 基金项目：上海市重点学科建设项目(Y0103)资助；上海市自然科学基金(06ZRl4037) 作者简介：刘林超(1979一)，男，讲师，博士研究生。E—mail：11c109@sohu．com 万方数据 应用基础与工程科学学报 料的复模量与频率的分数次幂成正比，后来经过力学工作者和数学家的共同努力，逐步建 立起来分数阶导数的数学概念，并用它来描述这种粘弹性本构关系，建立了广义的分数导 数粘弹性本构关系[5] (1) or(t)+∑6。Do”盯(t)=Eos(t)+∑E。D4“8(0 表达式为∞J 州巾)]=币％瓤潜曲 (2) 式中，0a1，厂(u)=I．广一eldt为Gamma函数． 对于分