道路结构中瑞利波特性.pdfVIP

第五届全国土动力学学术会议论文集 大连·1998．6 道路结构中瑞利波特性 夏唐代 丁狄刚朱少杰 (浙扛大学土木工程系，坑螂310027)(杭州市市政工程建设处，坑州310001 摘要车文介绍丁几种利用有限单元法及解祈法相结合的方法计算道路结构瑞利渡弥散曲线。并用该 方法详细地讨论丁道路结构中瑞利波弥散特性以及位移分布规律。 l引言 随着交通事业的发展，高等级公路及机场的建设日趋增多，常规检测，路面质量的弯 沉法已不能满足要求，需一套能快速可靠检测道路质量的测试分析系统。表面波法是一种 能快速无损测试土参数的方法，因此也可利用表面波来无损检验道路结构体系。’。在表面 渡测试分析道路质量时，需事先了解道路结构体系瑞利波的弥散特性，反分析才能进行。 以往对道路结构瑞利波弥散特性只局限于上硬下软地基(即结构的剪切波速随深度减小的 情况，参见本文算例1情况)o“)，这类道路结构在土木工程中很难遇到。一般，道路结构 中结构的剪切波速起初随深度减小(如混凝土结构、垫层、填土等)，而后剪切波速叉随深度 增加，这类道路结构在土木工程中经常遇到。因此研究这类道路结构瑞利波弥散特性更具 有工程意义，可进一步完善表面渡法在无损检验道路结构质量中的应用以及在动力响应计 算中的应用， 乖j用瑞利波特征方程计算弥散曲线时，一般是假设相速度c已知，再利用矩阵特征根 的计算方法计算相应模态的波数t“’。对予上软下硬地基，c与各模态的k是一一对映的“，； 但对软夹层地基，c与各模态的^可能是一多对映”，即对于速度C某一模态存在几个波数 ^，此时由特征方程计算出的波数^很难判定是第几模态的波数。本文将另外介绍两种计算 瑞利波弥散曲线的方法——迭代法和柔度系数法。 2方法简介 这里介绍几种利用有限元法及解析法相结合的方法，计算瑞利波弥散曲线。 2．1方法1 ·般方法o’利用有限元法及解析法相结合的方法建立了成层地基中瑞利波特征方程， 该特征方程具有求解容易的特点。特征方程可分解成渡数k的二次特征方程 CA3k2+CB]^+CC3l=0(1) 层的剪切波速、压缩波速、层厚度及褶速度‘构成。因此，假设相速度c为已知，由式(1)就 可以求得波数k，从而得瑞利波弥散曲线(波长^一2“／五)．对于上软下硬地基瑞利波弥散曲 线有如图1(a)所示形式“’，此时相连度c与波数^有一一对映关系，由c可求得多个模态 波的弥散曲线．文6讨论过三层夹层地基的瑞利波弥散曲线有如图1(b)所示形式，此时由 相速度c可求出两个第一模态波数A，C与{不是一一对映关系．对于更复杂情况(见算例分 析)，由c很难判断哪些波数k是第一模态和第二模态的，因此式(1)的计算很难适应图1 (b)所示情况． 图1璃利弥散醢缓 2．2方法2：迭代法 对于某一固定频率m，瑞利渡各模态的波数l是唯一的，同理也可以得到唯一的速度c (“=k)。因此可假设振动频率m已知，求解出渡数k或相速度c．但由式(1)可知，在无限 r—————一r—————一 层提供的边界力中存在不可解项n一√1一c2／诈和6一√1一一／瞻，弭和v。分别为半无限 层的压缩波速和剪切波速，因此由式(1)只能采用迭代逼近求解方法得^或c．此外，由于 是迭代求解，且由式(1)可求得多个模态的坡数^，因此初始值的选取极为重要，需判断所 求得的k是第几模态的波效，计算效率低。 2．3方法3：柔度系数法 文献7利用解析分析方法建立了地基动力响应的动刚度法(或传递矩阵法、柔度系数 法)，文献8利用有限元法建立了地基动力响应的动刚度法，后者可用于分析非匀质地基响 应。动刚度法中力与位移有如下关系 [国{；)一商) (2) 式中心]为动刚度矩阵，实际上矩阵(s)即为式(】)的特征方程矩阵，{Q)为外荷载，{“)为 节点位移。当有表面竖向外荷载单位力作用时，即令式(2)的{面)中第二项为单位1，其它 为零，由式(2)可求出表面质点的竖直向位移．进一步可得柔度系数为