纺织复合材料和结构多尺度耦合的数值分析

研究了纺织复合材料和结构多尺度耦合的数值分析模型。建立了微、细观单胞,给出了纺织复合材料平均弹性常数的逐级分析方法,着重研究了由宏观结构、到细观纤维束、再到微观纤维三个尺度耦合的应力分析方案。对于常用的板壳状纺织复合材料结构,在面内载荷下,假设每层细观单胞的平均面内应变是一致的,在弯曲、横向剪切及扭曲等非面内载荷下,在内力等效条件下将沿厚度方向连续分布的宏观应力简化为阶梯状分布,忽略了每层细观单胞范围内宏观应力沿厚度方向的梯度变化,由此利用细观单胞模型实现宏观应力与细观应力之间的传递,再利用微观单胞可得到纤维尺度的微观应力。最后以一种三维机织复合材料为例,用上述多尺度耦合的模型逐级分析了材料的平均弹性常数,并沿相反方向,由宏观结构分析逐级计算出纤维束尺度和纤维尺度的细、微观应力的局部波动。     

基于多尺度方法的三维编织C/C复合材料弯曲失效特性分析

为了预测三维编织C/C复合材料的弯曲失效行为,基于多尺度渐进展开理论,结合细观渐进损伤模型,建立了三维编织C/C复合材料宏细观多尺度分析模型。通过商业有限元软件ABAQUS用户子程序UMAT的二次开发,在宏观结构有限元分析中实时调用细观单胞模型进行细观渐进损伤分析,实现了宏细观尺度之间交互式信息传递和多尺度损伤模拟。利用上述模型对三点弯曲载荷下三维编织C/C复合材料梁的渐进损伤和失效过程进行了模拟,预测了梁的载荷-挠度曲线和弯曲强度,并与实验结果进行了对比分析,验证了基于多尺度方法的三维编织C/C复合材料弯曲强度预测模型的有效性,为此类材料及结构失效分析提供了一种手段。     

Biot-Cosserat连续体-离散颗粒集合体模型的非饱和土连接尺度方法

耦合了非饱和多孔多相介质有限元模型和颗粒介质离散元(DEM)模型,提出了以宏、细观尺度分别耦合Biot-Cosserat连续体模型和离散颗粒集合体模型的连接尺度方法来分析非饱和含液颗粒材料的力学渗流耦合问题。根据被动空气压力假定和对空间离散孔隙水质量守恒方程的约化,从非饱和土有限元控制方程的基本未知量中消去了孔隙水压力,而将其取作有限元积分点上定义取值的内状态变量,进而建立了节点未知量仅包含固相线位移和转角的非饱和Cosserat多孔连续体约化有限元数值模型。基于连接尺度方法(BSM)宏、细观尺度数值过程的解耦计算的特点,对宏、细观两尺度数值模型的时域积分分别采用隐式Newmark方法和显式中心差分法,且取不同时间步长以提高计算效率。与全域采用DEM的精细分析方法相比,本文BSM在保证计算精度的前提下可大幅节省计算时间。在不考虑湿化效应的边坡稳定算例中,在得到类似计算精度条件下它比全域采用DEM节省计算时间高达86.65%。二维边坡稳定算例结果验证了本文连接尺度方法的有效性,以及在揭示含液颗粒结构细观破坏机理上的优点。数值算例结果显示,边坡承载能力因降雨大幅下降约50%,这表明本文发展的计及伴随湿化过程的颗粒材料结构中饱和度及吸力分布演变及其对结构破坏失效影响的非饱和颗粒材料多尺度计算模型是很有必要的。     

含损伤演化的TM耦合数值模型及其应用研究

从岩石材料的细观结构层次出发，应用损伤力学和热弹性理论，对热力耦合作用下岩石破裂 过程中热－应力相互作用关系进行了分析. 初步建立了细观岩石热－力(TM)耦合数值模型， 探讨了TM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制，把岩石 热固耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中. 运用该数值模型对某 硬岩实验室开展的原位尺度实验中的废料处理井间柱稳定性进行了模拟分析，其应力场、岩 石剥离破坏形态及破坏诱发的AE特性等均与实验监测结果表现出了较好的一致性，证 明了该数值模型的合理性和有效性.     

周期性复合材料湿热力耦合问题的二阶渐近分析

