软土亚塑性模型与一维非线性固结理论研究

With the great development of Tianjing New Coastal District economy, people need more land to build and live. Land subsidence, which is caused by its special engineering geological conditions, has restricted the further development in the district. Soft soil consolidation is main factor of land subsidence ;thus , on the basis of consolidation theory, the paper make further study on soft soils one-dimension nonlinear consolidation which contains two parts:(1) the nonlinear consolidation of permeability coefficient and compressibility coefficient changing with time and depth, which means real one-dimension nonlinear consolidation;(2) the non-homogeneous consolidation of permeability coefficient and compressibility coefficient only changing with depth. Firstly, nonlinear characteristics of soft soils are elaborated. Hypoplastic theory is introduced to establish a modified soft soils nonlinear constitutive model; the nonlinear governing equation of compressibility coefficient is built, and the nonlinear characteristics of compressibility coefficient are analyzed. Secondly, Considering Load Fluctuation and soil thickness changing ,the consolidation characteristics of single layer is discussed in the paper; meanwhile, on the basis of the Davis and Raymond’s hypothesis and single layer nonlinear consolidation equation, the doubled-layer one-dimension nonlinear consolidation equation is also derived. The solution of the equation is obtained by analytical method, and the consolidation characteristics of doubled-layer soft soil nonlinear theory is also analyzed. Finally, based on assumption that permeability coefficient and compressibility coefficient is varying along depth, single layer soil one-dimension non-homogeneous consolidation differential equation is derived; and the approximate solution is obtained. Furthermore, the single layer non-homogeneous consolidation is extended to double layer non-homogeneous consolidation theory. By using parabolic differential scheme, the matrix equation is established; and the solution of the matrix equation is obtained by chase method. Consolidation characteristics of soil soft single (double) layer non-homogeneous consolidation theory and Terzaghi’s theory are also discussed.

铂族元素分析方法研究及其在有关地质问题中的应用

Characterization of Platinum Group Elements (PGE) has been applied to earth, space and environmental sciences. However, all these applications are based on a basic prerequisite, i.e. their concentration or ratio in the research objects can be accurately and precisely determined. In fact, development in these related studies is a great challenge to the analytical chemistry of the PGE because their content in the geological sample (non-mineralized) is often extremely low, range from ppt (10~(-12)g/g) to ppt (10~(-9)g/g). Their distribution is highly heterogeneous, usually concentrating in single particle or phase. Therefore, the accurate determination of these elements remains a problem in analytical chemistry and it obstructs the research on geochemistry of PGE. A great effort has been made in scientific community to reliable determining of very low amounts of PGE, which has been focused on to reduce the level of background in used reagents and to solve probable heterogeneity of PGE in samples. Undoubtedly, the fire-assay method is one of the best ways for solving the heterogeneity, as a large amount of sample weight (10-50g) can be hold. This page is mainly aimed at development of the methodology on separation, concentration and determination of the ultra-trace PGE in the rock and peat samples, and then they are applied to study the trace of PGE in ophiolite suite, in Kudi, West Kunlun and Tunguska explosion in 1908. The achievements of the study are summarized as follows: 1. A PGE lab is established in the Laboratory of Lithosphere Tectonic Evolution, IGG, CAS. 2. A modified method of determination of PGE in geological samples using NiS Fire-Assay with inductively coupled plasma-mass spectrometry (ICP-MS) is set up. The technical improvements are made as following: (1) investigating the level of background in used reagents, and finding the contents of Au, Pt and Pd in carbonyl nickel powder are 30, 0.6 and 0.6ng/g, respectively and 0.35, 7.5 and 6.4ng, respectively in other flux, and the contents of Ru, Rh, Os in whole reagents used are very low (below or near the detection limits of ICP-MS); (2) measuring the recoveries of PGE using different collector (Ni+S) and finding 1.5g of carbonyl nickel is effective for recovering the PGE for 15g samples (recoveries are more than 90%), reducing the inherent blank value due to impurities reagents; (3) direct dissolving nickel button in Teflon bomb and using Te-precipitation, so reducing the loss of PGE during preconcentration process and improving the recoveries of PGE (above 60% for Os and 93.6-106.3% for other PGE, using 2g carbonyl nickel); (4) simplifying the procedure of analyzing Osmium; (5)method detection limits are 8.6, 4.8, 43, 2.4, 82pg/g for 15g sample size ofRu, Rh, Pd, Ir, Pt, respectively. 3. An analytical method is set up to determine the content of ultra-trace PGE in peat samples. The method detection limits are 0.06, 0.1, 0.001, 0.001 and 0.002ng/mL for Ru, Rh, Pd, Ir and Pt, respectively. 4. Distinct anomaly of Pd and Os are firstly found in the peat sampling near the Tunguska explosion site, using the analytical method. 5. Applying the method to the study on the origin of Tunguska explosion and making the following conclusions: (1) these excess elements were likely resulted from the Tunguska Cosmic Body (TCB) explosion of 1908. (2) The Tunguska explosive body was composed of materials (solid components) similar to C1 chondrite, and, most probably, a cometary object, which weighed more than 10~7 tons and had a radius of more than 126 m. 6. The analysis method about ultra-trace PGE in rock samples is successfully used in the study on the characteristic of PGE in Kudi ophiolite suite and the following conclusions are made: (1) The difference of the mantle normalization of PGE patterns between dunite, harzburgite and lherzolite in Kudi indicates that they are residual of multi-stage partial melt of the mantle. Their depletion of Ir at a similar degree probably indicates the existence of an upper mantle depleted Ir. (2) With the evolution of the magma produced by the partial melt of the mantle, strong differentiation has been shown between IPGE and PPGE; and the differentiation from pyroxenite to basalt would have been more and more distinct. (3) The magma forming ophiolite in Kudi probably suffered S-saturation process.

