• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 21
    • 下载费用:30 金币  

    6月22日重庆时时彩: 改进双向渐进法的等效静载荷法动态响应拓扑优化方法.pdf

    关 键 词:
    改进 双向 渐进法 等效 载荷 动态 响应 拓扑 优化 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610807923.7

    申请日:

    2016.09.08

    公开号:

    CN106372347A

    公开日:

    2017.02.01

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20160908|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 厦门大学嘉庚学院
    发明人: 贾檀; 蓝萌; 吴晓明
    地址: 363105 福建省漳州市招商局经济技术开发区厦门大学漳州校区
    优先权:
    专利代理机构: 福州元创专利商标代理有限公司 35100 代理人: 蔡学俊
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610807923.7

    授权公告号:

    |||

    法律状态公告日:

    2017.03.01|||2017.02.01

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明涉及一种改进双向渐进法的等效静载荷法动态响应拓扑优化方法。该方法将改进的BESO方法应用于等效静载荷法中作为优化迭代方法,替换掉原始的变密度法,并建立新的收敛条件,从而改进动态响应拓扑优化流程。本发明提出的方法,能有效减少k大于1次设计循环阶段的优化时间,两个算例中优化所需CPU总时间比原始方法分别减少了21.8%、25.3%;该方法具有计算速度快、优化流程简单的特点;随着外载荷作用点的增多、优化规模的增大,优化所需的设计循环次数也会增加,该方法将会带来更大的优势。

    权利要求书

    1.一种改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其特征在于:将改进
    的BESO方法应用于等效静载荷法中作为优化迭代方法,替换掉原始的变密度法,并建立新
    的收敛条件,从而改进动态响应拓扑优化流程。
    2.根据权利要求1所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:所述改进的BESO方法具体如下,
    在对单元处理之前,应事先给出每次迭代的单步目标体积Vk+1;初始体积Vk大于或小于
    整体的体积约束V*;因此,体积进化迭代式可表示为:
    Vk+1=Vk(1±ER),(k=1,2,3...) 式(1)
    式(1)中k为当前迭代次数,ER为体积进化率;
    若目标体积Vk+1小于当前体积Vk,则必须将单元按灵敏度大小排列并删除灵敏度低的单
    元,在删除单元之前,空单元的灵敏度有可能大于待删单元灵敏度,因此要先考虑增添单
    元,这里引入增删阀值αth,αth定义为待删单元灵敏度中的最大值,有:
    αi>αth 式(2)
    式(2)中αi为增添单元的灵敏度;当增添单元的体积大于体积增添率AR时,需要重新计
    算αth直到满足增添单元的体积小于或等于AR为止;然后删除灵敏度小于该值的单元;显然,
    删除单元的体积Vdel可以表示为:
    Vdel=Vk-Vk+1+Vadd 式(3)
    式(3)中Vadd为增添单元的体积;
    若目标体积Vk+1大于当前体积Vk,则必须将单元按灵敏度大小排列并增添灵敏度高的单
    元,在增添单元之前,实单元的灵敏度有可能小于待增单元灵敏度,因此要先考虑删除单
    元,这里引入增删阀值βth,βth定义为待增单元灵敏度中的最小值,有:
    βi<βth 式(4)
    其中βi为删除单元的灵敏度;当删除单元的体积大于体积增添率AR时,需要重新计算
    βth直到满足删除单元的体积小于或等于AR为止;然后增添灵敏度大于该值的单元;显然,增
    添单元的体积Vadd可以表示为:
    Vadd=Vk-Vk+1+Vdel 式(5)
    重复进行有限元分析与单元增删的过程,直至迭代循环满足两个终止条件:一是体积
    约束,即当前进化体积已满足整体约束体积V*;二是收敛准则,考虑目标函数的变化量,收
    敛判别式规定为:

    式(6)中k为当前迭代步数,τ为内部许用收敛因子;M是一整数,用以限制并稳定柔顺度
    在多次连续迭代中的平均变化量;计入局部多个M值,能够使目标函数波动较小,曲线更为
    平滑,有效防止局部收敛。
    3.根据权利要求2所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:所述M取值为5。
    4.根据权利要求2所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:以最小化柔顺度为目标函数,以体积为约束条件,以单元本身为设计变量,基于
    灵敏度分析改进的BESO刚度优化数学表达式可描述为:

    式(7)中F和U分别为外载荷向量和结点位移向量,C为结构柔顺度,Vi是个体单元的体
    积,V*为体积约束,N为设计域中单元总数量;xi为双向赋值的非连续设计变量,即材料单元
    的密度;K为结构总体刚度矩阵,xmin为材料单元的密度的一个小值。
    5.根据权利要求4所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:所述改进的BESO方法中,删除单元的方法为软杀法,即将灵敏度低的待删单元的
    密度人为减至一个小值xmin,以防止刚度矩阵奇异,可令xmin=0.001,但后续仍参与有限元
    分析,而不是直接删除;因此,采用单元软杀法的BESO方法的单元灵敏度计算式可归结为:

    式(8)中Ki0表示实单元的单元刚度矩阵,ui为单元节点位移向量。
    6.根据权利要求5所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:将权值法引入到BESO法中,多载荷工况的刚度拓扑优化问题可以描述为:

    式(9)中,m是载荷工况的总数,ωs、Cs分别是工况s的预设权重系数和平均柔顺度,并有

    由于各载荷工况的位移场相互独立,带权重的目标函数设计变量xi求偏导数可得:

    于是灵敏度计算式可以归结为:

    7.根据权利要求6所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:在等效静态载荷法中,等效静态载荷被定义为能够产生与动态载荷相同位移场
    的一组载荷集;可以看出,第i个等效静载荷作用下产生的位移响应与第i个时间步上的动
    态响应完全相等;
    已知动态分析控制方程:

    式(12)中,M、C、K分别为结构质量、阻尼和刚度矩阵,均为设计变量x的函数;
    和u(t)分别为加速度、速度和位移向量,f(t)为外部施加的动载荷;
    对式(12)式进行变换有:

    式(13)中,t=1,2,3…m;
    再用静态变量s替换feq(t)中的时间变量t,即可算出等效静载荷:
    feq(s)=K(x)u(s),s=1,2,3...m 式(14)。
    8.根据权利要求7所述的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,其
    特征在于:该方法具体实现如下,
    S1:进行有限元建模,设定初始设计变量x(k)=x(0),初始设计循环次数k=0;
    S2:对以x(k)作为密度的有限元模型进行动力学分析,求其时间域上的节点位移响应u
    (t);
    S3:将时间域上的连续响应u(t)离散为m个时间点上的响应u(s),并计算m个时刻的等
    效静载荷feq(s),计算公式如下:
    feq(s)=K(x)u(s);s=1,2,3...m 式(15)
    S4:将等效静载荷输入静态多目标拓扑优化数学模型中,每一个时刻的等效载荷为一
    种工况,每个工况赋予ωk=1/l的权重因子,进行优化迭代直至收敛;

    S5:当k=0时,令k=k+1并进入步骤S2;当k≠0,如果满足收敛条件,优化结束;如果不
    满足收敛条件,依然令k=k+1并进入步骤S2;至此完成一次设计循环;收敛条件可以表示
    为:

    Convergent if||con||≤n×ε2
    该方法将渐进结构法替换变密度法作为内循环方法,第一次外循环时,初始设计域为
    满设计域,此时按照Vk+1=Vk(1±ER),其中k=1,2,3…,计算每次迭代的目标体积,从第二
    次外循环开始,均满足Vk+1=V*,即设计域同时包含0、1单元,此时有Vdel=Vadd=AR,这种情
    况下,渐进结构法的优化速度要比满设计域时快。

    说明书

    改进双向渐进法的等效静载荷法动态响应拓扑优化方法

    技术领域

    本发明属于结构动态响应拓扑优化领域,具体涉及一种改进双向渐进法的等效静
    载荷法动态响应拓扑优化方法。

    背景技术

    拓扑优化是指在满足给定约束下寻求材料的最优分布形式[1]。常见拓扑优化方法
    有变密度法、渐进结构法以及均匀化法。Bend-soe M P和Kikuchi N[2]于1988年提出了均匀
    化法,Bend-soe M P[3]于1989年提出了变密度法。渐进结构优化法(ESO)是近年来兴起的一
    种解决各类结构优化问题的数值方法,由Steven等人[4]于1993年提出,该方法通过删除结
    构中无效或低效的材料逐步去掉使剩下的结构渐渐趋于优化。2000年Querin等[5]采用双向
    渐进法(BESO)利用Von Mises应力准则对“满应力设计”进行了研究,BESO是对ESO的改进和
    延伸,优化过程中可以删除和增添材料。无论是ESO还是BESO都采用两个独立的参数(单元
    删除率RR和单元添加率IR)来控制所删除的和添加的单元数量。但是Rozvany等[6]于2009年
    指出,采用两个独立的参数RR和IR是不合理的,如果没有选择好这两个参数将得不到最优
    的结果。鉴于此,Huang和Xie等人[7]针对刚度优化,提出一种基于灵敏度分析的改进BESO
    法,该方法解决了连续体拓扑优化的很多问题,如优化问题的恰当描述、棋盘格现象、网格
    依赖性以及解的收敛性,使其优化结果更加趋于稳定。PARK J Y等人[8-9]先后将改进的BESO
    法用于汽车发动机机罩的多目标优化,取得了理想的结果。

    目前拓扑优化多用于静态载荷优化,而工程实际中结构受到的一般是动态载荷,
    直接进行动态响应优化分析[10-11]虽然在理论上可行,但是计算很难收敛,同时计算过程繁
    冗,耗时长,未能在工程中得到应用。针对这些问题,Park于2002年提出了等效静载荷法[12]
    (ESLM),将动态载荷转化为等效静态载荷,并作为边界条件对结构进行动态优化,该方法最
    初只是用于尺寸优化。在随后的研究中,PARK证明了该方法获得的结果满足Karush-Kuhn-
    Tucker必要条件[13],并通过相应的算例验证了该方法对线性结构和非线性结构进行拓扑
    优化的可行性[14-15]。该方法在国内亦得到广泛应用[16-20]。相较于传统动力学优化法,等效
    静载荷方法在计算效率上已经得到很大提升,但是由于优化时需要若干个设计循环才能得
    到稳定结果,且每次设计循环都要进行动力学分析以及静态拓扑优化迭代,依然会耗费大
    量时间。

    为了进一步提升等效静载荷法的计算效率,Kim E等人[21]将系统约简法与等效静
    载荷法结合,能够更有效地求解大规模动力学优化问题。Lee H A等人[15]针对非线性结构
    提出一种结点近似位移法来减少动态响应分析阶段所花费的时间,以此提高整体优化效
    率。

    发明内容

    本发明的目的在于解决当前等效静载荷法中单次设计循环耗时长,优化流程复杂
    等问题,提供一种改进双向渐进法的等效静载荷法动态响应拓扑优化方法,该方法将改进
    的BESO方法应用到等效静载荷法中作为优化迭代方法,替换掉原始的变密度法,从而提升
    计算效率,简化优化流程。

    为实现上述目的,本发明的技术方案是:一种改进双向渐进法的等效静载荷法动
    态响应拓扑优化方法,将改进的BESO方法应用于等效静载荷法中作为优化迭代方法,替换
    掉原始的变密度法,并建立新的收敛条件,从而改进动态响应拓扑优化流程。

    在本发明一实施例中,所述改进的BESO方法具体如下,

    在对单元处理之前,应事先给出每次迭代的单步目标体积Vk+1;初始体积Vk大于或
    小于整体的体积约束V*;因此,体积进化迭代式可表示为:

    Vk+1=Vk(1±ER),(k=1,2,3...) 式(1)

    式(1)中k为当前迭代次数,ER为体积进化率;

    若目标体积Vk+1小于当前体积Vk,则必须将单元按灵敏度大小排列并删除灵敏度
    低的单元,在删除单元之前,空单元的灵敏度有可能大于待删单元灵敏度,因此要先考虑增
    添单元,这里引入增删阀值αth,αth定义为待删单元灵敏度中的最大值,有:

    αi>αth 式(2)

    式(2)中αi为增添单元的灵敏度;当增添单元的体积大于体积增添率AR时,需要重
    新计算αth直到满足增添单元的体积小于或等于AR为止;然后删除灵敏度小于该值的单元;
    显然,删除单元的体积Vdel可以表示为:

    Vdel=Vk-Vk+1+Vadd 式(3)

    式(3)中Vadd为增添单元的体积;

    若目标体积Vk+1大于当前体积Vk,则必须将单元按灵敏度大小排列并增添灵敏度
    高的单元,在增添单元之前,实单元的灵敏度有可能小于待增单元灵敏度,因此要先考虑删
    除单元,这里引入增删阀值βth,βth定义为待增单元灵敏度中的最小值,有:

    βi<βth 式(4)

    其中βi为删除单元的灵敏度;当删除单元的体积大于体积增添率AR时,需要重新
    计算βth直到满足删除单元的体积小于或等于AR为止;然后增添灵敏度大于该值的单元;显
    然,增添单元的体积Vadd可以表示为:

    Vadd=Vk-Vk+1+Vdel 式(5)

    重复进行有限元分析与单元增删的过程,直至迭代循环满足两个终止条件:一是
    体积约束,即当前进化体积已满足整体约束体积V*;二是收敛准则,考虑目标函数的变化
    量,收敛判别式规定为:


    式(6)中k为当前迭代步数,τ为内部许用收敛因子;M是一整数,用以限制并稳定柔
    顺度在多次连续迭代中的平均变化量;计入局部多个M值,能够使目标函数波动较小,曲线
    更为平滑,有效防止局部收敛。

    在本发明一实施例中,所述M取值为5。

    在本发明一实施例中,以最小化柔顺度为目标函数,以体积为约束条件,以单元本
    身为设计变量,基于灵敏度分析改进的BESO刚度优化数学表达式可描述为:


    式(7)中F和U分别为外载荷向量和结点位移向量,C为结构柔顺度,Vi是个体单元
    的体积,V*为体积约束,N为设计域中单元总数量;xi为双向赋值的非连续设计变量,即材料
    单元的密度;K为结构总体刚度矩阵,xmin为材料单元的密度的一个小值。

    在本发明一实施例中,所述改进的BESO方法中,删除单元的方法为软杀法,即将灵
    敏度低的待删单元的密度人为减至一个小值xmin,以防止刚度矩阵奇异,可令xmin=0.001,
    但后续仍参与有限元分析,而不是直接删除;因此,采用单元软杀法的BESO方法的单元灵敏
    度计算式可归结为:


    式(8)中Ki0表示实单元的单元刚度矩阵,ui为单元节点位移向量,P为惩罚因子,可
    令P=3。

    在本发明一实施例中,将权值法引入到BESO法中,多载荷工况的刚度拓扑优化问
    题可以描述为:


    式(9)中,m是载荷工况的总数,ωs、Cs分别是工况s的预设权重系数和平均柔顺度,
    并有

    由于各载荷工况的位移场相互独立,带权重的目标函数设计变量xi求偏导数可
    得:


    于是灵敏度计算式可以归结为:


    在本发明一实施例中,在等效静态载荷法中,等效静态载荷被定义为能够产生与
    动态载荷相同位移场的一组载荷集;可以看出,第i个等效静载荷作用下产生的位移响应与
    第i个时间步上的动态响应完全相等;

    已知动态分析控制方程:


    式(12)中,M、C、K分别为结构质量、阻尼和刚度矩阵,均为设计变量x的函数;
    和u(t)分别为加速度、速度和位移向量,f(t)为外部施加的动载荷;

    对式(12)式进行变换有:


    式(13)中,t=1,2,3…m;

    再用静态变量s替换feq(t)中的时间变量t,即可算出等效静载荷:

    feq(s)=K(x)u(s),s=1,2,3...m 式(14)。

    在本发明一实施例中,该改进双向渐进法的等效静载荷法动态响应拓扑优化方
    法,具体实现如下,

    S1:进行有限元建模,设定初始设计变量x(k)=x(0),初始设计循环次数k=0;

    S2:对以x(k)作为密度的有限元模型进行动力学分析,求其时间域上的节点位移
    响应u(t);

    S3:将时间域上的连续响应u(t)离散为m个时间点上的响应u(s),并计算m个时刻
    的等效静载荷feq(s),计算公式如下:

    feq(s)=K(x)u(s);s=1,2,3...m 式(15)

    S4:将等效静载荷输入静态多目标拓扑优化数学模型中,每一个时刻的等效载荷
    为一种工况,每个工况赋予ωk=1/l的权重因子,进行优化迭代直至收敛;


    S5:当k=0时,令k=k+1并进入步骤S2;当k≠0,如果满足收敛条件,优化结束;如
    果不满足收敛条件,依然令k=k+1并进入步骤S2;至此完成一次设计循环;收敛条件可以表
    示为:


    该方法将渐进结构法替换变密度法作为内循环方法,第一次外循环时,初始设计
    域为满设计域,此时按照Vk+1=Vk(1±ER),其中k=1,2,3…,计算每次迭代的目标体积,从
    第二次外循环开始,均满足Vk+1=V*,即设计域同时包含0、1单元,此时有Vdel=Vadd=AR,这
    种情况下,渐进结构法的优化速度要比满设计域时快。

    相较于现有技术,本发明具有以下有益效果:本发明方法将改进的BESO法引入到
    等效静载荷优化中,改进了动态响应拓扑优化流程,并建立了新的收敛条件;通过matlab软
    件对2个算例进行了优化分析,经比较可以看出,本发明提出的方法,能有效减少k大于1次
    设计循环阶段的优化时间,两个算例中优化所需CPU总时间比原始方法分别减少了21.8%、
    25.3%。该方法具有计算速度快、优化流程简单的特点;随着外载荷作用点的增多、优化规
    模的增大,优化所需的设计循环次数也会增加,该方法将会带来更大的优势。

    附图说明

    图1为本发明静态载荷等效过程图。

    图2为本发明改进双向渐进法的等效静载荷法动态响应拓扑优化方法流程图。

    图3为算例一。

    图4(a)为针对算例一的原始方法优化结果图。

    图4(b)为针对算例一的本申请方法优化结果图。

    图5为算例二。

    图6(a)为针对算例二的原始方法优化结果图。

    图6(b)为针对算例二的本申请方法优化结果图。

    具体实施方式

    下面结合附图,对本发明的技术方案进行具体说明。

    本发明的一种改进双向渐进法的等效静载荷法动态响应拓扑优化方法,将改进的
    BESO方法应用于等效静载荷法中作为优化迭代方法,替换掉原始的变密度法,并建立新的
    收敛条件,从而改进动态响应拓扑优化流程。

    所述改进的BESO方法具体如下,

    在对单元处理之前,应事先给出每次迭代的单步目标体积Vk+1;初始体积Vk大于或
    小于整体的体积约束V*;因此,体积进化迭代式可表示为:

    Vk+1=Vk(1±ER),(k=1,2,3...) 式(1)

    式(1)中k为当前迭代次数,ER为体积进化率;

    若目标体积Vk+1小于当前体积Vk,则必须将单元按灵敏度大小排列并删除灵敏度
    低的单元,在删除单元之前,空单元的灵敏度有可能大于待删单元灵敏度,因此要先考虑增
    添单元,这里引入增删阀值αth,αth定义为待删单元灵敏度中的最大值,有:

    αi>αth 式(2)

    式(2)中αi为增添单元的灵敏度;当增添单元的体积大于体积增添率AR时,需要重
    新计算αth直到满足增添单元的体积小于或等于AR为止;然后删除灵敏度小于该值的单元;
    显然,删除单元的体积Vdel可以表示为:

    Vdel=Vk-Vk+1+Vadd 式(3)

    式(3)中Vadd为增添单元的体积;

    若目标体积Vk+1大于当前体积Vk,则必须将单元按灵敏度大小排列并增添灵敏度
    高的单元,在增添单元之前,实单元的灵敏度有可能小于待增单元灵敏度,因此要先考虑删
    除单元,这里引入增删阀值βth,βth定义为待增单元灵敏度中的最小值,有:

    βi<βth 式(4)

    其中βi为删除单元的灵敏度;当删除单元的体积大于体积增添率AR时,需要重新
    计算βth直到满足删除单元的体积小于或等于AR为止;然后增添灵敏度大于该值的单元;显
    然,增添单元的体积Vadd可以表示为:

    Vadd=Vk-Vk+1+Vdel 式(5)

    重复进行有限元分析与单元增删的过程,直至迭代循环满足两个终止条件:一是
    体积约束,即当前进化体积已满足整体约束体积V*;二是收敛准则,考虑目标函数的变化
    量,收敛判别式规定为:


    式(6)中k为当前迭代步数,τ为内部许用收敛因子;M是一整数,用以限制并稳定柔
    顺度在多次连续迭代中的平均变化量,通常取5;计入局部多个M值,能够使目标函数波动较
    小,曲线更为平滑,有效防止局部收敛。

    以最小化柔顺度为目标函数,以体积为约束条件,以单元本身为设计变量,基于灵
    敏度分析改进的BESO刚度优化数学表达式可描述为:


    式(7)中F和U分别为外载荷向量和结点位移向量,C为结构柔顺度,Vi是个体单元
    的体积,V*为体积约束,N为设计域中单元总数量;xi为双向赋值的非连续设计变量,即材料
    单元的密度;K为结构总体刚度矩阵,xmin为材料单元的密度的一个小值。

    所述改进的BESO方法中,删除单元的方法为软杀法,即将灵敏度低的待删单元的
    密度人为减至一个小值xmin,以防止刚度矩阵奇异,可令xmin=0.001,但后续仍参与有限元
    分析,而不是直接删除;因此,采用单元软杀法的BESO方法的单元灵敏度计算式可归结为:


    式(8)中Ki0表示实单元的单元刚度矩阵,ui为单元节点位移向量,P为惩罚因子,可
    令P=3。

    将权值法引入到BESO法中,多载荷工况的刚度拓扑优化问题可以描述为:


    式(9)中,m是载荷工况的总数,ωs、Cs分别是工况s的预设权重系数和平均柔顺度,
    并有

    由于各载荷工况的位移场相互独立,带权重的目标函数设计变量xi求偏导数可
    得:


    于是灵敏度计算式可以归结为:


    如图1所示,在等效静态载荷法中,等效静态载荷被定义为能够产生与动态载荷相
    同位移场的一组载荷集;可以看出,第i个等效静载荷作用下产生的位移响应与第i个时间
    步上的动态响应完全相等;

    已知动态分析控制方程:


    式(12)中,M、C、K分别为结构质量、阻尼和刚度矩阵,均为设计变量x的函数;
    和u(t)分别为加速度、速度和位移向量,f(t)为外部施加的动载荷;

    对式(12)式进行变换有:


    式(13)中,t=1,2,3…m;

    再用静态变量s替换feq(t)中的时间变量t,即可算出等效静载荷:

    feq(s)=K(x)u(s),s=1,2,3...m 式(14)。

    如图2所示,本发明的改进双向渐进法的等效静载荷法动态响应拓扑优化方法,具
    体实现如下,

    S1:进行有限元建模,设定初始设计变量x(k)=x(0),初始设计循环次数k=0;

    S2:对以x(k)作为密度的有限元模型进行动力学分析,求其时间域上的节点位移
    响应u(t);

    S3:将时间域上的连续响应u(t)离散为m个时间点上的响应u(s),并计算m个时刻
    的等效静载荷feq(s),计算公式如下:

    feq(s)=K(x)u(s);s=1,2,3...m 式(15)

    S4:将等效静载荷输入静态多目标拓扑优化数学模型中,每一个时刻的等效载荷
    为一种工况,每个工况赋予ωk=1/l的权重因子,进行优化迭代直至收敛;


    S5:当k=0时,令k=k+1并进入步骤S2;当k≠0,如果满足收敛条件,优化结束;如
    果不满足收敛条件,依然令k=k+1并进入步骤S2;至此完成一次设计循环;收敛条件可以表
    示为:


    该方法将渐进结构法替换变密度法作为内循环方法,第一次外循环时,初始设计
    域为满设计域,此时按照Vk+1=Vk(1±ER),其中k=1,2,3…,计算每次迭代的目标体积,从
    第二次外循环开始,均满足Vk+1=V*,即设计域同时包含0、1单元,此时有Vdel=Vadd=AR,这
    种情况下,渐进结构法的优化速度要比满设计域时快。

    传统变密度法如果直接对含有中间密度单元的模型进行动力学分析会引起的网
    格畸变[22],因此需要在下一次外循环开始前运用EEM(eliminate element method)法更新
    设计变量以避免该现象,而BESO法由于其单元密度只有1和Xmin两种,可以直接对其优化模
    型进行动力学分析,简化了优化流程。需要注意的是,由于BESO法单元密度为1或Xmin,不是
    一个区间内的连续数值,因此收敛条件也不同于变密度法。

    下面根据原始等效静载荷法和本发明提出的方法计算两个算例。

    单个动载荷条件下的优化:

    如图3,设计域为160mm*100mm的二维悬臂梁结构,将其划分为40*25个单元的有限
    元模型。模型右边中点受到一个周期的正弦动载荷f1=1000sin(5πt)[N]作用。每个时间点
    的载荷工况权值为ωs=1/l。算例其他优化参数如表1所示。

    表1拓扑优化参数



    两种方法的优化结果如图4所示。

    结果数据对比如表2所示。

    表2结果对比


    从优化结果图4(a)和图4(b)可看出,对于受到单个动载荷的悬臂梁结构,本发明
    研究所提出的方法得到的优化结果与原始方法相似,并且目标函数值也非常接近,可以认
    为该方法是可行的。同时本发明研究提出的方法在第二次设计循环阶段的优化时间为
    28.08s,明显快于原始方法在第二次设计循环阶段的56.20s,两种方法的最终优化时间分
    别为90.48s、115.74s。因此从整体的优化效率上看,本发明提出的方法要高于原始方法。

    多个动载荷条件下的优化:

    如图5,设计域为200mm*100mm的二维悬臂梁结构,将其划分为50*25个单元的有限
    元模型。模型右上方受到一个周期的正弦动载荷f1=10000sin(5πt)[N]作用,下方中点受
    到一个周期的余弦动载荷f2=20000cos(5πt)[N]作用。

    每个时间点的载荷工况权值为ωs=1/l。算例其他优化参数如表3所示。

    表3拓扑优化参数


    两种方法的优化结果如图6所示。

    结果数据对比如表4所示。

    表4结果对比


    从优化结果图6(a)和图6(b)可看出,对于受到两个动载荷的悬臂梁结构,得到的
    结论与算例1类似,两种方法下优化结果与目标函数非常接近,可以认为BESO法作为优化迭
    代方法可行。同时随着载荷作用点数的增加,设计循环次数增加到3次,由表4的优化结果中
    可以看出,本发明提出的方法在第二、三次设计循环阶段的优化时间分别为43.14s和
    40.60s,均明显快于原始方法在第二、三次设计循环阶段的71.14s和75.87s,两种方法的最
    终优化时间分别为171.27s、229.16s。因此从整体的优化效率上看,本发明提出的方法要高
    于原始方法。

    参考文献:

    [1]BENDSOE M P,SIGMUND O.Topology optimization:Theory,Methods,and
    Applications[M].New York:Methods and Applications,2003.

    [2]BENDSOE M P,KIKUCHI N.Generating optimal t-opologies in structural
    design using a homogenization method[J].Computer Methods in Applied Mechanics
    and Engineering,1988,71(2):197-224.

    [3]BENDSOE M P.Optimal Shape Design as a Materi-al Distribution
    Problem[J].Structural Optimization,1989,1(4):193-202.

    [4]Xie Y M,and StevenG P.A Simple Evolutionary Procedure for
    Structural Optimization[J].Computers and Structures,1993,49(5):885-896.

    [5]Querin O M,Young V,Steven G P,Xie Y M.Computational efficiency and
    validation of bi-directional evolutionary structural optimization[J]
    .Comput.Meth.Appl.Engng,2000(189):559-573.

    [6]Rozvany G I.A critical review of established metho-ds of
    structural topology optimization[J].Structural and Multidisciplinary
    Optimization,2009,37(3):217-237.

    [7]Huang X,Xie M.Evolutionary topology optimization of continuum
    structures:methods and applications[M].John Wiley&Sons,2010.

    [8]Choi S H,Kim S R,Park J Y,Han S Y.Multi-Objective Optimization Of
    The Inner Reinforcement For A Vehicle’S Hood Considering Static Stiffness And
    Natural Frequency[J].International Journal Of Automotive Technology,2007,8
    (3):337-342.

    [9]Choi K H,Park J Y,Ryu S P,Han S Y.Reliabili ty-Based Topology
    Optimization Based On Bidirectional Evolutionary Structural Optimization
    Using Multi-Objective Sensitivity Numbers[J].International Journal Of
    Automotive Technology,2011,12(6):849-856.

    [10]Afimiwala KA,Mayne RW.Optimal Design of an Impact Absorber[J]
    .Journal of Engineering for Industry,Transactions of the ASME 96:124-130.

    [11]Fox RL,Kapoor MP.Structural Optimization in the Dynamic Regime
    [J].A Computational Approach.AIAA Journal 8:1798-1804.

    [12]Choi W S,Park G J.Structural Optimization Using Equivalent Static
    Loads at All the Time Intervals[J].Computer Methods in Applied Mechanics and
    Engin-eering,2002,191(19-20):2105-2122.

    [13]PARK G J,KANG B S.Validation of a Structural Optimization
    Algorithm Transforming Dynamic Loadsinto Equivalent Static Loads[J].Journal
    Of Optimizat-ion Theory And Applications,2003,118(1):191-200.

    [14]Jang H H,Lee H A,Park G J.Preliminary Study on Linear Dynamic
    Response Topology OptimizationUsing Equivalent Static Loads[J].Transactions
    of theKorean Society of Mechanical Engineers,2009,33(12):1357-1493.

    [15]Lee H A,Park G J.Nonlinear dynamic response top-ology
    optimization using the equivalent static loads method[J].Computer Methods in
    Applied Mechanics and Engineering,2015,283(1):956-970.

    [16]苏义脑,唐雪平.初弯曲纵横弯曲梁的等效载荷法及其应用[J].力学与实践,
    20040.

    Su Yinao,Tang Xueping.The equivalent load method and its application
    to the bending and transverse beam of the initial bending beam[J].Mechanics
    and Practice,20040.

    [17]黄武龙.基于等效静态载荷方法的大型复杂结构的轻量化设计[D].广州:广
    东工业大学,2013(5)17.

    Huang Wulong.Based on lightweight design of equivalent static load
    method of large complex structures[D].Guangzhou:Guangdong University of
    technology,2013(5)17.

    [18]赵礼辉.ESL法在汽车结构优化设计中的应用[D].上海:上海交通大学,2009
    (2).

    Zhao Lihui.Application of ESL method in automobile structure
    optimization design[D].Shanghai:Shanghai Jiao Tong University,2009(2).

    [19]杨志军.基于等效静态载荷原理的高速机构结构拓扑优化方法[J].机械工程
    学报,2011,9(47):119-126.

    Yang Zhijun.Topology optimization method of high speed mechanism
    based on equivalent static load principle[J].Journal of mechanical
    engineering,2011,9(47):119-126.

    [20]芮强,王红岩.基于等效静态载荷法的结构动态优化[J].汽车工程,2014(1).

    Rui Qiang,Wang Hongyan.Dynamic optimization of structure based on
    equivalent static load method[J].Automotive Engineering,2014(1).

    [21]Kim E,Kim H,Baek S,Cho M.Effective structural o-ptimization based
    on equivalent static loads combined with system reduction method[J]
    .Structural and Mul-tidisciplinary Optimization,2014,50:775-786.

    [22]Lee J Y,Park G J.Dynamic Response Topology Optimization in the
    Time Domain Using Equivalent Static Loads[J].AIAA Journal,2012,50(1):226-
    234.。

    以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作
    用未超出本发明技术方案的范围时,均属于本发明的?;し段?。

    关于本文
    本文标题:改进双向渐进法的等效静载荷法动态响应拓扑优化方法.pdf
    链接地址://www.4mum.com.cn/p-6014057.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 看牌神器 扎金花 pk10一天稳赚5000图片 三肖6码 看极速时时开奖结果 11选5任二稳赚玩法 一分快三大小单双稳赚买法 上海时时五星走势图 北京pk拾人工精准计划 双色球走势图表近100期 后三包胆技巧 nba投注量 澳洲pk10计划软件 十期 倍投稳赚方案图片 8码倍投 稳赚三期必出六肖 江苏时时