Three-dimensional finite element analysis on effects of tunnel construction on nearby pile foundation

A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

Evaluation value of three-dimensional finite element model analysis for bone mineral density and bone metabolism activity in patients with osteoporosis

Objective: To study the evaluation value of three-dimensional finite element model analysis for bone mineral density (BMD) and bone metabolism activity in patients with osteoporosis. Methods: A total of 218 patients who were diagnosed with osteoporosis in the hospital between February 2014 and January 2017 were collected as observation group, and 100 healthy volunteers who received physical examination in the hospital during the same period were selected as normal control group. The femoral head of the two groups was analyzed by three-dimensional finite element model, and the femoral head BMD levels and serum bone metabolism index contents were measured. Pearson test was used to evaluate the evaluation value of femoral head three-dimensional finite element model for osteoporosis. Results: The cancellous bone and cortical bone Von Mises stress value of observation group were lower than those of normal control group, and femoral neck BMD value of observation group was lower than that of normal control group;serum bone metabolism index BGP content was lower than that of normal control group while NBAP, TRACP-5b and CTX-1 contents were higher than those of normal control group. Pearson test showed that the cancellous bone and cortical bone Von Mises stress value of patients with osteoporosis were directly correlated with BMD value and bone metabolism index contents. Conclusion: The three-dimensional finite element model analysis resultsof patients with osteoporosis can objectively reflect the femoral headBMD value and bone metabolism activity, and is a reliable way to evaluate the risk of long-term fractures.

Three-dimensional nonlinear analysis of creep in concrete filled steel tube columns

This paper proposes a based on 3D-VLE （three-dimensional nonlinear viscoelastic theory） three-parameters viscoelastic model for studying the time-dependent behaviour of concrete filled steel tube （CFT） columns. The method of 3D-VLE was developed to analyze the effects of concrete creep behavior on CFT structures. After the evaluation of the parameters in the proposed creep model, experimental measurements of two prestressed reinforced concrete beams were used to investigate the creep phenomenon of three CFT columns under long-term axial and eccentric load was investigated. The experimentally obtained time-dependent creep behaviour accorded well with the cu~＇es obtained from the proposed method. Many factors （such as ratio of long-term load to strength, slenderness ratio, steel ratio, and eccentricity ratio） were considered to obtain the regularity of influence of concrete creep on CFT structures. The analytical results can be consulted in the engineering practice and design.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Importance of Three-Dimensional Piezoelectric Coupling Modeling in Quantitative Analysis of Piezoelectric Actuators

This paper demonstrates the importance of three-dimensional(3-D)piezoelectric coupling in the electromechanical behavior of piezoelectric devices using three-dimensional finite element analyses based on weak and strong coupling models for a thin cantilevered piezoelectric bimorph actuator.It is found that there is a significant difference between the strong and weak coupling solutions given by coupling direct and inverse piezoelectric effects(i.e.,piezoelectric coupling effect).In addition,there is significant longitudinal bending caused by the constraint of the inverse piezoelectric effect in the width direction at the fixed end(i.e.,3-D effect).Hence,modeling of these effects or 3-D piezoelectric coupling modeling is an electromechanical basis for the piezoelectric devices,which contributes to the accurate prediction of their behavior.

Dynamic Analysis of Submarine Pipelines on an Elastic / Viscoelastic Foundation

This paper presents a new finite element method for solving static and dynamic problems in laying operation of pipelines. The effect of the viscoelastic soil behavior is considered by using the Pasternak foundation model. Some examples are also presented.

Viscoelastic relaxation of the upper mantle and afterslip following the 2014 M_(W)8.1 Iquique earthquake

An improved understanding of postseismic crustal deformation following large subduction earthquakes may help to better understand the rheological properties of upper mantle and the slip behavior of subduction interface.Here we construct a three-dimensional viscoelastic finite element model to study the postseismic deformation of the 2014 M_(W)8.1 Iquique,Chile earthquake.Elastic units in the model include the subducting slab,continental and oceanic lithospheres.Rheological units include the mantle wedge,the oceanic asthenosphere and upper mantle.We use a 2 km thick weak shear zone attached to the subduction fault to simulate the time-dependent stress-driven afterslip.The viscoelastic relaxation in the rheological units is represented by the Burgers rheology.We carry out grid-searches on the shear zone viscosity,thickness and viscosity of the asthenosphere,and they are determined to be 10^(17)Pa s,110 km and 2×10^(18)Pa s,respectively.The stress-driven afterlsip within the first two years is up to~47 cm and becomes negligible after two years(no more than 5 cm/yr).Our results suggest that a thin,low-viscosity oceanic asthenosphere together with a weak shear zone attached to the fault are required to better reproduce the observed postseismic deformation.

