考虑道砟侵切路基的有砟铁路离散元-有限元耦合方法

本文考虑有砟铁路道床结构中底层道砟颗粒会逐渐侵切路基上表层的实际情况,引入基于双曲线函数的土阻力公式,发展了离散元颗粒与有限元网格体内耦合的接触力算法,建立了考虑道砟颗粒侵切路基表层的有砟铁路离散元-有限元耦合模型。其中道床中的非规则形状道砟采用离散元镶嵌单元构造,路基采用有限元20节点六面体实体单元构造。分析了道砟侵切路基下列车荷载及路基软硬程度对道床整体沉降的影响,绘制了有砟道床-路基的横截面应力及位移云图。研究表明提出的耦合模型不仅能反映道床自身的弹性变形和累积塑性变形,还能真实反映底层道砟侵切软土路基表层所导致的侵切沉降,为深入了解高速重载下有砟道床的沉降劣化机理提供重要借鉴作用。     

针对有砟道床动力特性分析的嵌入式离散元-有限元耦合方法

在离散元-有限元耦合方法中,离散元和有限元交界面处的耦合方式对整体有砟道床的力学行为影响显著.采用基于球形单元的镶嵌单元或粘结单元模拟有砟道床时,由于球形单元和有限单元表面的自锁能力较差,使道砟层在列车载荷作用下容易产生侧向滑移,导致数值模型不稳定.此外,在实际铁路道床中,底部道砟均不同程度地嵌入路堤.为此,发展了一种嵌入式离散元-有限元耦合方法,通过设置一层嵌入地基有限元模型中的球形颗粒传递离散元域和有限元域间的力学参数,实现离散元和有限元方法的耦合.数值结果表明,嵌入式离散元-有限元耦合模型能够有效降低有砟道床的侧向位移,数值结果更加稳定,在处理与有砟道床类似的连续介质与散体介质的耦合问题时推荐采用嵌入式耦合算法.     

考虑道砟侵切土质路基的有砟轨道离散元-有限元耦合方法

论文考虑有砟轨道道床结构中底层道砟颗粒会逐渐侵切土质路基上表层的实际情况,引入基于双曲线函数的土阻力公式,发展了离散元颗粒与有限元网格体内耦合的接触力算法,建立了考虑道砟颗粒侵切路基表层的有砟轨道离散元-有限元耦合模型.其中道床中的非规则形状道砟采用离散元镶嵌单元构造,路基采用有限元20节点六面体实体单元构造.分析了道...     

高速铁路碎石道砟振动的离散元模拟

道砟振动对其磨损、破碎和道床累积变形有显著影响,为揭示高速车辆移动荷载作用下道砟动态响应特性,建立有砟道床离散元模型,开展车辆-轨道耦合动力学计算得到离散元模型输入荷载,模拟分析高速车辆以不同速度通过时有砟道床的振动响应,并与车辆-轨道耦合动力学计算结果进行对比分析。结果表明,轨枕、道砟和道床块振动位移波形相似,位移幅值沿道床深度方向减小,道床块振动位移与轨枕底面以下0.3m处道砟的振动位移相当;轨枕、道床块振动速度与加速度随行车速度提高而增大;受道砟颗粒间复杂相互作用的影响,道砟振动加速度会出现突变。道床离散元模型能合理反映道砟颗粒的振动响应特性,道床块模型体现了道床层在有砟轨道结构中的动力传递与减振特性,两种道床模型的计算结果具有一定的相似性。     

不同载荷频率下有砟铁路道床累积沉降特性的离散元分析

为研究有砟道床在高频载荷作用下的累积沉降规律,建立工程尺度下枕木-有砟道床的离散元计算模型。考虑道砟的非规则真实几何形态,利用不可破碎的簇颗粒模型近似构造。在10 Hz～100 Hz载荷频率区间内,对比分析枕木位移以及道砟颗粒在垂向和横向上的运动分布规律,研究分析道床在循环正弦载荷作用下的动力特性。数值结果表明,在载荷频率不超过60 Hz时,道床累积沉降随着载荷频率的提高而逐渐增大;当载荷频率超过60 Hz时,道床累积沉降随着载荷频率的提高而急剧增大,并在载荷频率为80 Hz时达到最大值,而后道床累积沉降量随着频率的提高反而逐渐减小。此外,分析对比不同载荷频率下道床速度场分布及典型位置道砟的加速度等物理量,发现80 Hz的载荷频率下道床离散元模型的振动最强烈。     

土工格栅加强下有砟铁路道床动力特性的离散元-有限元耦合分析

