三维各向异性TI介质中的P/S波快速解耦技术及在弹性波逆时偏移中的应用

随着多分量采集技术的发展,弹性波逆时偏移技术在三维各向异性介质复杂地质构造成像中得到了广泛的应用.然而耦合的P波场和S波场,会在传播过程中产生串扰噪声,降低弹性波逆时偏移的成像精度.为了解决这一问题,本研究针对具有倾斜各向异性对称轴的三维横向各向同性(Transverse Isotropy,TI)介质,提出了一种矢量弹性波场快速解耦方法,可以有效提高偏移剖面的成像质量.该方法首先通过坐标转换,将观测系统坐标系的垂直轴旋转到TI介质的对称轴方向,在新坐标系下,根据具有垂直对称轴的三维横向各向同性(Vertical Transverse Isotropy,VTI)介质中的分解算子,推导出三维TI介质解耦算子表达式.接着引入一种在空间域快速计算分解波场的方法,来实现空间域矢量P波场和S波场分离,极大地提高了计算效率.最后,通过点积成像条件,将提出的P/S波分解方法引入到三维TI介质弹性波逆时偏移中,得到高精度的PP和PS成像.与以往的波场分解方法相比,本文方法具有数值稳定和计算效率高的特点.数值算例表明,应用上述三维TI分解算子得到的偏移剖面有效压制了噪声,提高了成像质量.

Women’s empowerment and food consumption:Evidence from female-headed households in Tanzania

Despite the growing recognition of women’s increasing role in the household and corresponding empowerment programs in sub-Saharan Africa,intensive research on the relationship between women’s influence and household food consumption is minimal.Using the most recent(2017-2018)national household survey data from Tanzania,this study examined the influence of women’s empowerment on household food consumption.First,we compared the monthly consumption of eight food categories between female-headed households(FHHs)and male-headed households(MHHs)using both descriptive statistics and the propensity score matching(PSM)method.Furthermore,we adopted the two-stage Linear Expenditure System and Almost Ideal Demand System model(LES-AIDS)to estimate income and price elasticities for the two household types.The results show that FHHs consume bread and cereals,fish,oils and fats,vegetables,and confectionery(sugar,jam,honey,chocolate,etc.)more than MHHs.Moreover,FHHs have a significantly higher income elasticity of demand for all food groups than MHHs.They are also more price elastic than MHHs in meat,fish,oils,fats,sugar,jam,honey,chocolate,etc.

SIMULATION AND STUDY OF THE MODULUS OF ELASTICITY OF NANOCRYSTALLINE MATERIALS

In this paper, a molecular dynamics simulations are provided for atomic structure of nanocrystals(1～3nm)by which t he lattice parameter of X_ray diffraction, cohesive energy and modulus of elas ticity were computed. The results show that the structure of grain and grain bou ndaries in the same in both nanocrystal and coarse grain materials. The decrease of grain size and the increase volume fraction of grain boundaries lead to a se ries of different features, the modulus of elasticity of nanocrystalline materia ls have been found to be much reduced.

Assessment of the Elastic-Wave Well Treatment in Oil-Bearing Clastic and Carbonate Reservoirs

A set of techniques for well treatment aimed to enhance oil recovery are considered in the present study.These are based on the application of elastic waves of various types(dilation-wave,vibro-wave,or other acoustically induced effects).In such a context,a new technique is proposed to predict the effectiveness of the elastic-wave well treatment using the rank distribution according to Zipf’s law.It is revealed that,when the results of elastic wave well treatments are analyzed,groups of wells exploiting various geological deposits can differ in terms of their slope coefficients and free members.As the slope coefficient increases,the average increase in the well oil production rate(after the well treatment)becomes larger.An equation is obtained accordingly for estimating the slope coefficient in the Zipf’s equation from the frequency of the elastic wave.The obtained results demonstrate the applicability of the Zipf’s law in the analysis of the technological efficiency of elastic-wave well treatment methods.

城市轨道交通S型曲线夹直线与长直线地段振动对比试验研究

夹直线在一定程度上可以保证线型连续和行车平顺,但过短的夹直线长度会导致振动叠加,影响旅客舒适度和列车运行的稳定性,增加运营维护成本。本文通过对城市轨道交通S型曲线夹直线与长直线地段振动对比试验研究,从时域、时频域、振动加速度级及地面Z振级等方面对获得的数据进行分析,解释了夹直线段振动叠加的产生机理,验证了车体垂向加速度峰值在缓直点处会出现峰值提前现象,并阐述了夹直线段的振动叠加现象对地面环境振动的影响。结论如下:车辆在通过S型曲线线形变化点时产生的冲击振动会通过钢轨以弹性波的形式向前传播,与车辆在夹直线段行驶产生的振动叠加;由于弹性波的形式传播的速度远高于车辆行驶速度,因此会导致夹直线段出现振动加速度峰值提前现象;S型曲线夹直线的垂向振动加速度的优势频率要高于长直线段,且更向高频分散;S型曲线夹直线产生的振动叠加现象不会对地面环境振动产生影响。

