数值分析成层地基中应力波传播规律.pdfVIP

数值分析成层地基中应力波传播规律 李本平刘军章克凌陈燕顺 提要夺义编制；+峨变动态弹嗵性冉限凡i芏序，针对“I．软下破”及“I叫下轼” 载，将半无限介质作轴对符有限离融．边界采用多重适射边界采模拟，聚用，负校正f}三 A的N“ma出方法避{了教使租分。通过埘成甚地基进行算捌分析．褂Ⅲr战甚地基中 戍力波传播的一牡州律。 ≥ 关键词成层地皋．瞬蹙动恋，J五力波 1 前言 用数值方法分析地摹中戍力波}|勺传播规律这一研究领域十分活跃，而成层地 基中应力波的传播和均匀介质’}麻力波的传播有许多不同之处。本文编制了瞬变 动态弹塑性有限元程序，针对版种常见的分层地基}目式进行j-计算分析t得{}5了 相应的一些规律。 2 瞬变动态有限元基本方程 动力有限元基本公式为： 【Ml{／i)+lc】{“}-【Km}_{Q(r)} piN]r如为总体质量矩阵 其中：[MI：∑f[A7】7 IK】：∑￡【B1■D】l引㈦v为总体刚墁矩阵 {Q(，)}=∑I【|、，】。∥／s为等效结点衙载向封 “’ f(1为阻尼矩阵，利划瑞售【Raylei曲)阻尼 【C]=“【M】+醴朗 q——边界外力父皱 i Afl——形函数矩阵 【剀——鹰变矩阵 【D】——应力应变本构关系矩阵 口、6为常数，一般i，通过振型阻尼参数蔷柳振动狮率q求H 。。：迎二三塑 婚)’ 甜』一肛，。 为I厂获得离散系统的础jJ响应，必；i!二R斛1F线形动力方氍(1)式t目前普 遍采刚直接积分法．．经典的疗法有两种：隐，t浊承j显示法，但备有利弊。本文采 l{jHug。j等人提出的显隐结台法2，经过计搏}：t埘：蛙鼗积分参数选取适当，稳 定性能得到保障。 3 计算模型 3 1 介质本构模型 弹塑性本构方稚采瑚F述形式： dajl=D蔫d￡。 ㈨ f 式中，』聪j为弹塑胜拽苗张蜒．，廷救表J厶式为： ∥》以。一些一錾喜∑ ㈤ ¨蠡：：‰一≯jj 其中：g——塑性势函数 F——硝服函数 A——硬化参数 Dm，——弹性模鲈张量 在众多的屈服准则巾，由_rDrucker—pHg¨屈服准!j{『j比较适台岩r介质． f}f}般准则，其表达式为： 闪此，本文采H{bruckt,r--0rager 盯，1+√-，：=K (6) 式中，J，为应力偏鼙的第一不变醢： ，：=f(仃。一盯，)1+(盯：一仃，)2+(c，．一盯、)：1，6 c，) 2§in d l 8) 。2焉_万蓑需了 订c c。s∞