用ICM法拓扑优化静位移及频率约束下连续体结构

用ICM方法建立了静位移及频率约束下、重量最小为目标的连续体结构拓扑优化模型。采用独立于截面及形状参数的连续拓扑变量,借助于过滤函数,位移约束用莫尔定理显式化,频率约束用瑞利商求导数借助模态动能及模态应变能近似显式化。用图形过滤处理的方法解决了棋盘格及网格依赖问题。通过构造适当的过滤函数有效地防止了局部模态问题。动态引入防止模态交换的频率约束条件,使迭代过程不发生振荡。算例表明:用ICM方法建立的模型在处理多工况静位移约束、多频率约束及解决局部模态及模态交换等问题上有优势。     

结构拓扑优化ICM方法的改善

对结构拓扑优化的ICM(独立、连续、映射)方法进行了深入探讨,通 过选取不同的过滤函数可以不进行每步删除而得到清晰的拓扑图形. 以位移约束为例阐述了 ICM方法建模及求解过程. 对位移约束、频率约束、位移及频率约束、简谐载荷激励下动位 移幅值约束等拓扑优化进行了研究,计算算例表明ICM方法在处理静力问题及动力问题的拓 扑优化都是可行的. 程序算法都在MSC.Nastran及MSC.Patran的二次开发环境下实现,与 原软件有机结合在一起.     

有频率禁区的连续体结构拓扑优化

用ICM(独立连续映射)方法建立了以重量为目标、频率为约束的连续体结构拓扑优化模型,并转化为对偶模型利用序列二次规划求解.解决了“棋盘格”现象、网格依赖性、局部模态及模态交换等问题.以重量为目标建立的优化模型同截面及形状优化形成了统一模型,在处理多频率约束及解决模态交换问题上有优势.     

多约束作用下连续体结构的拓扑优化

基于ICM(独立、连续、映射)方法建立了以结构重量最小为目标,以屈曲临界力、位移及应力三种约束同时作用的连续体拓扑优化模型:采用独立的连续拓扑变量,借助泰勒展式、过滤函数将目标函数作二阶近似展开。借助瑞利商、泰勒展式、过滤函数将屈曲约束化为近似显函数,借助于过滤函数,将位移约束用莫尔定理显式化;将应力这种局部性约束采用全局化策略进行处理,即借助第四强度理论、过滤函数将应力局部性约束转化为应变能约束,大大减少了灵敏度分析的计算量;将优化模型转化为对偶规划,减少了设计变量的数目,并利用序列二次规划求解,缩小了模型的求解规模。数值算例表明,ICM方法在解决屈曲、位移及应力三种约束共同作用的连续体拓扑优化问题上有优势。     

具有动应力和动位移约束的离散变量结构形状优化设计方法

讨论了动应力、动位移约束下离散变量状优化设计问题。首先用拟静力算法,将结构惯性力极值作为静载荷施加到结构上,求得结构的动位移和动内力,然后将考虑动应力约束和动作移约束的离散变量结构优化设计问题化为静应力和静位移约束的优化问题。在求解过程中,将单元内力作了一阶近似,并将多约束问题转化为单约束问题,然后利用两类变量统一考虑的离散变量结构形状优化设计的综合算法进行求解。     

结构有限元分析中阶跃型弹性约束的一种有效处理方法

本文通过引入弹性约束刚度矩阵和结构位移约束列阵，提出了结构有限元分析中处理阶跃型弹性约束的一种有效方法。该法通过改变弹性约束系数及位移非约束量大小，可方便有效地处理结点常弹性约束，刚性约束，阶跃型弹性的约束及阶跃型刚性约束等问题。     

基于可行域调整策略的多约束材料优化设计

研究基于可行域调整策略的多约束结构材料优化问题。建立了以结构柔顺度为目标函数,考虑体积、位移、结构自振频率为约束条件的结构材料优化模型和目标函数及约束函数近似式,并采用可行域调整策略和对偶算法求解。典型结构材料优化设计结果表明:可行域调整策略使优化过程平稳,获得的序列细观结构拓扑构型清晰;相比于单一体积约束结构拓扑优化求解,考虑多约束结构拓扑优化能获得更好的黑白占优解。     

自旋挠性航天器非约束模态动力学特性研究