水库滑坡复活机理及其发展趋势预测研究

Various numbers of ancient landslides of various scales are frequently distributed on both banks of reservoirs, especially large reservoirs, both in China and abroad. During inundation and operation of theses reservoirs, some of the landslides are reactivated, which caused losses of people's lives and properties to various extents, some even disasters. Systematic studies are, however, very few on the reservoir-induced reactivation mechanism and development tendency prediction. Based on investigation of reservoir-induced reactivation phenomena of ancient landslides and relevant existing research problems, a systematic study is carried out on the field identification, induced reactivation mechanism, development tendency prediction, risk decision-making and treatment of reservoir-related ancient landslides, through analysis of large numbers of engineering geological investigation results, scientific experimental and research results, in combination with prevention and treatment practices of reservoir-related landslides both in China and abroad, and a series of research results have been obtained. 1. On the basis of study of the distribution features, genesis mechanism of ancient landslides on river banks, a set of scientific methods are summarized on field identification of ancient landslides, and a significant method named "lithologic sequence method" or "indicator layer method", is proposed, which is proved to be very useful. 2. A detail study is made on the reservoir-induced hydraulic effects and material mechanic effects (or softening effects) on the ancient landslide through model and case studies, which concludes that the magnitude and properties of reservoir-induced hydraulic effects are related to the shapes of sliding planes, water content and permeability of landslide materials and variation rate and magnitude of reservoir levels; the magnitude of material mechanic effects are related to the material composition (including mineral composition and grain size), natural water content and saturation state of sliding zones. Also a sensitive analysis is made on the factors that are related to the stabilities of the landslides, which indicate that the stability of a landslide is more sensitive to the groundwater head h_w in the slides and the inner friction angleψof sliding zones than others. 3. The joint inducing mechanism of rainfall and reservoir is also discussed in the paper through model analysis and case study, which proves that reservoir inundation increases firstly the sensitivity of a landslides to rainfall through reduction of its stability or cracking deformation which will increase the rainfall infiltration to the slide body, and then rainfall triggers reactivation or intensifies the reservoir-induced deformation of a landslide. 4. Based on rheologic test results of sliding zones of several reservoir-related ancient landslides, the rheologic characteristics of sliding zones have been discussed in detail and several typical rheologic models have been set up, which well explains the dynamic process of slide deformation. The response types to reservoir inundation and development tendency of reservoir -related ancient landslide are discussed in the paper based on field investigation results. And prediction methods for reservoir-related landslides have been studied based on the Mate-Synthetic principle of quantitative and qualitative analysis, as well as combination of computation and internal mechanism analysis, and a rheologic analytical method is proposed which is proved very useful for prediction of the landslide development tendency. 6. In disaster-prevention and treatment of reservoir-related landslides, risk decision-making has been proved very significant both in engineering and economics. Based on the practices in disaster-prevention and treatment of reservoir-related landslides both in China and abroad, the disaster-prevention risk decision-making for reservoir-related landslides has been proposed in terms of philosophy, methods and procedures, and well put into practice. A summary is also made through case study of the experiences of treatment of reservoir-related landslides both in China and abroad in terms of principle, methods and technical lines. 7 A detail study is made as a case study of the reactivated Maoping ancient landslide on the left bank of the Geheyan Reservoir on Qingjiang river in Hubei province, China, including its field identification features, reservoir-induced reactivation characteristics and mechanism, development tendency prediction and proposed counter measures based on risk analysis.

