An explicit representation of steady state response of a beam on an elastic foundation to moving harmonic line loads

A closed‐form deflection response of a beam rest is presented in this paper using the integral transform method. The theory of linear partial differential equations is used to represent the deflection of beam subjected to a moving harmonic line load in integration form. The solution is finally carried out using the inverse Fourier transform. To evaluate the integration analytically, poles of the integrand are identified with the help of algebraic equation theory. Residue theorem is then utilized to represent the integration as a contour integral in the complex plane. Closed‐form deflections and numerical results are provided for different combinations of load frequency and velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

Multistep matched-field inversion for broad-band data from ASIAEX2001

This paper discusses the results of geoacoustic inversion carried out using explosive charge data from the Asian Seas International Acoustic Experiment (ASIAEX) East China Sea (ECS) Experiment. A multifrequency incoherent matched-field inversion processor and a genetic algorithm (GA) are used for the inversion. A multistep matched field inversion approach is presented, which makes use of the varying sensitivities of wave fields at various frequencies to reduce the inversion problem into a sequence of smaller inversions with fewer unknowns to estimate at each stage. Different parameters are estimated using data at different frequencies according to their sensitivities. Inversion results for different areas in the ECS region are summarized and compared with core data.     

The global joule heat production rate and the AE index

The accuracy of the AE index as a measure of the joule heat production rate is examined for a typical substorm event on 18 March 1978 by estimating the global joule heat production rate, as a function of time, using data obtained from the IMS six meridian chains. In spite of the fact that the AE index had had an initial slow growth which was followed by a rapid growth, the joule heat production rate attained a high level during the slow growth and thus their variations were considerably different from those of the AE index. Therefore, although the AE index is statistically linearly related to the global joule heat production rate, one should be cautious in assuming that details of time variations of the AE index during individual events represent those of the joule heat production rate.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Himalayan magmatism and porphyry copper–molybdenum mineralization in the Yulong ore belt, East Tibet

Summary ?The NW–SE-trending Yulong porphyry Cu–Mo ore belt, situated in the Sanjiang0 area of eastern Tibet, is approximately 400 km long and 35 to 70 km wide. Complex tectonic and magmatic processes during the Himalayan epoch have given rise to favorable conditions for porphyry-type Cu–Mo mineralization. Porphyry masses of the Himalayan epoch in the Yulong ore belt are distributed in groups along regional NW–SE striking tectonic lineaments. They were emplaced mainly into Triassic and Lower Permian sedimentary-volcanic rocks. K–Ar und U–Pb isotopic datings give an intrusion age range of 57–26 Ma. The porphyries are mainly of biotite monzogranitic and biotite syenogranitic compositions. Geological and geochemical data indicate that the various porphyritic intrusions in the belt had a common or similar magma source, are metaluminous to peraluminous, Nb–Y–Ba-depleted, I-type granitoids, and belong to the high-K calc-alkaline series. Within the Yulong subvolcanic belt a number of porphyry stocks bear typical porphyry type Cu–Mo alteration and mineralization. The most prominent porphyry Co–Mo deposits include Yulong, Malasongduo, Duoxiasongduo, Mangzong and Zhanaga, of which Yulong is one of the largest porphyry Cu (Mo) deposits in China with approximately 8 × 10⁶ tons of contained Cu metal. Hydrothermal alteration at Yulong developed around a biotite–monzogranitic porphyry stock that was emplaced within Upper Triassic limestone, siltstone and mudstone. The earliest alteration was due to the effects of contact metamorphism of the country rocks and alkali metasomatism (potassic alteration) within and around the porphyry body. The alteration of this stage was accompanied by a small amount of disseminated and veinlet Cu–Mo sulfide mineralization. Later alteration–mineralization zones form more or less concentric shells around the potassic zone, around which are distributed a phyllic or quartz–sericite–pyrite zone, a silicification and argillic zone, and a propylitic zone. Fluid inclusion data indicate that three types of fluids were involved in the alteration–mineralization processes: (1) early high temperature (660–420 °C) and high salinity (30–51 wt% NaCl equiv) fluids responsible for the potassic alteration and the earliest disseminated and/or veinlet Cu–Mo sulfide mineralization; (2) intermediate unmixed fluids corresponding to phyllic alteration and most Cu–Mo sulfide mineralization, with salinities of 30–50 wt% NaCl equiv and homogenization temperatures of 460–280 °C; and (3) late low to moderate temperature (300–160 °C) and low salinity (6–13 wt% NaCl equiv) fluids responsible for argillic and propylitic alteration. Hydrogen and oxygen isotopic studies show that the early hydrothermal fluids are of magmatic origin and were succeeded by increasing amounts of meteoric-derived convective waters. Sulfur isotopes also indicate a magmatic source for the sulfur in the early sulfide mineralization, with the increasing addition of sedimentary sulfur outward from the porphyry stock. Received August 29, 2001; revised version accepted May 1, 2002 Published online: November 29, 2002     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

The distribution of QSO absorption lines and its implications

The distribution of L absorption lines has been investigated in the fractal scheme. It is found that (1) the L absorption clouds distribute completely different from that of galaxies; (2) the L absorption clouds are anti-associated with galaxies and quasars. These results may imply that there are two kinds of objects formed by different processes of clustering. This is favourable for the cosmic-string theory on the formation of large-scale structure of the Universe. In the string model, the objects can be divided into two kinds according to their clustering with or without string loops as their initial density perturbation.On leave from the Center of Astrophysics, University of Science and Technology of China, Hefei, P.R. China. Y. Chu is a research fellow of Alexander von Humboldt-Foundation.     

Timing of the Nihewan formation and faunas

Magnetostratigraphic dating of the fluvio-lacustrine sequence in the Nihewan Basin, North China, has permitted the precise timing of the basin infilling and associated Nihewan mammalian faunas. The combined evidence of new paleomagnetic findings from the Hongya and Huabaogou sections of the eastern Nihewan Basin and previously published magnetochronological data suggests that the Nihewan Formation records the tectono-sedimentary processes of the Plio-Pleistocene Nihewan Basin and that the Nihewan faunas can be placed between the Matuyama-Brunhes geomagnetic reversal and the onset of the Olduvai subchron (0.78-1.95 Ma). The onset and termination of the basin deposition occurred just prior to the Gauss-Matuyama geomagnetic reversal and during the period from the last interglaciation to the late last glaciation, respectively, suggesting that the Nihewan Formation is of Late Pliocene to late Pleistocene age. The Nihewan faunas, comprising a series of mammalian faunas (such as Maliang, Donggutuo, Xiaochangliang, Banshan, Majuangou, Huabaogou, Xiashagou, Danangou and Dongyaozitou), are suggested to span a time range of about 0.8-2.0 Ma. The combination of our new and previously published magnetostratigraphy has significantly refined the chronology of the terrestrial Nihewan Formation and faunas.