铺设角度与铺层顺序对层合板稳定性的影响

以层合板结构的临界屈曲载荷系数最大化为优化目标,基于改进型模拟退火算法对层合板结构铺设角度和铺层顺序进行优化。由于层合板结构的铺层角度是离散变量,模拟退火算法适合求解离散变量的优化问题。利用模拟退火算法优化层合板铺层,在算法内采用并行计算、引入记忆功能同时设置双阈值终止准则,有效地提高了优化过程的收敛速度,同时避免优化过程中出现局部最优解。以临界屈曲载荷系数作为目标函数,选取复合材料层合板的铺设角度顺序为设计变量,采用改进的模拟退火算法得出复合材料层合板的最优铺设角度以及铺层顺序。     

基于双值参数化方法且考虑制造约束的炭纤维复合材料铺层优化

基于梯度的优化方法对炭纤维复合材料层合板的铺层数量和顺序进行优化。优化问题中以铺层质量为目标,并以刚度和制造约束为约束。采用改进双值参数化方法对铺层的材料性能进行插值,并基于凸规划对偶算法对优化问题进行求解。为了适应凸规划对偶算法的特点,将关于铺层角度的制造约束表述为少量非线性约束。同时引入离散度约束和惩罚指数以消除优化结果中的中间变量。算例结果验证了该优化方法的有效性。     

基于可靠性的复合材料结构稳定性约束优化设计

基于结构的可靠性, 研究了复合材料结构的稳定性约束优化设计方法。考虑材料及载荷的不确定性, 通过结构可靠性分析的响应面法和有限元法的结合, 对复合材料结构稳定性进行可靠性分析; 利用优化软件iSIGHT集成可靠性分析程序, 实现了以铺层层数及铺层角度为设计变量的复合材料结构稳定性约束问题的可靠性优化方法。对层合板及层合圆柱进行算例分析, 验证了本文中可靠性优化方法的有效性, 为工程实际中的复合材料结构稳定性约束优化设计问题提供借鉴。      

基于改进模拟退火算法的复合材料层合板频率优化

针对复合材料层合板频率优化问题,结合可行规则法和直接搜索模拟退化算法,提出了一种自适应模拟退火(SA)改进算法。层合板优化目标是基频、频率带隙以及给定基频和带隙约束的层合板厚度。设计变量包括铺层角度和铺层数两种离散变量。改进算法的自适应新点产生模块采用依赖温度的动态调整搜索半径,改善了直接搜索模拟退化(DSA)算法易陷入局部极值的缺陷,而可行规则法的引入提高了SA算法求解约束问题的效率和简易性。采用Ritz法进行频率响应分析以考虑弯扭耦合影响。不同铺层数、角度增量和长宽比时的层合板3类算例结果显示:改进算法能有效求解层合板频率优化,可获得更多或更好的铺层顺序全局优化解。     

基于几何因子的复合材料气动弹性剪裁设计

实现了基于几何因子的复合材料层合板建模,解决了几何因子与Natran的参数输入问题,并根据工艺约束中的最小铺层比例对几何因子可行空间进行了推导补充。在此基础上,提出了一种基于几何因子和Nastran的复合材料气动弹性剪裁优化设计方法。首先以总厚度和几何因子作为设计变量以及以Nastran作为求解器,以强度、刚度、颤振和发散速度以及几何因子相关性约束作为约束条件进行结构寻优,得到最优的铺层总厚度和几何因子。其次,以最优几何因子作为目标,进行铺层结构逆问题求解,约束条件为复合材料铺层工艺约束。因几何因子为铺层厚度和铺层顺序的表达式,与传统的多级优化相比,以几何因子作为设计变量可以避免铺层厚度和铺层顺序的解耦,进而获得更大的设计空间,且得到的铺层结构可以满足工艺约束。最后,对一矩形悬臂复合材料层合板进行剪裁设计,使得铺层结构满足气动弹性约束且质量最小。结果显示,运用该优化方法可以得到质量更小且满足工艺约束的铺层结构。     

基于遗传算法的复合材料层合板削层结构铺层优化

针对应力变化较大的碳纤维增强复合材料层合板,提出削层结构铺层分级优化模式。通过将结构分解为若干子铺层并对各子铺层的位置、尺寸、铺层数以及铺层顺序进行优化,得到了满足强度和可制造性要求且质量最小的结构设计方案。该模式的第1、2级优化利用参考层对各子铺层位置及尺寸进行优化,第3级优化通过引入3次样条插值参数化方法对各子铺层层数和铺层顺序进行优化。参考层的引入可减少设计变量的数量,3次样条插值参数化方法可解决以铺层角为设计变量时设计变量数目不确定的问题。利用有限元方法对结构进行力学分析计算,并依据Tsai-Wu准则确定结构强度。在第2、3级优化中利用遗传算法对优化问题进行求解。算例计算表明:削层结构铺层分级优化模式结果合理可信。与均匀铺层方法结果比较可知:削层结构可有效减少结构质量。     

