Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

A machine learning approach for predicting computational intensity and domain decomposition in parallel geoprocessing

ABSTRACT

High performance computing is required for fast geoprocessing of geospatial big data. Using spatial domains to represent computational intensity (CIT) and domain decomposition for parallelism are prominent strategies when designing parallel geoprocessing applications. Traditional domain decomposition is limited in evaluating the computational intensity, which often results in load imbalance and poor parallel performance. From the data science perspective, machine learning from Artificial Intelligence (AI) shows promise for better CIT evaluation. This paper proposes a machine learning approach for predicting computational intensity, followed by an optimized domain decomposition, which divides the spatial domain into balanced subdivisions based on the predicted CIT to achieve better parallel performance. The approach provides a reference framework on how various machine learning methods including feature selection and model training can be used in predicting computational intensity and optimizing parallel geoprocessing against different cases. Some comparative experiments between the approach and traditional methods were performed using the two cases, DEM generation from point clouds and spatial intersection on vector data. The results not only demonstrate the advantage of the approach, but also provide hints on how traditional GIS computation can be improved by the AI machine learning.     

Fast genetic mapping in barley: case studies of cuticle mutants using RNA-sequencing

Barley(Hordeum vulgare L.) is one of the earliest domesticated crop species and ranked as the fourth largest cereal production worldwide. Forward genetic studies in barley have greatly advanced plant genetics during the last century; however, most genes are identified by the conventional mapping method. Array genotyping and exome-capture sequencing have also been successfully used to target the causal mutation in barley populations, but these techniques are not widely adopted because of associated costs and partly due to the huge genome size of barley. This review summarizes three mapping cases of barley cuticle mutants in our laboratory with the help of RNA-sequencing. The causal mutations have been successfully identified for two of them and the target genes are located in the pericentromeric regions. Detailed information on the mapping-by-sequencing, mapping-and-sequencing, and RNA-sequencing assisted linkage mapping are presented and some limitations and challenges on the mapping assisted by RNA sequencing are also discussed. The alternative and elegant methods presented in this review may greatly accelerate forward genetics of barley mapping, especially for laboratories without large funding.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

THE CHARACTERISTIC OF MARINE ENVIRONMENT IN LINGDINGYANG ESTUARY

1 HYDROLOGIC FEATURES Lingdingyang Estuary, located at the middle south of Guangdong Province, is a bell-shaped estuary with a north-south direction. Its area is about 2100km2. The north of Qi′ao Island and Inner-Lingding Island, and the south of Humen are grouped as Neilingdingyang Estuary, having an area of 1041km2. Affected by topography, runoff and tide, its dynamic condition is very complicated. Different water areas have different hydrologic features. The topography under …     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

基于统一强度理论的地基临界荷载公式

基于统一强度理论,按照临界荷载公式的推导方法,得到地基临界荷载的统一解公式。Mohr-Coulomb强度理论计算的结果为其特例。把计算结果与《建筑地基基础设计规范》中的结果进行了对比分析,认为当不同程度地考虑中间主应力σ2的影响时,可以有效地发挥地基的强度潜能。并通过实例计算了统一强度理论参数b取不同值时地基承载力特征值,并进行了对比分析。     

The Expansion of the Distance Between Two Arbitrary Gravitating Points as Applied to Specific Models of Celestial Bodies

For the complete system of biharmonic functions a suitable representation in spheroidal coordinates is found. It is used for expanding the distance between two arbitrary gravitating points and its inverse quantity in appropriate series. Such expansions can be of interest and useful in a number of problems of celestial mechanics and stellar dynamics.     

基于负反馈的网络负载调度算法

提出一种基于负反馈权值的动态网络负载调度算法。算法主要用在多出口路由器上，其执行效率高，开销小。基于算法的多出口路由器不仅能够很好的保持出口间的负载均衡。还具有良好的出口容错性能。