具有频率约束的桁架结构可靠性拓扑优化

基于结构的可靠性,建立具有频率约束的桁架结构拓扑优化模型,即以桁架的横截面积为设计变量、重量最小为优化目标,位移、应力等可靠性及基频为约束,然后从优化策略--拓扑组方法出发,先考虑桁架结构可能的优化布局,以避免优化过程中节点和杆件的增加或减少引起解的奇异性,且可在一定程度上减少计算工作量.从工程实际出发,对结构系统的可靠性隐形约束进行等价显化处理,最终转化为常规的横截面积优化问题.算例表明了文中所提方法的简单性、有效性和合理性.     

连续体结构拓扑优化

对连续体结构的拓扑优化，给出一种工程实用方法：将拓扑优化分两步进行，首先解决在弹性体内哪些区域需要删除的问题，然后再确定删除区的边界，这种方法适用于各种约束条件的问题，而且拓扑清晰。     

SIMULTANEOUS SHAPE AND TOPOLOGY OPTIMIZATION OF TRUSS UNDER LOCAL AND GLOBAL STABILITY CONSTRAINTS

A new approach for the solution of truss shape and topology optimization problem sunder local and global stability constraints is proposed. By employing the cross sectional areas of each bar and some shape parameters as topology design variables, the difficulty arising from the jumping of buckling length phenomenon can be easily overcome without the necessity of introducing the overlapping bars into the initial ground structure. Therefore computational efforts can be saved for the solution of this kind of problem. By modifying the elements of the stiffness matrix using Sigmoid function, the continuity of the objective and constraint functions with respect to shape design parameters can be restored to some extent. Some numerical examples demonstrate the effectiveness of the proposed method.     

基于拓扑描述函数的特定性能材料设计方法

研究了拓扑描述函数在材料设计中的应用,给出了一种基于拓扑描述函数的特定性能材料设计问题的提法和求解方法.将拓扑描述函数表示成含参数的基函数之和,将材料微结构拓扑优化问题转化为设计基函数描述参数的尺寸优化问题,使问题求解更方便.基于拓扑描述函数的方法可以准确确定设计域上任意点的材料分布,避免了变密度法常出现的棋盘效应、设计变量和有限元单元相关的缺点;与传统的水平集方法相比,其优化模型可以利用现有的优化方法求解,避免了差分法求偏微分方程.具有正泊松比和负泊松比的特定弹性性能材料的设计算例,说明了基于拓扑描述函数的材料设计方法的有效性.     

改进的遗传算法在连续体结构多目标拓扑优化中的应用

连续体结掏多目标拓扑优化是结构优化领域中的一个较难的研究课题.本文提出了一种改进的SPEA2多目标优化算法.该算法中采用数学形态学中的四方向链码的编码方式进行结构拓扑表达,使产生的拓扑结构清晰且无异议,完全消除了模糊的拓扑边界和棋盘格现象.提出了"质量向量"的概念,用以量化两个结构拓扑之间的相似性,从而使得不同结构拓扑之间的相似性比较成为可能.同时还将机器学习中的范例学习的思想融入到新算法中,利用前面的有限元分析结果对后来的结构分析进行指导,剔除了很多不必要的重复的计算,从而使算法的总计算量大幅度地减少.将新算法应用于悬臂和简支两类深梁的多目标拓扑优化,获得了高质量的Pareto最优解,且具有很好的分布.     

油船货舱结构拓扑优化工程应用

基于双向渐进结构优化方法,采用自编程的形式,结合成熟的有限元计算处理程序,对某型油船货舱船体横向强框架和水平桁进行了拓扑优化设计;并提出一套适用于复杂船体结构拓扑优化设计流程,包括船体结构的拓扑优化基结构合理选取以及多工况下油船货舱结构拓扑优化的加权系数的选取方式等。最终提取得到了合理且原创的拓扑构型,并将其成功地应用到新船型结构设计之中,使新船型的横向强框相较传统油船重量减轻28.6%,单位长度内货舱区结构的重量减轻5.4%。     

桁架结构拓扑优化设计的线性规划方法

本文提出了一种新的桁架结构拓扑优化设计方法,在该方法中,以杆件内力为设计变量,以由结构力学的基本方程构成的位移、应力等物理量为约束,构成了拓扑优化的线性规划模型。它克服了目前桁架结构拓扑优化的两大困难——预定设计位移场与在拓扑优化过程中无法考虑位移、应力等性态约束。文章最后给出了两个考题,说明了本方法的可行性与有效性。     

带有预应力的连续体组合结构拓扑优化

考虑梁结构的预应力，导出了结构的应力及其灵敏度公式。同时，结合结构单元应力水平和相对差商，对于带有尺寸和拓扑变量的连续体组合结构优化问题，建立了一套优化准则，在双方向渐进结构优化方法思想的基础上，形成了一种新的拓扑优化算法。最后，算例表明了该方法的正确性和有效性。     

