某混凝土面板堆石坝施工期变形特性分析

面板堆石坝具有地形适应力强、可利用当地材料筑坝、经济性好等特点,堆石体作为坝体重要组成部分,其变形对坝体结构稳定、面板防裂等具有较大影响。文章以某混凝土面板堆石坝为实例,建立三维有限元计算模型,采用有限元数值计算方法,对该堆石坝施工期变形特性进行模拟分析。结果表明:当填筑完成后,最大沉降值约为1.8m,出现在Ⅱ-1分期填筑的分界面;顺河向和横向位移均呈现出一定的对称性;受到地形的不对称影响,靠近较陡峭的左岸的堆石体变形更大。

Three-dimensional simulation of pore scale fluid flow in granular ore media with realistic geometry

The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.

Finite element analysis of contact fatigue and bending fatigue of a theoretical assembling straight bevel gear pair

The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulative fatigue criterion and the stress-life equation, the key meshing states of the gear pair were investigated for the contact fatigue and the bending fatigue. Then, the reliability of the proposed model was proved by comparing the calculation result with the simulation result. Further study was performed to analyze the variation of the contact fatigue stress and the bending fatigue stress under different loads. Furthermore, the roles of the driving pinion and the driven gear pair were evaluated in the fatigue life of the straight bevel gear pair and the main fatigue failure mode was determined for the significant gear. The results show that the fatigue failure of the driving pinion is the main fatigue failure for the straight bevel gear pair and the bending fatigue failure is the main fatigue failure for the driving pinion.

Numerical simulation of temperature distribution and heat transfer during solidification of titanium alloy ingots in vacuum arc remelting process

In order to get a better understanding of the vacuum consumable arc remelting(VAR) processes and thus to optimize them,a 3D finite element model was developed for the temperature fields and heat transfer of titanium alloy ingots during VAR process.The results show that the temperature fields obtained by the simulation are well validated through the experiment results.The temperature distribution is different during the whole VAR process and the steady-state molten pool forms at 329 s for d100 mm × 180 mm ingots.At the initial stage of remelting,the heat dissipation of crucible bottom plays an important role in the whole heat dissipation system.At the middle of remelting,the crucible wall becomes a major heat dissipation way.The effect of cooling velocity on the solidification structure of ingots was investigated based on the temperature fields and the results can well explain the macrostructure of titanium alloy ingots.

Cross-sectional distortion behaviors of thin-walled rectangular tube in rotary-draw bending process

The cross-sectional distortion usually appears during rotary-draw bending process of thin-walled rectangular tube with small bending radius.To study the cross-sectional distortion of the tube,a three-dimensional finite-element model of the process was developed based on ABAQUS/Explicit code and its reliability was validated by experiment.Then,the cross-sectional distortion behaviors of the tube were investigated.The results show that a zone of larger circumferential stress appears on the tube when bending angle reaches 30°.And in the larger circumferential stress zone,the sagging phenomenon is produced obviously.The maximum cross-sectional distortion is located in the larger circumferential stress zone and the angle between the plane of maximum cross-sectional distortion and the bending reference plane is about 50°.The position of the maximum cross-sectional distortion keeps almost unchanged with the variation of the clearances between dies and tube.

Numerical simulation of tire/soil interaction using a verified 3D finite element model

The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.

Nonlinear finite element analysis of effect of seismic waves on dynamic response of Shiziping dam

Many high earth-rockfill dams are constructed in the west of China. The seismic intensity at the dam site is usually very high, thus it is of great importance to ensure the safety of the dam in meizoseismal area. A 3D FEM model is established to analyze the seismic responses of Shiziping earth-rockfill dam. The nonlinear elastic Duncan-Chang constitutive model and the equivalent viscoelastic constitutive model are used to simulate the static and dynamic stress strain relationships of the dam materials, respectively. Four groups of seismic waves are inputted from the top of the bedrock to analyze the dynamic responses of the dam. The numerical results show that the calculated dynamic magnification factors display a good consistency with the specification values. The site spectrum results in larger acceleration response than the specification spectrum. The analysis of relative dynamic displacement indicates that the displacement at the downstream side of the dam is larger than that at the upstream side. The displacement response reduces from the center of river valley to two banks. The displacement responses corresponding to the specification spectrum are a little smaller than those corresponding to the site spectrum. The analysis of shear stress indicates that a large shear stress area appears in the upstream overburden layer, where the shear stress caused by site waves is larger than that caused by specification waves. The analysis of dynamic principal stress indicates that the minimum dynamic stresses in corridor caused by specification and site waves have little difference. The maximum and minimum dynamic stresses are relatively large at two sides. The largest tensile stress occurs at two sides of the floor of grouting corridor, which may result in the crack near the corridor side. The numerical results present good consistency with the observation data of the grouting corridor in Wenchuan earthquake.

