一种基于动态骨架曲线的土层时域积分方法及其验证 A Direct Time-domain Integral Method Based on Dynamic Skeleton Curve Constitutive Model and Its Verification.pdfVIP

第2l卷1期 应用基础与工程科学学报 V01．21。No．1 2013年2月 2013 JOURNALOFBASICSCIENCEANDENGINEERINGFebruary 文章编号：100541930(2013)01-0078-012中图分类号：P315文献标识码：A dd：10．3969／j．issn．1005-0930．2013．01．009 一种基于动态骨架曲线的土层时域 积分方法及其验证 荣棉水1”， 卢 滔3， 李小军2， 吕悦军1， 赵亚敏1 (1．中国地震局地壳应力研究所，北京100085；2．中国地震局地球物理研究所，北京100081；3．防灾科技学院， 河北燕郊065201) 摘要：由于土层等效线性化方法在强地震动输入下会低估地表峰值加速度且不 适用于处理具有较大非线性的软弱土层．本文提出了基于动态骨架曲线的土层 直接时域积分方法．该方法在动力方程的数值积分中采用中心差分格式，土体的 剪切动力本构关系采用了基于动态骨架曲线的简单函数形式以考虑其非线性， 在土层底部引入人工透射边界以考虑能量向基底的辐射，并基于MATLAB平台 编制了DYMAlfI．AB计算程序．以唐山响堂三维场地影响观测台阵的强震记录为 例对本文方法进行验证．结果表明，本文提出的方法能较好地考虑土体的非线 性，是一种合理可行的土层地震反应分析方法． 关t词：等效线性化方法；动态骨架曲线；直接时域积分方法 在工程抗震设防工作中，场地设计地震动参数的确定是首当其冲的重要问题．目前国 内工程界最为通用的土层地震反应分析方法为等效线性化方法，其基本思想就是在总体 动力学效应大致相当的意义上，用一个等效剪切模量及等效阻尼比替代所有不同应变幅 值下的剪切模量和阻尼比．由于剪切模量和阻尼比与应变幅值无关，整个问题化为线性问 题…．在强地震动输入下等效线性化方法会低估地表峰值加速度水平，实际台阵的观测 记录与计算结果的比较已证实了这一现象忙j．此外，由于等效线性化方法基于叠加原理， 仅适用于小的非线性情形，处理具有较大非线性的软弱土层时也存在局限性，发展时域土 层反应分析方法是工程实际的需求．．非线性土层地震反应分析中非线性动力方程数值积 分方法的选择以及土体应力应变关系的确定是两个最为核心的问题，另外，计算基底边界 的处理也是较为关键的环节．许多专家和学者在时域积分方法方面都曾做过探讨和研究， 如栾茂田pJ、GuangYuHl、Jianye 方法提出了各自的时域方法，但其对底边界的考虑均为粘滞边界，这就不可避免地在计算 中要形成总刚度矩阵和求解大型方程组．李小军【61基于土体动态骨架曲线和显式差分解 收疆日期：2012-01—19；修订日期：2012-05-23 项目)(加1lc舶13酏1)；国家自然科学基金资助 作者筒介：荣棉水(1粥2一)。男，博士生，助研．E—mail：wahrong@126．com 遗信作者：卢滔(1979一)，男。博士后，副教授．E-mail．grit_lu@yahoo．咖．cn 万方数据 No．1 荣棉水等：一种基于动态骨架曲线的土层时域积分方法及其验证 79 法提出了直接时域积分方法，其对基底能量辐射效应的考虑也并非直接采用透射人工边 界，且由于早期强震记录资料极其有限，难以利用强震记录予以检验．随着人工边界条件 研究的进展，廖振鹏…等提出的多次透射边界由于具有时空解耦特性，便于计算机实现， 也可通过多次透射使其精度与有限元精度一致，因而许多学者¨剖将其应用于时域求解， 提出了基于UE本构模型、弹塑性本构模型和新的动力方程积分格式的时域积分算法．另 外，强震观测资料的积累，特别是场地影响台阵的建设为各种时域算法提供了重要的验证 资料．卢滔【9o等引入李小军基于动态骨架曲线且能同时拟合动剪切模量比和阻尼比试验 曲线的本构关系与中心差分及多次透射人工边界条件结合，提出了一种新的土层时域积 分方