耦合地震动下地铁车站的动力响应研究

耦合地震动对地铁车站结构的响应较为复杂,本文运用有限单元法进行动力时程分析,研究地铁车站结构承受水平、竖向和双向耦合地震动作用下的动力响应规律.结果表明:(1)耦合地震动对结构相对水平位移影响较小,但对结构的相对竖向位移影响较大.(2)耦合地震动会减小中柱主应力值,但并未较大程度上改变中柱主应力值大小分布情况.(3)车...     

Machine learning for pore-water pressure time-series prediction: Application of recurrent neural networks

Knowledge of pore-water pressure(PWP)variation is fundamental for slope stability.A precise prediction of PWP is difficult due to complex physical mechanisms and in situ natural variability.To explore the applicability and advantages of recurrent neural networks(RNNs)on PWP prediction,three variants of RNNs,i.e.,standard RNN,long short-term memory(LSTM)and gated recurrent unit(GRU)are adopted and compared with a traditional static artificial neural network(ANN),i.e.,multi-layer perceptron(MLP).Measurements of rainfall and PWP of representative piezometers from a fully instrumented natural slope in Hong Kong are used to establish the prediction models.The coefficient of determination(R^2)and root mean square error(RMSE)are used for model evaluations.The influence of input time series length on the model performance is investigated.The results reveal that MLP can provide acceptable performance but is not robust.The uncertainty bounds of RMSE of the MLP model range from 0.24 kPa to 1.12 k Pa for the selected two piezometers.The standard RNN can perform better but the robustness is slightly affected when there are significant time lags between PWP changes and rainfall.The GRU and LSTM models can provide more precise and robust predictions than the standard RNN.The effects of the hidden layer structure and the dropout technique are investigated.The single-layer GRU is accurate enough for PWP prediction,whereas a double-layer GRU brings extra time cost with little accuracy improvement.The dropout technique is essential to overfitting prevention and improvement of accuracy.     

In situ stress field of eastern Australia

New in situ data based on hydraulic fracturing and overcoring have been compiled for eastern Australia, increasing from 23 to 110 the number of in situ stress analyses available for the area between and including the Bowen and Sydney Basins. The Bowen Basin displays a consistent north‐northeast maximum horizontal stress (σ_H) orientation over some 500 km. Stress orientations in the Sydney Basin are more variable than in the Bowen Basin, with areas of the Sydney Basin exhibiting north‐northeast, northeast, east‐west and bimodal σ_H orientations. Most new data indicate that the overburden stress (σ_V) is the minimum principal stress in both the Bowen and Sydney Basins. The Sydney Basin is relatively seismically active, whereas the Bowen Basin is relatively aseismic. Despite the fact that in situ stress measurements sample the stress field at shallower depth than the seismogenic zone, there is a correlation between the stress measurements and seismicity in the two areas. Mohr‐Coulomb analysis of the propensity for failure in the Sydney Basin suggests 41% of the new in situ stress data are indicative of failure, as opposed to 13% in the Bowen Basin. The multiple pre‐existing structural grains in the Sydney Basin further emphasise the difference between propensity for failure in the two areas. Previous modelling of intraplate stresses due to plate boundary forces has been less successful at predicting stress orientations in eastern than in western and central Australia. Nonetheless, stress orientation in the Bowen Basin is consistent with that predicted by modelling of stresses due to plate boundary forces. Variable stress orientations in the Sydney Basin suggest that more local sources of stress, such as those associated with the continental margin and with local structure, significantly influence stress orientation. The effect of local sources of stress may be relatively pronounced because stresses due to plate boundary forces result in low horizontal stress anisotropy in the Sydney Basin.     

A 3D discrete FEM iterative algorithm for solving the water pipe cooling problems of massive concrete structures

Water pipe cooling has been widely used for the temperature control and crack prevention of massive concrete structures such as high dams. Because both under‐cooling and over‐cooling may reduce the efficiency of crack prevention, or even lead to great harm to structures, we need an accurate and robust numerical tool for the prediction of cooling effect. Here, a 3D discrete FEM Iterative Algorithm is introduced, which can simulate the concrete temperature gradient near the pipes, as well as the water temperature rising along the pipes. On the basis of the heat balance between water and concrete, the whole temperature field of the problem can be computed exactly within a few iteration steps. Providing the pipe meshing tool for building the FE model, this algorithm can take account of the water pipe distribution, the variation of water flow, water temperature, and other factors, while the traditional equivalent algorithm based on semi‐theoretical solutions can only solve problems with constant water flow and water temperature. The validation and convergence are proved by comparing the simulated results and analytical solutions of two standard second‐stage cooling problems. Then, a practical concrete block with different cooling schemes is analyzed and the influences of cooling factors are investigated. In the end, detailed guidance for pipe system optimization is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.     

