Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane. 相似文献
In this study, a directional interpolation infinite element suited to a saturated porous medium is presented to account for dynamic problems with semi-infi 相似文献
This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features. 相似文献
We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AICC] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component. 相似文献
The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m3. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage. 相似文献
Southern right whales—Eubalaena australis (Desmoulins, 1822)—migrate seasonally from high‐latitude feeding grounds to coastal breeding and calving grounds at lower latitudes such as the southern coast of Brazil. Understanding how these whales are distributed along the coast is important for monitoring their postwhaling recovery and defining management strategies. In this study, we applied Kernel density estimators to aerial survey data to determine main occurrence and concentration areas of right whales in southern Brazil and investigate inter‐ and intra‐annual distribution patterns between 2003 and 2012. Our results show considerable variation in area usage within and among years, and changes in the general distribution pattern of right whales in the last years of the study. Intra‐annually, higher concentration area tended to expand from July to September and decrease in November. Some areas stood out as high‐density areas for right whales: Ribanceira/Ibiraquera, Itapirubá Sul/Sol, and from Arroio to Gaivota. Some evidences also suggest preferential areas for mother–calf pairs. The higher concentration area of right whales in southern Brazil was estimated at 52,541 km2 and the occurrence area was 682.69 km2, which is the whole study area. As right whale distribution in the region is likely expanding due to this population's current recovery, our study provides essential information for management plan of the Right Whale Environmental Protection Area. 相似文献
A series of fully three-dimensional(3 D) numerical simulations of flow past a free-to-oscillate curved flexible riser in shear flow were conducted at Reynolds number of 185–1015. The numerical results obtained by the two-way fluid–structure interaction(FSI) simulations are in good agreement with the experimental results reported in the earlier study. It is further found that the frequency transition is out of phase not only in the inline(IL) and crossflow(CF) directions but also along the span direction. The mode competition leads to the non-zero nodes of the rootmean-square(RMS) amplitude and the relatively chaotic trajectories. The fluid–structure interaction is to some extent reflected by the transverse velocity of the ambient fluid, which reaches the maximum value when the riser reaches the equilibrium position. Moreover, the local maximum transverse velocities occur at the peak CF amplitudes, and the values are relatively large when the vibration is in the resonance regions. The 3 D vortex columns are shed nearly parallel to the axis of the curved flexible riser. As the local Reynolds number increases from 0 at the bottom of the riser to the maximum value at the top, the wake undergoes a transition from a two-dimensional structure to a 3 D one. More irregular small-scale vortices appeared at the wake region of the riser, undergoing large amplitude responses. 相似文献
The transfer and evolution of stress among rock blocks directly change the void ratios of crushed rock masses and affect the flow of methane in coal mine gobs. In this study, a Lagrange framework and a discrete element method, along with the soft-sphere model and EDEM numerical software, were used. The compaction processes of rock blocks with diameters of 0.6, 0.8, and 1.0 m were simulated with the degrees of compression set at 0%, 5%, 10%, 15%, 20%, and 25%. This study examines the influence of stress on void ratios of compacted crushed rock masses in coal mine gobs. The results showed that stress was mainly transmitted downward through strong force chains. As the degree of compression increased, the strong force chains extended downward, which resulted in the stress at the upper rock mass to become significantly higher than that at the lower rock mass. It was determined that under different degrees of compression, the rock mass of coal mine gobs could be divided, from the bottom to the top, into a lower insufficient compression zone (ICZ) and an upper sufficient compression zone (SCZ). From bottom to top, the void ratios in the ICZ sharply decreased and those in the SCZ slowly decreased. Void ratios in the ICZ were 1.2–1.7 times higher than those in the SCZ.