Prediction of GC retention values under various column temperature conditions from temperature programmed data

A numerical approach has been developed for the correlation of retention limes (total retention lime) with temperature in gas chromatography, which allows the calculation of retention parameters including retention index from data acquired under two or more different temperature program conditions. By using this procedure the optimization of temperature condition can be further achieved, especially when a temperature-programmed run is the most suitable mode in the preliminary development of an analytical method for the analysis of an unknown sample.

Analytical study of the layered XY model with VCE method

Motivated by experiments on liquid-crystal films, we study the development of specific heat anomaly of finite layer system. With the VCE method, we introduce the strong surface interaction into the layered XY model and get the results of the forth-order analytical expansion. The results show that when the strong surface interaction becomes strong enough, the order trend defeats the quantum noise and the specific heat peak moves abnormally to the high temperature with the number of layers decreasing.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

Analytical Methods of Residual Strength in Injecting Water Pipeline

In this paper, common criterions about residual strength evaluation at home and abroad are generalized and seven methods are acquired, namely ASME-B31G, DM, Wes-2805-97, CVDA-84, Burdekin, Irwin and J integral methods. BP neural network are Combined with Genetic Algorithm (GA) named by modified BP-GA methods to successfully predict residual strength and critical pressure of injecting water, corrosion pipelines. Examples are shown that calculation results of every kind of method have great difference and calculating values of Wes-2805-97 criterion, ASME-B31G criterion, CVDA-84 criterion and Irwin fracture mechanics model are conservative and higher than, those of J integral methods while calculating values of Burdiken model and DM fracture mechanics model are dangerous and less than those of J integral methods and calculating values of modified BP-GA methods are close and moderate to those of J integral methods. Therefore modified BP-GA methods and J integral methods are considered better methods to calculate residual strength and critical pressure of injecting water corrosion pipelines

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

Approximate analytical description for fundamental-mode fields of graded-index fibers: Beyond the Gaussian approximation

An approximate analytical description for fundamental-mode fields of graded-index fibers is explicitly presented by use of the power-series expansion method, the maximum-value condition at the fiber axis, the decay properties of fundamental-mode fields at large distance from the fiber axis, and the approximate modal parameters U obtained from the Gaussian approximation. This analytical description is much more accurate than the Gaussian approximation and at the same time keep the simplicity of the latter. As two special examples, we present the approximate analytical formulas for the fundamental-mode fields of a step profile fiber and a Gaussian profile fiber, and we find that they are both highly accurate in the single-mode range by comparing them with the corresponding exact solutions.