海洋立管-海床土体接触作用数值分析

钢悬链线立管(SCR)与海床土体的接触问题会使立管疲劳寿命急剧减少因而受到广泛关注。引入包含土体吸力效应的简化帽盖型P-y曲线,据此利用等效弹簧单元模拟土体并建立了立管与弹簧的有限元接触模型,计算并分析了不同立管竖向位移荷载和土体吸力最值下管道形态和弯矩的变化规律,讨论了着陆点(TDP)位置和着陆区(TDZ)深度与土体受力变化之间的关系。结果表明,立管弯矩在着陆点处变化最剧烈,在其两侧分别达到正负弯矩峰值,同时土体在相应位置出现吸力和承载力最大值,着陆区以后管土各变量基本无变化。

Longitudinal forces of continuously welded track on high-speed railway cable-stayed bridge considering impact of adjacent bridges

A proven beam-track contact model was used to analyze the track-structure interaction of CWR (continuously welded track) on bridge. Considering the impact of adjacent bridges, the tower-cable-track-beam-pier-pile finite element model of the cable-stayed bridge was established. Taking a bridge group including 40-32m simply-supported beam and (32+80+112)m single-tower cable-stayed bridge and 17-32m simply-supported beam on the Kunming-Shanghai high-speed railway as an example, the characteristics of CWR longitudinal force on the cable-stayed bridge were studied. It is shown that adjacent bridges must be considered in the calculation of the track expansion force and bending force on cable-stayed bridge. When the span amount of adjacent bridges is too numerous, it can be simplified as six spans; the fixed bearing of adjacent simply-supported beams should be placed on the side near the cable-stayed bridge; the track expansion device should be set at the bridge tower to reduce the track force near the bridge abutment.

Using FEM to predict tree motion in a wind field

In this paper we propose a finite element(FE) simulation method to predict tree motion in a wind field. Two FE tree models were investigated:One model was generated based on a realistic nature-looking geometric tree model,and the other was a symmetric model to investigate the influence of asymmetric material properties on tree motion. The vortex-induced vibration(VIV) theory is introduced to estimate the fluctuating wind force being exerted on tree stems and the fluid-structure interaction(FSI) analysis is also included in the simulation. The results indicate that asymmetric material properties result in the crosswind displacement of the investigated node and the main swaying direction deviation. The simulation reveals that under wind loading,a tree with leaves has much larger swaying amplitude along the wind direction and longer swaying period than a tree without leaves. However,the crosswind swaying amplitude is mainly due to branch interaction. The numerical simulation proved that the inter-action of tree branches can prevent dangerous swaying motion developing.

3D finite element analysis on pile-soil interaction of passive pile group

The interaction between pile and soft soil of the passive pile group subjected to soil movement was analyzed with three-dimensional finite element model by using ANSYS software. The soil was assumed to be elastic-plastic complying with the Drucker-Prager yield criterion in the analysis. The large displacement of soil was considered and contact elements were used to evaluate the interaction between pile and soil. The influences of soil depth of layer and number of piles on the lateral pressure of the pile were investigated, and the lateral pressure distributions on the (2×1) pile group and on the (2×2) pile group were compared. The results show that the adjacent surcharge may result in significant lateral movement of the soft soil and considerable pressure on the pile. The pressure acting on the row near the surcharge is higher than that on the other row, due to the "barrier" and arching effects in pile groups. The passive load and its distribution should be taken into account in the design of the passive piles.

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.

Numerical Simulation of Geotechnical Problems by Coupled Finite and Infinite Elements

