Unified two-phase nonlocal formulation for vibration of functionally graded beams resting on nonlocal viscoelastic Winkler-Pasternak foundation

A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation,which were not considered in most literature on this subject,and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven(ε-D)and stress-driven(σ-D)two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered,which can address both the stiffness softening and toughing effects due to scale reduction.The generalized differential quadrature method(GDQM)is used to solve the complex eigenvalue problem.After verifying the solution procedure,a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained.Subsequently,the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.

Simulation of gas–liquid two-phase flow in a flow-focusing microchannel with the lattice Boltzmann method

A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.

Free Vibration Analysis of Rectangular Plate with Cutouts under Elastic Boundary Conditions in Independent Coordinate Coupling Method

Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.

Research on propagation properties of elastic waves in two-phase anisotropic media

With the development of seismic engineering and seismic exploration of energy, the underground media that westudy are more and more complicated. Conventional anisotropy theory or two-phase isotropy theory is difficult todescribe anisotropic media containing fluid, such as fractures containing gas, shales containing water Based onBlot theory about two-phase anisotropy, with the use of elastic plane wave equations, we get Christoffel equations.We calculate and analyze the effects of frequency on phase velocity, attenuation, amplitude ratio and polarizationdirection of elastic waves of two-phase, transversely isotropic media. Results show that frequency affects slow Pwave the greatest among the four kinds of waves, i.e., fast P wave, slow P wave, fast S wave and slow S wave.Fluid phase amplitude to solid phase amplitude ratio of fast P wave, fast S wave and slow S wave approaches unitfor large dissipation coefficients. Polarization analysis shows that polarization direction of fluid phase displacement is different from, not parallel to or reverse to, that of solid phase displacement in two-phase anisotropic media.

Computational Complex Method for Finite Elastic Plane Problem with Multiple Elliptical Holes

把多复变量应力函数和边界配点最小二乘法相结合，计算含有任意多个椭圆孔有限大小弹性板的应力场．每个孔的大小、位置为任意．板的边界由下列四种边界条件以及它们的组合所组成：１）已知Ｐｘ和Ｐｙ的应力边界，２）已知Ｖｘ和Ｖｙ的位移边界，３）已知法向位移ＶＮ和切向力ＰＴ的混合边界，４）已知法向力ＰＮ和切向位移ＶＴ的混合边界．编制了相应的ＦＯＲＴＲＡＮ程序，进行考题和算例分析，结果表明本方法行之有效．

Study of elastic wave propagation in two phase anisotropic media by numerical modeling of pseudospectral method

When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.

BOUNDARY ELEMENT METHOD FOR MOVING AND ROLLING CONTACT OF 2D ELASTIC BODIES WITH DEFECTS

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented.Both thedisplacement and the traction boundary conditions are satisfied on the contactingregion in the sense of discretization.An algorithm to deal with the moving of thecontact boundary on a larger possible contact region is presented.The algorithm isgeneralized to rolling contact problem as well.Some numerical examples of movingand rolling contact of 2D elastic bodies with or without friction,including the bodieswith a hole-type defect,are given to show the effectiveness and the accuracy of thepresented schemes.

Application of Two-dimensional Viscous CE/SE Method in Calculation of Two-phase Detonation

The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.

Buckling analysis of functionally graded plates partially resting on elastic foundation using the differential quadrature element method

We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.

Reconstruction of electrical capacitance tomography images based on fast linearized alternating direction method of multipliers for two-phase flow system

Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed and nonlinear inverse problem of ECT image reconstruction,a new ECT image reconstruction method based on fast linearized alternating direction method of multipliers(FLADMM)is proposed in this paper.On the basis of theoretical analysis of compressed sensing(CS),the data acquisition of ECT is regarded as a linear measurement process of permittivity distribution signal of pipe section.A new measurement matrix is designed and L1 regularization method is used to convert ECT inverse problem to a convex relaxation problem which contains prior knowledge.A new fast alternating direction method of multipliers which contained linearized idea is employed to minimize the objective function.Simulation data and experimental results indicate that compared with other methods,the quality and speed of reconstructed images are markedly improved.Also,the dynamic experimental results indicate that the proposed algorithm can ful fill the real-time requirement of ECT systems in the application.

STOCHASTIC BOUNDARY ELEMENT METHOD IN ELASTICITY

The stochastic boundary element method is developed to analyzeelasticity problems with random material and/or geometrical parameters and ran-domly perturbed boundaries. Based on the first-order Taylor series expansion, theboundary integration equations concerning the mean and deviation of the displace-ments are derived, respectively. It is found that the randomness of material param-eters is equivalent to a random body force, so the mean and covariance matricesof unknown boundary displacements and tractions can be obtained. Furthermore,the mean and covariance of displacements and stresses at inner points can also beobtained. Numerical examples show that the proposed stochastic boundary elementmethod gives satisfactory solutions, as compared with those obtained by theoreticalanalysis or other numerical methods.

ANALYSIS OF BENDING, VIBRATION AND STABILITY FOR THIN PLATE ON ELASTIC FOUNDATION BY THE MULTIVARIABLE SPLINE ELEMENT METHOD

In this paper, the bicubic splines in product form are used to construct the multi-field functions for bending moments, twisting moment and transverse displacement ofthe plate on elastic foundation. The multivariable spline element equations are derived,based on the mixed variational principle. The analysis and calculations of bending,vibration and stability of the plates on elastic foundation are presented in the paper.Because the field functions of plate on elastic foundation are assumed independently,the precision of the field variables of bending moments and displacement is high.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

NEW STUDYING OF LATTICE BOLTZMANN METHOD FOR TWO-PHASE DRIVEN IN POROUS MEDIA

By using the interaction of particles,such as the physical principle of the same attract each other and the different repulse each other,a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore,the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.

Eigenfunction expansion method and its application to two-dimensional elasticity problems based on stress formulation

This paper proposes an eigenfunction expansion method to solve twodimensional (2D) elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.

Determination of dynamic capillary effect on two-phase flow in porous media: A perspective from various methods

The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.

Transverse Vibration Analysis of Single-Walled Carbon Nanotubes Embedded in an Elastic Medium Using Bernoulli-Fourier Method

Based on the Timoshenko beam theory and Bernoulli-Fourier method, a single-elastic beam model is developed for transverse vibrations of single-walled carbon nanotubes under additional axial load, which includes the effects of the elastic medium around them. Explicit expressions are derived for the natural frequencies and transversal responses of simply supported single-walled carbon nanotubes. The influence of addition axial load and the properties of elastic medium on the vibrations are discussed. The results showed that the effects of addition axial load on the lower natural frequencies of single-walled carbon nanotubes are sensitive to the lower vibration modes and the stiff elastic medium. The lower natural frequencies depend on the axial load;they become smaller with increasing axial load and vary with the vibration modes. In addition, except for the first mode, the effects of the axial load on the stiff elastic medium are considerably greater than those on the flexible one. However, the constants of the elastic medium have little effect on the first mode. The critical axial buckling stress and strain for simply-supported single-walled carbon nanotubes are also obtained.

Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

ELASTIC MODULUS REDUCTION METHOD FOR LIMIT LOAD EVALUATION OF FRAME STRUCTURES

A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio (EBR),and the elastic modulus reduction method (EMRM) is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.

Discrete element method of improved performance of railway ballast bed using elastic sleeper

With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.