在研究挠性航天器动力学问题时,关注的问题是挠性航天器系统的刚柔耦合作用问题,即航天器挠性附件的振动可能会造成航天器运动失稳。针对中心刚体-双侧大挠性结构的自旋航天器,提出了航天器帆板结构的梁式简化模型,建立了一种非约束模态动力学模型。本研究考虑受到万有引力作用,探讨自旋挠性航天器非约束模态的动力学建模及动态特性。首先利用欧拉方程和哈密顿原理建立了自旋挠性航天器动力学方程,方程解释了刚性模态和弹性模态之间的耦合;然后进行了模态离散化,分别在约束模态和非约束模态下对特征值问题开展研究,对频率和相关振型进行了定量比较;最后进行了数值仿真,求解了自旋挠性航天器非约束模态特征值问题,比较约束模态与非约束模态之间的差异,并用有限元进行验证,得到了随着梁长度的增加,即刚柔惯量比、质量比的减小,非约束模态比约束模态更加准确的结论。     

考虑回转性能的三维结构拓扑优化设计

以考虑运动学性能的回转运动体材料最优分布问题为背景,提出并求解了转动惯量约束下的三维结构拓扑优化问题.常规的结构拓扑优化通常构造为给定材料体积约束下的结构最小柔顺性问题.本文指出,引入结构转动惯量约束对于回转结构的拓扑优化具有重要意义.本文将重量约束下的优化准则法计算格式拓广到单个转动惯量约束下的最小柔顺性问题,而对于包含转动惯量约束的多约束问题采用数学规划法求解.本文给出的数值算例验证了所提出算法的正确性,并且表明引入转动惯量约束后可得到与常规拓扑优化问题不同的结构布局形式.     

模型不确定压电柔性结构的多目标振动控制

以压电柔性结构为对象,考虑其被控模态参数的不确定性和因剔除残余模态所引起的动态不确定性,建立了结构的不确定线性分式模型;根据不确定模型设计了一个对结构进行振动控制的动态输出反馈控制器,使闭环系统满足鲁棒稳定、扰动抑制、极点配置约束及控制输入约束等性能要求;并利用线性矩阵不等式,将具有多个性能要求的振动控制问题转化为具有线性矩阵不等式约束和线性目标函数的凸优化问题。最后以简支压电柔性梁为例,设计了一个融合多性能指标的动态输出反馈非同位控制器,并给出了仿真结果。     

RESEARCH PROGRESS OF FRACTURE TOUGHNESS TEST METHOD FOR LOW CONSTRAINT SPECIMENS 1)

低约束试件断裂韧性测试对油气管道安全运营具有重要意义。本文回顾了低约束试件断裂韧性测试方法及发展过程,介绍了裂纹尖端张开位 移(crack tip opening displacement, CTOD)和J积分等常用断裂韧性表征参数,并对断裂韧性测试中应力强度因子、J积分塑性因子、J积分与CTOD转换因子、裂纹尺寸测量方法、数 字图像相关方法等关键问题进行对比分析,总结需要深入研究的问题,为低约束试件断裂韧性测试发展提供一定依据。     

频率约束板结构拓扑优化

本文利用ICM（独立、连续、映射）方法建立了频率约束下平板结构重量最轻的拓扑优化模型。采用指数函数作为单元重量、单元质量及单元刚度的过滤函数。通过瑞利商对刚度过滤函数倒变量的泰勒一阶展式,将频率约束近似显式化。利用对偶理论将含有大量设计变量的约束优化模型转化为易于求解的少设计变量拟无约束优化模型,通过序列二次规划将转化模型进行求解,提高了求解的效率。本文选择MSC.Patran&Nastran软件及PCL二次开发语言构架了平板结构频率约束拓扑优化问题的软件。数值算例表明：本文的方法具有迭代稳定性和收敛高效性。     

用拓扑优化设计多体分比蜂窝状结构（英文）

提出用拓扑优化及增材制造法设计同一材料制作的蜂窝状结构,并把设计可归结为在不同的设计域内使用不同的体分比的拓扑优化问题。首先推导了多体分比约束条件下的蜂窝状结构拓扑优化问题;然后提出了基于多阈值的等值面法及与商用有限元分析软件对接的算法。给出了有三个体积约束的最小柔顺性问题的数值算例。计算结果表明,对于具有三个体积约束的拓扑优化设计问题,存在一个最佳的体积约束组合;使用等量的材料,与两个体积约束的最佳拓扑设计相比,用最佳体积约束组合所得到的优化结构的柔顺性还要小10%。     

用拓扑优化设计多体分比蜂窝状结构