基于拓扑描述函数的连续体结构拓扑优化方法

提出了一种利用拓扑描述函数(TDF)作为拓扑设计变量求解连续体结构拓扑优化问题 的新方法. 优化问题的目标函数是结构的整体柔顺性，约束条件为对于可利用材料的体积限 制. 这种方法不仅可以消除拓扑优化中经常出现的棋盘格式等数值不稳定现象，而且能够有 效地抑制传统算法处理此类优化问题时所引发的边界扩散效应. 与其它的基于水平集描述函 数的拓扑优化方法相比，所提出的算法不仅无需求解控制水平集函数演化的双曲守恒方 程，而且合理地考虑了目标函数的拓扑导数信息，因而使得算法的计算效率有了显著的提高.     

结构优化技术在高速飞行器上的应用与面临的挑战

高速飞行器是航空航天领域发展的重要方向,具有重大战略意义和重要经济社会价值。首先,对结构优化技术及其在飞行器零部件和全机结构概念设计中的应用进行综述;然后,根据典型使用环境和设计要求,总结高速飞行器结构设计中存在的主要问题和应用需求,着重介绍近年来在高速飞行器零部件和全机结构概念设计中应用结构优化技术的典型案例;最后,提出近期亟待突破或完善的结构优化技术。结构优化技术与设计经验相结合带来设计理念的变革,将成为先进飞行器设计的必备工具和标准流程。面向重大需求的结构优化理论研究与工程实践对提升我国航空航天核心竞争力具有重要价值。     

复合材料周期性线弹性微结构的拓扑优化设计

提出复合材料周期性线弹性微结构拓扑优化设计的模型，模型1设计具有极值弹性特性的复合材料，模型2设计工况最刚微结构单胞。通过该模型和均匀化技术可以获得优化的微结构单胞，进而改善或者得到最优宏观特性的复合材料。为了便于制造和应用，用胞体材料而不是多相材料来得到复合材料的极值弹性特性和最大刚度。优化结果表明，该模型与数值方法相结合可以有效地实现微结构的拓扑优化设计。     

复合材料扭转轴截面微结构拓扑优化设计

提出复合材料扭转轴截面微结构拓扑优化设计新模型，模型的优化目标是获得具有最大宏观剪切特性加权和的单胞形式．通过模型和均匀化方法及优化技术可以获得优化的微结构单胞，进而改善或者得到最优宏观弹性特性的复合材料．为了便于制造和应用，胞体材料用来获得复合材料的极值剪切模量．最后的优化结果表明，该模型连同数值处理技巧可以非常有效地实现微结构的拓扑优化设计．     

A Universal Expression for Bounds on the Torsional Rigidity of Cylindrical Shafts Containing Fibers with Arbitrary Transverse Cross-Sections

We derive upper and lower bounds for the torsional rigidity of host shafts containing a number of cylindrical fibers. The transverse cross-sections of the host shaft and the fibers are simply connected, but could be arbitrary in shape. Utilizing the fact that the torsion solution of a homogeneous host shaft with simply connected cross-section can be known, we propose a method to construct statically and kinematically admissible fields interior to the matrix and to the fibers. Previous developments on bounding the torsional rigidity of composite shaft so far are confined to circular fibers. Here we try to simulate fibers with non-circular cross-section and incorporate the interactions of the cross-sectional shapes of the host shaft and the fibers at the same time. Proceeding from extremal principles of elasticity, together with propositions of some domain integration procedures, we provide a universal expression for bounds on the torsional rigidity of the composite shaft. The exact expressions depend on the constituent information of the fibers and the host shaft, which could offer useful information to tailor the shape and the arrangement of the constituents to achieve an optimal value.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

Rigorous Bounds on the Torsional Rigidity of Composite Shafts with Imperfect Interfaces

We derive upper and lower bounds for the torsional rigidity of cylindrical shafts with arbitrary cross-section containing a number of fibers with circular cross-section. Each fiber may have different constituent materials with different radius. At the interfaces between the fibers and the host matrix two kinds of imperfect interfaces are considered: one which models a thin interphase of low shear modulus and one which models a thin interphase of high shear modulus. Both types of interface will be characterized by an interface parameter which measures the stiffness of the interface. The exact expressions for the upper and lower bounds of the composite shaft depend on the constituent shear moduli, the absolute sizes and locations of the fibers, interface parameters, and the cross-sectional shape of the host shaft. Simplified expressions are also deduced for shafts with perfect bonding interfaces and for shafts with circular cross-section. The effects of the imperfect bonding are illustrated for a circular shaft containing a non-centered fiber. We find that when an additional constraint between the constituent properties of the phases is fulfilled for circular shafts, the upper and lower bounds will coincide. In the latter situation, the fibers are neutral inclusions under torsion and the bounds recover the previously known exact torsional rigidity.      