Influences of the heterogeneity of viscoelastic medium on postseismic deformation of the 2008 M_(W) 7.9 Wenchuan earthquake

The study of postseismic deformation is important for constraining the viscoelastic properties of the Earth and inverting the post-earthquake process.The levelling survey revealed that the area near Bei-chuan elevated 5.3 cm about two years after the M_(W) 7.9 Wenchuan earthquake(05/12/2008),during which the area underwent significant downward movement.The GPS horizontal displacements showed a non-monotonic variation after the Wenchuan earthquake.In this study,a 3-D viscoelastic finite element model is employed to simulate the coseismic and postseismic deformation of the Wenchuan earthquake.The numerical simulations show that the lateral heterogeneity across the Longmenshan fault plays an important role in the postseismic displacements.The results reveal that the coseismic defor-mation is not sensitive to the horizontal heterogeneity,but the postseismic deformation is sensitive to it.The postseismic deformation of the horizontally heterogeneous model is generally consistent with the observations of all geodetic surveys,such as GPS,InSAR and levelling,but not for the horizontally homogenous model.We also find that the non-monotonous variation of the postseismic deformation of the Wenchuan earthquake could be explained by a viscoelastic relaxation model with lateral heterogeneous medium across the Longmenshan fault.

VISCOELASTIC CONSTITUTIVE MODEL RELATED TO DEFORMATION OF INSECT WING UNDER LOADING IN FLAPPING MOTION

Flexible insect wings deform passively under the periodic loading during flapping flight. The wing flexibility is considered as one of the specific mechanisms on improving insect flight performance. The constitutive relation of the insect wing material plays a key role on the wing deformation, but has not been clearly understood yet. A viscoelastic constitutive relation model was established based on the stress relaxation ex- periment of a dragonfly wing （in vitro）. This model was examined by the finite element analysis of the dynamic deformation response for a model insect wing under the action of the periodical inertial force in flapping. It is revealed that the viscoelastic constitutive relation is rational to characterize the biomaterial property of insect wings in contrast to the elastic one. The amplitude and form of the passive viscoelastic deformation of the wing is evidently dependent on the viscous parameters in the constitutive relation.

Stress Distribution in Maxillary Central Incisors without Ferrules: A Finite Element Analysis of Post and Core Systems

Purpose: The purpose of this study was to identify optimal post and core materials for central incisors without ferrules using three-dimensional finite element analysis and three-point bending tests. Methods: Stress analyses were performed with six models: cast metal post and core (MP), composite resin core alone, straight fiber-reinforced post-composite resin core (FSR), tapered fiber-reinforced post-composite resin core, straight titanium post-composite resin core (TSR), and tapered titanium post-composite resin core (TTR). A 100-N load was applied to the lingual surface at a 45° angle to the long axis of the tooth. Maximum von Mises stress distributions were calculated with finite element analysis software. Five samples each of composite resin, straight fiber-reinforced post, straight titanium post, straight fiber-reinforced post and composite resin, and straight titanium post and composite resin were subjected to three-point bending tests, followed by analysis of variance and Tukey’s multiple comparison test. Results: Stress distribution was optimal on TTR. Maximum von Mises stress on the cervical side of the post was greatest in TSR (693 MPa) and TTR (556 MPa). Maximum stress on the apical side of the post was greatest in MP (110 MPa). Maximum stress in surrounding dentin was lowest in MP (203 MPa) and TTR (250 MPa). Gap distance was smallest in MP (0.09 mm) and largest in FSR (0.26 mm). Mean maximum three-point bending force was lowest in composite resin (26.9 N/mm) and highest in titanium post and composite resin (97.1 N/mm). Titanium post bending strength was consistently greater than that of the fiber-reinforced post (p < 0.01). Conclusion: These results revealed optimal stress distribution and high bending strength with the tapered titanium post and resin combination, suggesting that this combination can most effectively prevent root or post fracture in an anterior tooth without a ferrule.

