动力分析的刚体－界面元与有限元耦合法

提出了刚体－界面元、有限元与半无限元相耦合的算法，研究了挡水结构系统的动力特性及在地震激励下的动力响应。对于结构的主要部分用刚体－界面元离散，对结构的迎水面部分用有限元离散，在两种单元的交接面上设置了过渡单元，采用有限元和半无限元相结合模拟了水体的宽广性。文中导出了相关公式，编制了相应的计算程序，数值算例结果与有关文献比较，证实了提供方法和编制程序的正确性，并说明该方法具有简便，工作量少和计算精度高的优点。     

高炉内炉料流动模式及力链分布的离散元模拟

基于离散元数值计算方法,建立了高炉内炉料颗粒尺度运动行为的数学模型,主要研究固体炉料的运动模式和颗粒间相互作用力链的分布.结果表明:建立的离散元模型计算获得了炉内颗粒间的介观力链结构,炉底中心部位存在强力链结构支撑高炉料柱,最强力链结构对应于死料柱区,而且离散元模拟也给出炉内固体料运动模式由四个区域构成,分别为死料柱区、活塞流区、准静态滑流区和沟流区,而沟流区的力链最弱.     

基于ABAQUS的水平成层场地地震反应无限元分析

地下结构的抗震研究是当今地震工程界的重要研究方向之一,其研究方法与地面结构相比有较大的不同.本文基于动力无限元分析理论,利用ABAQUS软件的动力时程分析结合无限元边界实现了一维岩体的地震动响应分析,进而同基于SHAKE程序进行的一维土体地震动反应分析结果作了对比,说明在ABAQUS软件中实现的一维岩体动力反应分析方法是具有较好精度和易用性的.     

动力分析的刚体—界面元与有限元耦合法

提出了刚体－界面元、有限元与半无限元相耦合的算法，研究了挡水结构系统的动力特性及在地震激励下的动力响应，对于结构的主要部分用刚体－界面元离散，对结构的迎水面部分用有限元离散，在两种单元的交换面上设置了过渡单元，采用有限元和半元限元相结合模拟了水体的宽广性，文中导出了相关公式，编制了相应的计算程序，数值算例结果与有关文献比较，证实了提供方法和编程序的正确性，并说明该方法具有简便，工作量少和计算精度高     

非饱和土三维离散元计算中宏-细观参数关系探究

提出了一种建立非饱和土体宏-细观参数之间关系的方法,来建立在不同孔隙比和含水率等初始条件下非饱和土在不同应力路径下的离散元计算模型;通过编制离散元程序,基于接触粘结模型,对现有的PFC3D(Particle Flow Code in three dimensions)离散元程序进行改进,从而对在不同颗粒间粘结强度下的离散元试样进行一维固结的数值模拟试验来确定其结构屈服应力,并以结构屈服应力为桥梁建立颗粒间粘结强度随含水率变化的函数关系,最后建立能够反映真实非饱和土试样颗粒级配和在不同含水量下的离散元数值模型,为通过PFC3D等三维离散元软件研究非饱和土的基本力学特性提供思路.     

非饱和土力学行为的三维离散元分析

基于离散单元法颗粒流理论,土体颗粒单元间采用Hill接触模型来考虑非饱和土中吸力的作用,建立了非饱和土的颗粒流模型。根据室内试验结果对模型细观参数进行标定,开展了固结排水三轴剪切试验离散元数值模拟,揭示了在不同围压和吸力下,模型试样的应力-应变关系与强度特性等宏观力学性质以及细观信息平均配位数的变化规律。根据不同吸力下的莫尔-库伦破坏包线,建立了吸力与抗剪强度参数之间的关系。此外,通过标定室内试验结果中不同含水率下的应力-应变曲线,建立了Hill模型中吸力与含水率的关系。研究结果表明：Hill接触模型可以较好的模拟非饱和土的力学行为,可为非饱和土的离散元模拟提供参考。     

基于离散元的南海软粘土剪切变形模拟

