An approach to the modeling of viscoelastic damage. Application to the long‐term creep of gypsum rock materials     

B. Nedjar  R. Le Roy 《国际地质力学数值与分析法杂志》2013,37(9):1066-1078

A three‐dimensional phenomenological model is developed to describe the long‐term creep of gypsum rock materials. The approach is based on the framework of continuum damage mechanics where coupling with viscoelasticity is adopted. Specifically, a local damage model based on the concept of yield surface is proposed and deeply investigated. Among the many possibilities, we choose in this work its coupling with a generalized Kelvin–Voigt rheological model to formulate the whole behavior. Long‐term as well as short‐term relaxation processes can be integrated in the model by means of as many as necessary viscoelastic processes. The numerical discretization is described for an easy integration within a finite element procedure. Finally, a set of numerical simulations is given to show the possibilities of the presented model. It shows good agreement with some experimental results found in the literature. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Water waves generated by instantaneous disturbances at the bed of a sloping beach     

Rumpa Chakraborty 《地球物理与天体物理流体动力学》2013,107(5):481-496

The two-dimensional problem of the generation of water waves due to instantaneous disturbances prescribed at the bed of a beach sloping at an arbitrary angle is studied here. It is formulated in terms of an initial-boundary-value problem for the velocity potential describing the motion in the fluid region assuming the linear theory. Using the Laplace transform in time and the Mellin transform in distance, the problem is reduced to solving a difference equation whose method of solution is of considerable importance in the literature. The form of the free surface is obtained in terms of a multiple infinite integral that is evaluated by the method of steepest-descent. For some prescribed forms of the disturbance at the bed of the beach, the free surface is depicted in a number of figures for different beach angles. It is observed that as the beach angle decreases, the maximum wave height increases, which is plausible.  相似文献   

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

Songhao Shang 《水文研究》2012,26(22):3338-3343

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.  相似文献   

Particle manifold method (PMM): A new continuum‐discontinuum numerical model for geomechanics     

L. Sun  G.F. Zhao  J. Zhao 《国际地质力学数值与分析法杂志》2013,37(12):1711-1736

Particle manifold method (PMM) is a new extension of the numerical manifold method (NMM). PMM uses a mathematical cover system to describe the motion and deformation of a particle‐based physical domain. By introducing the concept of particle into NMM, PMM takes the advantages of easy topological and contact operations with particles. In this article, the methodology, formulations and implementation of the method are presented, together with modelling examples for validation. It is found that good solutions for both continuous and discontinuous problems are obtained by the new developed PMM. Due to the underlying coupled continuum‐discontinuum property of PMM, it has great potential for modelling of geomechanical problems. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Spatial and temporal distributions of nitrogen, phosphorus and heavy metals in the intertidal sediment of the Chang jiang River Estuary in China   总被引：2，自引：0，他引：2  

QUAN Weimin  SHI Liyan  HAN Jindi  PING Xianyin  SHEN Anglv  CHEN Yaqu 《海洋学报(英文版)》2010,29(1):108-115

Spatial and temporal distributions of nutrients (N and P) and heavy metals (Cu, Zn, Pb and Cd) in the intertidal sediments at Dongtan wetland of the Changjiang River (Yangtze River) Estuary were examined by laying one transect running seaward to landward during January, April, June, August and October of 2005. Total nitrogen (TN) and heavy metal concentrations in the intertidal sediment showed an increasing trend from open bare °at to the high marsh, with the maximum concentrations in the high marsh and the minimum values in open bare °at. In contrast, sediment total phosphorus (TP) concentrations changed little among sampling sites. One-way ANOVA found that there were signiˉcantly temporal variations of nutrients and heavy metals concentrations in the intertidal sediment (P <0.05). The concentrations of heavy metals generally showed an increasing trend from January to October. As compared with other large estuaries, heavy metal contamination in the intertidal sediment of the Changjiang River Estuary was relatively low.  相似文献   

The form of the asymptotic depth-limited wind-wave spectrum: Part III — Directional spreading     

I.R. Young 《Coastal Engineering》2010

The directional spreading of both the wavenumber and frequency spectra of finite-depth wind generated waves at the asymptotic depth limit are examined. The analysis uses the Wavelet Directional Method, removing the need to assume a form for the dispersion relationship. The paper shows that both the wavenumber and frequency forms are narrowest at the spectral peak and broaden at wavenumbers (frequencies) both above and below the peak. The directional spreading of the wavenumber spectrum is bi-modal above the spectral peak. In contrast, the frequency spectrum is uni-modal. This difference is shown to be the result of energy in the wind direction being displaced from the linear dispersion shell. A full parametric relationship for the directional spreading of the wavenumber spectrum is developed. The analysis clearly shows that typical dispersion relationships are questionable at high frequencies and that such effects can be significant. This result supports greater attention being focussed on the routine recording of wavenumber spectra, rather than frequency spectra.  相似文献   

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir     

Bin Xu  Ron C.K. Wong 《国际地质力学数值与分析法杂志》2013,37(18):3131-3149

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Analytical and computational procedure for solving the consolidation problem of layered soils     

S. Sadiku 《国际地质力学数值与分析法杂志》2013,37(16):2789-2800

One‐dimensional consolidation analysis of layered soils conventionally entails solving a system of differential equations subject to the flow conditions at the bounding upper and lower surfaces, as well as the continuity conditions at the interface of every pair of contiguous layers. Formidable computational efforts are required to solve the ensuing transcendental equations expressing the matching conditions at the interfaces, using this method. In this paper, the jump discontinuities in the flow parameters upon crossing from one layer to the other have been systematically built into a single partial differential equation governing the space–time variation of the excess pore pressure in the entire composite medium, by the use of the Heaviside distribution. Despite the presence of the discontinuities in the coefficients of the differential equation, a closed‐form solution in the sense of an infinite generalized Fourier series is obtained, in addition to which is the development of a Green's function for the differential problem. The eigenfunctions of the composite medium are the coordinate functions of the series, obtained computationally through the application of the extended equations of Galerkin. The analysis has been illustrated by solving the consolidation problem of a four‐layer composite, and the results obtained agree very well with the results obtained by previous researchers. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

胶州湾水体镉的迁移过程和环境本底值结构     

杨东方  陈豫  王虹  杨春茹  郭军辉 《海岸工程》2010,29(4):73-82

根据1979年的胶州湾水域调查资料,分析重金属Cd在胶州湾水域的水平分布、垂直分布和季节变化。研究结果表明:在整个胶州湾水域,水质没有受到任何Cd的污染。在胶州湾和湾口水域只有两个Cd的来源。一个是来源于地壳岩石风化,通过陆地径流把重金属Cd注入胶州湾水域;另一个是来源于海底,经过海洋水流的作用把重金属Cd注入湾口水域。作者提出了重金属在水域的环境本底值结构:基础本底值、陆地径流的输入量和海洋水流的输入量,并且应用于胶州湾水域。作者认为在来源的迁移过程中,有陆地来源迁移和海底来源迁移。根据杨东方的水域迁移过程研究,陆地径流决定了Cd的表、底层质量浓度变化。  相似文献   

Curve‐Number/Green–Ampt mixed procedure for streamflow predictions in ungauged basins: Parameter sensitivity analysis     

S. Grimaldi  A. Petroselli  N. Romano 《水文研究》2013,27(8):1265-1275

A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献