in geotechnical engineering, numerical simulation of problems is of great importance. This work proposes a new formulation of coupled finite-infinite elements which can be used in numerical simulation ofgeotechnical problems in both static and dynamic conditions. Formulation and various implementation aspects of the proposed coupled finite-infinite elements are carefully discussed. To the authors＇ knowledge, this approach that considers coupled finite-infinite elements is more efficient in the sense that appropriate and accurate results are obtained by using less elements. The accuracy and efficiency of the proposed approach is considered by comparing the obtained results with analytical and numerical results. In a static case, the problem of circular domain ol infinite length is considered. In a dynamic case, one dimensional wave propagation problems arising from the Heaviside step fimction and impulse functions are considered. In order to get a more complete picture, two dimensional wave propagation in a circular qtmrter space is considered and the results are presented. Finally, a soil-structure interaction system subjected to seismic excitation is analyzed. In the analysis of soil-structure interaction phenomenon, frames with different number of storeys and soil media with various stiffness characteristics have been taken into consideration. In the analysis, the finite element software ANSYS has been used. For the newly developed infinite element, the programming has been done by the help of the User Programmable Features of the ANSYS software, which enable creating new elements in the ANSYS software.

Effect of crack face contact on dynamic stress intensity factors for a hole-edge crack

In order to determine the dynamic stress intensity factors(DSIFs)for a single edge crack at the center hole of a finite plate under a compressive step loading parallel to the crack,the finite element method was employed to solve the cracked plate problem.The square-root stress singularity around the crack tip was simulated by quarter point singular elements collapsed by 8-node two-dimensional isoparametric elements.The DSIFs with and without considering crack face contact situations were evaluated by using the displacement correlation technique,and the influence of contact interaction between crack surfaces on DSIFs was investigated.The numerical results show that if the contact interaction between crack surfaces is ignored,the negative mode I DSIFs may be obtained and a physically impossible interpenetration or overlap of the crack surfaces will occur.Thus the crack face contact has a significant influence on the mode I DSIFs.

The Influence of Surface Building's Relative Position on the Seismic Response of Subway Structures

By using the finite element method and viscoelastic artificial boundary, a soil-structure interaction system is established to simulate the influence of surface buildings on the seismic response of subway structures. The conditions of different relative positions between ground building and subway structure are analyzed. The result~ indicate that when considering the existence of surface buildings, the relative story displacements and internal forces of subway structures are changed greatly. Further the influence of surface buildings on subway structure changes as the distance increases.

Comparative analysis of seismic response characteristics of pile-supported structure

In order to analyze the seismic response characteristics of pile-supported structure,a computational model considering pile-soil-structure interaction effect was established by finite element method.Then,numerical implementation was made in time domain.At the same time,a simplified approximation for seismic response analysis of pile-soil-structure system was briefly presented.Furthermore,comparative study was performed for an engineering example.Through comparative analysis,it is shown that the results obtained by the simplified method well agree with those achieved by the finite element method.These results show that spectrum characteristics and intensity of input earthquakes are two important factors that can notablely influence the seismic response characteristics of superstructure.When the input ground motion acceleration amplitude gradually increases from 1 to 4 m/s2,the acceleration of pier top will increase,but it will not be simply proportional to the increase of input acceleration amplitude.

Static Analysis of Buried Pipes Using Coupling between Layerwise Finite Element and Boundary Element Method

This work deals with the analysis of soil-structure interaction modeling of pipeline problems in static behavior using the coupling between FEM （finite element method） and BEM （boundary element method）. The representation of the pipe is made by MEF using one fmite element in the cylindrical panel formulated from the theory of equivalent discrete layers （Layerwise theory）, proposed by J. N. Reddy. The soil is represented by elastic continum infimite or semi-infinite and modeled using boundary elements with special curved surface, associated with cylindrical panel used to represent the soil-structure interaction within the soil, especially at the contact surface with the pipe.

On the Lateral Behavior of the Backfill of a Skew Abutment