软粘土剪切变形是决定海洋构筑物、管线安全存在的重要因素,离散元能有效模拟粘土剪切变形。通过对南海神狐海域原状软粘土开展固结不排水三轴试验,得到离散元模拟所需宏观剪切强度参数指标;采用试错法调整细观参数使得双轴压缩模拟结果与室内试验相吻合,标定出一组适用于南海软粘土的细观参数组合。离散元模拟结果得到,在应变达到5.3%时,土体发生破坏,形成一条与X轴反方向为49°的剪切带;应变增加到10%时,倾角增大到57°。测量圆对剪切带内外的应变进行监测得到,剪切带内土体初始应变较大,在达到4%应变后不再继续发生变形;剪切带边缘区域,土体应变持续增大;剪切带外应变较小。一般的室内三轴试验无法得到剪切带内外土体应变变化定量数据,离散元模拟是对室内三轴试验的重要补充。     

用无穷元模拟半无限平面弹性地基

本文提出一种七节点超参无穷元用于模拟半无限平面弹性地基的计算方法。文中阐述了它的基本原理与公式，对其中衰减函数与无穷无边中节点位置的选择进行了比较研究．对均布及非均布荷载作用下的半无限弹性地基给出了数值结果；并与精确解作了比较。最后用这类模型计算了重力坝的地基应力分布，表明了它具有较高的计算精度与显著的经济性。     

极区船舶及海洋结构冰荷载的离散元分析

 冰荷载是影响极区船舶及海洋结构振动响应和疲劳寿命的重要因素。针对海冰的离散特性,分别采用具有黏结-破碎效应的球体单元、扩展圆盘单元以及基于闵可夫斯基和方法的扩展多面体单元构造海冰离散单元模型,分析了船体、固定式、自升式以及浮式海洋结构的冰荷载特性及分布规律。此外,为计算海冰作用下海洋结构的振动响应,建立了相应的冰激结构振动的离散元-有限元模型。在上述基础上自主研发了基于图像处理单位（GPU）并行算法的海冰离散元软件（IceDEM）,实现了对船舶以及海洋平台结构冰荷载离散元分析的大规模计算。     

基于地理空间实体的矢量元胞自动机研究

元胞自动机是能模拟复杂行为的离散动力学模型，GIS与元胞自动机的集成对于模拟地理时空过程具有很大的优势。但是标准CA的一些限制因素如空间尺度、精度及数据量问题、邻居定义等制约了其模拟真实世界的能力。本文通过分析地理现象和元胞自动机之间的关系，从元胞及状态、元胞空间、邻居、转换规则、离散时间五个方面对标准元胞自动机的组成部分进行了扩展，构建了在不规则的地理空间中基于地理实体的元胞自动机模型，同时提出了一个通用的基于地理实体的元胞自动机模型和地理信息系统集成的框架。     

可变形体离散元模型

为解决岩土工程中的大变形和不连续问题，根据变形体动力学和离散元法的基本原理，建立了适用于可变形体的单元模型．采用弹塑性边-边接触模型及动态松弛法，推导出相应的理论公式并编制了计算程序．通过选取适当的计算参数对典型算例进行分析计算；得到了静力平面问题的应力、位移分布及静态解的收敛过程；验证了计算程序和计算参数选取的正确性；证明了该方法用于岩土工程静力计算的可行性．     

适用于可变形体的离散元法

由变形体动力学原理和模态分解法,推导出可变形块体单元模型,并根据离散元法的基本原理,选择合理接触模型,依照动态松弛法,编制相应的计算程序．由静力问题为例,验证了离散元应用于可变形体的可行性．     

离散元法土-地下结构动力相互作用分析

基于变形体动力学原理,建立了新的可变形块体单元模型.根据离散元法原理,采用边-边接触关系及动态松弛法,推导出其理论公式并编制了计算程序;由静力问题计算结果的收敛性,验证了计算程序和计算参数选取的正确性,求出了地下结构在人工地震波作用下的动力反应.     

节理岩体的动力响应分析及工程应用

节理岩体经常受到地震、爆破等动载荷的作用,在动载荷作用下的稳定性经常是结构安全的关键。只有对节理岩体的动力响应有本质的理解,才有可能对岩体工程的动力稳定性进行合理的评价。在常规离散元法的基础上提出了动力离散元法的思想及计算方法,并编制了计算程序DDEM,最后给出了该方法在岩体动力学中的应用实例。     

风暴状态下自升式平台非线性动力分析模型