渐近均匀化方法是一种很好的分析复合材料力学性能的方法。通过建立考虑材料湿热效应的双尺度渐近展开模型,研究了周期性复合材料在湿、热、力耦合作用下的结构响应。将带有小参数的位移渐近展开式代入湿热力耦合控制方程,利用等效积分"弱"变换推导出复合材料等效参数的单胞控制方程以及均匀化方程。为了提高计算精度,本文对均匀化解进行了一阶和二阶修正,利用FreeFem++软件进行了数值模拟。结果显示:双尺度二阶渐近分析能够高效、准确地计算出一系列多场耦合问题,为研究周期性复合材料的力学性能提供了新的方法。     

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

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

Nb3Sn高场复合超导体临界性能力学变形效应的多尺度模拟

A15型Nb₃Sn超导体是制造高场( > 10 T)超导磁体线圈的主要材料, 被广泛应用于磁约束可控核聚变、高能物理等强磁场超导磁体装备制造领域. 力学变形诱导的Nb₃Sn超导临界性能退化给高场超导磁体装备的电磁性能指标和安全稳定运行造成了极其不利的影响. 鉴于Nb₃Sn超导体具有复杂的多尺度结构特征, 不同尺度下变形与超导电性能耦合行为是相互关联的, 本文建立了考虑微/细/宏观关联的非线性力电磁耦合本构模型, 提出了从原子尺度A15晶体结构到超导体微结构到宏观非均质Nb₃Sn复合超导体的多尺度模拟模型. 基于多晶体有限元方法, 对静水压加载条件下Nb₃Sn多晶体超导临界温度衰退和单轴拉压加载条件下Nb₃Sn复合多晶体临界性能衰退行为进行了模拟预测, 预测结果与实验观测结果定性吻合. 该模型揭示了Nb₃Sn复合超导体变形-超导电性能多尺度耦合机理, 实现对高场超导体力、电、磁、热耦合行为的预测, 有助于提高对A15型金属间化合物高场超导复合材料力、电、磁、热多尺度耦合行为的认识和描述能力, 为强磁场超导磁体的设计与制造提供有力的理论支撑.      

混凝土受拉性能多尺度均匀化数值模拟研究

为研究混凝土的细观结构对其受拉性能的影响,本文首先采用蒙特卡罗方法生成了多组混凝土随机骨料模型;然后划分有限元网格并映射到所建立的混凝土随机骨料数值模型上,采用复合材料均匀化方法建立了混凝土多尺度均匀化数值模型。考虑骨料随机分布和骨料形状,通过数值分析得到了在单轴拉伸情况下混凝土损伤分布特性及宏观应力-应变关系,在考虑随机性与非均质性的同时得到了较为理想的结果。运用多尺度均匀化的建模方法不仅能反应混凝土细观结构的影响,而且等效化模型的模拟计算时间大大减少,节约计算资源,提高了计算效率,具有优越性。     

复合材料弹塑性多尺度分析模型与算法

对材料非线性多尺度分析的计算模型与算法进行研究.在构建周期分布单胞分析算法的基础上,发展针对复合材料结构材料非线性多尺度分析的一般有限元方法.方法的特点是将所建立的单胞分析过程作为有限元分析的子程序嵌入到总程序系统当中,完成对应的高斯点应力计算,因而使所发展的方法具有实现方便的特点.给出数值计算结果,验证了方法与所发展的多尺度有限元分析程序的正确与有效性.     

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

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

Particle methods for complex flows in chemical engineering — the pseudo-particle approach