Stress‐dependent hardening and failure surfaces of dry sand

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.     

Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.     

A multiscale finite element formulation for axisymmetric elastoplasticity with volumetric locking

Similar to plane strain, axisymmetric stress problem is also highly kinematics constrained. Standard displacement‐based finite element exhibits volumetric locking issue in simulating nearly/fully incompressible material or isochoric plasticity under axisymmetric loading conditions, which severely underestimates the deformation and overestimates the bearing capacity for structural/geotechnical engineering problems. The aim of this paper is to apply variational multiscale method to produce a stabilized mixed displacement–pressure formulation, which can effectively alleviate the volumetric locking issue for axisymmetric stress problem. Both nearly incompressible elasticity and isochoric J₂ elastoplasticity are investigated. First‐order 3‐node triangular and 4‐node quadrilateral elements are tested for locking issues. Several representative simulations are provided to demonstrate the performance of the linear elements, which include the convergence study and comparison with closed‐form solutions. A comparative study with pressure Laplacian stabilized formulation is also presented. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of stress mapping in cross‐section to understanding ore geometry,predicting ore zones and development of drilling strategies

Stress mapping is a numerical modelling technique used to determine the distribution and relative magnitude of stress during deformation in a mineralised terrane. It is based on the general principle that fluid flow in the Earth's crust is primarily related to pressure gradients. It is best applied to epigenetic hydrothermal mineral deposits, where fluid flow and fluid flux are enhanced in dilational sections of structures and in sites of enhanced rock permeability due to high fracture density. These are defined by sites of low minimum principal stress (σ₃). Most stress mapping is carried out in two dimensions in plan view using geological maps. This is suitable for terranes with steeply dipping lithostratigraphy and structures in which the distribution of mineral deposits is largely controlled by fault structures portrayed on the maps. However, for terranes with gently dipping sequences and structures, and for situations where deposits are sited in and near the hinges of complex fold structures, stress mapping in cross‐section is preferable. The effectiveness of stress mapping is maximised if mineralisation was late in the evolutionary history of the host terrane, and hence the structural geometry of the terrane and contained deposits were essentially that expressed today. The orientation of syn‐mineralisation far‐field stresses must also be inferred. Two examples of orogenic gold deposits, which meet the above criteria, are used to illustrate the potential of stress mapping in cross‐section. Sunrise Dam, located in the Archaean Yilgarn Craton, is a lode‐gold deposit sited in a thrust‐fold belt. Stress mapping illustrates the heterogeneity of stress distribution in the complex structural geometry of the deposit, and predicts the preferential siting of ore zones around the intersections of more steeply dipping, linking thrusts and banded iron‐formation units, and below the controlling more gently dipping basal thrust, the Sunrise Shear. The Howley Anticline in the Pine Creek block hosts several Palaeoproterozoic gold deposits, sited in complex anticlinal structures in greywacke sequences. Stress mapping indicates that gold ores should develop in the hinge zones of symmetrical anticlines, in the hinge zones and more steeply dipping to overturned limbs of asymmetric anticlines, and in and around thrusts in both anticlines and parasitic synclines. The strong correlation between the predictions of the stress mapping, based on the distribution of low σ₃, and the location of gold ores emphasises the potential of stress mapping in cross‐section, not only as an exploration tool for the discovery of additional resources or deposits, but also as a test of geological models. Knowledge of the potential siting of gold ores and their probable orientations also provides a guide to drilling strategies in both mine‐ and regional‐scale exploration.     

Long waves due to a symmetrical wind stress on the surface of a rotating ocean

Abstract

The problem of unsteady long waves generated by any horizontal and symmetrically distributed, time-periodic surface wind on a rotating ocean is analysed for large times and distances. Uniform asymptotic estimates of the surface displacement in the unsteady state are obtained. The steady-state wave and velocity fields at any distance are also determined. Some characteristics of the unsteady and steady motions are described. Also noted are the features that distinguish the motion from its one-dimensional analogue for which a non-uniform analysis in the unsteady state along with a large-distance form of the surface elevation are already known.