硅和镁质量分数对Al-Si-Mg合金机械力学性能参数影响的计算分析

基于第一性原理,利用CASTEP模块,通过虚拟晶体近似的方法建模,采用GGA+PBE泛函关联函数,研究了Al-Si7-Mg X(X=0.1~0.5)、Al-Si X-Mg0.3(X=6~9)合金中的Si、Mg质量分数的变化对机械力学性能参数的影响。研究表明,随着研究组元质量分数的增加,晶格常数有所减小,最大值与最小值相差约0.265%;弹性常数C_(11)、C_(12)、C_(44)满足结构稳定性条件;体积模量B、剪切模量G、杨氏模量E、泊松比u随着A356中Si、Mg质量分数变化而变化;Mg质量分数的变化对剪切模量G、杨氏模量E的影响比Si质量分数的变化更敏感。计算结果与现有实验数据对比表明,计算结果接近实验数据,第一性原理计算Al-Si-Mg合金材料参数可行。

经会阴三维超声联合SWE评估分娩方式对初产妇肛提肌的影响

目的探讨经会阴三维超声联合剪切波弹性成像技术(shear wave elastography,SWE)在评估分娩方式对初产妇肛提肌形态及功能影响中的应用价值,为早期诊断盆底功能障碍性疾病(pelvic floor disorders,PFD)提供新方法。方法选取在笔者医院分娩后6~8周复查的92例初产妇为产后组,根据其分娩方式分为经阴道分娩组(n=45)和剖宫产组(n=47),选取同期健康未育女性作为对照组(n=43)。采用经会阴三维超声联合SWE进行检查,比较各组间不同状态下肛提肌裂孔形态参数及弹性参数的差异。结果经阴道分娩组在静息、最大缩肛及最大Valsalva状态下肛提肌裂孔面积(levator ani hiatus area,LHA)均大于剖宫产组和对照组,剖宫产组大于对照组,差异有统计学意义(P<0.01);经阴道分娩组在静息及最大缩肛状态下弹性参数较对照组、剖宫产组均减小,差异有统计学意义(P<0.01);而剖宫产组与对照组弹性参数比较,差异无统计学意义(P>0.05);最大Valsalva状态与静息状态下LHA差值(ΔA′)与双侧PR杨氏模量值收缩前后的差值(ΔE)呈显著负相关(左侧r=-0.444,P=0.008;右侧r=-0.488,P=0.002)。结论经会阴三维超声联合SWE可定量评估肛提肌的形态及功能,为PFD的早期诊断提供依据。

Three-Dimensional Static Analysis of Nanoplates and Graphene Sheets by Using Eringen’s Nonlocal Elasticity Theory and the Perturbation Method

A three-dimensional(3D)asymptotic theory is reformulated for the static analysis of simply-supported,isotropic and orthotropic single-layered nanoplates and graphene sheets(GSs),in which Eringen’s nonlocal elasticity theory is used to capture the small length scale effect on the static behaviors of these.The perturbation method is used to expand the 3D nonlocal elasticity problems as a series of two-dimensional(2D)nonlocal plate problems,the governing equations of which for various order problems retain the same differential operators as those of the nonlocal classical plate theory(CST),although with different nonhomogeneous terms.Expanding the primary field variables of each order as the double Fourier series functions in the in-plane directions,we can obtain the Navier solutions of the leading-order problem,and the higher-order modifications can then be determined in a hierarchic and consistent manner.Some benchmark solutions for the static analysis of isotropic and orthotropic nanoplates and GSs subjected to sinusoidally and uniformly distributed loads are given to demonstrate the performance of the 3D nonlocal asymptotic theory.

NEW SECOND ORDER NONCONFORMING TRIANGULAR ELEMENT FOR PLANAR ELASTICITY PROBLEMS

In the use of finite element methods to the planar elasticity problems,one diffculty is to overcome locking when elasticity constant λ→∞.In the case of traction boundary condition,another diffculty is to make the discrete Korn＇s second inequality valid.In this paper,a triangular element is presented.We prove that this element is locking-free,the discrete Korn＇s second inequality holds and the convergence order is two.