The multi-scale structures of complex flows in chemical engineering have been great challenges to the design and scaling of such systems, and multi-scale modeling is the natural way in response. Particle methods (PMs) are ideal constituents and powerful tools of multi-scale models, owing to their physical fidelity and computational simplicity. Especially, pseudo-particle modeling (PPM, Ge & Li, 1996; Ge & Li, 2003) is most suitable for molecular scale flow prediction and exploration of the origin of multi-scale structures; macro-scale PPM (MaPPM, Ge & Li, 2001) and similar models are advantageous for meso-scale simulations of flows with complex and dynamic discontinuity, while the lattice Boltzmann model is more competent for homogeneous media in complex geometries; and meso-scale methods such as dissipative particle dynamics are unique tools for complex fluids of uncertain properties or flows with strong thermal fluctuations. All these methods are favorable for seamless interconnection of models for different scales.However, as PMs are not originally designed as either tools for complexity or constituents of multi-scale models, further improvements are expected. PPM is proposed for microscopic simulation of particle-fluid systems as a combination of molecular dynamics (MD) and direct simulation Monte-Carlo (DSMC). The collision dynamics in PPM is identical to that of hard-sphere MD, so that mass, momentum and energy are conserved to machine accuracy. However, the collision detection procedure, which is most time-consuming and difficult to be parallelized for hard-sphere MD, has been greatly simplified to a procedure identical to that of soft-sphere MD. Actually, the physical model behind such a treatment is essentially different from MD and is more similar to DSMC, but an intrinsic difference is that in DSMC the collisions follow designed statistical rules that are reflection of the real physical processes only in very limited cases such as dilute gas.PPM is ideal for exploring the mechanism of complex flows ab initio. In final analysis, the complexity of flow behavior is shaped by two components on the micro-scale: the relative displacements and interactions of the numerous molecules. Adding to the generality of the characteristics of complex system as described by Li and Kwauk (2003), we notice that complex structures or behaviors are most probably observed when these two components are competitive and hence they must compromise, as in the case of emulsions and the so-called soft-matter that includes most bio-systems. When either displacement or interaction is dominant, as in the case of dilute gas or solid crystals, respectively, complexity is much less spectacular. Most PMs consist explicitly of these two components, which is operator splitting in a numerical sense, but it is physically more meaningful and concise in PPM.The properties of the pseudo-particle fluid are in good conformance to typical simple gas (Ge et al., 2003; Ge et al., 2005). The ability of PPM to describe the dynamic transport process on the micro-scale in heterogeneous particle-fluid systems has been demonstrated in recent simulations (Ge et al., 2005). Especially, the method has been employed to study the temporal evolution of the stability criterion in the energy minimization multi-scale model (Li & Kwauk, 1994), which confirms its monotonously asymptotic decreasing as the model has assumed (Zhang et al., 2005). Massive parallel processing is also practiced for simulating particle-fluid systems in PPM, indicating an optimistic prospect to elevate the computational limitations on their wider applications, and exploring deeper underlying mechanism in complex particle-fluid systems.     

A multi-scale architecture for multi-scale simulation and its application to gas–solid flows

A multi-scale hardware and software architecture implementing the EMMS (energy-minimization multi-scale) paradigm is proven to be effective in the simulation of a two-dimensional gas–solid suspension. General purpose CPUs are employed for macro-scale control and optimization, and many integrated cores (MICs) operating in multiple-instruction multiple-data mode are used for a molecular dynamics simulation of the solid particles at the meso-scale. Many cores operating in single-instruction multiple-data mode, such as general purpose graphics processing units (GPGPUs), are employed for direct numerical simulation of the fluid flow at the micro-scale using the lattice Boltzmann method. This architecture is also expected to be efficient for the multi-scale simulation of other complex systems.     

A multi-scale architecture for multi-scale simulation and its application to gas-solid flows

A multi-scale hardware and software architecture implementing the EMMS （energy-minimization multi-scale） paradigm is proven to be effective in the simulation of a two-dimensional gas-solid suspension. General purpose CPUs are employed for macro-scale control and optimization, and many integrated cores （MlCs） operating in multiple-instruction multiple-data mode are used for a molecular dynamics simulation of the solid particles at the meso-scale. Many cores operating in single-instruction multiple- data mode, such as general purpose graphics processing units （GPGPUs）, are employed for direct numerical simulation of the fluid flow at the micro-scale using the lattice Boltzmann method. This architecture is also expected to be efficient for the multi-scale simulation of other comolex systems.     

Multi-scale FE computation for the structures of composite materials with small periodic configuration under condition of coupled thermoelasticity

In this paper, the multi-scale computational method for a structure of composite materials with a small periodic configuration under the coupled thermoelasticity condition is presented. The two-scale asymptotic (TSA) expression of the displacement and the increment of temperature for composite materials with a small periodic configuration under the condition of thermoelasticity are briefly shown at first, then the multi-scale finite element algorithms based on TSA are discussed. Finally the numerical results evaluated by the multi-scale computational method are shown. It demonstrates that the basic configuration and the increment of temperature strongly influence the local strains and local stresses inside a basic cell. The project supported by the National Natural Science Foundation of China (19932030) and Special Funds for Major State Basic Research Projects     

延性金属层裂自由面速度曲线特征多尺度模拟研究