Flexural response of reinforced concrete beams strengthened with post-poured ultra high toughness cementitious composites layer

Four-point bending tests were conducted up to failure on eleven reinforced concrete (RC) beams and strengthening beams to study the effectiveness of externally pouring ultra high toughness cementitious composites (UHTCC) on improving the flexural behavior of existing RC beams.The strengthening materials included UHTCC and high strength grade concrete.The parameters,such as thickness and length of strengthening layer and reinforcement in post-poured layer,were analyzed.The flexural behavior,failure mode and crack propagation of composite beams were investigated.The test results show that the strengthening layer improves the cracking and ultimate load by increasing the cross section area.Introducing UHTCC material into strengthening not only improves the bearing capacity of the original specimens,but also disperses larger cracks in upper concrete into multiple tightly-spaced fine cracks,thus prolonging the appearance of harm surface cracks and increasing the durability of existing structures.Compared with post-poured concrete,UHTCC is more suitable for working together with reinforcement.The load?deflection plots obtained from three-dimensional finite-element model (FEM) analyses are compared with those obtained from the experimental results,and show close correlation.

Establishment of a three-dimensional finite element model for gunshot wounds to the human mandible

Objective:To investigate the feasibility of a finite element model as an ideal research tool for human maxillofacial gunshot wounds.Methods:Mandible CT scan data on the Chinese Visible Human were imported into MIMICS software to obtain the surface mesh of the mandible.Then,these surface-meshed models were imported into ANSA software for automatic net generation.Elements and nodes were partitioned on the basis of the mesh to obtain a three-dimensional finite element model for the mandible with every internal parameter consistent with those of our previously developed model in the pig mandible.The finite element model was imported into LS-DYNA for computation.Finally,the LS-POST was used to complete the simulation and the measurements.Results:A three-dimensional finite element model was successfully established for gunshot wounds in the human mandible.The stress distribution and the degree of injury were simulated dynamically for shots from two types of projectiles in the mandible at one entry angle and at three impact velocities.Conclusion:Three-dimensional finite element models will become ideal research tools for treatment of ballistic wounds of the human maxillofacial region.Using this human mandibular model as a foundation,we will be able to successfully develop three-dimensional finite element models for human maxillofacial gunshot wounds.

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.

Modelling and simulation of contact force in cold rotary forging

Cold rotary forging is an advanced and complex metal forming technology with continuous local plastic deformation.Investigating the contact force between the dies and the workpiece has a great significance to improve the life of the dies in cold rotary forging.The purpose of this work is to reveal the contact force responses in cold rotary forging through the modelling and simulation.For this purpose,a 3D elastic-plastic dynamic explicit FE model of cold rotary forging is developed using the FE code ABAQUS/Explicit.Through the modelling and simulation,the distribution and evolution of the contact force in cold rotary forging is investigated in detail.The experiment has been conducted and the validity of the 3D FE model of cold rotary forging has been verified.The results show that: 1) The contact force distribution is complex and exhibits an obvious non-uniform characteristic in the radial and circumferential directions; 2) The maximum contact force between the upper die and the workpiece is much larger than that between the lower die and the workpiece; 3) The contact force evolution history is periodic and every period experiences three different stages; 4) The total normal contact force is much larger than the total shear contact force at any given time.

Effect of cooling structure on thermal behavior of copper plates of slab continuous casting mold

A three-dimensional finite-element model of slab continuous casting mold was conducted to clarify the effect of cooling structure on thermal behavior of copper plates. The results show that temperature distribution of hot surface is mainly governed by cooling structure and heat-transfer conditions. For hot surface centricity, maximum surface temperature promotions are 30 ℃and 15 ℃ with thickness increments of copper plates of 5 mm and nickel layers of 1 ram, respectively. The surface temperature without nickel layers is depressed by 10 ℃ when the depth increment of water slots is 2 mm and that with nickel layers adjacent to and away from mold outlet is depressed by 7℃ and 5 ℃, respectively. The specific trend of temperature distribution of transverse sections of copper plates is nearly free of cooling structure, but temperature is changed and its law is similar to the corresponding surface temperature.

Numerical analysis on Joule heating of double-loop channel induction furnaces

