平面膜结构拓扑优化的有无复合体方法

将作者对桁架在应力约束下结构拓扑优化的有无复合体模型发展到平面膜结构在应力、位移约束下结构拓扑优化的建模与求解。同时提出了该模型的有效解法,获得了令人满意的数值结果。本文工作表明独立连续拓扑变量的提出对于结构拓扑优化的研究是有意义的。     

位移约束集成化处理的连续体结构拓扑优化

为解决多工况下多位移约束的连续体结构拓扑优化问题,引入了K-S函数对位移约束进行集成化处理.在建立优化模型时,基于莫尔定理按ICM方法导出约束点位移与设计变量之间的近似显函数关系,然后采用Lagrange乘子法进行求解.给出三个算例对该方法进行验证,并与ESO法、BESO法、MD法以及均匀化方法的结果进行比较.结果表明:该方法计算效率较高,并且能够计算出更合理的结构拓扑.     

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)     

非齐次动力方程Duhamel项的精细积分

提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性.     

包含两类变量的离散变量桁架结构拓扑优化设计

建立了包含截面和拓扑两类变量的离散变量结构拓扑优化设计的数学模型,该模型考虑了截面变量与拓扑变量间的耦合关系,反映了拓扑优化问题的组合优化本质,可以较好地解决"极限应力"、"最优解的奇异性"等困扰结构拓外优化设计的问题．同时采用相对差商法进行离散变量桁架结构拓扑优化,直接求解包含两类变量的离散变量结构拓扑优化设计数学模型,收到了比较满意的效果．     

含三次耦合项的两自由度Duffing系统的共振及混沌行为

研究了一类含三次耦合项的两自由度Duffing系统的动力学行为。首先应用多尺度方法近似求解系统的一阶稳态响应。通过讨论系统的主共振和1∶1内共振,分析了三次耦合项对系统响应的影响。随后研究系统随外加周期力强度变化的分岔过程,发现除了常见的倍周期分岔通向混沌外,还存在一种直接由周期运动进入混沌的突发路径。结合对系统的最大Lyapunov指数,相轨图及Poincar啨映射的分析验证了上述结论。     

连续体结构拓扑优化应力约束凝聚化的ICM方法

为克服应力约束下拓扑优化问题约束数目多、应力敏度计算量大的困难，提出 了应力约束化凝聚化的ICM方法. 在利用Mises强度理论将应力约束转换成应变能约束后， 提出了应力约束凝聚化的两条途径：其一为应力全局化的方法，其二为应力约束集成化的方 法. 由此建立了多工况下以重量为目标、以凝聚化应变能为约束的连续体结构优化模型，并 利用对偶理论对优化模型进行了求解. 4个数值算例表明：该方法具有较高的计算效率，得 到的拓扑结构比较合理，不仅适用于二维连续体结构，也适用于三维连续体结构.     

类升力体外形高超声速粘性绕流拓扑结构分析

采用交替方向隐式分解的隐式NND格式求解全N-S方程模拟“类升力体”外形在高超声速下的大攻角流动,给出了“类升力体”外形表面极限流线随攻角变化的拓扑结构及40°攻角下垂直于体轴的横截面流线拓扑结构。结果表明:类升力体外形三维流场结构十分复杂,攻角从0°～40°变化时,背风面表面极限流线依次由不分离、开式分离向起始于鞍、结点组合的高阶奇点的分离方式转化,翼面横向分离亦随攻角增大而增大;垂直于体轴的横截面流动拓扑结构与文献[2]给出的理论分析一致。     

多工况多约束下离散变量桁架结构的拓扑优化设计

提出了一个多工况下受应力、位移约束的离散变量桁架结构的拓扑优化方法，给出了结构拓扑形式变更时的约束处理方法及杆件删除策略，使基结构设计空间的维数不断降低，达到最优拓扑，避免了奇异解的出现，在算例中指出了截面离散集和位移约束对最优拓扑的影响；算例给出了满意的拓扑优化解。     

CONTINUUM TOPOLOGY OPTIMIZATION FOR MONOLITHIC COMPLIANT MECHANISMS OF MICRO-ACTUATORS

A multi-objective scheme for structural topology optimization of distributed compliant mechanisms of micro-actuators in MEMS condition is presented in this work, in which mechanical flexibility and structural stiffness are both considered as objective functions. The compliant micro-mechanism developed in this way can not only provide sufficient output work but also have sufficient rigidity to resist reaction forces and maintain its shape when holding the work-piece. A density filtering approach is also proposed to eliminate numerical instabilities such as checkerboards, mesh-dependency and one-node connected hinges occurring in resulting mechanisms. SIMP is used as the interpolation scheme to indicate the dependence of material modulus on element-regularized densities. The sequential convex programming method, such as the method of moving asymptotes (MMA), is used to solve the optimization problem. The validation of the presented methodologies is demonstrated by a typical numerical example.