以风暴状态下的自升式平台为研究对象 ,建立了用于分析平台结构的非线性动力模型。该模型采用改进的NewWave模型描述随机波浪 ,借助非线性弹簧模拟平台结构的桩—土之间的相互作用 ,综合考虑了影响平台结构动力响应的动力敏感性、非线性以及波浪载荷的随机性因素 ;同时以波面极值出现前后的若干载荷循环作用下的结构响应结果作为对平台结构响应的估计 ,大大节省了计算时间。建立了简化的三维有限元结构模型。通过具体实例分析 ,给出了使用该模型进行风暴状态下自升式平台动力响应的计算步骤 ,并对分别采用非线性弹簧模型和传统的铰接模型时的平台结构动力响应进行了对比。计算结果表明 ,风暴状态下自升式平台甲板位移响应极值服从对数正态分布。使用该模型进行自升式平台动力响应分析时 ,能够兼顾模型复杂程度和计算效率 ,这为今后进行自升式平台结构动力可靠性分析奠定了基础。     

大型发电机冷却器和机座的动态分析

首先概述了计算复杂系统模态参数的动态予结构法的基本原理。然后采用有限元计算与试验模态分析相结合的研究方法,分析了发电机主要部件的有机玻璃结构模型的动态特性,以确定发电机结构的合理的离散化数学模型。最后据此按照设计图纸研究了发电机结构原型的动态特性,并进而论证了发电机结构设计的可靠性。     

风力发电结构有限元建模及其地震响应分析

为研究不同精细程度的风力发电结构有限元模型在不同分析目的下的适用性,基于2MW三桨叶水平轴风力发电机,采用有限元软件ANSYS建立七种不同精细程度的风力发电结构有限元模型,以风力发电结构地震响应分析为例,得到各模型的响应结果,分别从变形、内力、应力以及应力集中四个方面分析各个模型的计算结果.研究结果表明:在计算风力发电结构的变形时,由于上部机舱及叶片对其影响较小,采用低阶单元即可得到较好的模拟结果;在计算塔底截面剪力和截面弯矩时,考虑机舱及叶片的有限元模型得到的结果较为精确,单元选择时需采用高阶实体单元;不考虑叶片、轮毂以及机舱得到的塔筒截面应力离散性较大,采用壳单元可较好地模拟塔筒截面应力;叶片、机舱等上部结构对门洞局部产生的应力集中影响较大,在计算时建议应用塔筒采用壳单元、门洞局部采用实体单元的多尺度模型进行模拟.     

基于DEM-FEM耦合的滑坡作用下输气管道动力响应分析

油气长输管道输送距离长,途径地形地貌复杂多样,沿途地质灾害频发对油气管道的安全输送造成严重威胁与破坏。本文以联络线输气管道横向滑坡为例,通过离散元与有限元单向耦合方式,对管道进行滑坡作用下的动力响应分析。首先通过EDEM离散元软件模拟滑坡过程中土体颗粒对管道的动态冲击过程,分析在滑坡过程中土体颗粒与管道的接触变化以及受力变化,然后将管道所受最大滑坡推力作为荷载导入ANSYS中,模拟得到在滑坡及内压协同作用下管道的应力和位移分布,最后将模拟结果与线上监测数据进行对比。结果表明：由离散元模拟得到滑坡作用下输气管道受力时程变化曲线,在第3秒管道受到的滑坡推力最大,将其导入有限元中模拟得到滑坡段管道等效应力、剪切应力与位移分布云图,将等效应力与实际监测数据进行对比,计算得到相对误差在7.30%以内,因此表明该模型具有一定的可靠性。     

筒型基础负压沉贯中的动态物性参数

为研究筒型基础内海积土的物性参数，结合筒内沉贯土层性质及沉贯环境，应用J2流动理论和各向同性硬化原理计算沉贯过程中土颗粒的应变，考虑到土中出现贯通裂隙，应用离散元方法计算土体骨架应变,结合土颗粒和土体骨架应变，从土的孔隙率和渗透率的定义出发，推导出沉贯过程中筒内海积土物性参数的动态模型-最后，通过实例计算验证了模型和计算过程的正确性．     

渐变刚度钢板弹簧建模方法对比

为更好地分析钢板弹簧刚度特性,对钢板弹簧建模方法进行了对比研究。依据钢板弹簧国家标准GB/T 19844—2005,针对某型渐变刚度钢板弹簧分别在有限元软件Hypermesh和多体动力学软件Adams中建立其有限元模型和离散梁模型,并进行了仿真,得到刚度特性。将仿真结果与试验值进行了比较。结果表明,在针对该钢板弹簧的刚度分析中,运用离散梁法得到的刚度特性更加接近试验值。