针对柔性机构的稳定性要求,基于Heaviside三场映射方案和柔顺度变化率约束措施,研究了屈曲约束下无铰链多输入多输出柔性机构柔度最小化的拓扑优化问题。首先,为解决低密度单元引起的伪屈曲模态和计算效率等问题,构建了单元刚度矩阵和应力矩阵的光滑惩罚函数,提出了解决低阶伪屈曲模态问题的综合措施。其次,结合Pian混合应力单元公式、凝聚函数法、Heaviside三场映射方案和柔顺度变化率约束措施,构建了考虑屈曲的无铰链多输入多输出柔性机构拓扑优化模型。而后,推导了目标函数和屈曲约束的灵敏度,并利用MMA(Method of Moving Asymptotes)算法进行优化求解。最后,给出的数值算例,说明了此方法的可行性和有效性。
Minimum length scale control on real and void material phases in topology optimization is an important topic of research with direct implications on numerical stability and solution manufacturability.And it also is a challenge area of research due to serious conflicts of both the solid and the void phase element densities in phase mixing domains of the topologies obtained by existing methods.Moreover,there is few work dealing with controlling distinct minimum feature length scales of real and void phase materials used in topology designs.A new method for solving the minimum length scale controlling problem of real and void material phases,is proposed.Firstly,we introduce two sets of coordinating design variable filters for these two material phases,and two distinct smooth Heaviside projection functions to destroy the serious conflicts in the existing methods(e.g.Guest Comput Methods Appl Mech Eng 199(14):123-135,2009).Then,by introducing an adaptive weighted 2-norm aggregation constraint function,we construct a coordinating topology optimization model to ensure distinct minimum length scale controls of real and void phase materials for the minimum compliance problem.By adopting a varied volume constraint limit scheme,this coordinating topology optimization model is transferred into a series of coordinating topology optimization sub-models so that the structural topology configuration can stably and smoothly changes during an optimization process.The structural topology optimization sub-models are solved by the method of moving asymptotes(MMA).Then,the proposed method is extended to the compliant mechanism design problem.Numerical examples are given to demonstrate that the proposed method is effective and can obtain a good 0/1 distribution final topology.
A methodology for topology optimization based on element independent nodal density(EIND) is developed.Nodal densities are implemented as the design variables and interpolated onto element space to determine the density of any point with Shepard interpolation function.The influence of the diameter of interpolation is discussed which shows good robustness.The new approach is demonstrated on the minimum volume problem subjected to a displacement constraint.The rational approximation for material properties(RAMP) method and a dual programming optimization algorithm are used to penalize the intermediate density point to achieve nearly 0-1 solutions.Solutions are shown to meet stability,mesh dependence or non-checkerboard patterns of topology optimization without additional constraints.Finally,the computational efficiency is greatly improved by multithread parallel computing with OpenMP.
