基于无单元Galerkin方法的受迫振动下的连续体结构拓扑优化

基于无单元Galerkin法(EFG)对受迫振动下的连续体结构进行拓扑优化设计。选取节点的相对密度作为设计变量,以动柔度最小化为目标函数,基于带惩罚的各向同性固体材料模型（SIMP）建立了受迫振动下的连续体结构拓扑优化模型,采用伴随法求解得到目标函数的敏度分析公式,利用优化准则法对优化模型进行求解。通过经典的二维连续体结构拓扑优化算例证明了该方法的可行性和有效性。     

基于柔顺性与热变形双目标的多孔材料与结构几何多尺度优化设计

论文采用拓扑优化技术针对由具有均一微结构的多孔陶瓷材料构成的结构,从结构和材料两个方面进行优化设计,在轻量化的同时使得结构具有较高的刚度和特定区域较小的热变形.以机械载荷下的结构柔顺性和热载荷下结构特定区域的热变形为目标,以宏观密度和微观密度为两类独立的设计变量,对两个目标进行归一化处理,建立多目标优化模型.采用体积守...     

基于热固耦合问题的多尺度拓扑优化设计

基于均质材料的拓扑优化逐渐难以适应现代化生产对工业产品高品质、轻量化的需求,同时很多工业设备经常面临着高温和高负载的工作环境.为了提高结构在温度载荷和机械载荷共同作用下的力学性能,本文提出了一种在稳态热源作用下的热固耦合连续体结构并行化拓扑优化方法.以结构刚度作为目标函数,材料的体分比为约束,利用能量均匀化方法预测微结...     

考虑自重载荷作用的连续体结构拓扑优化

针对自重载荷作用下连续体结构拓扑优化中常遇到的柔顺度随单元密度变化的非单调性、非有效的材料体积约束和材料低密度区出现的``附属'效应等主要困难,提出了RAMP材料插值模型和平均敏度过滤技术相结合的求解策略,并采用MMA方法进行了优化求解. 研究了惩罚因子的选择对最优拓扑以及材料体积约束的影响. 结果表明,合理增大RAMP材料模型中惩罚因子(例如取值达到20.0), 可使材料体积约束达到有效约束. 此外,通过与Sigmund提出的敏度过滤技术比较表明,当引入平均敏度过滤技术时,结果0/1化程度较Sigmund的敏度技术高,可获得清晰的黑白拓扑.      

连续体结构拓扑优化的高精度逼近ICM方法

采用指数类函数为快滤函数的高精度逼近ICM (independent continuous and mapping)方法, 建立了以结构重量为目 标, 应力和位移共同约束下的连续体结构拓扑优化模型. 利用结构畸变比能的方法全局化应 力约束, 单位虚载荷法显式化位移约束, 归一化约束以解决约束限数量级不一致的问题. 针 对不同性态的过滤函数, 给出了指数类快滤函数参数的取值方法. 单工况和多工况的算例表 明了高精度逼近的ICM方法处理多种约束下连续体结构拓扑优化的可行性与有效性.     

基于类桁架连续体的结构拓扑优化方法与应用

以各向异性连续体为基结构,采用类桁架连续体材料模型进行结构拓扑优化。以材料在结点位置的密度和方向作为优化设计变量,使材料在设计域内连续分布。并以此建立材料的弹性矩阵和刚度矩阵。优化过程没有抑制中间密度,这从根本上避免了许多拓扑优化方法普遍存在的单元铰接、棋盘格现象以及单元依赖性等数值不稳定问题。采用满应力准则法,借助有限元结构分析,经过少量迭代,建立优化的材料连续分布场,即类桁架连续体结构。由于首先建立的拓扑优化结构是各向异性连续体,从而得到更大优化空间。然后可以结合工程实际需要将其转化为离散的拓扑优化杆系结构。最后,以1个经典Michell桁架和3种形式的拱桥为数值算例,演示了其结构拓扑优化过程。     

