Numerical modeling for the coupled thermo-mechanical processes and spalling phenomena in sp Pillar Stability Experiment (APSE)

In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.

A rigid-flexible coupling finite element model of coupler for analyzing train instability behavior during collision

Rail vehicles generate huge longitudinal impact loads in collisions.If unreasonable matching exists between the compressive strength of the intermediate coupler and the structural strength of the car body,the risk of car body structure damage and train derailment will increase.Herein,a four-stage rigid-flexible coupling finite element model of the coupler is established considering the coupler buckling load.The influence of the coupler buckling load on the train longitudinal-vertical-hori-zontal buckling behavior was studied,and the mechanism of the train horizontal buckling instability in train collisions was revealed.Analysis results show that an intermediate coupler should be designed to ensure that the actual buckling load is less than the compressive load when the car body structure begins to deform plastically.The actual buckling load of the coupler and the asymmetry of the structural strength of the car body in the lateral direction are two important influencing factors for the lateral buckling of a train collision.If the strength of the two sides of the car body structure in the lateral direction is asymmetrical,the deformation on the weaker side will be larger,and the end of the car body will begin to deflect under the action of the coupler force,which in turn causes the train to undergo sawtooth buckling.

Influence of radial forging process on strain inhomogeneity of hollow gear shaft using finite element method and orthogonal design

Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.

Two-dimensional simulation on the rolling process of semi-solid 60Si2Mn by finite element method

The effect of various process variables on the law of metal flow for semi-solid rolling 60Si2Mn was studied by finite element method. Semi-solid 60Si2Mn can be described as compressible rigid visco-plastic porous material saturated with liquid. In terms of ther-mo-mechanical coupling condition, the distributions of stress, velocity and temperature were studied using software MARC. The simulation results show that the rigid visco-plastic model can accurately describe the semi-solid 60Si2Mn rolling process. The great deformation can achieve completely in view of low flow stress of semi-solid slurry.

Finite element simulation of chip formation in peripheral milling mild carbon steel

A two-dimensional coupled thermo-mechanical model is used to simulate the progress of milling mild carbon steel with continuous chip formation. Deformation of the workpiece is treated as elastic-plastic with isotropic strain-hardening. An analysis of the properties of temperature-dependent materials is employed. The chip separation is achieved by the adaptive remeshing strategy module found in commercial finite element code Marc. In order to increase the computational efficiency and accuracy, we implement a customized subroutine into the Marc code to define local refinement at the tool tip. Measurements are taken of the shape of the chip, the temperature, stress, strain and strain-rate fields. The values of the cutting forces obtained from the simulations agree well with those obtained from the cutting experiments.

Wave-Tide-Surge Coupled Simulation for Typhoon Maemi

The main task of this study focuses on studying the effect of wave-current interaction on currents, storm surge and wind wave as well as effects of current induced wave refraction and current on waves by using numerical models which consider the bottom boundary layer and sea surface roughness parameter for shallow and smooth bed area around Korean Peninsula. The coupled system （unstructured-mesh SWAN wave and ADCIRC） run on the same unstructured mesh. This identical and homogeneous mesh allows the physics of wave-circulation interactions to be correctly resolved in both models. The unstructured mesh can be applied to a large domain allowing all energy from deep to shallow waters to be seamlessly followed. There is no nesting or overlapping of structured wave meshes, and no interpolation is required. In response to typhoon Maemi （2003）, all model components were validated independently, and shown to provide a faithful representation of the system＇s response to this storm. The waves and storm surge were allowed to develop on the continental shelf and interact with the complex nearshore environment. The resulting modeling system can be used extensively for prediction of the typhoon surge. The result show that it is important to incorporate the wave-current interaction effect into coastal area in the wave-tide-surge coupled model. At the same time, it should consider effects of depth-induced wave breaking, wind field, currents and sea surface elevation in prediction of waves. Specially, we found that： （1） wave radiation stress enhanced the current and surge elevation otherwise wave enhanced nonlinear bottom boundary layer decreased that, （2） wind wave was significantly controlled by sea surface roughness thus we cautiously took the experimental expression. The resulting modeling system can be used for hindcasting （prediction） the wave-tide-surge coupled environments at complex coastline, shallow water and fine sediment area like areas around Korean Peninsula.

MECHANISM OF WATER-SOIL COUPLED ACTION DURING MINING SUBSIDENCE