From Gradient Elasticity to Gradient Interatomic Potentials: The Case-Study of Gradient London Potential

Motivated by the special theory of gradient elasticity (GradEla), a proposal is advanced for extending it to construct gradient models for interatomic potentials, commonly used in atomistic simulations. Our focus is on London’s quantum mechanical potential which is an analytical expression valid until a certain characteristic distance where “attractive” molecular interactions change character and become “repulsive” and cannot be described by the classical form of London’s potential. It turns out that the suggested internal length gradient (ILG) generalization of London’s potential generates both an “attractive” and a “repulsive” branch, and by adjusting the corresponding gradient parameters, the behavior of the empirical Lennard-Jones potentials is theoretically captured.

A Simple Theory of Asymmetric Linear Elasticity

Rotation is antisymmetric and therefore is not a coherent element of the classical elastic theory, which is characterized by symmetry. A new theory of linear elasticity is developed from the concept of asymmetric strain, which is defined as the transpose of the deformation gradient tensor to involve rotation as well as symmetric strain. The new theory basically differs from the prevailing micropolar theory or couple stress theory in that it maintains the same basis as the classical theory of linear elasticity and does not need extra concepts, such as “microrotation” and “couple stresses”. The constitutive relation of the new theory, the three-parameter Hooke’s law, comes from the theorem about isotropic asymmetric linear elastic materials. Concise differential equations of translational motion are derived consequently giving the same velocity formula for P-wave and a different one for S-wave. Differential equations of rotational motion are derived with the introduction of spin, which has an intrinsic connection with rotation. According to the new theory, S-wave essentially has rotation as large as deviatoric strain and should be referred to as “shear wave” in the context of asymmetric strain. There are nine partial differential equations for the deformation harmony condition in the new theory;these are given with the first spatial differentiations of asymmetric strain. Formulas for rotation energy, in addition to those for (symmetric) strain energy, are derived to form a complete set of formulas for the total mechanical energy.

Chemical bonding and elastic properties of quaternary arsenide oxides YZnAsO and LaZnAsO investigated by first principles

The structural parameters, chemical bonding and elastic properties of the tetragonal phase quaternary arsenide oxides YZnAsO and LaZnAsO were investigated by using density-functional theory （DFT） within generalized gradient approximation （GGA）. The GGA calculated structural parameters are in agreement with the experimental results. Population analysis suggests that the chemical bonding in YZnAsO and LaZnAsO can be classified as a mixture of ionic and covalent characteristic. Single-crystal elastic constants were calculated and the polycrystalline elastic modules were estimated according to Voigt, Reuss and Hill＇s approximations （VRH）. The result shows that both YZnAsO and LaZnAsO are relatively soft materials exhibiting ductile behavior. The calculated polycrystalline elastic anisotropy result shows that LaZnAsO is more anisotropy in compressibility and YZnAsO is more anisotropy in shear.

Effects of porosity on seismic velocities, elastic moduli and Poisson's ratios of solid materials and rocks

The generalized mixture rule(GMR) is used to provide a unified framework for describing Young’s(E),shear(G) and bulk(K) moduli, Lame parameter(l), and P- and S-wave velocities(Vpand Vs) as a function of porosity in various isotropic materials such as metals, ceramics and rocks. The characteristic J values of the GMR for E, G, K and l of each material are systematically different and display consistent correlations with the Poisson’s ratio of the nonporous material(v0). For the materials dominated by corner-shaped pores, the fixed point at which the effective Poisson’s ratio(n) remains constant is at v0=0.2, and J(G) > J(E) > J(K) > J(l) and J(G) < J(E) < J(K) < J(l) for materials with v0> 0.2 and v0< 0.2, respectively.J(Vs) > J(Vp) and J(Vs) < J(Vp) for the materials with v0> 0.2 and v0< 0.2, respectively. The effective n increases, decreases and remains unchanged with increasing porosity for the materials with v0< 0.2,v0> 0.2 and v0=0.2, respectively. For natural rocks containing thin-disk-shaped pores parallel to mineral cleavages, grain boundaries and foliation, however, the n fixed point decreases nonlinearly with decreasing pore aspect ratio(a: width/length). With increasing depth or pressure, cracks with smaller a values are progressively closed, making the n fixed point rise and finally reach to the point at v0=0.2.

A methodology to determine the elastic properties of anisotropic rocks from a single uniaxial compression test