多工况下连续体结构的多刚度拓扑优化设计和二重敏度过滤技术

针对多工况载荷条件下连续体结构拓扑优化设计的研究现状,以各载荷下对应结构的最小柔度作为目标函数,基于带权重的折衷规划法和SIMP密度函数插值模型,提出并建立了多载荷条件下线弹性连续体结构拓扑优化设计的数学模型．对优化过程中出现的棋盘格式和网格依赖性等数值计算问题进行了研究,在此基础上提出一种二重敏度过滤技术,该方法能有效消除优化结构中的棋盘格式,并使优化结构体现出较好的网格无关性．通过几个典型的算例证明了该文方法的有效性．     

模糊-区间参数连续体结构的拓扑优化设计

研究位移可靠性约束下的模糊-区间参数连续体结构的拓扑优化问题.利用3σ准则和信息熵分别将区间变量和模糊变量近似转换为随机变量,构建了随机参数平面连续体结构的拓扑优化模型,以结构的形状拓扑信息为设计变量,结构总质量均值极小化为目标函数.以满足位移可靠性为约束条件;利用代数综合法导出了位移响应的数字特征的计算表达式;通过单位虚载荷求得位移灵敏度然后采用双方向渐进结构优化方法求解.通过两个算例验证了文中模型的合理性和求解策略的有效性.     

各向异性薄板的无网格法结构拓扑优化研究

将无网格Galerkin法引入正交各向异性薄板的结构拓扑优化中,建立了以节点相对密度为设计变量、以结构柔度为目标函数的结构拓扑优化模型;采用罚函数施加本质边界条件,结合固体各向同性惩罚插值模型和OC优化准则,推导了目标函数的灵敏度分析算法;利用数值算例验证了所建模型及算法的可行性,完成了单载荷工况、多载荷工况下各向异性材料的拓扑优化设计,探讨了材料性能与铺设角度对各向异性薄板结构拓扑优化结果的影响。结果表明,各向异性薄板在弹性模量较小的方向上,材料分布较多,且拓扑结构呈现周期性变化。     

多工况应力和位移约束下连续体结构拓扑优化

将文「１」所提出的对拓扑变量的独立连续映射（ＩＧＭ）的拓扑估化方法应用于连续体结构,从而建立了统一的以重量为目标,考虑应力和位移约束的连续体结构拓扑优化模型。通过对位移－应力拓扑解和各工况下应力拓扑解的综合协调,进而对于协调拓扑解按照阈值完成从离散到连续的反演,并且采用分层与加权策略克服了“荷载病态”困难。给出的经典的二维平面问题和三维连续体结构拓扑优化算例表明,这种统一的模型由骨架结构发展到连续     

Optimal internal architectures of femoral bone based on relaxation by homogenization and isotropic material design

The influence of the loading conditions on the trabecular architecture of a femur is investigated by using topology optimization methods. The response of the bone to physiological loads results in changes of the internal architecture of bone, reflected by a modification of internal effective density and mechanical properties. The homogenization based optimization model is developed for predicting optimal bone density distribution, wherein bone tissue is assumed to be a composite material consisting of a mixture of material and void. The homogenization scheme treats the geometric parameters of the microstructures and their orientation as design variables and homogenizes the properties in that microstructure, which is generally anisotropic. The penalization of the optimal material density then leads to a classical optimal structure which consists of regions with bone material and regions without bone material. The IMD (Isotropic Material Design) approach is next applied to determine the optimal elasticity tensor in terms of the bulk and shear moduli for the present loading applied to the femoral bone sample. IMD is able to provide both the external shape and topology together with the optimal layout of the isotropic moduli. Both topology optimization methods appear to be complementary. Simulations of the internal bone architecture of the human proximal femur results in a density distribution pattern with good consistency with that of the real bone.     

蜂窝夹芯圆环的拓扑优化设计及尺度效应研究