In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with experimental data from low load trials for a 400 kW inductor. The results,such as power factor and Joule heating power,show reasonable correlation with experimental data,and Joule heating rate reaches the maximum at the corners and the minimum at the centre of the cross-section area. With increasing relative permeability of iron core,length of coils,current frequency and resistivity of metal melt,the power factor and Joule heating power change. It is concluded that current frequency,the resistivity and length of the coil play a critical role in determining the power factor and Joule heating power,whereas relative permeability of the magnetic core shows no significant influence on them.

Simulation of hydrogen diffusion in welded joint of X80 pipeline steel

Hydrogen diffusion coefficients of different regions in the welded joint of X80 pipeline steel were measured using the electro-chemical permeation technique. Using ABAQUS software, hydrogen diffusion in X80 pipeline steel welded joint was studied in consideration of the inhomogeneity of the welding zone, and temperature-dependent thermo-physical and mechanical properties of the metals. A three dimensional finite element model was developed and a coupled thermo-mechanical-diffusion analysis was performed. Hydrogen concentration distribution across the welded joint was obtained. It is found that the postweld residual hydrogen exhibits a non-uniform distribution across the welded joint. A maximum equivalent stress occurs in the immediate vicinity of the weld metal. The heat affected zone has the highest hydrogen concentration level, followed by the weld zone and the base metal.Simulation results are well consistent with theoretical analysis.

Application of 1D/3D finite elements coupling for structural nonlinear analysis

An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.

A New 3D Meso-mechanical Modeling Method of Coral Aggregate Concrete Considering Interface Characteristics

On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.

Warpage prediction of the injection-molded strip-like plastic parts

For most strip-like plastic injection molded parts, whose cross section size is much smaller than their length, the traditional Hele-Shaw model and three-dimensional model do not work well in the prediction of the warpage be- cause of their special shape. A new solution was suggested in this work. The strip-like plastic part was regarded as a little-curved beam macrnscopically, and was divided into a few one-dimensional elements. On the section of each elemental node location, two-dimensional thermal finite element analysis was made to obtain the non- uniform thermal stress caused by the time difference of the solidification of the plastic melt in the mold. The stress relaxation, or equivalently, strain creep was dealt with by using a special computing model. On the bases of in-mold elastic stress, the final bending moment to the beam was obtained and the warpage was predict- ed in good a^reement with practical cases.

THE EFFECT OF POSTERIOR ELEMENT RESECTION ON THE STRESS DISTRIBUTION IN THE LUMBAR SPINEA

A three-dimensional finite element model of the lumbar motion segment was used to predict the stress distribution in lumbar spine with posterior element resection. It was shown that the stress level in all parts of the lumbar spine was elevated although the stress distribution remained unchanged. The authors concluded that the posterior element resection in lumbar surgery should be avoided as much as possible.

Numerical analysis on thermal regime in double-loop channel inductor

In order to investigate the temperature distribution, a three-dimensional finite element model （FEM） was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated results were compared with experimental data from low load trials of a 400 kW inductor. The results of numerical simulations, such as the temperature and Joule heating rate, show reasonable correlation with experimental data. The results indicate that Joule heating rate and the temperature reach the maximum at the comers and the minimum at the centre of the cross-section area. The temperature difference between the inlet and outlet is in an inverse proportion to mass transport. Joule heating rate and the temperature are directly proportional to power frequency. It is concluded that mass transport and power frequency play a critical role in determining the temperature regime and Joule heating rate, the relative permeability of the magnetic core shows no significant influence on temperature regime and Joule heating rate, when the relative permeability varies from 5 000 to 10 000.

Combined back-analysis method of ground stress based on refined geological modeling

A new back-analysis method of ground stress is proposed with comprehensive consideration of influence of topography, geology and nonlinear physical mechanical properties of rock on ground stress. This method based on non-uniform rational B-spline (NURBS) technology provides the means to build a refined three-dimensional finite element model with more accurate meshing under complex terrain and geological conditions. Meanwhile, this method is a back-analysis of ground stress with combination of multivariable linear regression model and neural network (ANN) model. Firstly, the regression model is used to fit approximately boundary loads. Regarding the regressed loads as mean value, some sets of boundary loads with the same interval are constructed according to the principle of orthogonal design, to calculate the corresponding ground stress at the observation positions using finite element method. The results (boundary loads and the corresponding ground stress) are added to the samples for ANN training. And on this basis, an ANN model is established to implement higher precise back-analysis of initial ground stress. A practical application case shows that the relative error between the inversed ground stress and observed value is mostly less than 10 %, which can meet the need of engineering design and construction requirements.