Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Rayleigh 阻尼作用下斜拉桥扭转振动

在侧向风力或侧向水平地源力作用下,斜拉桥容易发生扭转振动。本文采用非线性Rayleigh阴尼,它能反映在低速振动时有激发而加速振动,至一定高速后又衰减的完整振动过程。对非线性的统辖方程及定解条件,先使之无量纲化,再按非线性项的系数这个小参数展开未知函数,得到线性化的各级近似方程。把待求函数展成梁轴坐标的Fourier级数,其系数为待定的时间函数,应用Fourier级数的正交性,得到求解未知时间函数的二阶常微分方程,完全类似单自由度质点在线性粘滞性阻尼作用下的统辖方程,其解为已知。     

涡旋诱发重联模型（Ⅰ）——动力学特性

本文首先从理论上系统地研究了涡旋诱发重联的基本动力学特性.包括:涡旋诱发重联的物理机制;Kelvin-Helmholtz（K-H）不稳定性和撕裂模不稳定性耦合过程;磁流体涡管的基本特性;不同尺度涡旋间的相互作用,动能和磁能间的转化及准稳态时流线、磁力线、等涡度线和等电流密度线的相似性等.用二维MHD数值模拟方法研究了上述一些特性,并与理论分析结论进行了比较.结果表明,模拟结果证实了理论推论.涡旋诱发重联是产生局部磁场重联的重要机制之一,该模型在磁顶区的应用将另文发表.     

新丰江水库区1989年11月26日4.5级地震的主要特征

本文总结分析了发生在新丰江水库老震区，距1962年6.1级主震27年的1次4.5级地震的活动和特点，认为此次4.5级地震及其伴生的小震群仍系新丰江水库诱发地震序列的中晚期强余震。     

地震反应分析中的主时段

在计算反应谱及土层反应时,发现反应谱的形态和大小只受时程的主要部分控制,由此给出主时段的定义和确定主时段的方法。主时段的应用不仅可以大大缩短计算时间,而且能够成为持时的一种有物理意义的度量。     

Time effects on dynamic soil properties under random excitation conditions

Two types of clayey soils, a kaolinite and a bentonite, were tested using a resonant column apparatus under random excitation conditions. The concept of root mean square (rms) strain was utilized for the purpose of strain calculations during random loading. The conventional estimator of the transfer function was used for random vibration analysis. The effect of confinement duration (at a constant pressure) on dynamic soil properties, namely damping and shear modulus, was evaluated. The results indicate that for both cohesive soils, the effect of time was less pronounced during random vibration than sinusoidal loading at the same rms strain. This effect is however more pronounced when the peak shearing strain of sinusoidal loading is considered. Furthermore, time effects were more pronounced at low strain levels than at high strain levels.     

Application of stress-pore pressure coupling theory for porous media to the Xinfengjiang reservoir earthquakes

Theory of the coupling of stress-pore pressure in the saturated, elastic porous media is used in the study of the formation mechanism of the Xinfengjiang reservoir-induced earthquakes. Based on the results, it is believed that compared with the mechanism of additional stress in the vicinity of the reservoir, the mechanism of the coupling of additional stress and pore pressure may be more well-founded for the occurrence of reservoir-induced earthquakes.     

Earthquakes in the Czech Republic and Surrounding Regions in 1995–1999

In the time span from January 1995 to December 1999 the Czech National Seismological Network (CNSN), consisting of ten permanent digital broadband stations, several local networks and two data centers, detected and recorded 9530 regional natural seismic events, 27 greater than magnitude 2. Most of these events were located by the Czech Seismological Service (CSS), and the most prominent of them were analyzed in detail. A large number of quarry blasts were recorded as well but were not included in the analysis. We provide basic information on the configuration of the CNSN and on the way of routine data processing employed by the CSS in this paper. The over-all regional seismicity monitored by the CNSN in 1995–1999 is briefly reviewed. The main results of observations and evaluation of the local (NW-Bohemia/Vogtland, South Bohemia, Sudeten) and induced (Kladno, Píbram, Upper Silesia, Lubin/Poland) seismic activity within this period are presented in a condensed form. Finally, a summary on macroseismic observations on the territory of the Czech Republic in 1995–1999 is also presented.