采用尺度关联的一体化设计方法开展了旋转周期圆环结构的拓扑优化设计研究,以宏观 结构的最大刚度为目标,研究了材料表征体胞尺度、构型以及不同载荷作用形式对蜂窝夹芯 圆环结构优化结果的影响. 所提出的无量纲结构构型因子实现了优化结构的结构效率量化评 估. 结合SIMP材料模型和周长控制方法,实现了宏观结构和细观表征体胞的优化设计,获得 清晰的材料分布. 数值算例表明,尺度关联的一体化设计方法能有效地完成圆环结构的拓扑 优化设计,设计结果充分反映体胞尺度效应对旋转周期圆环结构夹芯构型的影响.     

Principal Compliance and Robust Optimal Design

The paper addresses a problem of robust optimal design of elastic structures when the loading is unknown and only an integral constraint for the loading is given. We propose to minimize the principal compliance of the domain equal to the maximum of the stored energy over all admissible loadings. The principal compliance is the maximal compliance under the extreme, worst possible loading. The robust optimal design is formulated as a min-max problem for the energy stored in the structure. The maximum of the energy is chosen over the constrained class of loadings, while the minimum is taken over the design parameters. It is shown that the problem for the extreme loading can be reduced to an elasticity problem with mixed nonlinear boundary conditions; the last problem may have multiple solutions. The optimization with respect to the designed structure takes into account the possible multiplicity of extreme loadings and divides resources (reinforced material) to equally resist all of them. Continuous change of the loading constraint causes bifurcation of the solution of the optimization problem. It is shown that an invariance of the constraints under a symmetry transformation leads to a symmetry of the optimal design. Examples of optimal design are investigated; symmetries and bifurcations of the solutions are revealed.     

CRACK TIP PLASTICITY OF A THERMALLY LOADED PENNY-SHAPED CRACK IN AN INFINITE SPACE OF 1D QC

The present work is concerned with a penny-shaped Dugdale crack embedded in an infinite space of one-dimensional(1D) hexagonal quasicrystals and subjected to two identical axisymmetric temperature loadings on the upper and lower crack surfaces. Applying Dugdale hypothesis to thermo-elastic results, the extent of the plastic zone at the crack tip is determined.The normal stress outside the plastic zone and crack surface displacement are derived in terms of special functions. For a uniform loading case, the corresponding results are presented by simplifying the preceding results. Numerical calculations are carried out to show the influence of some parameters.     

A MODIFIED INTERPOLATION APPROACH FOR TOPOLOGY OPTIMIZATION

In view of the fact that the follow-up search for an optimal topology is affected by deleting a large number of high-relative-density elements. When the typical density interpolation approach, namely, solid isotropic microstructures with penalization(SIMP), is employed in the continuum structural topology optimization, a new density interpolation approach based on the logistic function is proposed in this paper. This method can weaken low-relative-density elements while enhancing high-relative-density elements by polarization, and then rationally realize polarization of the intermediate density elements. It can reduce the number of gray-scale elements as much as possible to get the optimal topology with distinct boundaries in conjunction with the sensitivity filtering method based on particle swarm optimization(PSO). Several typical numerical examples are given to demonstrate this method.     

传热结构的多目标拓扑优化设计研究

深入分析了传热结构多目标拓扑优化设计中的几个关键问题。提出了基于结构柔度最小化和结构散热弱度最小化的多目标拓扑优化设计方法,建立了传热结构的多目标拓扑优化设计模型,推导了传热结构多目标拓扑优化中用于迭代分析求解的优化准则算法和敏度分析方程。通过数值计算验证了理论和算法的有效性。     

Thermo-mechanical measurement of elasto-plastic transitions during cyclic loading

The surface temperature of stainless steel SS304 low cycle fatigue specimens subjected to cyclic loading was studied using infrared thermography technique. The thermal data mapped onto the various stages of cyclic stress-strain curve shows the ability of these measurements to identify the yield points in both the compression and tension loading. Based on the results of this study, it is possible to identify the state of stress for materials such as elastic tension, plastic tension, elastic compression, plastic compression during cyclic loading using infrared thermographic data. The thermo-elastic slope and thermo-plastic slope was observed to be dependent on the prior loading cycles.     

