A hybrid contact approach for modeling soil-structure interaction using the material point method

The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.

Fluid-solid Interaction Model for Hydraulic Reciprocating O-ring Seals

Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. Only elastic deformations of hydraulic reciprocating seals were discussed, and hydrodynamic effects were neglected in many studies. The physical process of the fluid-solid interaction effect did not be clearly presented in the existing fluid-solid interaction models for hydraulic reciprocating O-ring seals, and few of these models had been simultaneously validated through experiments. By exploring the physical process of the fluid-solid interaction effect of the hydraulic reciprocating O-ring seal, a numerical fluid-solid interaction model consisting of fluid lubrication, contact mechanics, asperity contact and elastic deformation analyses is constructed with an iterative procedure. With the SRV friction and wear tester, the experiments are performed to investigate the elastohydrodynamic lubrication characteristics of the O-ring seal. The regularity of the friction coefficient varying with the speed of reciprocating motion is obtained in the mixed lubrication condition. The experimental result is used to validate the fluid-solid interaction model. Based on the model, The elastohydrodynamic lubrication characteristics of the hydraulic reciprocating O-ring seal are presented respectively in the dry friction, mixed lubrication and full film lubrication conditions, including of the contact pressure, film thickness, friction coefficient, liquid film pressure and viscous shear stress in the sealing zone. The proposed numerical fluid-solid interaction model can be effectively used to analyze the operation characteristics of the hydraulic reciprocating O-ring seal, and can also be widely used to study other hydraulic reciprocating seals.

Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries

The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.

Land-Air Interaction over Arid/Semi-arid Areas in China and Its Impact on the East Asian Summer Monsoon. Part I:Calibration of the Land Surface Model (BATS)Using Multicriteria Methods

To improve the land surface simulation in the arid and semi-arid areas of northern China, the observational data from two field experiments in Dunhuang and Tongyu are used to optimize the parameters in the land surface model, BATS, through calibration with the multicriteria method. Sensitivity analysis to the parameters in Dunhuang and Tongyu indicates that different parameters need to be calibrated in two sites with different environmental and climate regimes. Comparison of observed sensible heat flux, latent heat flux, and ground surface temperature with the simulated ones shows the simulations with the optimized parameters have been substantially improved. Especially, the holistic simulations with the calibration of the parameter values are much closer to the observations in the arid region （Dunhuang）, and the energy partition with the calibrated parameters can also be simulated well in the semi-arid region （Tongyu）. Whole results demonstrate that the parameter calibration of the land surface model is important when the model is to be used to investigate the land-air interaction.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.

Numerical simulation of unsteady interaction of centrifugal impeller with its diffuser using Lagrangian discrete vortex method

In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.

A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

Fluid-solid interaction of resistance loss of flexible hose in deep ocean mining

The resistance loss of transportation was studied and the influences of buoyancy layout,mineral content and elastic modulus of flexible hose were investigated based on three-dimensional finite element model of fluid-solid interaction by MSC.MARC/MENTAT software.The numerical results show that the resistance losses increase with the increase of mineral content Cv and velocity of internal fluid v and decrease with the increase of elastic modulus E of flexible hose.The buoyancy layout and the velocity of internal fluid have greater impacts on the resistance losses than the elastic modulus of flexible hose.In order to reduce the resistance losses and improve the efficiency of the deep-ocean mining,Cv and v must be restricted in a suitable range (e.g.10%-25% and 2.5-4 m/s).Effective buoyancy layout (such as Scheme C and D) should be adopted and the suitable material of moderate E should be used for the flexible hose in deep-ocean mining.

MPS-FEM Coupled Method for Study of Wave-Structure Interaction

Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.

Alterable-element Method for Vehicle-bridge Interaction Considering the Transient Jump of Wheel

The so called "alterable-element method" (AEM) was introduced to deal with the coupling interac-tion of vehicle and sub-structure considering the actual transient jump of wheel, while the classical "contact all along" assumption based on which wheels and lower structure are always contact was abandoned. The alterable element used in this method is a conceptional element, which is used to calculate the coupling interaction of upper and lower structures and has some typical characteristics: firstly it flows along with the moving of contact point; secondly whether it is used for calculation depends on the contact state; thirdly its sizes could change according to specific problems and so on. VISUAL FORTRAN program was coded, and different moving vehicle models were presented taking into consideration the effects of random corrugation in the numerical study. The numerical solutions are favored comparing with the results obtained by alternative methods when there is no jump phenomenon existed. With abrupt irregularity, the transient jump of wheel was studied using the present method.

OpenIFEM:A High Performance Modular Open-Source Software of the Immersed Finite Element Method for Fluid-Structure Interactions

We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

