有限元计算细观力学对复合材料力学行为的数值分析

有限元计算细观力学的发展是近十年来细观计算力学发展的主要特征和推动力．本文综述了有限元计算细观力学近十年来应用于复合材料力学行为分析研究方面的进展．介绍了基本的数值模型和计算方法，重点评述了强度和损伤等协同效应问题上的最新研究成果．最后对有限元计算细观力学应用于材料设计的前景做了展望     

THE APPLICATION OF VARIATIONAL MULTISCALE METHOD ON A ONE-DIMENSIONAL MECHANICAL MODEL

将变分多尺度方法应用于一维缆索模型,导出受力缆索的宏观有限元模型并求得细观位移解析解,总结出变分多尺度方法应用于具体模型的关键点和缺陷. 假定刚度为常值,数值模拟一定边界和受力下的缆索,得到宏观和细观位移. 将细观与宏观位移叠加,相比于精确位移得出：细观位移可视为常规有限元模型的后验误差. 变分多尺度方法在一维力学模型中的成功应用,推进了其实用性,为其在更多力学及工程问题中的运用和发展提供了参考.     

具有不同功能梯度分布函数的板件的三维分析

发展了一种新颖的功能梯度结构分析的细观元法。细观元法对结构的常规有限单元内部设置密集细观单元以反映材料特性梯度变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算。这种细观元法既能充分反映材料功能梯度变化特性;而其计算单元又和常规有限元一样,是一种针对功能梯度结构分析的有效数值方法。文中通过细观元技术进行计算、分析,给出了具有不同功能梯度分布函数板件的力学量三维分布形态。     

基于降阶均匀化理论的混凝土双尺度力学性能分析

混凝土材料细观特性对宏观力学性能有着重要影响。为进一步分析混凝土细观特性对宏观力学行为的影响规律,将混凝土材料简化为由骨料、砂浆和界面三相组成,编制了随机凸多面体骨料生成、投放和网格剖分算法,建立可用于有限元计算的满足级配要求的随机细观模型。针对直接使用细观力学模型计算量较大的问题,采用降阶均匀化理论,对混凝土细观胞元模型进行预处理并编制了相应的双尺度计算程序。对不同强度混凝土进行了单轴静态压缩双尺度计算,与实验数据和细观力学模拟结果符合较好。研究表明,降阶均匀化理论在加快求解速度的同时具有较高的精度,可以用于混凝土的多尺度力学性能分析。     

沿板平面材料梯度变化功能材料板件分析

本文发展一种功能梯度构件分析的细观元法.细观元法在构件的常规有限元内部设置密集观细单元以反映材料特性变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算.这种细观元法既能充分反映材料功能梯度沿各方向任意变化特性;而其计算单元又和常规有限元一样,是一种针对功能梯度结构分析的有效数值方法.现有功能梯度板件分析中无论对不同形状还是不同边界的功能梯度构件,其材料特性均沿板厚度方向梯度变化,本文用细观元法进行计算与分析,给出了目前尚未得到的沿板平面方向功能梯度变化构件的力学量三维分布形态.     

纤维增强复合材料风机叶片宏细观一体化有限元分析An unified FEM for macroscopic and mesoscopic analysis of fiber reinforced composite fan blades

为了建立纤维增强复合材料风机叶片宏观性能和细观组分的直接关联,得到一般有限元分析时无法获得的细观参量值,利用FORTRAN程序把细观力学的失效／损伤分析模块,嵌入到有限元软件ABAQUS中的USD‐FLD 用户子程序中,建立了风机叶片宏细观一体化模型。该模型能够实现基于细观组分级损伤／失效判据的宏细观渐进损伤分析和强度预报功能。该模型计算结果与文献中的试验结果有较好的一致性。     

中国力学学会2004年国际、国内学术活动计划

序号 会议名称第六届全国细观力学学术会议 学术内容及议题交流近两年来在细观力学领域取得的新成果,促进各方面对细观力学的了解和重视,进一步推动细观力学的发展第六届全国实验流体力学会议流体力学与空气动力学实验技术研究;流体力学与空气动力学各学科课题实验研究4月18、23日801太原全国非线性有限元学习班邀请美国西北大学Ted Belytsehko教授讲学,内容:连续体和结构的非线性有限元等6月15、17日待定!北京复合材料力学研讨会复合材料力学行为的表征、测试;复合材料细观力学;复合材料结构力学;复合材料应用中的力学问题6月301大连350一…     

纤维增强复合材料风机叶片宏细观一体化有限元分析

为了建立纤维增强复合材料风机叶片宏观性能和细观组分的直接关联,得到一般有限元分析时无法获得的细观参量值,利用FORTRAN程序把细观力学的失效/损伤分析模块,嵌入到有限元软件ABAQUS中的USDFLD用户子程序中,建立了风机叶片宏细观一体化模型。该模型能够实现基于细观组分级损伤/失效判据的宏细观渐进损伤分析和强度预报功能。该模型计算结果与文献中的试验结果有较好的一致性。     