级联碰撞对层状珠光体力学行为的影响

铁素体合金钢是目前在核能工程界应用最为广泛的一种金属结构材料,以渗碳体和铁素体基体构成的层状珠光体是铁素体合金钢中常见的金相结构。深入理解辐照效应对层状珠光体力学性能的影响对高辐照条件下铁素体钢的材料设计与寿命评估有着重要的理论参考意义。基于以上考虑,本文采用分子动力学（MD）模拟,研究了连续低能铁原子级联碰撞对渗碳体/铁素体两相界面的破坏情况,探讨了经历不同程度级联碰撞的两相结构在单向拉伸以及压缩荷载下的初始屈服情况。通过对MD模拟结果的深入分析,得到了以下主要结论：a.辐照会破坏渗碳体/铁素体两相界面的失配位错结构,引起渗碳体的分解,并促进碳原子向铁素体的扩散;b.在单轴拉伸荷载作用下,级联碰撞会使初始屈服机制由{112}<111>位错滑移系的开动转变为间隙原子团簇附近位错环的形核与长大;c.在单轴压缩荷载作用下,级联碰撞会使初始塑性变形机制由{110}<111>滑移系的开动转变为{112}<111>滑移系的开动;d)无论在单轴拉伸还是压缩情况下,级联碰撞（及辐照效应）都会导致位错初始形核应力的提升。本文的研究结果为铁素体合金钢的辐照硬化和辐照脆化行为提供了新的微观解释,对于辐照条件下铁素体合金钢材料的优化设计有一定的参考意义。     

基于NURBS插值的三维渐进结构优化方法

提出一种基于等几何控制点密度变量的三维双向渐进结构拓扑优化方法。在当前列式下,高阶NURBS基函数被同时用于CAD模型中NURBS实体片的几何场、位移场和温度场以及密度场插值,实现了几何模型、分析模型和优化模型的有效统一,确保了位移场、温度场及密度场的高阶连续性;详细推导了基于等几何控制点密度变量的三维渐进结构法模型及其灵敏度分析列式;最后几个典型的数值算例,包括最小柔顺性、热传导优化问题及三维结构自由振动的基频最大化问题,验证了本文方法的有效性。     

级联碰撞对层状珠光体力学行为的影响

铁素体合金钢是目前在核能工程界应用最为广泛的一种金属结构材料,以渗碳体和铁素体基体构成的层状珠光体是铁素体合金钢中常见的金相结构。深入理解辐照效应对层状珠光体力学性能的影响对高辐照条件下铁素体钢的材料设计与寿命评估有着重要的理论参考意义。基于以上考虑,本文采用分子动力学（MD）模拟,研究了连续低能铁原子级联碰撞对渗碳体/铁素体两相界面的破坏情况,探讨了经历不同程度级联碰撞的两相结构在单向拉伸以及压缩荷载下的初始屈服情况。通过对MD模拟结果的深入分析,得到了以下主要结论：a.辐照会破坏渗碳体/铁素体两相界面的失配位错结构,引起渗碳体的分解,并促进碳原子向铁素体的扩散;b.在单轴拉伸荷载作用下,级联碰撞会使初始屈服机制由{112}<111>位错滑移系的开动转变为间隙原子团簇附近位错环的形核与长大;c.在单轴压缩荷载作用下,级联碰撞会使初始塑性变形机制由{110}<111>滑移系的开动转变为{112}<111>滑移系的开动;d)无论在单轴拉伸还是压缩情况下,级联碰撞（及辐照效应）都会导致位错初始形核应力的提升。本文的研究结果为铁素体合金钢的辐照硬化和辐照脆化行为提供了新的微观解释,对于辐照条件下铁素体合金钢材料的优化设计有一定的参考意义。