传递矩阵法分析复合材料层合板的传声损失

为了得到不同频率下正交各向异性复合材料层合板的传声损失,基于传递矩阵的方法,推导出层合板的传声损失计算公式。通过建立复合材料层合板的传声计算模型,研究了层合板铺设角度、板厚度和板密度等结构参数对层合板的传声损失影响。计算结果表明：复合材料的密度与传声损失之间没有明显的线性关系,而是随着频率的增加而上升;层合板的总厚度越大,传声损失也越大,而且各层之间厚度不同,也会引起传声损失的较大改变;层合板铺层角度越大,传声损失也越大。采用传递矩阵法能充分考虑复合材料层合板的铺设方式和铺层角度等因素的影响,利用层合板层间的速度和应力连续边界条件,准确的反应复合材料层合板隔声性能。     

基于铺层设计特征的碳纤维增强复合材料悬架控制臂结构优化

基于铺层设计特征,提出一种使用碳纤维复合材料对承载结构件进行结构优化设计的方法和流程.该方法综合考虑结构几何特征、材料铺层方式、铺层厚度及铺层角度在设计环节中的序列关系,通过几何设计空间构建、离散变量多目标优化、基于工艺可行性的最优决策等方法实现结构设计.以碳纤维增强复合材料悬架控制臂的轻量化设计为例:首先,以钢质控制臂结构为参考建立复合材料控制臂的几何设计空间;然后,以复合材料铺层便利性为原则对其进行结构设计,采用准各向同性铺层对控制臂的铺层厚度进行设计;进而,以提高控制臂刚度和1阶固有频率为目标,使用优化算法对铺层角度进行多目标优化设计;最后,以工艺可行性为约束对优化结果进行筛选并最终完成结构设计.结果表明,所设计复合材料结构具有更大的刚度和1阶固有频率,并且与钢质结构相比减重47.9%.所提出的方法能够较好地兼顾结构特征和复合材料设计要求之间的关系,为复合材料结构优化设计理论与方法的发展提供有益参考.     

面内线性变化载荷作用复合材料层合板的振动与屈曲优化

采用广义微分求积（GDQ）法开展了不同边界条件下承受面内线性变化载荷作用复合材料层合板振动与屈曲的分析与优化。针对GDQ法求解面内线性变化载荷工况复合材料层合板屈曲问题存在计算振荡、不收敛现象,提出载荷扰动策略实现了GDQ法对复合材料层合板屈曲问题的稳定高效求解。基于基础圆频率和临界屈曲载荷系数的归一化指标,分析了铺层角度对复合材料层合板综合性能的影响,并结合直接搜索模拟退火算法开展了复合材料层合板的铺层顺序优化。结果表明:铺层角度变化对屈曲性能的影响明显强于频率特性;面内线性变化载荷中,以弯曲载荷作用下复合材料层合板的优化综合性能受边界条件变化的影响最小,而优化铺层角度受边界条件变化的影响最大。研究结果为复杂载荷作用下复合材料层合板的设计提供了参考。      

复合材料层合板及机翼格栅结构的动态特性优化设计

以复合材料层合板各单层连续变化的铺层角度为设计变量, 在有限元软件中对层合板结构的基频进行优化分析, 在四边简支和固支两种不同的边界条件下, 结构的基频分别提高了4.9%和16.2%, 并对优化前后结构的静力失效强度进行了对比分析。随后将这种优化方法应用到某无人机复合材料机翼格栅结构中, 针对格栅结构蒙皮和肋板共计24个纤维铺层角度进行了优化设计, 使结构基频提高了10.6%, 同时结构的承载能力也有了一定程度的提高。     

Integrated layout design of multi‐component system

A new integrated layout optimization method is proposed here for the design of multi‐component systems. By introducing movable components into the design domain, the components layout and the supporting structural topology are optimized simultaneously. The developed design procedure mainly consists of three parts: (i) Introduction of non‐overlap constraints between components. The finite circle method (FCM) is used to avoid the components overlaps and also overlaps between components and the design domain boundaries. (ii) Layout optimization of the components and supporting structure. Locations and orientations of the components are assumed as geometrical design variables for the optimal placement while topology design variables of the supporting structure are defined by the density points. Meanwhile, embedded meshing techniques are developed to take into account the finite element mesh change caused by the component movements. (iii) Consistent material interpolation scheme between element stiffness and inertial load. The commonly used solid isotropic material with penalization model is improved to avoid the singularity of localized deformation in the presence of design dependent loading when the element stiffness and the involved inertial load are weakened by the element material removal. Finally, to validate the proposed design procedure, a variety of multi‐component system layout design problems are tested and solved on account of inertia loads and gravity center position constraint. Solutions are compared with traditional topology designs without component. Copyright © 2008 John Wiley & Sons, Ltd.     