提出用拓扑优化及增材制造法设计同一材料制作的蜂窝状结构,并把设计可归结为在不同的设计域内使用不同的体分比的拓扑优化问题。首先推导了多体分比约束条件下的蜂窝状结构拓扑优化问题;然后提出了基于多阈值的等值面法及与商用有限元分析软件对接的算法。给出了有三个体积约束的最小柔顺性问题的数值算例。计算结果表明,对于具有三个体积约束的拓扑优化设计问题,存在一个最佳的体积约束组合;使用等量的材料,与两个体积约束的最佳拓扑设计相比,用最佳体积约束组合所得到的优化结构的柔顺性还要小10%。     

遗传算法求解可行域分离的结构优化问题

应用遗传算法求解了两类可行域分离的结构优化问题：局部屈曲约束的桁架拓扑优化问题和动力响应约束优化问题．对第一类问题，提出了新的数学表达式，适合于遗传算法求解．采用了改进的适应度函数及约束处理方法、约束凝聚选择、交叉操作改进和竞争最优保留，提高了遗传算法的效率和可靠性．算例说明，该方法能够克服可行域分离给传统优化算法带来的困难，有效地在多连通可行域中搜索全局最优解．     

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。     

频率约束的三维连续体结构动力拓扑优化设计

针对频率约束和结构重量最小的动力拓扑优化问题, 基于(independent continuous mapping, ICM)独立、连续、映射方法,建立了频率约束下的三维连续体拓扑优化模型. 利用瑞利商和一阶泰勒展式对频率约束进行了显式化处理,并采用幂函数与复合指数函数作为过滤函数,将优化模型进行了标准化转换, 利用对偶理论及数学规划法进行了求解. 另外,利用质量矩阵和刚度矩阵过滤函数比值与动态约束处理了局部模态和模态交换等数值问题. 最后,通过应用两类不同过滤函数的数值算例表明了文中模型及方法在处理动力拓扑优化问题上的合理性与有效性.      

考虑轮胎力耦合约束的智能汽车轨迹跟踪控制算法

针对在低附着变速工况下,忽略纵、侧向轮胎力耦合约束可能导致轨迹跟踪时车辆失稳问题,提出了一种模型预测控制框架下的考虑轮胎力耦合约束的车辆轨迹跟踪控制方法。首先,通过摩擦圆假设建立纵、侧向轮胎力的耦合关系,并推导与之等效的输入量边界约束,将该问题转化为约束二次规划问题;其次,提出了一种基于交叉方向乘子法的数值求解构型,降低了求解约束优化问题时Karush-Kuhn-Tucker方程的维数,实现了求解加速。仿真结果表明,在低附着变速工况下,所提出的算法能够实现最大0.166 m误差的稳定跟踪;同时数值求解过程最多仅需8次迭代,增强了控制过程的实时性。     

应力约束处理为应变能集成的连续体结构拓扑优化

为克服应力约束下拓扑优化问题约束多、应力敏度计算量大的困难,本文提出了将应力约束处理为应变能集成的结构拓扑优化ICM方法。文中用几个典型算例对该方法进行了检验,结果表明:本方法优化效率较高,而且得到的结构最优拓扑较为合理。     

多体系统摩擦接触问题的DAE-LCP方法

当多体系统的约束全部是摩擦接触时, 其动力学问题可归结为一个常微分方程(ordinarydifferential equation, ODE)与线性互补问题(linear complementarityproblem, LCP)的混合动力学问题. 如果除了摩擦接触之外还增加了光滑的双边约束, 则需要将ODE-LCP混合动力学模型推广为微分代数方程(differential algebra equation, DAE)与LCP的混合动力学模型. 该文采用DAE与LCP混合动力学方法求解不考虑碰撞但同时含有持续摩擦接触及光滑等式约束的多体系统动力学问题. 在建立系统动力学模型时,首先将含摩擦的约束从系统中移去得到基本动力学系统. 由于基本系统中带有等式约束, 所以基本系统的动力学方程为一组DAE. 结合基本系统的DAE与约束的互补条件便可以得到DAE-LCP混合动力学模型. 数值计算采用基于DAE与LCP的步进(time-stepping)算法, 将系统动力学方程及其约束离散化并转化为一个混合LCP进行求解. 该算法无需进行滞-滑状态检测, 避免了事件检测导致的繁复计算. 利用所提方法对典型机构的非光滑非线性特征进行了数值分析,验证了该文方法的正确有效性.