以延性金属钽为研究对象,对钽在平板撞击下的层裂行为进行了多尺度下的数值模拟研究,从微观视角对自由面速度曲线上的典型特征进行了新的解读。在宏观尺度,对比分析了光滑粒子流体动力学法（smootfied particle hydrodynamics, SPH）与Lagrange网格法以及几种本构模型的模拟结果及其适用性。通过与实验数据的对比表明,Steinberg-Cochran-Guinan本构模型在层裂模拟中与实验数据吻合较好,通过改变加载条件获得了不同应变率下的自由面速度曲线,分析了不同应变率下的自由面速度曲线中的典型特征。在微观尺度,采用分子动力学方法获得层裂区域内损伤演化情况,揭示了宏观尺度自由面速度曲线典型特征所蕴含的物理内涵。分析表明,层裂表现为材料内部微孔洞形核、长大和聚集的损伤演化过程,自由面速度曲线上的典型特征与层裂区域的损伤演化过程存在密切关联。Pullback信号是层裂区域内微孔洞形核的宏观表征;自由面速度曲线的下降幅值在一定程度上反映了微孔洞的形核条件,由此计算得到的层裂强度实际上是微孔洞的形核强度。此外,Pullback信号后的速度回跳速率反映了微损伤演化的速率。     

连接结构接触界面非线性力学建模研究

连接界面上存在的跨尺度、多物理场和非线性行为是引起结构复杂非线性动力学的主要原因。由于连接界面力学行为的复杂性,以及对连接界面进行直接试验观测的困难,连接界面的力学建模一直是非常具有挑战性的科学问题。本文首先从分析结合面的跨尺度物理机理入手,将名义的光滑平面视作凹凸不平的粗糙面,考虑单个微凸体的黏滑摩擦行为,建立接触载荷与变形的非线性关系,然后采用GW(Greenwood和Williamson, GW)模型数理统计方法建立整个粗糙界面的跨尺度力学模型,并与公开文献中试验结果进行对比。然后考虑连接界面典型非线性特征,提出一种改进的Iwan唯象模型,利用精细有限元方法获得非线性特征结果,采用系统辨识理论建立连接结构的降阶力学模型,并利用有限元结果进行模型验证。结果表明,本文提出的粗糙界面跨尺度模型在法向载荷较小时与试验结果吻合较好,改进的Iwan模型能够较好描述连接界面非线性特征,并与有限元结果吻合较好。      

箱型结构损伤分析的空间网格多尺度模拟方法

为高效模拟空间效应显著的结构在关键局部的损伤，本文研究建立了以损伤分析为目标的空间网格多尺度模拟方法．首先基于变形协调法和内力平衡法，研究了空间网格多尺度建模中的跨尺度界面连接方法，对比分析了两种界面连接合理性．以三跨连续刚构混凝土箱梁在地震载荷下的损伤为分析案例，验证了空间网格多尺度模型在结构损伤分析中的可行性及其在计算效率上的优越性．分析结果表明：空间网格多尺度模型可以精确模拟结构的静力效应和动力特性；空间网格多尺度模型既考虑了结构空间效应，又可以高效分析箱梁结构局部易损部位的损伤演化过程，从而为空间效应显著的结构损伤分析提供了更为实用精细的计算模拟方法．     

基于渐进均匀化的平纹编织复合材料低速冲击多尺度方法

针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明:宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.      

Multi-scale numerical simulation of fluidized beds: Model applicability assessment

In the past few decades, multi-scale numerical methods have been developed to model dense gas-solid flow in fluidized beds with different resolutions, accuracies, and efficiencies. However, ambiguity needs to be clarified in the multi-scale numerical simulation of fluidized beds: (i) the selection of the sub-models, parameters, and numerical resolution; (ii) the multivariate coupling of operating conditions, bed configurations, polydispersity, and additional forces. Accordingly, a state-of-the-art review is performed to assess the applicability of multi-scale numerical methods in predicting dense gas-solid flow in fluidized beds at specific fluidization regimes (e.g., bubbling fluidization region, fast fluidization regime), with a focus on the inter-particle collision models, inter-phase interaction models, collision parameters, and polydispersity effect. A mutual restriction exists between resolution and efficiency. Higher-resolution methods need more computational resources and thus are suitable for smaller-scale simulations to provide a database for closure development. Lower-resolution methods require fewer computational resources and thus underpin large-scale simulations to explore macro-scale phenomena. Model validations need to be further conducted under multiple flow conditions and comprehensive metrics (e.g., velocity profiles at different heights, bubbles, or cluster characteristics) for further improvement of the applicability of each numerical method.