Thermomechanical analysis of triangular zone cracks in vertical continuous casting slabs based on viscoelastic-plastic model

The triangular zone cracks in 2101 duplex stainless steel produced by the vertical continuous caster have troubled company A for a long time. To simulate the temperature and thermal stress distributions in the solidification process of 2101 duplex stainless steel produced by the vertical continuous caster, a two-dimensional viscoelastic-plastic thermomechanically coupled finite element model was established by the secondary development of the commercial nonlinear finite element analysis software MSC Marc. The results show that the thermal stress on the surface reaches a maximum at the exit of the mould, and the highest thermal stresses at the centre of the wide face and the narrow face are 75 and 115 MPa, respectively. Meanwhile, the internal temperature of slab is still higher than the solidus temperature, resulting in no thermal stress. The slab shows different high-temperature strengths and suffers from different stresses at different positions; thus, the risk of cracking also varies. At a location of 6-8 m from the meniscus, the temperature of the triangular zone is 1270-1360℃ and the corresponding permissible high-temperature strength is about 10-30 MPa, while the thermal stress at this time is 60 MPa, which is higher than the high-temperature strength. As a result, triangular zone cracks form easily.

Cell membrane tensile strain under cyclic compression:A viscoelastic myoblast finite element model

Cyclic mechanical stimulation could lead to subsequent biomechanical and biological effects on cells.Viscoelastic cells could deform or show energy dissipation with hysteresis behavior in response to external cyclic compression.The aim of this study was to investigate the effect of cyclic compression on a single viscoelastic myoblast through a confocal–based cell–specific finite element model,including cell membrane tensile strain and damaged elements.Sinusoidal compression was applied to the apical surface of the myoblast(cell membrane)with compressive stress of 500
500 Pa(with stress offset at 500 Pa and amplitude of 500 Pa)at 0 Hz(static compression of 500 Pa),0.25 Hz,0.5 Hz,0.75 Hz,1 Hz,5 Hz,and 10 Hz.Results showed that the ratio of average tensile strain integral in all cell membrane elements over a certain period of time(T)to that duration(T)(MAS index)decreased under cyclic compression compared to that of static compression in the short term(within 4 s).Furthermore,compared to static compression,the percentage of damaged elements of cell membrane under cyclic compression decreased assuming a 3%cell membrane tensile strain damage threshold.The optimal cyclic compression frequency 0.25 Hz led to the largest difference of MAS index under cyclic compression and static compression.These results may provide support for the application of cyclic compressive stimulation in the prevention of cell damage.

A New Method for 3D Finite Element Modeling of Human Mandible Based on CT Data

This study presents a reliable method for the semi-automatic generation of an FE model, which determines both geometrical data and bone properties from patient CT scans.3D FE analysis is one of the best approaches to predict the stress and strain distribution in complex bone structures, but its accuracy strongly depends on the precision of input information. In geometric reconstruction, various methods of image processing, geometric modeling and finite element analysis are combined and extended. Emphasis is given to the assignment of the material properties based on the density values computed from CT data. Through this technique, the model with high geometric and material similarities were generated in an easy way. Consequently, the patient-specific FE model from mandible CT data is realized also.

MATHEMATICAL MODEL OF MILLING TEMPERATURE AND TEMPERATURE FIELD ANALYSIS FOR THREE-DIMEN-SIONAL COMPLEX GROOVE MILLING INSERT

The mathematical model on the temperature of the waved-edge is constructedaccording to Jaeger's theory of moving solid and based on the used temperature model of the flatinsert. It is possible to forecast the milling temperature through programming. The comparableexperiments have been done between the two new three-dimension groove inserts (waved-edge insert,great edge insert) and flat fake insert. The theoretic forecast is in good agreement with theexperimental result. According to the cutting conditions, the boundary condition of finite elementanalysis on cutting temperature field is established, and the three-dimensional temperature fieldsof inserts with grooves are analyzed by FEM, so as to offer a reference basis for the design andoptimization of insert grooves.

Mechanical analysis of the femoral neck dynamic intersection system with different nail angles and clinical applications

