Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

计算机模拟引擎及其在铀成矿过程模拟中的应用研究

本文介绍了计算机模拟引擎的设计与开发 ,铀成矿过程模拟的实施步骤 ,并对其在伊犁盆地南缘的应用作了探讨。     

长江口深水航道30年空间特征变化研究

根据1973-2003年6幅长江口海图资料,建立不同时期长江口深水航道水下数字高程模型, 对航道治理前后30 a的冲淤状况、纵剖面水深变化以及空间模拟等进行研究. 结果表明,30 a来深水航道冲淤交替出现,总变化趋势是冲刷.航道冲淤除与 水动力条件有关外,与治理工程也有密切联系.长江口深水航道治理一期工程结束后,水深明显加深,坡度比变大,主航道呈现平顺的“U”形河道景观.     

Discretization error for the superposition of Gaussian beams

Summary. A high-frequency asymptotic integral expansion of a time-harmonic wavefield into Gaussian beams was derived in a previous paper by Klimeš. The discretization error caused by replacing this integral superposition by a discrete summation of Gaussian beams is estimated in this paper.     

Wind tunnel simulation experiment and investigation on the electrification of sandstorms

~~Wind tunnel simulation experiment and investigation on the electrification of sandstorms~~     

Granulometric and chemical composition of the Sava River sediments upstream and downstream of the Krsko nuclear power plant

Total concentrations of 13 elements (K, Ca, Ti, Cr, Mn, Fe, Cu, Zn, Rb, Sr, Y, Zr, Pb) in the size-fractionated Sava River sediments upstream and downstream of the Krsko nuclear power plant together with metal speciation within bulk sediment have been investigated. Trace metals generally increase with decreasing particle size, however, because of entrapment of organic matter in the 0.63–1 mm fraction, concentrations in the coarser sediment fraction are higher than expected. Exchangeable Pb, Zn, Cu, Mn, Cr and Fe are generally found to represent a negligible fraction of the total metal concentration of the bulk sediment. Seasonal variations of the Pb, Zn and Cu concentrations in the <0.5 mm fraction reflect decreased values during the spring period. Heavy metal concentrations in the 2003 waste water discharges from the Krsko nuclear power plant released into the Sava River were much lower than their maximum allowed values. Combined rubidium and organic matter normalization of the Zn, Pb and Cu concentrations, which was applied on the minus 0.063 mm fraction, indicated three potential sources of contaminants.     

锦屏电站坝厂区卸荷岩体力学参数研究

选择具有代表性的区块,根据结构面的性状和规模拟定模型层次和尺寸,建立岩体数值仿真模型,通过计算机数值仿真,采用卸荷非线性力学理论对锦屏电站坝厂区岩体进行卸荷岩体宏观力学特性及参数的数值分析;同时,用RMR法和BP神经网络法对岩体宏观力学特性及参数进行分析,综合确定岩体力学参数,为设计部门提供设计参考依据.     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.