一种基于结构分解的大型构件动态优化设计方法
【技术领域】
[0001] 本发明设及大型复杂构件的动态优化设计领域，具体地说是一种基于结构分解的 大型构件动态优化设计方法。
【背景技术】
[0002] 随着我国现阶段工业化的蓬勃发展，一些更加大型和复杂的重载设备逐渐出现在 我们的视线中，例如:大型重载的数控机床、大型掘进机和采煤机、大吨位的锻造机W及大 型高伸缩行程的吊车和起重机等。运些设备都是由一些大型复杂的零部件组合而成，运些 大型结构件的动态特性将会直接影响整个设备的动态性能。因此，针对于我国目前在工程 应用中大型复杂结构零件的动态特性优化设计问题一直是目前的研究热点。
[0003] 结构优化设计是随着最优化方法的不断发展和改善逐渐得W发展。近些年来，在 结构优化算法方面，结构优化设计趋向于采用接近实际的复杂结构模型模拟大型结构系 统，由于设计变量数目大，研究新的有效的准则优化方法受到重视，但仍有如何去解决针对 各种特殊的结构优化问题建立相应的公式，解决解析推导和数值计算的实现问题;再是使 用大型系统的分解优化方法，对于大型结构优化，可W按子结构分解或者进行多级分解优 化，对于多学科的复杂系统可W按学科分解优化。分解算法的关键在于建立各个子问题之 间的禪合关系，比如通过使用最优解对参数的灵敏度和采用线性分解等法建立起禪合关 系，使得子问题的解相容，从而保证迭代收敛。
[0004] 特别是针对于复杂重载大型的机械结构零件，采用基于结构分解的多目标动态结 构优化技术对于准确找到整个结构的最优设计方案是一种快速而有效的方法。
[0005] 现阶段我国大型重载设备在作业过程中由于自身的动态性能(包括刚度、强度、阻 尼和振动特性等)不足或不匹配而导致整个设备自身出现结构破坏、结构变形、振动剧烈或 控制效果不理想，最终导致整个作业任务的没有达到既定目标和失败。
【发明内容】
[0006] 针对现有技术的不足，本发明提供一种基于结构分解的大型构件动态优化设计方 法，根据已知的边界载荷和约束的前提下，建立起综合结构动态优化技术流程，深入研究= 维实体结构的拓扑优化和尺寸优化技术，完善优化算法的鲁棒性。
[0007] 本发明为实现上述目的所采用的技术方案是：
[000引一种基于结构分解的大型构件动态优化设计方法，
[0009] 步骤1:在=维建模软件中，建立大型构件的几何模型，识别构成结构件的外部框 架尺寸变量，识别内部组成单元的基本样式和尺寸变量；
[0010] 步骤2:对结构件的外部框架进行拓扑优化，同时对内部的组成单元进行尺寸优 化，得到尺寸数据；
[0011] 步骤3:将优化后的尺寸数据导入到Isight-FD软件中，根据约束条件和目标函数 分别对外部框架和内部组成单元的不同组合W及不同尺寸配置进行优化迭代分析，输出优 化结果。
[0012] 所述识别构成结构件的外部框架尺寸变量的过程为将结构件的外部轮廓尺寸定 义为外部框架尺寸变量；
[0013] 所述识别内部组成单元的基本样式和尺寸变量过程为:将组成大型零件内部的结 构用不同样式的基本单元体来进行填充，再针对于每一种基本单元体不同部位的尺寸进行 变量定义，基本单元体的样式即为内部组成单元的基本样式，基本单元体的尺寸变量即为 内部组成单元的尺寸变量。
[0014] 所述对结构件的外部框架进行拓扑优化采用变密度法：
[001引 Find:P=(Pl,P2，...，Pn)TeRn
[0016] Min:C(p)=F^ = X\X
[0017] S.T:F=KX 引 n抑=艺二 i 二 1
[0019] 0<Pmin<f?i< l，i = l，2,---,n
[0020] 其中，P为设计变量，表示材料的相对密度;n为大型构件结构离散成的单元总数;C (P)是优化目标函数，表示大型构件结构的柔度，K为结构的总体刚度矩阵，X为结构的总体 位移向量，F为结构所受的载荷向量，Vi为结构单元的体积，V为大型构件优化后的体积，Vo为 大型构件初始状态的体积，Vmax为体积上限，f为体分比，Pmin为最小相对密度。
[0021] 所述对内部的组成单元进行尺寸优化采用基于Kriging插值的目标响应面的方 法：
[0022] Find:X=[xiX2. . .Xn]T
[0023] Min:F(X)=W/fi'
[0024] S.T:g(xiX2. . .Xn) =fi-f/ <0
[0025] XER
[0026] 其中，W为被优化结构件的质量，Xn为被优化结构的尺寸设计变量，fi为优化前结构 的第一阶固有频率，fi'为优化后结构的第一阶固有频率，F(X)为目标函数，表示被优化结 构件的质量与优化后结构的第一阶固有频率之比。
[0027] 所述优化迭代分析过程为：
[0028] 步骤1:根据设计要求将设计变量采用拉下方超立方的方法进行实验设计，得到设 计空间内的实验矩阵；
[0029] 步骤2:对得到的实验矩阵进行基于Kriging插值和目标响应面方法的优化算法， 得到近似模型；
[0030] 步骤3:将近似模型随机采样，分别通过计算机仿真和近似计算，得到真实响应值 和近似响应值；
[0031 ]步骤4:将真实响应值与近似响应值进行精度检验，如果精度大于被优化对象许用 精度，执行步骤5;否则需要新增插值点作为样本点，并更新近似模型样本空间，返回步骤3;
[0032] 步骤5:使用近似模型代替仿真模型，并输出优化结果。
[0033] 所述约束条件为：
[0034] 结构件优化后的体积不大于给定的体积上限数值，即V(P)含Vmax,最小相对密度 Pmin在(0, I)之间；
[0035] 结构件优化后的基频不小于原始结构的基频，即fV > fi。
[0036] 本发明具有W下有益效果及优点：
[0037] 本发明可W广泛适用于类似大型重载复杂结构零件的轻量化和减重设计，具有由 简入深，快速可靠的优点，由于算法的鲁棒性只需要少量的采样点即可W找到满足力学性 能要求的几何结构，为现阶段大型重载设备的结构设计工作提供了可靠的技术支持。
【附图说明】
[0038] 图1是本发明的方法流程图；
[0039] 图2是本发明的大型悬臂构件针对其外部框架拓扑优化的最优材料分布图；
[0040] 图3是本发明的大型悬臂构件针对不同内部单元样式和设计变量的尺寸优化筋格 模型图；
[0041 ]图4是本发明的有限元分析位移云图；
[0042] 图5是本发明的优化迭代分析流程图。
【具体实施方式】
[0043] 下面结合附图及实施例对本发明做进一步的详细说明。
[0044] 如图1所示为本发明的方法流程图。
[0045] 步骤一，首先采用=维造型软件建立起复杂零件的几何模型；识别出构成结构件 的主要框架尺寸变量和组成内部单元的基本样式和尺寸变量。具体的识别方法为:将结构 件的外部轮廓尺寸定义为框架尺寸变量，例如长、宽、高和横截面积的大小等。将组成大型 零件内部的结构用不同样式的筋格(基本单元体)来进行填充，例如:十字筋、米子筋、斜线 筋和太阳筋等，详见图1。此外，针对于每一种筋格不同部位的尺寸进行变量定义，例如:壁 厚、直径和筋格大小等。运些筋格的样式和尺寸变量就被识别为结构件内部单元的基本样 式和尺寸变量。
[0046] 步骤二，然后，根据结构的质量最轻和基频最高为目标，在Hypermesh/optistruct 软件中分别对结构的外部框架进行拓扑优化，对内部的组成单元进行尺寸优化。从而找到 最优的材料分布路径W及最佳的结构单元尺寸配置。运里面的拓扑优化采用的是变密度的 方法，而尺寸优化采用的是基于Kriging插值的目标响应面的方法，具体的优化算法和优化 原理详见【具体实施方式】中内容。
[0047] 步骤S，将步骤二中优化后的所有尺寸数据导入到Isight-FD软件中，根据具体的 约束条件和目标函数分别对外部框架和内部单元的不同组合W及不同尺寸配置进行优化 迭代分析，从而找到整个结构设计问题的最优解。运里面采用商业优化软件进行多目标的 综合迭代优化。
[004引它将构成大型重载设备的大件结构分为外部和内部两个部分，然后采用有限元分 析方法分别对运两个部分进行动态优化设计。采用具体的设计变量和优化目标，将对影响 整个结构动态性能的重要因素量化，并采用性能优越的迭代算法来进行求解计算。运样，就 可W保证从组成结构的外形也就是宏观上的外部材料分布W及内部的具体样式也就是微 观上的内部细节所构成的结构具有优越的综合性能，运种方法易于操作，实用性强，计算时 采用S维建模软件、Hyperworks/Optistruct和Isight-FD软件进行S维建模和协同优化设 计工作。
[0049] 我们将复杂大件结构划分出外部框架部分和组成其整体的内部单元，然后分别选 取他们的不同组合来进行基于外部框架的有限元分析，从而形成数据库。
[0050] 在上述组成大型复杂零件的结构框架尺寸和基本单元样式基础上，通过选取最优 的组合进行特定约束和目标函数的动态特性分析，使得结构无论是在外部的具体拓扑形态 上，还是内部的详细构成方面都具有最优的力学性能。运样我们即在"微观"和"宏观"两个 方面均考虑到了最终设计出来的大件结构的动态性