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.

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.

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.

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

为进一步研究弯管内外表面腐蚀缺陷力学-电化学(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计算,并采用两种网格映射方式和四点插值法进行流场与结构的数据传递,完成了旋翼单向稳态及单向瞬态流固耦合分析。然后采用双向弱流固耦合方法计算结构动响应及非定常流场,分析比较单、双向流固耦合计算方法下桨叶振动响应的差异。最后,基于单、双向流固耦合计算得到桨叶表面的时变气动压强,计算旋翼旋转噪声,分析比较单、双向流固耦合计算方法的旋转噪声大小及分布特征。结果表明,文中基于高精度桨叶模型及双向流固耦合方法计算出的桨叶动响应及噪声分布符合悬停旋转噪声规律,可为未来桨叶结构设计、强度校核及噪声预估提供参考。

多胞子弹冲击泡沫夹芯梁的动力学响应分析

为了研究均匀/梯度多胞子弹冲击泡沫夹芯梁的耦合响应过程和多胞子弹对夹芯梁的加载效果,对该冲击过程开展了理论分析、数值模拟和试验研究:通过将泡沫夹芯梁等效为单梁以简化分析,基于多胞子弹的冲击波模型和泡沫夹芯梁的等效单梁响应模型,构建了多胞子弹冲击泡沫夹芯梁的耦合分析模型,给出了冲击过程中各响应阶段的控制方程,并结合龙格-库塔方法对方程进行了数值求解;基于三维Voronoi技术,开展了均匀/梯度多胞子弹冲击泡沫夹芯梁的细观有限元模拟;在多胞子弹的冲击测试平台上进行了试验研究,结合高速摄影技术获取了多胞子弹和泡沫夹芯梁的速度响应。结果表明:耦合分析模型可以准确地预测多胞子弹和泡沫夹芯梁的速度历程曲线以及多胞子弹产生的冲击压强;在初始动量相同但密度分布或初速度不同的多胞子弹冲击下,同一构型的泡沫夹芯梁展现出不同的力学响应,这说明多胞子弹的加载不能简单地等效为脉冲加载,多胞子弹与夹芯梁之间的耦合效应不可忽略;相较于均匀多胞子弹,梯度多胞子弹的冲击压力波形更加尖锐,在其衰减过程中展现出更强的非线性特征。

车辆服役环境对高速列车轴箱轴承温度的影响分析

构建考虑轴箱表面散热及轴承内部传热的功率损耗模型,分析轴箱轴承在不同服役环境下的载荷;综合考虑轴箱表面空气流场对流换热的影响,建立精细化轴箱轴承温度模型,分析不同服役环境对轴箱轴承温度分布和温度特性的影响,并通过轴承试验台验证模型的有效性。研究结果表明:当车辆速度由220 km/h增至300 km/h时,轴承的总摩擦力矩增大11.4%;当车轮多边形阶数由16阶增加到22阶时,摩擦力矩平均增大2.8%;轴箱轴承最高温度出现在内圈与滚动体接触的区域,最低温度出现在轴上且接近环境温度;当车速由220 km/h增加到300 km/h时,轴承的最高温度上升9.2℃,各节点处温度均有一定程度增加,当车轮多边形阶数由16阶增加到22阶时,最高温度平均升高1.1%;当多边形深度幅值由10 dB增加到18 dB时,最高温度平均升高1.4%。

不同起爆条件下双钢壁LNG储罐抗爆性能分析

由于双钢壁储罐具有建造周期短、投资小等优点,在LNG储存中应用越来越广泛,但其易在外罐受到外部冲击时引发内罐的破坏。为了分析储罐在不同外爆载荷条件下的破坏行为,以双钢壁LNG储罐为研究对象,考虑内、外罐及保冷材料之间的载荷传递与相互作用,建立了外罐-保冷材料-内罐多介质的多相耦合有限元模型,并对多影响因素下储罐在外爆载荷作用下的破坏进行了分析,探究其抗爆性能。结果表明:在外爆载荷作用下,内、外罐从迎爆面中心向内凹陷直至发生破坏,最大应力均在迎爆面中心,且越靠近中心应力值越大;内、外罐材料均先达到屈服强度后发生塑性变形,随后达到失效应变值并发生破坏;与外罐相比,内罐发生破坏的时间晚42.9%,变形范围和破坏范围分别小12.3%、6.7%;储罐的抗爆性能随着储存液位的升高和罐体容积的增大而增强,随着起爆位置从罐壁中心向上、下两侧移动而提高,在罐壁中心起爆时储罐的破坏程度最大。该研究成果可为双钢壁LNG储罐的设计和储存液位的控制等方面提供理论参考依据。

基于刚柔耦合模型的某车载高炮火控诸元修正方法研究

针对车载高炮射击时车体姿态对射击精度的影响问题,提出了一种基于刚柔耦合模型修正火炮射击诸元的方法。以某车载高炮为研究对象,通过建立其刚柔耦合模型对车体姿态进行预测,火控系统将获取的射击方向上的车体姿态响应实时分解到沿车头方向和垂直于车头方向,经过纵、横倾及车向反变换后由载车不稳定坐标系转换为载车地理系下,再经滤波与解相遇问题求解得到地理系下的目标诸元,随后将地理系下的目标诸元通过纵、横倾及车向变换后得到载车不稳定系下的火炮射击诸元,从而将修正后的瞄准诸元提供给火炮。为验证修正方法的可行性与适用性,以另一型车载炮为对象进行了试验验证,试验结果表明:经诸元修正后的火炮射击精度得到明显改善,火炮射击初期高低和方位弹丸脱靶量均值分别减小了53.1%和41.3%,修正效果显著,研究可为车载高炮火控射击诸元修正提供一种新思路。