孔隙水压力静力触探动态贯入过程的有限元模拟

孔隙水压力静力触探作为一种新型的土体勘测技术已被广泛应用,但是,其机理的研究尚存在很大的不足。笔者系统分析了孔压触探的贯入特征,提出了“孔隙水压力体积基准值”的概念,借此推导出用土体体积变形量计算超孔隙水压力的公式。在此基础上,用轴对称有限单元法模拟探头的动态贯入过程,揭示应力－应变的变化特征,寻找超孔隙水压力生成与消散的规律,对孔压触探的机理进行了研究。     

System identification of a concrete arch dam and calibration of its finite element model

A magnitude 4.3 earthquake occurred near Pacoima Dam on 13 January 2001. An accelerometer array that had been upgraded after the Northridge earthquake recorded the motion with 17 channels on the dam and the dam–foundation interface. Using this data, properties of the first two modes are found from a system identification study. Modal properties are also determined from a forced vibration experiment performed in 2002 and indicate a significantly stiffer system than is estimated from the 2001 earthquake records. The 2001 earthquake, although small, must have induced temporary nonlinearity. This has implications for structural health monitoring. The source of the nonlinear behaviour is believed to be loss of stiffness in the foundation rock. A finite element model of Pacoima Dam is constructed and calibrated to match modal properties determined from the system identification study. A dynamic simulation of the 2001 earthquake response produces computed motions that agree fairly well with the recorded ones. Copyright © 2006 John Wiley & Sons, Ltd.     

Transient analysis of wave propagation in non‐homogeneous elastic unbounded domains by using the scaled boundary finite‐element method

The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.     

巴姆地震变形场和应力场：Ⅱ．用FEPG有限元方法求解

本文利用 FEPG 有限元程序对2003年12月的巴姆地震变形场和应力场进行了模拟计算。2003年12月26日,巴姆地区发生6.6级地震。Nakamura、Suzuki 等人利用余震资料和地震破坏程度的分布图分析了余震震中的统计分布形式,得到了此次地震的发震断层是巴姆断层以西大约5km 的 Arg-e-Bam 隐伏断层的结论。凌勇等人从差分干涉雷达(D-InSAR)出发,得到了此次地震的同震干涉条纹和雷达视线方向的位移场,并从雷达干涉的相干图上确定了发震断层的痕迹。他们还利用 Okada 程序模拟计算了此次地震的位移场,模拟计算的结果与 D-InSAR 测量的位移场比较符合。Okada 程序从断层的位错为出发点进行计算,从运动学角度对此次地震变形场和应力场进行研究;本文从动力学的角度出发,研究远场区域边界的应力状态以及断层的物性参数的变化对巴姆地震形变场的影响,将各种计算结果进行对比的同时也和前人的结果相比较。本文计算中采用了凌勇等人的断层参数,即 Arg-e-Bam 隐伏断层为近南北走向,向东倾斜,近直立的倾角,右旋走滑,其(Strike,dip,slip)为(175°,80°,173°)。断层未延伸到地表,自地面以下1km 处往下延伸16km 的宽度,断层长度为20km。对断层的杨氏模量、泊松比、断层附近区域的边界受力大小、两个主应力的比例以及主应力的方向的改变对计算结果的影响做了讨论。最后通过上述几组计算结果建立了一个与前人结果较为符合的巴姆地震参数模型。通过计算,一方面揭示了应力状态和物性参数的影响,另一方面也说明 Arg-e-Bam 隐伏断层是此次地震的发震断层这一结论是可信的。     

井地直流电法三维数值模拟中若干问题研究

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。     

洞庭湖泥沙淤积数值模拟模式

在确保沙量守恒的非恒定流、非均匀沙的二维水沙数学模型的基础上,通过洞庭湖水沙输运和河床变形计算实践,比较了两种冲泻质与床沙质的转化模式,提出了包括冲泻质在内的泥沙淤积模式和水沙动边界计算模式.从洞庭湖泥沙输运和河床变形计算的结果看,洞庭湖泥沙淤积数值模式具有较好的稳定性和较高的模拟精度,说明所建湖泊泥沙模型的良好性能和具有推广运用的前景.     

长板短桩工法加固机理数值研究

长板短桩工法是采用长的塑料排水板（或砂井）与短的水泥土搅拌桩联合的新型地基处理方法。作者依据淮盐高速公路盐城段第十四合同段K97＋461-K97＋561的试验研究，在比较了各种数值模拟计算方法的基础上，选用复合地基单元体方法，利用原位测试和土工实验得到的指标作为计算所使用的参数，以现场量测结果作为数值模型检验的依据，对长板短桩工法复合地基数值模型的合理性作进一步的探索和验证。     

大欧拉角的空间直角坐标转换方法探讨

首先指出建立使用大旋转角的空间直角坐标转换模型的必要性，接着推导出以罗德里格矩阵为旋转矩阵的大欧拉角空间直角坐标转换的数学模型，然后以布尔沙模型为例、以Visual C#．NET2003编制了计算软件。该算法具有性能稳定、计算收敛速度快、精确度高、转换参数初始值可以为任意值的优点。     

密肋复合墙体非线性数值分析及抗剪承载力研究

结合课题组前期的研究成果,通过对墙体中不同材料的本构模型、墙体中裂缝、多种材料的联结以及墙体边界条件的研究与探讨,运用通用有限元程序ANSYS对密肋复合墙体进行非线性有限元分析;并在试验结果与有限元模拟的基础上,依据极限平衡理论,建立起来复合墙体的斜截面抗剪极限承载力公式。理论分析与试验研究表明:复合墙体非线性有限元模型能够较好地模拟墙体内力和变形发展的全过程,描述裂缝的形成和扩展,对于实际工程和科学研究具有足够的精度;墙体斜截面抗剪承载力公式与实测值吻合较好,具有一定的精度,满足工程设计需要。     

Fine upper crustal structure of Jiashi strong earthquake swarm region in Xinjiang inferred from high resolution seismic refraction profile data

Introduction Jiashi-Artux area in southwest Xinjiang is one of the most active earthquake provinces at pre-sent in Chinese mainland. In the last century, about 3/4 strong earthquakes in Chinese mainland hit this area, and especially from January 21 to April 16 in 1997, 7 earthquakes with the magnituderanging from 6.0 to 6.9 occurred in a very small area of 9 km×18 km near Jiashi (ZHU et al, 1998). It has never taken place before in Chinese mainland that a series of strong earthquakes shoo…