平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体元法的增量平衡方程。     

空间刚架结构大变形弹塑性屈曲及后屈曲分析

本文在当前Lagrange列式法的基础上,采用第二类P—K增量应力与Grecn增量应变之间的切线本构关系,并结合非线性杂交应力梁元进行空间刚架结构大变形弹塑性屈曲及后屈曲分析、计算结果表明了本文方法的有效性。     

TL法和UL法对刚架弹性大位移分析的适用性

本文分别基于Total Lagrange(TL)和Updated Lagrange(UL)描述分别导出了适用于刚架问题几何非线性分析的有限元公式系统，并研制了相应的计算程序，文中从量化的角度讨论了上述两种列式方法在通常情况下的适用性，并研究了诸如有限单元划分、杆件中轴向力的大小以及采用通常积分方法或积分平均化方法时等多种因素对于有限单元TL法和UL法适用性方面的影响，本文得出了一些相对量化的结论，对于非线性有限元方法的实际工程应用具有一定的指导意义。     

一般壳体结构大位移和大应变弹塑性有限元分析与损伤力学的某些应用

离散的Kirchhoff理论三角形单元对一般壳体结构的分析是十分可靠的和有效的，本文将推广此种单元去分析壳体结构的大位移大应变弹塑性响应，校正的拉格朗日增量的Jaumann应力公式和塑性损伤相耦合的塑性流动本构关系已应用于列式中，在对带环向裂纹和环向槽受四点弯曲载荷的圆醉壳的分析中，显示了和实验结果符合一致，且说明了与塑性损伤相耦合的大应变弹塑性分析无论在宏观水平上或微观水平上都能显示比断裂力学方法更符合实际的结果。     

弹塑性自黏结颗粒聚合体碰撞破损的离散元法模拟研究

利用颗粒离散元方法,并结合接触力学理论,对自黏连、弹塑性颗粒聚合体碰撞破损的细观力学机理进行了模拟研究。在颗粒间塑性变形存在的条件下,研究了颗粒物碰撞损伤。将弹性颗粒体碰撞与弹塑性颗粒体碰撞结果作了比较。结果表明:其他同等条件下,与弹性颗粒聚合体相比,弹塑性颗粒聚合体的损伤模式多为衰变或解体,而非断裂,其原由可以归结为塑性变形引起的额外能量损耗;其他同等条件下,与弹性颗粒聚合体比较而言,弹塑性颗粒聚合体碰撞加载阶段平板撞击力波动振幅较小,加载时间较长,撞击力峰值较大,且峰值出现的时间较晚,而损伤率却较大,其原因可以归结为颗粒间塑性变形对动能损耗、局部结构变化、内部力传播的影响。     

Galerkin边界元法用于弹塑性分析的准高次元方法

提出了一种适用于Ｇａｌｅｒｋｉｎ边界元法的高次单元插值方法和半解析半数值的积分方法 准高次元法，建立了有关数值模型和将该方法应用于结构弹塑性分析的算法模型．有关算例结果表明，本文建议的方法是切实可行的．     

大位移弹塑性样条有限条分析

本文将弹性样条有限条分析推广到大位移弹塑性分析过程．开发了大位移样条有限条分析的Ｃ语言程序包．可用于任意边界条件及变厚度平板及柱壳结构的分析．与有限元法相比可节约大量计算机内存和机时，并可望推广到塑性加工领域中．     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

高速冲击问题的离散元法数值模拟

以能量等效为原则,构建了新的求解弹塑性轴对称问题的离散元模型。通过引入断裂准则,使模型不仅可以应用于连续介质问题,还可以应用于从连续介质向非连续介质的动态转化过程。最后,通过模拟钢板受冲击载荷产生层裂的过程,说明模型在高速冲击问题中的适用性及其有效性。     

岩土弹塑性有限元分析的隐式积分弹性刚度法

