花岗岩风化土微观结构及其工程地质特性研究

Saprolite is the residual soil resulted from completely weathered or highly weathered granite and with corestones of parent rock. It is widely distributed in Hong Kong. Slope instability usually happens in this layer of residual soil and thus it is very important to study the engineering geological properties of Saprolite. Due to the relic granitic texture, the deformation and strength characteristics of Saprolite are very different from normal residual soils. In order to investigate the effects of the special microstructure on soil deformation and strength, a series of physical, chemical and mechanical tests were conducted on Saprolite at Kowloon, Hong Kong. The tests include chemical analysis, particle size analysis, mineral composition analysis, mercury injection, consolidation test, direct shear test, triaxial shear test, optical analysis, SEM & TEM analysis, and triaxial shear tests under real-time CT monitoring.Based on the testing results, intensity and degree of weathering were classified, factors affecting and controlling the deformation and strength of Saprolite were identified, and the interaction between those factors were analyzed.The major parameters describing soil microstructure were introduced mainly based on optical thin section analysis results. These parameters are of importance and physical meaning to describe particle shape, particle size distribution (PSD), and for numerical modeling of soil microstructure. A few parameters to depict particle geometry were proposed or improved. These parameters can be used to regenerate the particle shape and its distribution. Fractal dimension of particle shape was proposed to describe irregularity of particle shapes and capacity of space filling quantitatively. And the effect of fractal dimension of particle shape on soil strength was analyzed. At the same time, structural coefficient - a combined parameter which can quantify the overall microstructure of rock or soil was introduced to study Saprolite and the results are very positive. The study emphasized on the fractal characteristics of PSD and pore structure by applying fractal theory and method. With the results from thin section analysis and mercury injection, it was shown that at least two fractal dimensions Dfl(DB) and Df2 (Dw), exist for both PSD and pore structure. The reasons and physical meanings behind multi-fractal dimensions were analyzed. The fractal dimensions were used to calculate the formation depth and weathering rate of granite at Kowloon. As practical applications, correlations and mathematical models for fractal dimensions and engineering properties of soil were established. The correlation between fractal dimensions and mechanical properties of soil shows that the internal friction angle is mainly governed by Dfl 9 corresponding to coarse grain components, while the cohesion depends on Df2 , corresponding to fine grain components. The correlations between the fractal dimension, friction angle and cohesion are positive linear.Fractal models of PSD and pore size distribution were derived theoretically. Fragmentation mechanism of grains was also analyzed from the viewpoint of fractal. A simple function was derived to define the theoretical relationship between the water characteristic curve (WCC) and fractal dimension, based on a number of classical WCC models. This relationship provides a new analytical tool and research method for hydraulic properties in porous media and solute transportation. It also endues fractal dimensions with new physical meanings and facilitates applications of fractal dimensions in water retention characteristics, ground water movement, and environmental engineering.Based on the conclusions from the fractal characteristics of Saprolite, size effect on strength was expressed by fractal dimension. This function is in complete agreement with classical Weibull model and a simple function was derived to represent the relationship between them.In this thesis, the phenomenon of multi-fractal dimensions was theoretically analyzed and verified with WCC and saprolite PSD results, it was then concluded that multi-fractal can describe the characteristics of one object more accurately, compared to single fractal dimension. The multi-fractal of saprolite reflects its structural heterogeneity and changeable stress environment during the evolution history.

Mechanics of metallic foams: A fractal approach

A fractal approach was proposed to investigate the meso structures and size effect of metallic foams: For a series At foams of different relative densities, the information dimension method was applied to measure meso structures. The generalized sierpinski carpet was introduced to map the meso structures of the foam according to specific dimension. The results show that the fractal-based model can not only reveal the variation of yield strength with specimen size, but also bridge the meso structures and mechanical proper-ties of Al foams directly. Key words: metallic foams; fractal; size effect; meso structures

An extension of the Quasicontinuum Treatment of Multiscale Solid Systems to Nonzero Temperature