颗粒增强复合材料有效性能的三维数值分析

将细观力学和计算力学方法相结合用以确定复合材料中的局部和平均应力－应变场．对旋转体和非旋转体颗粒增强复合材料的有效模量进行了三维有限元数值计算，数值与实验结果对比表明，该方法是有效的、可靠的．分析了颗粒的排列分布、颗粒取向和颗粒的几何形状对有效模量的影响．数值结果表明，颗粒的排列对有效轴向弹性模量影响较大．颗粒的取向和颗粒的形状对有效性能的影响也是显著的     

混凝土材料宏观力学特性分析的细观单元等效化模型

提出了一种混凝土材料宏观力学特性分析的新方法—细观单元等效化模型。该方法从描述混凝土材料的细观尺度入手,采用Monte Carlo法生成由骨料颗粒及砂浆基质组成的混凝土试件的随机骨料模型;然后,依据混凝土材料特征单元尺度来剖分有限元网格并投影到建立的随机骨料模型上,各细观单元的有效力学特性则采用复合材料等效化方法来确定。本文方法体现了材料非线性宏观力学特性源于其内在的不均匀性这一认识,而对不均匀性的描述则是以网格剖分是否影响其宏观力学特性为准则。因此,本文方法较其他细观力学方法最大的优点在于极大地减小了体系自由度数目（特别是对于三维问题）,从而提高了计算效率。算例分析初步验证了本文方法的高效性。     

Nonlinear micromechanical formulation of the high fidelity generalized method of cells

The recent High Fidelity Generalized Method of Cells (HFGMC) micromechnical modeling framework of multiphase composites is formulated in a new form which facilitates its computational efficiency that allows an effective multiscale material–structural analysis. Towards this goal, incremental and total formulations of the governing equations are derived. A new stress update computational method is established to solve for the nonlinear material constituents along with the micromechanical equations. The method is well-suited for multiaxial finite increments of applied average stress or strain fields. Explicit matrix form of the HFGMC model is presented which allows an immediate and convenient computer implementation of the offered method. In particular, the offered derivations provide for the residual field vector (error) in its incremental and total forms along with an explicit expression for the Jacobian matrix. This enables the efficient iterative computational implementation of the HFGMC as a stand alone. Furthermore, the new formulation of the HFGMC is used to generate a nested local-global nonlinear finite element analysis of composite materials and structures. Applications are presented to demonstrate the efficiency of the proposed approach. These include the behavior of multiphase composites with nonlinearly elastic, elastoplastic and viscoplastic constituents.     

HOMOGENIZATION—BASED TOPOLOGY DESIGN FOR PURE TORSION OF COMPOSITE SHAFTS

In conjunction with the homogenization theory and the finite element method, the mathematical models for designing the corss-section of composite shafts by maximizing the torsion rigidity are developed in this paper. To obtain the extremal torsion rigidity, both the cross-section of the macro scale shaft and the representative microstructure of the composite material are optimized using the new models. The micro scale computational model addresses the problem of finding the periodic microstructures with extreme shear moduli. The optimal microstructure obtained with the new model and the homogenization method can be used to improve and optimize natural or artificial materials. In order to be more practical for engineering applications, cellular materials rather than ranked materials are used in the optimal process in the existence of optimal bounds for the elastic properties. Moreover, the macro scale model is proposed to optimize the cross-section of the torsional shaft based on the tailared composites. The validating optimal results show that the models are very effective in obtaining composites with extreme elastic properties, and the cross-section of the composite shaft with the extremal torsion rigidity. The project supported by the National Natural Science Foundation of China (10172078 and 10102018)     

Experimental and Theoretical Simulation of Heterogeneous Catalysis in Aerothermochemistry (a Review)

The design of aerospace vehicles has required the solution of radically new scientific and technological problems. One of the important problems has been to create reusable heat shield materials. In [1, 2] information concerning the methods and results of solving these problems, including the development of composites from ultrathin quartz fibers and carbon-carbon materials for the “Buran” orbital vehicle heat shield, was presented. The basic thermophysical characteristics of these materials include both the rate or probability coefficients of heterogeneous nitrogen and oxygen atom recombination and the accommodation coefficients of energy recombination at high surface temperatures. In the present paper the experimental and computational aspects of determining these parameters, which are also of interest for new heat shield materials for future space transport systems, are discussed.     

An investigation of intelligent tires using multiscale modeling of cord-rubber composites

A computational model based on the multiscale progressive failure analysis is employed to provide the theoretical predictions for damage development in the cord-rubber composites in tires. Vulcanized rubber, reinforcing belts, and carcass used in tire structures cause the anisotropic behavior under different loading conditions. Steel reinforcement layers made of steel wires combined with rubber complicate the macro-scale finite element modeling of tires. This paper presents a new three-dimensional model of the cord-rubber composite used in tires in order to predict the different types of damage including matrix cracking, delamination, and fiber failure based on the micro-scale analysis. Additionally, intelligent tires have the potential to be widely used to enhance the safety of road transportation systems, and this paper provides an estimation of the effects of void volume fraction, fiber volume fraction, and stacking sequence of the cord-rubber composites on the acceleration profile of the tire measured at the inner-liner.     