It is generally accepted that the interaction between a bridge and its abutment＇s backfill soil is highly nonlinear, especially under a strong earthquake loading that contains a velocity pulse. For bridges with skew abutments, the superstructure-abutment interaction can dominate the overall bridge performance. This study puts forth a new approach for predicting the lateral capacity of a skew abutment using verified high-fidelity three-dimensional continuum finite element （FE） models. The core idea is that the lateral capacity of a straight abutment is bounded from above and below by that of the abutment of a skew bridge that has the same deck-width, and that of another skew bridge （with the same angle） that has the same backwall length as the original/straight bridge, respectively. This postulation is then used in reverse to estimate the lateral capacity of a skew abutment, given the capacity of a straight but otherwise identical one with an arbitrary length. In prior research, the latter information had already been obtained in closed-form expressions that use physical parameters, such as backfill cohesion, internal friction angle and density, backwall height, and backwall-backfill friction angle. The approach presented here is constrained by the assumption that bridge deck will not rotate during loading. While this assumption is generally violated in a strong earthquake--because a skew bridge will tend to rotate, especially if its in-plane torsional rigidity is low, the model presented does serve as an anchor for parameterizing more advanced （e.g., macro-element plasticity） models that allow rotation, and also as fully parametric lateral response models for torsionally stiff （ile., multi-span, multi-bent） skew bridges.

Influence of Structure Plane Size on Seismic Response of Soil-Structure Interaction

The influence of the change of structure plane size on seismic response was studied for a soil-structure interaction system.Based on the finite element method,a soil-structure interaction calculation model was established to analyze the seismic response by changing the structure plane size and choosing different earthquake waves for different soil fields.The results show that when the natural periods of vibration for different structure plane sizes are close,under the same earthquake wave,the total displacement on the top layer of the structure and the foundation rotation displacement decrease with the increase of structure plane size,and the proportion of superstructure elastic selfdeformation displacement to the total displacement increases with the increase of structure plane size.While for different types of sites and seismic waves,under the horizontal and vertical seismic waves,the seismic responses of different plane sizes have a similar change rule.

A Brief Summary of Finite Element Method Applications to Nonlinear Wave-structure Interactions

We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.

Numerical Simulation of Multiple Collisions Between Ice Ridge and Polar Ship

Ships navigating in ice-covered regions will inevitably collide with ice ridges.Compared to other ice bodies,ice ridges exhibit more complicated mechanical behaviors due to the scale and structure characteristics.In this paper,nonlinear finite element method is used to investigate the interaction between a polar ship and an ice ridge.The ice ridge is modelled as elastic-plastic material based on Drucker-Prager yield function,with the consideration of the influence of cohesion,friction angle and material hardening.The material model is developed in LS-DYNA and solved using semi-implicit mapping algorithm.The stress distribution of ice ridge and ship,and the ice load history are evaluated through the simulation of multiple collisions.In addition,parametric analysis is performed to investigate the influence of ridge thickness and impact velocity on the ice load and energy absorption.

Investigation on the Forced Response of a Radial Turbine under Aerodynamic Excitations

Rotor blades in a radial turbine with nozzle guide vanes typically experience harmonic aerodynamic excitations due to the rotor stator interaction. Dynamic stresses induced by the harmonic excitations can result in high cycle fatigue(HCF) of the blades. A reliable prediction method for forced response issue is essential to avoid the HCF problem. In this work, the forced response mechanisms were investigated based on a fluid structure interaction(FSI) method. Aerodynamic excitations were obtained by three-dimensional unsteady computational fluid dynamics(CFD) simulation with phase shifted periodic boundary conditions. The first two harmonic pressures were determined as the primary components of the excitation and applied to finite element(FE) model to conduct the computational structural dynamics(CSD) simulation. The computed results from the harmonic forced response analysis show good agreement with the predictions of Singh's advanced frequency evaluation(SAFE) diagram. Moreover, the mode superposition method used in FE simulation offers an efficient way to provide quantitative assessments of mode response levels and resonant strength.

Study on interaction between induced and natural fractures by extended finite element method

Fracking is one of the kernel technologies in the remarkable shale gas revolution. The extended finite element method is used in this paper to numerically investigate the interaction between hydraulic and natural fractures, which is an important issue of the enigmatic fracture network formation in fracking. The criteria which control the opening of natural fracture and crossing of hydraulic fracture are tentatively presented. Influence factors on the interaction process are systematically analyzed, which include the approach angle, anisotropy of in-situ stress and fluid pressure profile.