This paper,on the basis of the scientific research of engineering geological exploration in a mining area,systematically studies the reasons and influence factors of consolidation and deformation of the saturated soil included in the thick loose water-bearing overburden due to mining subsidence,and analyses the dissipation of hyperstatic pore water pressure during the change of original stress and strain state of the soil. Again,by means of the coupled model based on Cambridge model and Biot's three-dimensional consolidation theory,adopting a great many physico-mechanical parameters measured in various soil layers,the paper analyses the consolidation and deformation of saturated soil affected by mining subsidence with elasto-plastic finite element method.Thus,the research not only reveals the regulation of stress,strain,displacement and hyperstatic pore water pressure dissipation in overlying soil mass,but also opens up a new direction and way for the research of mining subsidence.

Thermal-hydro-mechanical coupled analysis of unsaturated frostsusceptible soils

Damage caused by frost heave leads to costly maintenance in cold regions, like Hokkaido, Japan. Therefore, thestudy of the frost mechanism with experimental and numerical methods has been of great interest. Numerousmodels have been developed to describe the freezing process of saturated soil, which differs from the partiallysaturated conditions in the field. In fact, most subsurface soils are unsaturated. The freezing process of partiallysaturated soils is more complex than saturated soils, as the governing equations show strongly nonlinear characteristics. This study proposes a thermo-hydro-mechanical coupled model considering the heat transfer, waterinfiltration, and deformation of partially saturated soil to reproduce the freezing process of partially saturatedfrost susceptible soils distributed in Hokkaido. This model better considers the water-ice phase change and soilfreezing characteristic curve (SFCC) during freezing under field conditions. The results from the multiphysicssimulations agree well with the frost heave and water migration data from frost heave tests of Touryo soil andFujinomori soil. In addition, this study discussed the influence of the various factors on frost heave amount,including temperature gradients, overburden pressures, water supply conditions, cooling rates, and initial saturation. The simulation results indicate that the frost heave ratio is proportional to the initial degree of saturationand is inversely proportional to the cooling rate and overburden pressure.Moreover, simulation under the open system generates much more frost heave than under the closed system.Finally, the main features of the proposed model are revealed by simulating a closed-system frost heave test. Thesimulation results indicate that the proposed model adequately captures the coupling characteristics of water andice redistribution, temperature development, hydraulic conductivity, and suction in the freezing process. Togetherwith the decreased hydraulic conductivity, the increased suction controls the water flow in the freezing zone. Theinflow water driven by cryogenic suction gradient feeds the ice formation, leads to a rapid increase in total watercontent, expanding the voids that exceed the initial porosity and contributing to the frost heave.

Finite element implementation of coupled hydro-mechanical modeling of transversely isotropic porous media in DEAL.II

In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on homogenous and isotropic materials.However,soil and rock are often modeled in heterogeneous and anisotropic media because of the inherent anisotropy of sedimentary rock masses and their stratified structure.In recent decades,coupled hydro-mechanical(HM)interactions in isotropic porous media have been studied;however,the behavior of transversely isotropic porous media is rarely considered.In addition,it is difficult for commercial software such as Plaxis and Flac3D to express complex rock formation where the anisotropy of the material and the associated cracks and fractures could be assembled into a single model.In this study,a finite element implementation using Differential Equation Analysis Library(DEAL.II),an open-source library of finite element codes,was developed to model the fully coupled HM behavior of transversely isotropic porous media.The proposed implementation can be applied to both isotropic and transversely isotropic porous media based on Biot’s theory.The developed code can be used to model poroelastic media with(1)equations of linear elasticity for the solid matrix and(2)diffusion equations for fluid flow based on mass and linear-momentum conservation laws.We verified the performance and accuracy of the code through two examples,i.e.,Mandel’s problem with a compared analytical solution and a tunnel excavation process with the Flac3D software.On the basis of these numerical applications,we present the code to model the behavior of various geo-structures such as tunnels and pile–soil interactions with anisotropic materials.

A Coupling Model of the Discontinuous Deformation Analysis Method and the Finite Element Method

Neither the finite element method nor the discontinuous deformation analysis method can solve problems very well in rock mechanics and engineering due to their extreme complexities. A coupling method combining both of them should have wider applicability. Such a model coupling the discontinuous deforma- tion analysis method and the finite element method is proposed in this paper. In the model, so-called line blocks are introduced to deal with the interaction via the common interfacial boundary of the discontinuous deformation analysis domain with the finite element domain. The interfacial conditions during the incre- mental iteration process are satisfied by means of the line blocks. The requirement of gradual small dis- placements in each incremental step of this coupling method is met through a displacement control proce- dure. The model is simple in concept and is easy in numerical implementation. A numerical example is given. The displacement obtained by the coupling method agrees well with those obtained by the finite ele- ment method, which shows the rationality of this model and the validity of the implementation scheme.

Simulation of Coupling Process of Flexible Needle Insertion into Soft Tissue Based on ABAQUS

