含双椭圆形孔洞凸起地形对瞬态SH波的散射与动应力集中研究.pdfVIP

含双椭圆形孔洞凸起地形对瞬态SH波的散射与动应力集中 齐辉郭 晶姚东王兴涛 (哈尔滨工程大学建筑工程学院 黑龙江哈尔滨 150001) 摘要 本文研究在地震波传播过程中，当其接近地表面时．平面SH波对两个相邻椭圆形孔洞的散射 问题，并确定其对地震的影响。利用了有限元和软件ANSYS模拟了双椭圆形孔洞对SH波散射和动应 力集中问题，并利用SH波散射的对称性，构造了满足半无限空间的位移解和应力解以及孔心距不同时 的动应力集中系数。各种情况下波传到各个点时基本是同步的，只是位移幅值有一定的变化。在孔心 距较小时，位移的幅值较大。孔心距不同时右侧动应力集中系数明显大于左侧，且随着孔心距离的增 大动应力集中系数变小。 关键词 波动有限元法；sH波散射；动应力集中；ANSYS程序；瞬态分析 引言 研究不同形式地形对地面运动的影响是地震工程学研究的基本课题之一。而地形变化的模式 则可以分成凹陷地形和凸起地形两种基本形式，研究SH波对多个椭圆形孔洞的散射，不仅具有 理论价值，而且在工程中具有广泛的应用前景…，特别是在材料科学和结构抗震抗暴等研究领域。 随着计算力学的发展，对含椭圆孔凸起对弹性波散射问题的研究，在一般意义上讲，求解该问题 的理论解是有相当困难的。随着计算力学的发展，该问题可以用不同的数值方法进行求解，并可 得到数值结果[5~71。 本文使用有限元方法研究SH波对椭圆形孔洞凸起的散射阀题。其关键是在求解该问题时， 人工边界的设置问题，由于使用软件的限制，本文中将有限元模型取得较大，以减少边界反射对 结果的影响【2~4】。采用有限元法研究近场含双椭圆形孔洞对瞬态SH波的散射和动应力集中问题。 然后，用NEWMARK直接积分法进行求解，得到各节点和单元的位移、应力的时程解。本文对 要解决的问题建立了有限元模型，并用通用有限元分析软件ANSYS进行了计算，给出了部分节 点的时程解和孔边动应力集中系数，并对结果进行了讨论。 1 问题的表述 孑L洞这样的缺陷存在于很多地方，含椭圆形孔洞凸起的弹性半空间如图l所示，这种结构存 在于很多问题中，如地下通道和人防工程等都可以简化成该问题。介质的密度和横波速度分别为 P和C，椭圆的长轴和短轴分别为a和b。 ·-———141·--—— y 图l SH波作用于双椭圆形孔洞凸起的半无限空间模型 2算例分析与讨论’ 2．1有限元模型的建立 本文是模拟半无限空间中含双椭圆形孔洞问题，所以由地表产生的反射波的影响要考虑。这 样，计算时间就要取道使反射波达到并通过孔洞。考虑以上不同，用ANSYS软件建立的有限元 模型如图2所示，单元名称为PLANE83。本章中取模型长100m，宽50m。上部圆弧半径为5m， 满足横向宽度为结构尺寸的5-6倍的要求。两个椭圆孔洞的长轴和短轴分别为10m和5m，两个 椭圆问的距离为10m。 图2有限元模型 2．2材料参数的选取(表1) 表1物理参数 密度／(kg／m2) 弹性模量／GPa 泊松比 1650 458 O．28 2．3边界条件的处理 边界条件为： 上部水平边界：自由 ·--——142--_—— 两侧边界：UX=UY=UZ=0 底部水平边界：吸=UY=0 2．4荷载的施加 荷载的计算时间足够长，能够使地表反射波达到孔洞并通过孔洞。时间积分步长取O．0005s， 满足＆≤I}l／c0。荷载的周期为O．02s，频率为50Hz，小于截止频率皱=2c0／,5X所要求。 2．5 NEWMARK积分法 NEWMARK积分法实质上是线性加速法的一种推广。它采用下列假设 喀+Af=喀+【(1-8)／／,+艿茸+址]At (1) 1 at+山