考虑土一结构相互作用效应的结构地震响应时域子结构分析法.pdfVIP

第33卷第5期 北京工业大学学报 VbI．33No．5 ING 20()7年5月 JOURNAI。oFBEIJ UNIVERSITY0FTECHNOI。0GY Mav2007 考虑土一结构相互作用效应的结构地 震响应时域子结构分析法 杜修力，赵建锋 (北京工业大学城市与工程安全减灾省部共建教育部重点实验室，北京 100022) 摘要：为了在地基阻抗力的计算中能考虑基础动力阻抗函数频率相关性和避免反复进行傅里叶变换，采用地 基阻抗力的时域差分计算方法，进行了基于通用有限元软件的结构时程分析程序二次开发，增加了考虑土一结构 相互作用效应的结构动力响应时程分析计算模块．对圆盘基础上的多层结构在三维地震激励下的响应进行了 数值计算和分析．数值模型结果表明，该方法物理意义明确，能反映土一结构相互作用耦合体系的基础动力阻抗 对激振频率的依赖性． 关键词：土一结构相互作用；阻抗函数；频率相关；动力响应；子结构法 470 中图分类号：Tu 文献标识码：A 文章编号：0254一0037(2007)05～0517一07 随着震害经验的积累和抗震计算理论的发展，在结构动力响应分析中完整、准确地考虑土一结构相互 作用(soil—structure 时，为考虑土一结构动力相互作用和结构的非线性行为，必然需要计算基于频率相关基础阻抗函数的地基 阻抗力…．频域子结构法[2]能充分考虑基础阻抗函数频率相关特性，但无法处理结构自身的非线性问题， 须先将结构的非线性构件进行等效线性化．时一频域混合法[3】和频一时域混合法[4。5]均需反复进行傅里叶 变换和对每一计算步骤进行多次校核．而时域子结构法由于基础阻抗函数随频率变化复杂，使得通过基 础阻抗函数卷积表述的时域地基阻抗力的计算变得十分复杂和耗时． 时域地基阻抗力的实现目前主要采用集总参数模型法、脉冲反应法和集总参数模型耦合虚拟力法． 为了避免反复进行傅里叶变换并节省存储空间，众多国内外学者已提出了几种集总参数模型[6。11．集总 参数模型实际上是对阻抗函数进行等效线性化处理．然而，基础阻抗函数通常具有强烈的频率相关性，仅 采用个别频率动力阻抗函数值的集总参数模型难以充分包含所有频率成分，不能反映真实性态．根据傅 质是基于时频变换的积分算法，虽然避免了反复卷积，但需要较多的计算量和存储空间．wolf等提出了时 域递归方法来计算地基阻抗力【15|，其中递归方程的系数通过脉冲响应不变法在时域内计算获得．wolf等 又基于2变换给出了利用频域阻抗函数的无限地基阻抗力时域递归计算方法[16|，并用最小二乘法通过频 域2个频率多项式的比值近似拟合阻抗函数来获得递归参数，使得某时刻的阻抗力由当前时刻的交界面 位移以及前几时步的相互作用力与位移来计算．Safak提出将时域差分递归模型用于时域地基阻抗力的 在本质的缺陷，这就降低了对时域地基阻抗力的拟合精度．杜修力等将集总参数模型与时域差分递归模 型相结合，提出了一种可完整考虑基础阻抗函数中奇异项特征和正则项特征的简便、实用且精度可控的时 域地基阻抗力的计算方法[1J．该方法将当前时刻的时域地基阻抗力的分解表述为2项有明确物理意义的 收稿日期：2007—01—26． 北京市属市管高校人才强教计划资助项目(05004012200515)． 作者简介：杜修力(1962一)，男，四川广安人，教授，博士生导师． 万方数据 518 北京工业大学学报 2007年 阻抗力之代数和，即一部分是由弹簧阻尼器集总参数模型表示的反映瞬态影响的阻抗力，另一部分是由前 几个时步的虚拟力和位移值表示的反映延续影响的虚拟力．作者采用该方法，通过在有限元软件中引入 虚拟力的方法实现了考虑阻抗函数频率相关性的地基阻抗力的时域计算． 1时域地基阻抗力计算方法 在频域中，基础一结构交界面上的相互作用力一位移关系可写为 F(∞)=S(叫)U(∞) (1) 数S(叫)可