In order to get to the desired target inside the body,it is essential to investigate the needle-tissue coupling process and calculate the tissue deformation.A cantilever beam model is presented to predicting the deflection and bending angle of flexible needle by analyzing the distribution of the force on needle shaft during the procedure of needle insertion into soft tissue.Furthermore,a finite element(FE)coupling model is proposed to simulate the needle-tissue interactive process.The plane and spatial models are created to relate the needle and tissue nodes.Combined with the cantilever beam model and the finite element needle-tissue coupling model,the simulation of needle-tissue interaction was carried out by the ABAQUS software.The comparing experiments are designed to understand the needle-tissue interactions,by which the same points in the experiments and simulation are compared and analyzed.The results show that the displacements in x and z directions in the simulation can accord with the experiments,and the deformation inside the tissue mainly occurs in the axial direction.The study is beneficial to the robot-assisted and virtual needle insertion procedure,and to help the physicians to predict the inside tissue deformation during the treatments.

Multi-physics coupling field finite element analysis on giant magnetostrictive materials smart component

This study presents a new method to solve the difficult problem of precise machining a non-cylinder pinhole of a piston using embedded giant magnetostrictive material (GMM) in the component. We propose the finite element model of GMM smart component in electric, magnetic, and mechanical fields by step computation to optimize the design of GMM smart com-ponent. The proposed model is implemented by using COMSOL multi-physics V3.2a. The effects of the smart component on the deformation and the system resonance frequencies are studied. The results calculated by the model are in excellent agreement (relative errors are below 10%) with the experimental values.

Numerical analyses of coupled thermo-mechanical behaviors of cement sheath in geothermal wells

The stability of cement sheath under high temperature and high pressure is one of the most critical issues for the durability of geother-mal well systems.In this study,a two-dimensional plane-strain finite element code was developed to investigate the coupled thermo-mechanical behaviors of the casing-cement-formation system.Different from previous linear elastic analyses,a thermoelasto-plastic con-stitutive model based on the thermodynamic theory was adopted for the cement sheath.It is shown that the finite element simulations using the proposed model provide a more accurate and realistic prediction of stress–strain responses of the cement sheath under high temperature.The results demonstrate that the radial stress concentration and the tensile strain concentration occur at both the cement–casing interface and the cement–formation interface,where the cement sheath is most likely to fail.High strength and low stiff-ness in the cement sheath and the formation are preferred for the integrity of the system.Both large thermal cycles and large differences between the internal fluid pressure and the external pressure should be avoided during operation.The new code is an alternative tool for guiding the geothermal well design.The finite element framework described herein is universal for other thermo-mechanical applications,such as energy foundations and energy tunnels.

弯管内外表面腐蚀缺陷力学-电化学相互作用规律研究

为进一步研究弯管内外表面腐蚀缺陷力学-电化学(M-E)相互作用规律,建立基于多物理场耦合的弯管内外表面含腐蚀缺陷的三维模型,研究内压和缺陷几何形状对M-E相互作用规律及对弯管的影响。研究结果表明:当弯管内外表面存在腐蚀缺陷时,弯管内侧的M-E效应更强,增加应力和腐蚀的增长速率,使内侧更容易发生失效;当缺陷几何增长时,应力在缺陷长度和宽度方向上增长,腐蚀主要在缺陷长度方向上增长;外表面缺陷附近的应力和腐蚀增长规律基本一致,而内表面缺陷附近的应力增长较腐蚀增长更为复杂。研究结果可为含腐蚀缺陷弯管评估模型的完善以及管道系统完整性的管理提供理论参考。

不同制动工况下地铁车轮踏面磨耗研究

为研究不同制动方式对地铁车轮踏面磨耗的影响规律,联合ANSYS和SIMPACK软件,建立了车轮-闸瓦热机耦合有限元模型和含有柔性轮对的车辆动力学模型,同时考虑不同制动初速度和制动减速度对车轮踏面温度、踏面硬度、踏面与闸瓦间接触应力及制动距离的影响,分析了不同制动方式下车轮-闸瓦磨耗、车轮-钢轨磨耗和踏面总磨耗的变化规律。结果表明:踏面磨耗随着制动初速度的增大而增大;随着制动减速度的增大,车轮-闸瓦磨耗和总磨耗增大,车轮-钢轨磨耗减小;制动时踏面总磨耗主要受车轮-闸瓦磨耗影响。

不同季节影响下复合道面力学响应分析