机织复合材料的本构关系与成形性研究

对机织复合材料的成形性研究文献进行综述.首先介绍了纺织复合材料的种类、力学和物理特性及其工业背景.接着,从实验、理论、计算3个方面概述了目前对平纹机织复合材料的成形性以及工业冲压技术研究的进展.对于平纹机织复合材料,剪切变形是最主要的变形形式,其本构关系具有非线性的特征.通过对材料样品纯剪切行为的观测,发现了相应的皱曲现象并引入了``锁定角'的概念.大量的实验研究了材料在工业冲压中的可能出现的变形模式及其对皱曲的影响.随后,从理论上介绍了对材料本构关系(包含了拉、压、剪)的系统性模拟,以及冲压过程中有限元数值计算方面的进展.材料的成形性还受到冲压温度、模具对材料的压力、模具与材料表面的摩擦等多种因素的影响.文章对这些问题的研究一一作了简要的介绍,其中许多研究工作不仅得到了重要的结论,还开创了相关课题研究的先河,受到了人们的普遍关注.文章的最后,在总结前人成果的基础上讨论了今后的研究方向,并对该领域的发展作了展望.      

Material sensitivity analysis in homogenization of linear elastic composites

Summary  The main goal of the paper is to present theoretical aspects and the finite element method (FEM) implementation of the sensitivity analysis in homogenization of composite materials with linear elastic components, using effective modules approach. The deterministic sensitivity analysis of effective material properties is presented in a general form for an n-components periodic composite, and is illustrated by the examples of 1D as well as of 2D heterogeneous structures. The results of the sensitivity analysis presented in the paper confirm the usefulness of the homogenization method in computational analysis of composite materials the method may be applied to computational optimization of engineering composites, to the shape-sensitivity studies and, after some probabilistic extensions, to stochastic sensitivity analysis of random composites. Received 10 November 2000; accepted for publication 24 April 2001     

橡胶基复合靶抗贯穿特性的近似分析方法

采用体积含量加权的方法确定了复合靶等效状态方程参数,采用纤维含量平方外推的方法确定了复合靶的等效强度。然后利用球腔膨胀和柱腔贯穿相结合的方法提出了一种分析计算钢球贯穿橡胶基复合靶的近似分析模型。计算方法简明、高效,而且对剩余速度和临界穿透速度等的计算结果与已有实验有良好的符合。     

PVC/纳米Al2O3复合材料的力学性能研究

通过熔融共混方法制备了PVC/纳米Al2O3复合材料,研究了纳米粒子对PVC的增强、增韧效果。采用细观力学方法理论上求解了纳米复合材料的有效弹性模量,比较并分析了试验值和理论计算值的偏差。     

丝束变角度复合材料板壳结构静动力问题研究进展

近年来随着先进自动铺丝技术的发展，生产材料特性随空间位置连续变化的丝束变角度(Variable Angle Tow, VAT)复合材料已成为可能．由此制备的新型复合材料板壳结构不仅具有比传统结构更强的可设计性，而且在提升结构效率方面显示出极大的优越性，是实现航空航天工业装备高性能、轻质化发展的新的重要途径．然而，自动铺丝技术使VAT复合材料板壳结构具有了一般各向异性及面内变刚度特性，给其静动力学问题的分析带来了极大的困难．因此，发展针对该种新型变刚度复合材料板壳结构的力学模型和计算方法至关重要，这也是深入理解其复杂的力学响应机制并进一步促进其在航空航天工程中广泛应用的前提和基础．本文旨在总结近年来有关VAT复合材料板壳结构静动力问题的研究进展，着重从理论分析模型和数值计算方法等方面来简述其最新研究成果，最后讨论了目前VAT复合材料板壳结构静动力问题研究的局限性，并对未来的理论研究进行了展望．     

RELIABILITY-BASED DESIGN OF COMPOSITES UNDER THE MIXED UNCERTAINTIES AND THE OPTIMIZATION ALGORITHM

This paper proposed a reliability design model for composite materials under the mixture of random and interval variables. Together with the inverse reliability analysis technique, the sequential single-loop optimization method is applied to the reliability-based design of composites. In the sequential single-loop optimization, the optimization and the reliability analysis are decoupled to improve the computational efficiency. As shown in examples, the minimum weight problems under the constraint of structural reliability are solved for laminated composites. The Particle Swarm Optimization （PSO） algorithm is utilized to search for the optimal solutions. The design results indicate that, under the mixture of random and interval variables, the method that combines the sequential single-loop optimization and the PSO algorithm can deal effectively with the reliability-based design of composites.