本文提出了一种可用于分析土工格栅与道砟颗粒间相互作用的离散元-有限元耦合方法。道砟颗粒采用镶嵌单元模拟,土工格栅采用梁单元计算,并在接触面上实现了离散元和有限元域间力学参数的传递。通过循环载荷作用下的计算,对比分析了土工格栅对道砟颗粒的加固性能。通过加筋后道砟颗粒的直剪试验进一步揭示了土工格栅加固作用的内在机理。结果表明：土工格栅可有效地降低有砟道床横向位移以及顶部沉降量,加筋后道砟材料的内聚力和摩擦角都有所提高。由此可见,土工格栅能够增加道砟颗粒间的自锁作用,提高道砟材料的强度,并有效地阻止剪切带从上部道砟箱向下部道砟箱贯穿。     

振动效应对铁路道砟剪切性能影响的离散元数值模拟

铁路道砟在传递移动列车载荷时会产生明显振动及接触迁徙行为,从而诱发道床的破碎劣化及非均匀沉降变形,道砟材料的抗剪性能研究对于碎石道床动力学行为及维养评估非常重要。为进一步探究铁路道砟的剪切力学特性,基于Minkowski Sum理论构造了扩展多面体单元,采用三维Voronoi切割算法生成了非规则扩展多面体道砟颗粒。通过最小投影方法对单元进行粒径识别和筛选,建立了符合铁路道砟级配要求的扩展多面体道砟模型数据库。对比了不同法向应力下铁路道砟准静态直剪试验及离散元模拟结果,验证了扩展多面体道砟接触参数的合理性。在此基础上,开展了振动效应下的道砟动态直剪离散元模拟研究,分析了振动角度、幅值及频率对道砟材料抗剪性能的影响。结果表明：构造的扩展多面体道砟离散元模型能够较好地反映非规则道砟间的强咬合互锁接触作用,并有效模拟直剪过程中的剪切应力–应变的变化及剪缩–剪胀行为。正则化下不同法向应力对应的铁道道砟剪切强度更适合采用幂函数拟合表示。振动场的存在会明显降低道砟颗粒的抗剪性能。当振动方向与剪切方向一致时（振动角度为0°）,道砟的剪切强度被明显削弱;随着振幅及频率的增大,道砟颗粒的剪切强度也会显著降...     

无砟轨道约束对高速铁路列车-桥梁系统地震响应的影响

建立了38自由度高速列车（动车/拖车）空间振动模型,基于非线性Hertz接触理论和Kalker线性蠕滑理论,引进PEER-NGA强震记录,建立地震作用下的高速铁路列车-无砟轨道/钢轨-桥梁-支座-桥墩-桩土系统精细计算模型,考虑支座刚度和桩土作用及近场地震效应的影响,以高速铁路32 m跨简支箱梁桥和CRTSⅡ板式无砟轨道为研究对象,编制了MATLAB程序,计算了有/无砟轨道约束下高速列车-桥梁时系统的地震响应。计算结果表明,不考虑三层无砟轨道时列车-桥梁系统模型由于桥梁刚度过大,列车-无砟轨道-桥梁系统的各项地震响应较列车-钢轨-桥梁系统有较大增加;对于有/无地震作用时,采用无砟轨道均可有效改进列车走行性能,近断层地震脉冲效应对有/无砟轨道约束系统动力响应及列车走行安全指标影响均较大。     

铁路道砟破碎特性的离散元分析

为降低有砟铁路轨道维护成本并提高乘客舒适度,对铁路道砟材料破碎机理深入研究具有重要作用。本文设计并开展了测量铁路道砟材料特征强度的试验,对道砟受压破碎过程的力和位移等数据加以分析;同时,通过三维激光扫描系统获取道砟三角网格多面体边界,依此采用颗粒平行粘结方式构造具有不规则真实外观的铁路道砟,并采用离散单元方法(DEM)数值模拟道砟径向加载破碎过程,分析应力应变曲线和粘结断裂数等数据以分析颗粒破碎机制。通过对多个道砟颗粒加载试验及数值模拟计算发现,随着道砟尺寸的增加,其有效压缩强度逐渐减小。道砟颗粒的有效压缩强度符合威布尔概率分布函数,并确定了其残存概率的平均强度和威布尔参量,试验和数值模拟结果相吻合。     

关于采用粗粒化提高颗粒材料多尺度模拟守恒特性的研究

连续体-颗粒耦合方法常用来描述连续-非连续颗粒行为或解决颗粒材料与其他可变形构件间相互作用问题。粗粒化coarse-graining (CG)是基于统计力学的均匀化方法,由离散的颗粒运动定义连续的宏观物理场。本文利用粗粒化(CG)推导有限元-离散元(FEM-DEM)表面和体积耦合的一般性表达式。对于表面耦合,CG可以将耦合力分布到颗粒-单元接触点以外的位置,如相邻的积分点;对于体积耦合,CG可以将颗粒尺度的运动均匀化到耦合单元上。由粗粒化推导出的耦合项仅包含一个参数,即粗粒化宽度,为均匀化后的宏观场定义了一个可调整的空间尺度。当粗粒化宽度为零时,表面和体积耦合表达式简化为常规局部耦合。本文通过弹性立方体冲击颗粒床和离散-连续介质间波传播两个数值算例,展示使用粗粒化方法提高耦合系统能量守恒的优势,并结合其他耦合参数(如体积耦合深度)讨论了粗粒化参数对数值稳定性和计算效率的影响。     