多相材料微结构多目标拓扑优化设计

在采用多尺度均匀化方法求解微结构等效特性的基础上，提出了多相材料 微结构的多目标优化设计模型. 以组分材料用量为约束，采用周长控制消除棋盘格，结合有 限元方法和对偶凸规划求解技术，对两相和三相材料微结构多项等效模量的组合进行了优化 设计. 研究比较了微结构网格粗细、材料组分以及三相材料微结构优化中的两相实体材料弹 性模量相对比例不同对优化结果的影响. 数值算例验证了优化模型和优化算法的有效性，表 明了相关因素对优化结果的影响.     

Homogenization of rectangular cross-section fibre-reinforced materials: bending–torsion effects

We study the homogenization of an elastic material in contact with periodic parallel elastic rectangular cross-section fibres of higher rigidity. The interactions between the matrix and the fibres are described by a local adhesion contact law with interfacial adhesive stiffness parameter depending on the period. Assuming that the Lamé constants in the fibres and the stiffness parameter have appropriate orders of magnitude, we derive a class of energy functionals involving extension, flexure and torsion terms.     

A micromechanics-based approach for the derivation of constitutive elastic coefficients of strain-gradient media

A micromechanics-based approach for the derivation of the effective properties of periodic linear elastic composites which exhibit strain gradient effects at the macroscopic level is presented. At the local scale, all phases of the composite obey the classic equations of three-dimensional elasticity, but, since the assumption of strict separation of scales is not verified, the macroscopic behavior is described by the equations of strain gradient elasticity. The methodology uses the series expansions at the local scale, for which higher-order terms, (which are generally neglected in standard homogenization framework) are kept, in order to take into account the microstructural effects. An energy based micro–macro transition is then proposed for upscaling and constitutes, in fact, a generalization of the Hill–Mandel lemma to the case of higher-order homogenization problems. The constitutive relations and the definitions for higher-order elasticity tensors are retrieved by means of the “state law” associated to the derived macroscopic potential. As an illustration purpose, we derive the closed-form expressions for the components of the gradient elasticity tensors in the particular case of a stratified periodic composite. For handling the problems with an arbitrary microstructure, a FFT-based computational iterative scheme is proposed in the last part of the paper. Its efficiency is shown in the particular case of composites reinforced by long fibers.     

Multiscale homogenization of n-component composites with semi-elliptical random interface defects

Effective elastic characteristics of periodic multicomponent composite materials with random interface defects are studied in the paper. The defects are assumed to be semi-elliptical and lying with major semi axes along the interfaces, where minor and major semi-axes as well as the defects number are given as input random variables. The homogenization approach has a multiscale character—some algebraic approximation is used first to calculate effective elastic parameters of the interphase including all defects located at the same interface. Equations for interphase random elastic parameters are obtained using MAPLE symbolic mathematics in conjunction with probabilistic generalized perturbation method. A different homogenization method is applied at the micro scale, where the cell problem is solved numerically using the Finite Element Method (FEM) program. Since the composites considered exhibit random variations of both elastic properties and the interface defects, the overall homogenized characteristics must be obtained as random quantities, which is realized on the micro scale by the Monte-Carlo simulation. The proposed interface defects model obeys the porosity effects resulting from the nature of some matrices in engineering composites as well as the interface cracks appearing as a result of composites ageing during static or fatigue fracture.     

Objective evaluation of linearization procedures in nonlinear homogenization: A methodology and some implications on the accuracy of micromechanical schemes

A systematic methodology for an accurate evaluation of various existing linearization procedures sustaining mean fields theories for nonlinear composites is proposed and applied to recent homogenization methods. It relies on the analysis of a periodic composite for which an exact resolution of both the original nonlinear homogenization problem and the linear homogenization problems associated with the chosen linear comparison composite (LCC) with an identical microstructure is possible. The effects of the sole linearization scheme can then be evaluated without ambiguity. This methodology is applied to three different two-phase materials in which the constitutive behavior of at least one constituent is nonlinear elastic (or viscoplastic): a reinforced composite, a material in which both phases are nonlinear and a porous material. Comparisons performed on these three materials between the considered homogenization schemes and the reference solution bear out the relevance and the performances of the modified second-order procedure introduced by Ponte Castañeda in terms of prediction of the effective responses. However, under the assumption that the field statistics (first and second moments) are given by the local fields in the LCC, all the recent nonlinear homogenization procedures still fail to provide an accurate enough estimate of the strain statistics, especially for composites with high contrast.