裂缝介质渗流多尺度混合有限元数值模拟研究

针对裂缝介质具有多尺度特点,建立了Darcy/Stokes-Brinkman多尺度耦合模型,采用多尺度混合有限元方法,对裂缝介质渗流问题进行了研究.阐述了多尺度混合有限元方法的基本原理,并推导得到Darcy/Stokes-Brinkman方程的多尺度混合有限元计算格式.数值计算结果表明,大尺度Darcy模型能够捕捉到小尺度上裂缝网络渗流特征;与网格粗化、传统有限元方法相比,多尺度混合有限元方法的基函数具有能反映单元内参数变化的优点,在保证计算精度的同时能够减少计算量,对于裂缝油藏具有良好的适用性.     

可压缩流体流动多尺度混合有限元数值方法研究

可压缩流体是天然油藏中广泛存在的一种流体,研究其在多孔介质中的渗流规律对于油藏开发具有重要意义。本文采用多尺度混合有限元方法,对可压缩流体渗流问题进行了研究。考虑流体的可压缩性以及介质形变,推导得到了可压缩流体渗流问题的多尺度计算格式。数值计算结果表明,多尺度混合有限元适于求解非均质性和可压缩流问题,具有节省计算量、计算精度高等优势,对于实际大规模油藏模拟具有重要意义。     

微纳米尺度下材料性能多尺度模拟方法进展

微纳米材料的力学行为正日益引起研究者的关注.微纳米材料的性能取决于从微观、细观到宏观多个空间、时间尺度上不同物理过程非线性耦合演化的结果,发展相应的多尺度数值模拟方法已成为该领域研究工作的一个热点.本文对微纳米材料模拟中比较典型的几种协同多空间尺度和协同多时间尺度方法进行了介绍,着重介绍这些方法的的基本思想、应用情况, 以及各自的优缺点,并对微纳米材料多尺度方法的发展趋势进行总结和评述.      

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

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

多尺度模拟与计算研究进展

简要介绍了多尺问题与研究方法.重点论述了两类常见多尺度问题的模拟计算方法与研究进展,分析了各自的优缺点和使用范围.对现有研究的局限性和存在的问题进行分析,指出了进一步研究多尺度模拟与计算的必要性.介绍了求解含有孤立缺陷问题的非局部准连续体法,MAAD等方法以及求解基于微观模型本构模拟问题的局部连续体法、HMM等方法.文章对多尺度模拟与计算的前景进行展望,提出了一些亟待解决的问题.     

微/纳尺度粘着滑动接触的多尺度分析

采用分子动力学与有限元耦合的多尺度方法,求解二维刚性圆柱表面压头与弹性平面的微/纳尺度粘着滑动接触问题,通过与全分子动力学模拟结果的比较验证了多尺度方法的有效性。对压头半径、滑动速度、下压深度以及是否考虑粘着效应等对滑动接触性能的影响进行了全面研究,通过不同条件下摩擦力及接触力分布的比较,揭示了上述各参数对粘着滑动接触...     

大型土木结构的结构行为一致多尺度模拟——模拟方法与策略

针对大型土木结构以结构健康监测和损伤评估为目标的有限元分析的需要,从理论上探讨了复杂土木结构的结构行为一致多尺度模拟的方法和策略,指出对于大型复杂土木结构而言,适合考虑局部细节特性的结构行为多尺度建模的策略是:对于有限个需要重点关注的局部细节采用"小尺度"建模,而结构其他部分采用"大尺度"建模.针对目前大型复杂土木结构的建模都是借助于大型商用有限元软件,而现有软件都无法直接做到结构行为一致多尺度模拟,文中引入有限元方法中的子结构方法,从数学描述上对传统子结构方法加以改进,使之适合解决结构行为一致多尺度模拟问题.     

基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟

针对缝洞型油藏具有多尺度特征,基于离散缝洞网络模型,建立了Darcy/Stokes-Brinkman多尺度耦合数学模型,采用多尺度混合有限元方法,对缝洞型油藏的流体流动问题进行了研究。阐述了多尺度混合有限元方法的基本原理,推导得到了Darcy/Stokes-Brinkman多尺度模型的多尺度混合有限元计算格式。数值计算结果表明,在大尺度模型上进行计算能够捕捉到小尺度上的流动特征;相对于传统有限元,多尺度混合有限元能够捕捉小尺度上的非均质特征而具有更高的计算精度,在保证计算精度的同时能够减少计算量。     

基于拓展多尺度有限元的点阵材料结构最小柔顺性设计

论文应用拓展的多尺度有限元法(Extended Multiscale Finite Element Method),以微观构件的截面积为设计变量,研究了体积约束下点阵材料构成结构的最小柔顺性设计问题.建立了适应具有复杂几何形状和载荷边界的点阵材料结构的优化模型,应用序列二次规划算法对悬臂梁和L形梁算例在线性边界条件和周期性边界条件下进行了优化设计,讨论了点阵材料微结构尺寸效应对优化结果的影响,验证了优化模型和求解算法的可靠性,为点阵材料应用于复杂实际工程结构的优化设计提供了新的技术手段.     

计算材料科学中桥域多尺度方法的若干进展

材料科学中存在固有的多尺度特性,桥域多尺度方法是在宏观尺度(如连续介质力学)中引入不同的细微观尺度的计算区域,乃至纳米尺度的分子动力学、量子力学计算区域,将不同尺度的研究方法通过一定的数学模型耦合在一起。该方法既能节约计算成本,又能保证所研究问题的物理特性。本文对多尺度方法的基本概念、跨尺度桥域多尺度方法的发展、基本原理、耦合方法和离散方程进行了讨论,给出了几个应用算例,并在最后进行了总结,展望了今后的可能发展方向。     

非均质材料动力分析的广义多尺度有限元法

自然界和工程中的大部分材料都具有多尺度特征,当考察尺度小到一定程度后,都将表现出非均质性.针对非均质材料的动力问题,提出了一种广义多尺度有限元方法,其基本思想是利用静态凝聚法以及罚函数法构造能够反映单元内部材料非均质特性的多尺度位移基函数.与传统扩展多尺度有限元法中的基函数构造方式不同,广义多尺度有限元法的基函数无需通过在子网格域上多次求解椭圆问题得到,而可直接通过矩阵运算获得.其主要步骤如下:利用数值基函数将一个非均质单胞等效为一个宏观单元,进而形成整个结构的等效刚度矩阵,并得到宏观网格的节点位移,最后再次利用数值基函数得到微观尺度上的位移结果.该广义多尺度有限元法是扩展多尺度有限元法的一种新的拓展,可模拟具有更加复杂几何的非均质单胞的力学行为.通过数值算例,模拟了非均质材料的静力问题、广义特征值问题以及瞬态响应问题,计算结果表明:在边界条件一样的情况下,广义多尺度有限元法的计算结果与传统有限元的计算结果保持高度一致.与传统有限元相比,该方法在保证计算精度的同时极大地提高了计算效率.研究结果表明,广义多尺度有限元法能够很好地模拟非均质单胞的力学行为,具有良好的工程应用潜力.      

结构动力分析自适应有限元方法综述

结构动力分析自适应有限元方法主要研究有限元动力分析的误差估计理论,建立适用于复杂结构动力分析的有限元网格自适应过程．介绍了结构动力问题自适应有限元方法的重要发展,包括固有振动和动响应分析的误差估计及相应的自适应策略;且简要介绍了几种现有的网格生成技术及其特点．最后指出这种方法存在的问题和今后的研究方向．     

应变局部化分析的嵌入强间断多尺度有限元法

