基于有限体积直接平均方法(Finite-volume direct averaging micromechanics,FVDAM),建立了一种从复合材料细观到宏观的统一本构模型。根据均匀化方法和连续介质力学构建复合材料的宏细观相关矩阵,通过该矩阵将细观组分材料的损伤性能...基于有限体积直接平均方法(Finite-volume direct averaging micromechanics,FVDAM),建立了一种从复合材料细观到宏观的统一本构模型。根据均匀化方法和连续介质力学构建复合材料的宏细观相关矩阵,通过该矩阵将细观组分材料的损伤性能传递到宏观复合材料中,计算了湿热环境中复合材料的细观应力场。结果表明:FVDAM采用子胞边界平均位移作为未知量,使本构模型中的未知量总数大为减小,相对提高了模型的效率,但这些方程都是建立在平均意义上的,因此预测的应力场存在一定的不连续性;湿热环境下,前期的吸湿有缓解热残余应力的作用,随着时间的增加,吸湿的影响逐渐超过热残余应力的影响。展开更多
The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom...The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.展开更多
It is well known that the macroscopic Maxwell’s equations can be obtained from the corresponding microscopic or atomic equations by a proper averaging process. The purpose of this paper is to present the macroscopic ...It is well known that the macroscopic Maxwell’s equations can be obtained from the corresponding microscopic or atomic equations by a proper averaging process. The purpose of this paper is to present the macroscopic Maxwell’s equations which are valid in all regions of space, including an interface between two different media; and the boundary conditions can naturally emerge from the macroscopic equations as an effect of average of the microscopic Maxwell’s equations. In addition, the application of the unit step functions and the Dirac delta function to our discussion not only permits great mathematical simplicity but also gives rise to convenient physical concepts for the description and representation of the actual fields in the vicinity of the interface.展开更多
We derive exact near-wall and centerline constraints and apply them to improve a recently proposed LPR model for finite Reynolds number(Re) turbulent channel flows.The analysis defines two constants which are invarian...We derive exact near-wall and centerline constraints and apply them to improve a recently proposed LPR model for finite Reynolds number(Re) turbulent channel flows.The analysis defines two constants which are invariant with Re and suggests two more layers for incorporating boundary effects in the prediction of the mean velocity profile in the turbulent channel.These results provide corrections for the LPR mixing length model and incorrect predictions near the wall and the centerline.Moreover,we show that the analysis,together with a set of well-defined sensitive indicators,is useful for assessment of numerical simulation data.展开更多
A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain ...A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain smoothing operations. The strain smoothing is defined through strain averaging in the element sub-domains which are subsequently used for numerical integration of the Galerkin weak form. This formulation satisfies the orthogonality condition of the assumed strain methods. Meanwhile the present formulation totally avoids the gradient computation of the rational NURBS basis functions in the formulation of stiffness matrix. A transformation method is employed to accurately enforce the displacement boundary conditions. Numerical results demonstrate that the present formation gives very satisfactory solution accuracy simultaneously with improved computational efficiency.展开更多
文摘基于有限体积直接平均方法(Finite-volume direct averaging micromechanics,FVDAM),建立了一种从复合材料细观到宏观的统一本构模型。根据均匀化方法和连续介质力学构建复合材料的宏细观相关矩阵,通过该矩阵将细观组分材料的损伤性能传递到宏观复合材料中,计算了湿热环境中复合材料的细观应力场。结果表明:FVDAM采用子胞边界平均位移作为未知量,使本构模型中的未知量总数大为减小,相对提高了模型的效率,但这些方程都是建立在平均意义上的,因此预测的应力场存在一定的不连续性;湿热环境下,前期的吸湿有缓解热残余应力的作用,随着时间的增加,吸湿的影响逐渐超过热残余应力的影响。
基金supported by the High-Tech Research and Development Program of China(Grant No.2008AA093001)China Petroleum & Chemical Corporation(Grant No.YPH08006)
文摘The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.
文摘It is well known that the macroscopic Maxwell’s equations can be obtained from the corresponding microscopic or atomic equations by a proper averaging process. The purpose of this paper is to present the macroscopic Maxwell’s equations which are valid in all regions of space, including an interface between two different media; and the boundary conditions can naturally emerge from the macroscopic equations as an effect of average of the microscopic Maxwell’s equations. In addition, the application of the unit step functions and the Dirac delta function to our discussion not only permits great mathematical simplicity but also gives rise to convenient physical concepts for the description and representation of the actual fields in the vicinity of the interface.
基金supported by the National Natural Science Foundation of China (Grant Nos. 90716008 and 10921202)the National Basic Research Program of China (Grant No. 2009CB724100)
文摘We derive exact near-wall and centerline constraints and apply them to improve a recently proposed LPR model for finite Reynolds number(Re) turbulent channel flows.The analysis defines two constants which are invariant with Re and suggests two more layers for incorporating boundary effects in the prediction of the mean velocity profile in the turbulent channel.These results provide corrections for the LPR mixing length model and incorrect predictions near the wall and the centerline.Moreover,we show that the analysis,together with a set of well-defined sensitive indicators,is useful for assessment of numerical simulation data.
基金supported by the National Natural Science Foundation of China (Grant No. 10972188)the Program for New Century Excellent Talents in University from China Ministry of Education (Grant No. NCET-09-0678)the Fundamental Research Funds for the Central Universities of China (Grant No. 2010121073)
文摘A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain smoothing operations. The strain smoothing is defined through strain averaging in the element sub-domains which are subsequently used for numerical integration of the Galerkin weak form. This formulation satisfies the orthogonality condition of the assumed strain methods. Meanwhile the present formulation totally avoids the gradient computation of the rational NURBS basis functions in the formulation of stiffness matrix. A transformation method is employed to accurately enforce the displacement boundary conditions. Numerical results demonstrate that the present formation gives very satisfactory solution accuracy simultaneously with improved computational efficiency.