Automated design of threats and shields under hypervelocity impacts by using successive optimization methodology

In this study, predictive hydrocode simulations are coupled with approximate optimization (AO) methodology to achieve successive design automation for a projectile-Whipple shield (WS) system at hypervelocity impact (HVI) conditions. Successive design methodology is first applied to find the most dangerous threat for a given WS design by varying the shape and orientation of a projectile while imposing constraints on the total projectile mass and radar cross section (RCS). Subsequent optimization procedure is then carried on to improve the baseline WS design parameters. A parametric multi-layered stuffed WS model is considered with varying thicknesses of each layer and variable positions of the inter-layers while having a constraint on the areal density. HVI simulations are conducted by using a non-linear explicit dynamics numerical solver, LS-DYNA. Coupled finite element and smoothed particle hydrodynamics (SPH) parametric models are developed for the predictive numerical simulations. LS-OPT is employed to implement the design optimization process based on response surface methodology. It is found that the ideal spherical projectiles are not necessarily presenting the most dangerous threat compared to the ones with irregular shapes and random orientations, which have the same mass and RCS. Therefore, projectiles with different shapes and orientations should be considered while designing a WS. It is also shown that, successive AO methodology coupled with predictive hydrocode simulations can easily be utilized to enhance WS design.     

考虑结构稳定性的变密度拓扑优化方法

将稳定性问题引入传统变密度法中,可实现包含稳定性约束的平面模型结构拓扑优化。以单元相对密度为设计变量,结构柔度最小为目标函数,结构体积和失稳载荷因子为约束条件建立优化问题数学模型,提出了一种考虑结构稳定性的变密度拓扑优化方法。通过分析结构柔度、体积、失稳载荷因子对设计变量的灵敏度,并基于拉格朗日乘子法和Kuhn-Tucker条件,推导了优化问题的迭代准则。同时,利用基于约束条件的泰勒展开式求解优化准则中的拉格朗日乘子。通过推导平面四节点四边形单元几何刚度矩阵的显式表达式,得到了优化准则中的几何应变能。最后,通过算例对提出的方法进行了验证,并与不考虑稳定性的传统变密度拓扑优化方法进行对比,结果表明该方法能显著提高拓扑优化结果的稳定性。研究结果对细长受压结构的优化设计有重要指导意义,对结构的稳定性设计有一定参考价值。     

Implementation of Discrete Capability into the Enhanced Multipoint Approximation Method for Solving Mixed Integer-Continuous Optimization Problems

Multipoint approximation method (MAM) focuses on the development of metamodels for the objective and constraint functions in solving a mid-range optimization problem within a trust region. To develop an optimization technique applicable to mixed integer-continuous design optimization problems in which the objective and constraint functions are computationally expensive and could be impossible to evaluate at some combinations of design variables, a simple and efficient algorithm, coordinate search, is implemented in the MAM. This discrete optimization capability is examined by the well established benchmark problem and its effectiveness is also evaluated as the discreteness interval for discrete design variables is increased from 0.2 to 1. Furthermore, an application to the optimization of a lattice composite fuselage structure where one of design variables (number of helical ribs) is integer is also presented to demonstrate the efficiency of this capability.     

Solution to the topology optimization problem using a time-evolution equation

This article describes a numerical solution to the topology optimization problem using a time-evolution equation. The design variables of the topology optimization problem are defined as a mathematical scalar function in a given design domain. The scalar function is projected to the normalized density function. The adjoint variable method is used to determine the gradient defined as the ratio of the variation of the objective function or constraint function to the variation of the design variable. The variation of design variables is obtained using the solution of the time-evolution equation in which the source term and Neumann boundary condition are given as a negative gradient. The distribution of design variables yielding an optimal solution is obtained by time integration of the solution of the time-evolution equation. By solving the topology optimization problem using the proposed method, it is shown that the objective function decreases when the constraints are satisfied. Furthermore, we apply the proposed method to the thermal resistance minimization problem under the total volume constraint and the mean compliance minimization problem under the total volume constraint.     

基于遗传算法和拓扑优化的结构多孔洞损伤识别