Covering the solid lattice with a finite-element mesh produces a coarse-grained system of mesh nodes as pseudoatoms interacting through an effective potential energy that depends implicitly on the thermodynamic state. Use of the pseudoatomic Hamiltonian in a Monte Carlo simulation of the two-dimensional Lennard-Jones crystal yields equilibrium thermomechanical properties (e.g., isotropic stress) in excellent agreement with ``exact'' fully atomistic results.

Multiscale Coupling in Complex Mechanical Systems

Multiscale coupling attracts broad interests from mechanics, physics and chemistry to biology. The diversity and coupling of physics at different scales are two essential features of multiscale problems in far-from-equilibrium systems. The two features present fundamental difficulties and are great challenges to multiscale modeling and simulation. The theory of dynamical system and statistical mechanics provide fundamental tools for the multiscale coupling problems. The paper presents some closed multiscale formulations, e.g., the mapping closure approximation, multiscale large-eddy simulation and statistical mesoscopic damage mechanics, for two typical multiscale coupling problems in mechanics, that is, turbulence in fluids and failure in solids. It is pointed that developing a tractable, closed nonequilibrium statistical theory may be an effective approach to deal with the multiscale coupling problems. Some common characteristics of the statistical theory are discussed.

Permeability of Fractal Porous Media by Monte Carlo Simulations

The permeability of the fractal porous media is simulated by Monte Carlo technique in this work. Based oil the fractal character of pore size distribution in porous media, the probability models for pore diameter and for permeability are derived. Taking the bi-dispersed fractal porous media as examples, the permeability calculations are performed by the present Monte Carlo method. The results show that the present simulations present a good agreement compared with the existing fractal analytical solution in the general interested porosity range. The proposed simulation method may have the potential in prediction of other transport properties (such as thermal conductivity, dispersion conductivity and electrical conductivity) in fractal porous media, both saturated and unsaturated.

Multiscale Treatment of Thin-Film Lubrication

A multiscale technique that combines an atomistic description of the interfacial (near) region with a coarse-grained (continuum) description of the far regions of the solid substrates is proposed. The new hybrid technique, which represents an advance over a previously proposed dynamically-constrained hybrid atomistic-coarse-grained treatment (Wu et al.J. Chem. Phys., 120, 6744, 2004), is applied to a two-dimensional model tribological system comprising planar substrates sandwiching a monolayer film. Shear–stress profiles (shear stress versus strain) computed by the new hybrid technique are in excellent agreement with “exact” profiles (i.e. those computed treating the whole system at the atomic scale).

Multiscale Coupling: Challenges and Opportunities

Multiscale coupling is ubiquitous in nature and attracts broad interests of scientists from mathematicians, physicists, machinists, chemists to biologists. However, much less attention has been paid to its intrinsic implication. In this paper, multiscale coupling is introduced by studying two typical examples in classic mechanics: fluid turbulence and solid failure. The nature of multiscale coupling in the two examples lies in their physical diversities and strong coupling over wide-range scales. The theories of dynamical system and statistical mechanics provide fundamental methods for the multiscale coupling problems. The diverse multiscale couplings call for unified approaches and might expedite new concepts, theories and disciplines.

The Pressure Transient Analysis Of Deformation Of Fractal Medium

The assumption of constant rock properties in pressure-transient analysis of stress-sensitive reservoirs can cause significant errors in the estimation of temporal and spatial variation of pressure. In this article, the pressure transient response of the fractal medium in stress-sensitive reservoirs was studied by using the self-similarity solution method and the regular perturbation method. The dependence of permeability on pore pressure makes the flow equation strongly nonlinear. The nonlinearities associated with the governing equation become weaker by using the logarithm transformation. The perturbation solutions for a constant pressure production and a constant rate production of a linear-source well were obtained by using the self-similarity solution method and the regular perturbation method in an infinitely large system, and inquire into the changing rule of pressure when the fractal and deformation parameters change. The plots of typical pressure curves were given in a few cases, and the results can be applied to well test analysis.