固体材料的应变局部化行为是导致结构破坏失效的重要因素之一,开展相关数值模拟分析对于结构安全性评估具有重要意义.然而由于材料的非均质和多尺度特性,采用传统数值方法进行求解时通常需要从最小特征尺度离散求解的结构,这将大幅度增加计算规模和成本.针对这一问题,本文提出了一种基于嵌入强间断模型的多尺度有限元方法.该方法从粗细两个尺度离散求解模型,首先在细尺度单元上引入嵌入强间断模型来描述单元间断特性,所附加的跳跃位移自由度则通过凝聚技术进行消除,从而保持细尺度单元刚度阵维度不变.其次,提出了一种增强多节点粗单元技术,其可根据局部化带与粗单元边界相交情况自适应动态地增加粗节点,新构造的增强数值基函数可以捕捉细尺度间断特性,完成物理信息从细单元到粗单元的准确传递以及宏观响应的快速分析;再次,在细尺度解的计算中,将细尺度解分解为降尺度解与单胞局部摄动解,从而消除弹塑性分析时单胞内部的不平衡力.最后,通过两个典型算例分析,并与完全采用细单元的嵌入有限元结果进行对比,验证了所提出算法的正确性与有效性.     

A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

自适应有限元方法及其工程应用

自适应有限元方法是一种能通过自适应分析自动调整算法以改进求解过程的数值方法．它以常规有限元法为基础,以误差估计和自适应网格改进技术为核心,是一种效率高、可靠性高的计算方法．文中简要介绍并综述了自适应有限元方法的重要发展及应用情况．并对其发展前景作了概要的预测      

基于扩展多尺度有限元法的含液闭孔材料拓扑优化

基于扩展多尺度有限元方法,提出了含液闭孔结构多尺度拓扑优化方法.该多尺度优化方法旨在研究含液闭孔胞元布局对整体含液闭孔结构力学性能的影响.首先针对含液闭孔结构的整体结构柔顺性问题,采用类似SIMP模型对结构的宏观粗网格等效刚度阵进行插值,建立含液闭孔结构柔顺性的拓扑优化列式;其次,针对含液闭孔材料能够利用胞体内部液体腔体积增量产生变形的特性,提出含液闭孔材料柔性机构的概念,并以结构指定位置方向输出位移为目标,建立液体体积膨胀作用下的含液闭孔柔性机构多尺度拓扑优化数学模型.论文基于自主软件平台SiPESC完成了程序研发,并通过数值算例验证了所提出的拓扑优化方法的有效性.     

A FINITE ELEMENT METHOD FOR DIELECTRIC ELASTOMER TRANSDUCERS

We present a finite element method for dielectric elastomer(DE) transducers based on the nonlinear field theory of DE.The method is implemented in the commercial finite element software ABAQUS,which provides a large library functions to describe finite elasticity.This method can be used to solve electromechanical coupling problems of DE transducers with complex configurations and under inhomogeneous deformation.     

非线性方程组的解法:局部弧长法

描述了一个新的非线性方程组的求解方法——局部弧长法．该方法是在弧长法的基础上发展起来的适合于材料非线性有限元分析的数值解法．其约束方程充分利用了结构中破坏区域内的非线性变形信息，有效地解决了材料非线性分析中的稳定性与收敛性问题．数值计算表明，该方法不仅适合于求解结构的极限承载能力，也适合于求解结构达到极限承载参力以后的荷载－变形的全过程     

非均匀介质有限元法

提出适合非均匀介质应力分析的有限元法．文中在有限单元内部采用等参变换方法模拟材料特性的变化，算例表明该法计算效率高，计算精度好．     

FINITE ELEMENT SIMULATION OF GUST LOADING

In the structural design of civil aircraft the critical loads are often those encountered in a gust or atmospheric turbulence. The traditional ‘indicial’ solution is restricted to a simple plate. In this paper a finite element formulation is proposed for an aerofoil or arbitrary shape entering a uniform sharp-edged or sinusoidal gust. The thin rotational gust front and wake in an irrotational flow field are successfully modelled by a novel superposition technique. The finite element solutions are compared with the Kussner function and results by other numerical methods. The agreement is good.