BACKGROUND The femoral neck dynamic intersection system(FNS)is mechanically more stable than other internal fixation techniques.Current studies have confirmed that the structural design of FNS has good biomechanical properties in European and American populations.However,whether the suitability of the FNS's 130°main nail angle design for Asian populations has been thoroughly investigated remains unclear.AIM To compare the biomechanical stability differences among different main nail angles of the FNS in the treatment of femoral neck fractures in Asian populations.METHODS Computed tomography data of the femur of healthy adult male volunteers were imported into Mimics software to create a three-dimensional model of the femur.The model was adapted to the curve using Geomagic software and imported into Solidworks software to construct the Pauwels I femoral neck fracture model and design the FNS internal fixation model using different main nail angles.Afterward,the models were assembled with the FNS fracture model and meshed using the preprocessing Hypermesh software.Subsequently,they were imported into Abaqus software to analyze and evaluate the biomechanical effects of different angles of the FNS main nail on the treatment of femoral neck fractures.RESULTS The peak displacement of the proximal femur under different angles of FNS fixation under stress was 7.446 millimeters in the 120°group and 7.416 millimeters in the 125°group;in the 130°,135°,and 140°FNS fixation groups,the peak displacement was 7.324 millimeters,8.138 millimeters,and 8.246 millimeters,respectively.In the 120°and 125°FNS fixation groups,the maximum stresses were concentrated at the main nail and the anti-rotation screw,which intersected the fracture line of the femur neck,resulting in peak stresses of 200.7 MPa and 138.8 MPa,respectively.Peak stresses of 208.8 MPa,219.8 MPa,and 239.3 MPa were observed on the angular locking plate distal to the locking screw in the 130°,135°,and 140°fixation groups.CONCLUSION FNS has significant stress distribution properties,a minimal proximal femoral displacement,and an optimal stability for treating femoral neck fractures in Asian populations when performed with a 130°main nail angle.

Stochastic analysis of excavation-induced wall deflection and box culvert settlement considering spatial variability of soil stiffness

In this study, a three-dimensional (3D) finite element modelling (FEM) analysis is carried out to investigate the effects of soil spatial variability on the response of retaining walls and an adjacent box culvert due to a braced excavation. The spatial variability of soil stiffness is modelled using a variogram and calibrated by high-quality experimental data. Multiple random field samples (RFSs) of soil stiffness are generated using geostatistical analysis and mapped onto a finite element mesh for stochastic analysis of excavation-induced structural responses by Monte Carlo simulation. It is found that the spatial variability of soil stiffness can be described by an exponential variogram, and the associated vertical correlation length is varied from 1.3 m to 1.6 m. It also reveals that the spatial variability of soil stiffness has a significant effect on the variations of retaining wall deflections and box culvert settlements. The ignorance of spatial variability in 3D FEM can result in an underestimation of lateral wall deflections and culvert settlements. Thus, the stochastic structural responses obtained from the 3D analysis could serve as an effective aid for probabilistic design and analysis of excavations.

蜂窝夹层结构中粘结剂层的建模与仿真

为分析夹层结构制造过程中粘结剂层对表面质量的影响,以聚甲基丙烯酸甲脂(PMMA)为主要研究对象,在有限元分析软件ABAQUS中进行建模分析,通过加热和加压2个过程研究在夹层结构的热压制造过程中粘结剂对表面质量的影响,结果表明:当仅考虑加热过程时,由于粘结剂的粘弹性,随着温度的升高,PMMA变软,导致表面波纹度增大;当仅考虑加压过程时,粘弹性对表面波纹度几乎无影响;当热压耦合时,表面波纹度比单独考虑加热或加压的更大;粘结剂越软,表面波纹度越大。通过调整粘结剂的特性,可以有效控制夹层结构的表面质量,从而提高产品的整体质量。

移动荷载作用下正交各向异性地基—黏弹性路面结构动力响应的数值模拟

【目的】沥青路面结构在重载车辆的作用下易产生路面典型损害,通过ABAQUS软件进行数值模拟,研究常温状态下半刚性基层沥青路面结构在实际车辆荷载作用下的动力响应。【方法】将半刚性基层沥青路面道路结构简化为面层-基层-路基3层结构,沥青面层黏弹性采用Burgers模型描述,路基土体通过正交各向异性弹性介质模拟,路面车辆荷载采用移动简谐荷载进行模拟。【结果】路基土体正交各向异性,荷载移动速度、振动频率对道路表面竖向位移及竖向正应力有不同程度的影响。【结论】考虑土体正交各向异性能更准确地描述路基路面的动力响应特性,荷载振动频率、荷载移动速度会对路面竖向位移和竖向正应力产生较为明显的影响。