This paper introduces a new methodology to measure the elastic constants of transversely isotropic rocks from a single uniaxial compression test.We first give the mathematical proof that a uniaxial compression test provides only four independent strain equations.As a result,the exact determination of all five independent elastic constants from only one test is not possible.An approximate determination of the Young’s moduli and the Poisson’s ratios is however practical and efficient when adding the Saint-Venant relation as the fifth equation.Explicit formulae are then developed to calculate both secant and tangent definitions of the five elastic constants from a minimum of four strain measurements.The results of this new methodology applied on three granitic samples demonstrate a significant stress-induced nonlinear behavior,where the tangent moduli increase by a factor of three to four when the rock is loaded up to 20 MPa.The static elastic constants obtained from the uniaxial compression test are also found to be significantly smaller than the dynamic ones obtained from the ultrasonic measurements.

Dynamics Modeling and Numerical Analysis of Rotor with Elastic Support/Dry Friction Dampers

The elastic support/dry friction damper is a type of damper which is used for active vibration control in a rotor system.To establish the analytical model of this type of damper,a two-dimensional friction model-ball/plate model was proposed.By using this ball/plate model,a dynamics model of rotor with elastic support/dry friction dampers was established and experimentally verified.Moreover,the damping performance of the elastic support/dry friction damper was studied numerically with respect to some variable parameters.The numerical study shows that the damping performance of the elastic support/dry friction damper is closely related to the stiffness distribution of the rotor-support system,the damper location,the pressing force between the moving and stationary disk,the friction coefficient,the tangential contact stiffness of the contact interface,and the stiffness of the stationary disk.In general,the damper should be located on an elastic support which has a large vibration amplitude in order to achieve a better damping performance,and the more vibration energy in this elastic support concentrates,the better performance of the damper will be.The larger the tangential contact stiffness of the contact interface,and the stiffness of the stationary disk are,the better performance of the damper will be.There will be an optimal value of the friction force at which the damper performs best.

Effect of Axial Deformation on Elastic Properties of Irregular Honeycomb Structure

Irregular honeycomb structures occur abundantly in nature and in man-made products,and are an active area of research.In this paper,according to the optimization of regular honeycomb structures,two types of irregular honeycomb structures with both positive and negative Poisson’s ratios are presented.The elastic properties of irregular honeycombs with varying structure angles were investigated through a combination of material mechanics and structural mechanics methods,in which the axial deformation of the rods was considered.The numerical results show that axial deformation has a significant influence on the elastic properties of irregular honeycomb structures.The elastic properties of the structure can be considered by the enclosed area of the unit structure,the shape of the unit structure,and the elastic properties of the original materials.The elastic properties considering the axial deformation of rods studied in this study can provide a reference for other scholars.

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied. First, the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier expanding with respect to azimuth, and the state equation is established by Hankel integral transform method, furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation. Secondly, by the transfer matrixes, the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of ground as well as the contact conditions. Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundaryvalue problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily. At the end of this paper, a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.

Neumann＇s method for boundary problems of thin elastic shells

The possibility of using Neumann＇s method to solve the boundary problems for thin elastic shells is studied. The variational statement of the static problems for the shells allows for a problem examination within the distribution space. The convergence of Neumann＇s method is proven for the shells with holes when the boundary of the domain is not completely fixed. The numerical implementation of Neumann＇s method normally requires significant time before any reliable results can be achieved. This paper suggests a way to improve the convergence of the process, and allows for parallel computing and evaluation during the calculations.

Topology and Shape Optimization of 2-D and 3-D Micro-ArchitecturedThermoelastic Metamaterials Using a Parametric Level Set Method

2-D and 3-D micro-architectured multiphase thermoelastic metamaterials are designed and analyzed using a parametric level set method for topology optimization and the finite element method.An asymptotic homogenization approach is employed to obtain the effective thermoelastic properties of the multiphase metamaterials.Theε-constraint multi-objective optimization method is adopted in the formulation.The coefficient of thermal expansion(CTE)and Poisson’s ratio(PR)are chosen as two objective functions,with the CTE optimized and the PR treated as a constraint.The optimization problems are solved by using the method of moving asymptotes.Effective isotropic and anisotropic CTEs and stiffness constants are obtained for the topologically optimized metamaterials with prescribed values of PR under the constraints of specified effective bulk modulus,volume fractions and material symmetry.Two solid materials along with one additional void phase are involved in each of the 2-D and 3-D optimal design examples.The numerical results reveal that the newly proposed approach can integrate shape and topology optimizations and lead to optimal microstructures with distinct topological boundaries.The current method can topologically optimize metamaterials with a positive,negative or zero CTE and a positive,negative or zero Poisson’s ratio.

Forward modeling of ocean-bottom cable data and wave-mode separation in fluid–solid elastic media with irregular seabed

In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.