鉴于拓扑优化和遗传算法在结构损伤识别中各自的优点,本文将遗传算法、有限元和拓扑优化三种方法相结合,提出了一种用于二维结构多损伤识别的新方法。这种方法将拓扑优化的设计变量和遗传算法的参数统一化,将拓扑优化中的目标函数和约束方程与遗传算法的适应度函数联系起来,并以拓扑优化的约束方程作为控制条件参与整个遗传运算的控制。采用二进制编码遗传算法代替连续变量拓扑优化的方式对发生孔洞损伤形式的二维结构进行损伤识别,避免了利用连续变量拓扑优化进行损伤识别时参数阈值的确定可能给识别结果带来的不良影响。通过对两个二维结构模型的多损伤识别仿真计算,结果显示本方法能够很好地识别二维结构中多个位置的损伤,对于仅用拓扑优化法很难识别的轻微孔洞损伤情况,该方法也能得出与实际情况吻合良好的结果。     

宏观结构性能约束下的材料微结构拓扑优化

为了设计周期性多孔钢或钢/铝复合材料优化微结构,基于独立连续映射法,建立了以结构总质量最小化为目标,节点位移为约束的拓扑优化模型。假设宏观结构由多孔材料或复合材料组成,其等效特性采用均匀化理论计算得到。定义了微观材料拓扑变量,节点位移约束采用一阶泰勒展开近似。各种材料设计要求作为约束条件纳入到优化模型中。推导了节点位移和总质量的敏度表达式。采用基于求解偏微分的过滤方法消除了数值不稳定性。在二维数值算例中获得了各种满足设计要求的优化材料微结构。结果表明:提出的方法在材料微结构拓扑优化设计中具有可行性和有效性。      

Developing a fuzzy proportional–derivative controller optimization engine for engineering design optimization problems

In real world engineering design problems, decisions for design modifications are often based on engineering heuristics and knowledge. However, when solving an engineering design optimization problem using a numerical optimization algorithm, the engineering problem is basically viewed as purely mathematical. Design modifications in the iterative optimization process rely on numerical information. Engineering heuristics and knowledge are not utilized at all. In this article, the optimization process is analogous to a closed-loop control system, and a fuzzy proportional–derivative (PD) controller optimization engine is developed for engineering design optimization problems with monotonicity and implicit constraints. Monotonicity between design variables and the objective and constraint functions prevails in engineering design optimization problems. In this research, monotonicity of the design variables and activities of the constraints determined by the theory of monotonicity analysis are modelled in the fuzzy PD controller optimization engine using generic fuzzy rules. The designer only needs to define the initial values and move limits of the design variables to determine the parameters in the fuzzy PD controller optimization engine. In the optimization process using the fuzzy PD controller optimization engine, the function value of each constraint is evaluated once in each iteration. No sensitivity information is required. The fuzzy PD controller optimization engine appears to be robust in the various design examples tested.     

Continuous approximation of material distribution for topology optimization

In this paper, we propose a checkerboard‐free topology optimization method without introducing any additional constraint parameter. This aim is accomplished by the introduction of finite element approximation for continuous material distribution in a fixed design domain. That is, the continuous distribution of microstructures, or equivalently design variables, is realized in the whole design domain in the context of the homogenization design method (HDM), by the discretization with finite element interpolations. By virtue of this continuous FE approximation of design variables, discontinuous distribution like checkerboard patterns disappear without any filtering schemes. We call this proposed method the method of continuous approximation of material distribution (CAMD) to emphasize the continuity imposed on the ‘material field’. Two representative numerical examples are presented to demonstrate the capability and the efficiency of the proposed approach against some classes of numerical instabilities. Copyright © 2004 John Wiley & Sons, Ltd.     

Generalized Timoshenko modelling of composite beam structures: sensitivity analysis and optimal design

This article describes a new approach to design the cross-section layer orientations of composite laminated beam structures. The beams are modelled with realistic cross-sectional geometry and material properties instead of a simplified model. The VABS (the variational asymptotic beam section analysis) methodology is used to compute the cross-sectional model for a generalized Timoshenko model, which was embedded in the finite element solver FEAP. Optimal design is performed with respect to the layers’ orientation. The design sensitivity analysis is analytically formulated and implemented. The direct differentiation method is used to evaluate the response sensitivities with respect to the design variables. Thus, the design sensitivities of the Timoshenko stiffness computed by VABS methodology are imbedded into the modified VABS program and linked to the beam finite element solver. The modified method of feasible directions and sequential quadratic programming algorithms are used to seek the optimal continuous solution of a set of numerical examples. The buckling load associated with the twist–bend instability of cantilever composite beams, which may have several cross-section geometries, is improved in the optimization procedure.