Rapidly Infrared-Assisted Cooperatively Self-Assembled Highly Ordered Multiscale Porous Materials

In this paper, cooperative self-assembly (CSA) of colloidal spheres with different sizes was studied. It was found that a complicated jamming effect makes it difficult to achieve an optimal self-assembling condition for construction of a well-ordered stacking of colloidal spheres in a relatively short growth time by CSA. Through the use of a characteristic infrared (IR) technique to significantly accelerate local evaporation on the growing interface without changing the bulk growing environment, a concise three-parameter (temperature, pressure, and IR intensity) CSA method to effectively overcome the jamming effect has been developed. Mono- and multiscale inverse opals in a large range of lattice scales can be prepared within a growth time (15-30 min) that is remarkably shorter than the growth times of several hours for previous methods. Scanning electron microscopy images and transmittance spectra demonstrated the superior crystalline and optical qualities of the resulting materials. More importantly, the new method enables optimal conditions for CSA without limitations on sizes and materials of multiple colloids. This strategy not only makes a meaningful advance in the applicability and universality of colloidal crystals and ordered porous materials but also can be an inspiration to the self-assembly systems widely used in many other fields, such as nanotechnology and molecular bioengineering.

Fractal geometry model for effective thermal conductivity of three-phase porous media

An approximate model, a fractal geometry model, for the effective thermal conductivity of three-phase/unsaturated porous media is proposed based on the thermal-electrical analogy technique and on statistical self-similarity of porous media. The proposed thermal conductivity model is expressed as a function of porosity (related to stage n of Sierpinski carpet), ratio of areas, ratio of component thermal conductivities, and saturation. The recursive algorithm for the thermal conductivity by the proposed model is presented and found to be quite simple. The model predictions are compared with the existing measurements. Good agreement is found between the present model predictions and the existing experimental data. This verifies the validity of the proposed model. (C) 2004 American Institute of Physics.

Fractal features of oscillatory convection in the half-floating zone

Unsteady and two-dimensional numerical simulation is applied to study the transition process from steady convection to turbulence via subharmonic bifurcation in thermocapillary convection of a liquid bridge in the half-floating zone. The results of numerical tests show clearly the fractal structure of period-doubling bifurcations, and frequency-locking at f/4, f/8, f/16 with basic frequency f is observed with increasing temperature difference. The Feigenbaum universal constant is given by the present paper as delta(4) = 4.853, which can be compared with the theoretical value 4.6642016.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

The autocorrelation of speckles in deep Fresnel diffraction region and characterizations of random self-affine fractal surfaces

This paper studies the correlation properties of the speckles in the deep Fresnel diffraction region produced by the scattering of rough self-affine fractal surfaces. The autocorrelation function of the speckle intensities is formulated by the combination of the light scattering theory of Kirchhoff approximation and the principles of speckle statistics. We propose a method for extracting the three surface parameters, i.e. the roughness w, the lateral correlation length xi and the roughness exponent alpha, from the autocorrelation functions of speckles. This method is verified by simulating the speckle intensities and calculating the speckle autocorrelation function. We also find the phenomenon that for rough surfaces with alpha = 1, the structure of the speckles resembles that of the surface heights, which results from the effect of the peak and the valley parts of the surface, acting as micro-lenses converging and diverging the light waves.

The accuracy of Kirchhoff's approximation in describing the far field speckles produced by random self-affine fractal surfaces

Based on the rigorous formulation of integral equations for the propagations of light waves at the medium interface, we carry out the numerical solutions of the random light field scattered from self-affine fractal surface samples. The light intensities produced by the same surface samples are also calculated in Kirchhoff's approximation, and their comparisons with the corresponding rigorous results show directly the degree of the accuracy of the approximation. It is indicated that Kirchhoff's approximation is of good accuracy for random surfaces with small roughness value w and large roughness exponent alpha. For random surfaces with larger w and smaller alpha, the approximation results in considerable errors, and detailed calculations show that the inaccuracy comes from the simplification that the transmitted light field is proportional to the incident field and from the neglect of light field derivative at the interface.