基于离散元-有限元模型的冰激锥体海洋平台结构振动分析

冰荷载是影响海洋平台结构及其上部设备作业安全的重要因素。针对海洋平台与海冰相互作用的动力过程,采用具有粘接-破碎性能的离散元方法(DEM)对海冰的破碎特性进行分析,采用有限元方法(FEM)对海洋平台结构的动力响应进行计算;同时考虑海冰与海洋平台结构的耦合作用,将海冰的离散元方法与海洋平台结构的有限元方法相结合,建立了冰激海洋平台结构振动的DEM-FEM模型,并由此计算海冰作用下海洋平台结构的振动响应以及冰荷载特性。     

离散元与壳体有限元结合的多尺度方法及其应用

在深入研究复杂结构和非均质材料冲击响应和破坏机理的过程中,往往遇到多尺度计算问题。本文尝试建立三维离散元与壳体有限元结合的多尺度方法用于处理圆柱壳问题,该方法采用三维离散元对感兴趣的局域进行局部模拟,利用平板壳体有限元进行整体模拟,采用一种特殊的过渡层使离散元区和有限元区能很好的衔接。我们将这一方法应用于激光辐照下充压柱壳的热/力耦合冲击破坏响应,得到的模拟结果与文献报道有较好的吻合。     

基于DEM-FEM耦合方法的海上风机结构冰激振动分析

冰载荷是海上风机在寒区安全运行的重要影响因素之一,由其引发的冰激振动给风机结构带来了严重的危害. 本文通过离散元(discrete element method, DEM)-有限元(finiteelement method, FEM)耦合方法建立了寒区单桩式风机结构的 冰激振动模型.采用具有粘结-破碎性能的球体离散单元描述平整海冰损伤破坏行为,采用梁单元和三角形平板壳单元构造带有抗冰锥体的单桩式风机有限元模型.采用DEM-FEM耦合方法模拟不同冰速、冰厚条件下单桩式风机与平整冰相互作用过程,并且与IEC规范和ISO标准经验公式对比验证该耦合模型计算冰载荷的准确性.对比风机塔筒顶端和基础顶端的位移和加速度响 应时程,定性地给出风机结构不同部位振动响应行为差异性.风机不同部位动力特性差异原因为风机结构独特结构特点:下部为大刚度桩基和上部为高柔度塔筒,使其动力特征表现为主从式结构特性. “主-从式结构”特征使得结构在复杂的冰载荷作用下,风机塔筒(子结构)和桩基(主结构)表现为不同的响应行为,风机不同部位振动周期和加速度谱两者出现差异. 本文研究成果为海上风机抗冰设计和疲劳分析提供了有益参考.      

基于离散元法的刚性履带低含水软地面附着特性数值分析

针对低含水软地面低黏性、高摩擦的力学特性,以刚性履带在干燥壤土行驶这一典型工况为试验条件,基于离散单元法选取合理的颗粒细观参数,在三轴剪切试验的基础上建立并验证了土壤数值仿真模型。通过数值牵引试验,对刚性履带低含水软地面的附着特性进行了离散元细观分析。分析表明,土壤附着力随履刺的向后排列而递减;垂直于履带行驶方向的附着作用不可忽略;离散元法能更好地还原低含水软地面的力学特性,揭示了土壤颗粒在动态荷载下非连续的本质,适合于干燥、松散的低含水软地面数值研究。     

Interface damage and its effect on vibrations of slab track under temperature and vehicle dynamic loads

This paper presents a three-dimensional finite element model to investigate the interface damage occurred between prefabricated slab and CA (cement asphalt) mortar layer in the China Railway Track System (CRTS-II) slab track system. In the finite element model, a cohesive zone model with a non-linear constitutive law is introduced and utilized to model the damage, cracking and delamination at the interface. Combining with the temperature field database obtained from the three-dimensional transient heat transfer analysis, the interface damage evolution as a result of temperature change is analyzed. A three-dimensional coupled dynamic model of a vehicle and the slab track is then established to calculate the varying rail-supporting forces which are utilized as the inputs to the finite element model. The non-linearities of the wheel–rail contact geometry, the wheel–rail normal contact force and the wheel–rail tangential creep force are taken into account in the model. Setting the maximum interface damaged state calculated under temperature change as the initial condition, the interface damage evolution and its influence on the dynamic response of the slab track are investigated under the joint action of the temperature change and vehicle dynamic load. The analysis indicates that the proposed model is capable of predicting the initiation and propagation of cracks at the interface. The prefabricated slab presents lateral warping, resulting in severe interface damage on both the sides of the slab track along the longitudinal direction during temperature drop process, while the interface damage level does not change significantly under vehicle dynamic loads. The interface damage has great effects on the dynamic responses of the slab track.