Numerical investigation of the effects of soil-structure and granular material-structure interaction on the seismic response of a flat-bottom reinforced concrete silo

In this work,a numerical study of the effects of soil-structure interaction(SSI)and granular material-structure interaction(GSI)on the nonlinear response and seismic capacity of flat-bottomed storage silos is conducted.A series of incremental dynamic analyses(IDA)are performed on a case of large reinforced concrete silo using 10 seismic recordings.The IDA results are given by two average IDA capacity curves,which are represented,as well as the seismic capacity of the studied structure,with and without a consideration of the SSI while accounting for the effect of GSI.These curves are used to quantify and evaluate the damage of the studied silo by utilizing two damage indices,one based on dissipated energy and the other on displacement and dissipated energy.The cumulative energy dissipation curves obtained by the average IDA capacity curves with and without SSI are presented as a function of the base shear,and these curves allow one to obtain the two critical points and the different limit states of the structure.It is observed that the SSI and GSI significantly influence the seismic response and capacity of the studied structure,particularly at higher levels of PGA.Moreover,the effect of the SSI reduces the damage index of the studied structure by 4%.

The Boundary Element Method for Pile-Soil Interaction

In this paper, several mathmatical models for the pile- soil interaction are outlined. The Boundary Element Method is one of the very effective methods for the reasonable models of elasticity and elastoplasticity. The major of this paper is concerned with the Boundary Element Method for the pile-soil interaction, including general methods and calculating formulation of static and dynamic analysis of the pile and pile groups. Some results of analysis are also given.

Improved Simulation Method for Soil-Geogrid Interaction of Reinforced Earth Structure in FEM

The interaction between geogrid and soil is crucial for the stability of geogrid-reinforced earth structure. In finite element （FE） analysis, geogrids are usually assumed as beam or truss elements, and the interaction between geogrid and soil is considered as Coulomb friction resistance, which cannot reflect the true stress and displacement developed in the reinlbrcement. And the traditional Lagrangian elements used to simulate soil always become highly distorted and lose accuracy in high-stress blocks. An improved geogrid model that can produce shear resistance and passive resistance and a soil model using the Eulerian technique, in combination with the coupled Eulerian-Lagrangian （CEL） method, are used to analyze the interaction between geogrid and soil of reinforced foundation test in ABAQUS. The stress in the backfill, resistance of geogrid, and settlement of foundation were computed and the results of analysis agree well with the experimental results. This simulation method is of referential value for FE analysis of reinforced earth structure.

Opinion consensus incorporating higher-order interactions in individual-collective networks

In the current information society, the dissemination mechanisms and evolution laws of individual or collective opinions and their behaviors are the research hot topics in the field of opinion dynamics. First, in this paper, a two-layer network consisting of an individual-opinion layer and a collective-opinion layer is constructed, and a dissemination model of opinions incorporating higher-order interactions(i.e. OIHOI dissemination model) is proposed. Furthermore, the dynamic equations of opinion dissemination for both individuals and groups are presented. Using Lyapunov's first method,two equilibrium points, including the negative consensus point and positive consensus point, and the dynamic equations obtained for opinion dissemination, are analyzed theoretically. In addition, for individual opinions and collective opinions,some conditions for reaching negative consensus and positive consensus as well as the theoretical expression for the dissemination threshold are put forward. Numerical simulations are carried to verify the feasibility and effectiveness of the proposed theoretical results, as well as the influence of the intra-structure, inter-connections, and higher-order interactions on the dissemination and evolution of individual opinions. The main results are as follows.(i) When the intra-structure of the collective-opinion layer meets certain characteristics, then a negative or positive consensus is easier to reach for individuals.(ii) Both negative consensus and positive consensus perform best in mixed type of inter-connections in the two-layer network.(iii) Higher-order interactions can quickly eliminate differences in individual opinions, thereby enabling individuals to reach consensus faster.

ALE Fractional Step Finite Element Method for Fluid-Structure Nonlinear Interaction Problem

A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.

Numerical method of studying nonlinear interactions between long waves and multiple short waves

Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we present a numerical method of studying nonlinear interactions between a long wave and multiple short harmonic waves in infinitely deep water. Specifically, this method is applied to the calculation of the temporal and spatial evolutions of the surface elevations in which a given long wave interacts with several short harmonic waves. Another important application of our method is to quantitatively analyse the nonlinear interactions between an arbitrary short wave train and another short wave train. From simulation results, we obtain that the mechanism for the nonlinear interactions between one short wave train and another short wave train （expressed as wave train 2） leads to the energy focusing of the other short wave train （expressed as wave train 3）. This mechanism occurs on wave components with a narrow frequency bandwidth, whose frequencies are near that of wave train 3.