本文提出了用于岩土弹塑性有限元分析的隐式积分弹性刚度算法。该算法既具有隐式积分法精度好、效率高、无条件稳定等优点，也具有弹性刚度法中刚度矩阵正定、对称的特点，更重要的是它避免了传统切线刚度法在处理岩土非相关联塑性流动和屈服面“角点”所遇到的非对称性和奇异性计算问题。通过算例分析了该算法的精度、效率     

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

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

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

预应力混凝土平面杆系结构的有限元方法

建立了基于有限元方法的考虑材料和几何非线性的任意截面预应力混凝土平面杆系结构的数值分析模型,可用于模拟预应力混凝土大跨度梁、单向偏压细长柱等的非线性全过程结构响应。引入修正的Rodriguez截面模型确定截面切线刚度,其中混凝土的贡献通过截面边界顶点定义的梯形单元来实现;在此基础上利用传统的平面非线性杆单元导出了标准有限元公式。通过两个算例验证了该模型的可靠性和适用性。     

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

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

极地船舶冰区操纵性能的六自由度离散元分析

极地船舶操纵破冰性能是破冰船设计建造过程中需重点考虑的问题。为分析极地船舶操纵性能,本文发展了冰区船舶六自由度操纵破冰运动模型,采用扩展多面体离散元模拟海冰,舵桨操纵模型提供破冰力。开展平整冰区定常直航模拟计算,并与Lindqvist船体冰阻力经验公式展开对比验证;开展雪龙号敞水35°舵角操纵性仿真,并与实船试航结果进行对比。在此基础上,开展冰厚及舵角影响下的船体结构冰载荷及破冰轨迹计算,模拟得到操纵破冰航行中船体垂荡、横摇及纵摇运动时程;最后分析操纵破冰船体线载荷分布与直航破冰下的差异。     

Improved experimental tracking techniques for validating discrete element method simulations of tumbling mills

We report on further developments in the three-dimensional tracking of a particle deep within the tumbling ball charge of an experimental mill. The experimental X-ray program employing the use of bi-planar X-ray angiography now includes the tracking of a typical 6 mm bulk charge particle in three dimensions with a spatial resolution that is accurate to within 0.15 mm. The improved experimental tracking techniques presented were developed for the purpose of generating accurate three-dimensional particle trajectory data against which to validate a numerical method for the simulation of discrete media, namely the discrete element method (DEM). These improvements are complimented with techniques for comparing charge profiles between numerical DEM simulations and three-dimensional experimental trajectory data.     

不同重力场下颗粒冲击过程的离散元分析

开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟,并与地球重力下冲击的试验结果进行对比验证。在此基础上,进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明,在所有重力加速度下,冲击物的穿透深度d与冲击速度v₀的关系可以用Poncelet模型表达;d与冲击物下落的总高度H表现为d～Hⁿ的幂律关系,当H<10 m时,d与H的幂率标度为0.322,而H>10 m时,d与H的幂率标度下降到0.211。此外,穿透深度小于冲击物半径时,重力加速度对冲击物减速过程无影响。在所有的重力加速度下,当冲击速度大于5 m/s时,冲击物的持续碰撞时间t_c是恒定的,且与重力的-1/2次方呈正比。     

大规模颗粒系统的精确缩尺和粗粒化离散元方法

计算效率是制约工程尺度大规模颗粒系统离散元计算发展的重要因素,现有的粗粒化处理方法局限于特定应用并且缺少一般的理论依据。本文采用量纲分析方法,描述了在精确缩尺系统中各物理量应当满足的缩放定律;通过在粗粒化系统和原始系统的代表性体积单元之间建立质量、动量和能量的近似守恒关系,采用多尺度的描述方法得到了粗粒化系统与原始系统之间宏观和细观两种不同尺度的缩放关系,即双尺度粗粒化模型;精确缩尺系统中得到的缩放定律及离散元接触模型处理方法,完全适用于粗粒化系统中细观颗粒层面相关物理量的缩放,通过筒仓侧壁压力和休止角两个算例对精确缩尺模型在粗粒化系统中的有效性进行了验证。     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。