为了研究复合道面在不同季节时的力学响应,以中原地区某机场改造项目为项目背景,借助三维有限元软件,建立了道面有限元模型,分析了复合道面在温度影响下的力学响应,并与机场监测数据进行对比,验证了模型的准确性。研究发现:大气温度所能影响的道面深度范围在0.56 m以内,和大气温度相比,道面面层温度峰值远大于气温峰值;道面内部不同结构层之间温度变化趋势相同,随着道面深度的加深,温度峰值逐渐递减,并且出现时间具有延迟现象;在不同季节,道面在相同机型作用下所产生的变形不同,说明温度对于道面的受力是有影响的,其中混凝土层在冬季时所产生的应变值比夏季所产生的减少约74%,变化比较明显。

大型风力机叶片前缘的雨滴冲击动力学响应机制

随着近年海上风力机设计的大型化,叶片前缘的雨蚀失效问题变得愈发突出,不仅影响机组的风能转化效率,对结构的稳定运行也构成潜在威胁。采用光滑粒子流体动力学(smooth particle hydrodynamics,SPH)法研究雨滴内部的本构关系,以有限元法(finite element method,FEM)建立叶片前缘的代表性体积单元(representative volume element,RVE)模型,两者相互耦合,研究雨滴在叶片表面形成的冲击响应过程。考虑自然降雨的实际工况,建立与降雨强度相关联的雨滴尺寸模型,以及雨滴的空间分布模型;通过单雨滴的冲击仿真,研究叶片表面冲击载荷以及雨滴内部的速度分布,从而解构雨滴冲击的物理过程,并通过冲击的应力、应变场分析,为潜在损伤区提供评价;建立多雨滴冲击的仿真模型,研究冲击应力场之间的耦合作用以及涂层表面的塑性应变累积效果。研究结果表明,水锤冲击是造成塑性应变累积的关键因素,虽然横向喷射阶段的应力响应幅值小并呈现一定的无序特征,但如果存在多雨滴耦合冲击的情况,则会在耦合区内出现应力峰值,并对叶片变形和失效存在潜在影响。

考虑气体的渗流应力耦合模型在黄壁庄水库中的应用

目前建立的多孔介质耦合模型通常认为气体是不可压缩的,无法描述气体的输运过程,无法反映气体对水体入渗和固体变形的影响。为了描述多孔介质中固液气三相的渗流应力耦合过程,基于多孔介质理论和有效应力原理,提出了一种耦合了改进守恒水平集(ICLS)方法并考虑气体可压缩性的固液气三相渗流应力耦合模型。将该方法应用于黄壁庄水库坝前冒泡现象的渗流应力耦合分析中,模拟了水位骤升后,固体变形、水体入渗和气体运移过程,及其耦合效应。计算结果表明,该方法能够简单直观的描述水气两相的运移过程及渗流应力耦合过程,且模拟结果与实际观测结果基本一致。

C型钢缀板柱的有限元模拟方法研究

缀板柱以其制造简便、稳定性高且价格低廉的优势,在工程钢结构中得到了广泛应用,如何准确模拟缀板与立柱之间的连接情况,是工程结构有限元模拟时的重要问题。文章采用有限元软件ANSYS,研究C型钢缀板柱中缀板与立柱之间连接情况的有限元模拟问题,通过对比单元整体建模与节点耦合建模的模拟结果,分析了节点耦合连接方法的模拟原理并提出了使用时的注意事项。计算分析表明:整体建模模拟效果较好但不便于模型修改和精细化模拟,节点耦合便于修改但需要考虑耦合方法和连接面积。文章提出的简化统一耦合方法模拟效果好,适应性强,对于工程实践中的有限元模拟工作具有一定指导意义。

基于双向流固耦合的旋翼振动响应及噪声仿真分析

为精确预估流固耦合效应对旋翼桨叶振动响应及噪声计算的影响,本文建立了计算流体力学(Computational fluid dynamics,CFD)模型和高精度结构有限元模型,基于弱流固耦合方法完成了旋翼单向及双向流固耦合计算,并采用基于Kirchhoff方法对FW⁃H方程进行求解,分析了旋翼近场旋转噪声。首先,基于地面模态分析试验对桨叶结构有限元模型进行高精度建模,同时完成旋翼旋转状态的CFD计算,并采用两种网格映射方式和四点插值法进行流场与结构的数据传递,完成了旋翼单向稳态及单向瞬态流固耦合分析。然后采用双向弱流固耦合方法计算结构动响应及非定常流场,分析比较单、双向流固耦合计算方法下桨叶振动响应的差异。最后,基于单、双向流固耦合计算得到桨叶表面的时变气动压强,计算旋翼旋转噪声,分析比较单、双向流固耦合计算方法的旋转噪声大小及分布特征。结果表明,文中基于高精度桨叶模型及双向流固耦合方法计算出的桨叶动响应及噪声分布符合悬停旋转噪声规律,可为未来桨叶结构设计、强度校核及噪声预估提供参考。