CN110110413A - 一种基于材料场缩减级数展开的结构拓扑优化方法 - Google Patents
一种基于材料场缩减级数展开的结构拓扑优化方法 Download PDFInfo
- Publication number
- CN110110413A CN110110413A CN201910341472.6A CN201910341472A CN110110413A CN 110110413 A CN110110413 A CN 110110413A CN 201910341472 A CN201910341472 A CN 201910341472A CN 110110413 A CN110110413 A CN 110110413A
- Authority
- CN
- China
- Prior art keywords
- optimization
- material field
- correlation
- design
- function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000000463 material Substances 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000009467 reduction Effects 0.000 title claims abstract description 13
- 238000013461 design Methods 0.000 claims abstract description 57
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000005314 correlation function Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 239000006390 lc 2 Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010206 sensitivity analysis Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 abstract description 6
- 230000035945 sensitivity Effects 0.000 abstract description 5
- 238000012827 research and development Methods 0.000 abstract description 2
- 230000003595 spectral effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000109 continuous material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
一种基于材料场缩减级数展开的结构拓扑优化方法，解决传统密度法拓扑优化，由于设计变量过多、需要相对密度或灵敏度过滤措施等产生的计算效率低下问题。通过定义一种考虑相关性的有界材料场，采用谱分解方法变换为一系列待定系数的线性组合，并以这些待定系数作为设计变量，基于单元密度插值模型构建优化模型，采用梯度类或非梯度类优化算法对拓扑优化问题进行求解，进而高效率获得带有清晰边界的拓扑构型。该方法能够大幅度减少密度法拓扑优化中的设计变量个数，同时具有完全避免网格依赖性和棋盘格式问题的天然优势。该方法还继承了密度法形式简单，便于工程化推广等优点，优化求解速度快，将保证复杂装备结构创新拓扑设计的研发效率。
Description
技术领域
本发明属于机械、航空航天工程装备结构轻量化设计领域，涉及一种基于 材料场级数展开的结构拓扑优化方法。
背景技术
目前连续体结构拓扑优化的主流方法包括：变密度法，水平集法及(双向) 渐进优化方法。其中变密度法因模型简单、实施简便在机械、航空航天工程结 构的创新拓扑优化设计中被广泛采用，并集成于许多优化设计商业软件中。变 密度法通过引入介于0和1之间的中间材料相对密度并对其进行惩罚，将0-1 离散变量拓扑优化问题转化为设计变量取值连续的优化问题，并通过基于梯度 的优化算法对其进行高效求解。然而，在变密度法中，设计变量的数目依赖于 有限单元的数目。其次，变密度法本身并不能够解决拓扑优化问题固有的网格 依赖性问题及棋盘格式数值不稳定现象，需要通过密度过滤、灵敏度过滤或其他过滤方法施加最小长度尺度来控制上述问题，从而带来了额外的计算量。因 此，当处理网格离散较密的大规模拓扑优化问题时，对大量灵敏度或相对密度 的过滤及对相对密度的更新成为了除有限元分析外最为耗时的优化问题求解环 节。此外，传统密度法方法由于设计变量过多，只能采用基于梯度的算法求解， 无法适用于难以直接得到灵敏度的复杂问题。为有效减少迭代求解大规模设计 变量所带来的计算量，一种可行的选择是在变密度法框架下，基于材料场缩减 级数展开提出新的拓扑优化方法，大幅度减少设计变量数目，提高优化求解效 率，适用于非梯度类算法求解，同时有效消除网格依赖性和棋盘格式现象。
发明内容
针对传统密度法在处理大规模复杂结构拓扑优化问题中设计变量过多的缺 点，本发明提供一种用于大幅度减少设计变量规模的拓扑优化设计方法，该方 法能够大幅度减少密度法拓扑优化中的设计变量个数，提高优化求解效率，同 时具有完全避免网格依赖性和棋盘格式问题的天然优势。本发明适合于航空航 天、机械工程复杂装备的创新拓扑优化设计，适用于非梯度类优化方法求解， 有利于提高优化效率，尤其适用于大规模三维结构拓扑设计问题。
为了达到上述目的，本发明采用的技术方案为：
一种基于材料场缩减级数展开的结构拓扑优化方法，主要包括材料场缩减级 数展开和结构拓扑优化建模两部分，具体步骤如下：
第一步，对设计域材料场进行离散化和缩减级数展开
1.1)根据结构的实际情况和尺寸要求确定二维或三维设计域，定义具有空 间相关性的有界材料场函数，在设计域中均匀选取若干观察点对材料场进行离 散化；所述的观察点个数控制在10,000个以内；所述的材料场函数界限为[-1,1]， 材料场任意两点之间的相关性采用依赖于两点空间距离的相关函数进行定义， 即C(x1,x2)＝exp(-||x1-x2||2/lc 2)；其中，x1和x2为两点空间位置，lc为相关长度， || ||为2-norm范数。
1.2)确定相关长度，计算所有观察点之间的相关性，构造对角线为1的对 称正定相关矩阵；所述的相关长度不大于设计域长边尺寸的25％。
1.3)对步骤1.2)中的对称正定相关矩阵进行特征值分解，将特征值从大到 小排序，依照截断准则选取前几阶特征值；所述的截断准则为：选取的特征值 之和占所有特征值总和的99.9999％。
1.4)对材料场进行缩减级数展开，即其中，η为待定 的级数展开系数，Λ为1.3)中选取的特征值组成的对角矩阵，ψ为相应的特征 向量，C(x)为步骤1.1)中相关函数得到的相关向量。
第二步，对结构进行拓扑优化
2.1)首先，对设计域划分有限元单元网格，建立有限元单元弹性模量与材 料场的幂函数插值关系；其次，在设计域中施加荷载和约束边界，进行有限元 分析；最后建立结构拓扑优化模型，优化目标为使结构的整体刚度最大化或者 整体柔顺性最小化，约束条件和设计变量如下：
a)约束条件一：要求每个观察点的材料场函数值不超过1；
b)约束条件二：确定结构材料用量，作为材料体积约束上限；所述的材料 用量为设计域体积的5％-50％；
c)设计变量：材料场的缩减级数展开系数η，η取值为-100和100之间。
2.2)根据步骤2.1)建立的结构拓扑优化模型，对优化目标和约束条件进行 灵敏度分析；利用梯度类优化算法或者非梯度类优化算法进行迭代求解，在迭 代过程中采用紧约束处理策略，仅将当前观察点材料场函数值大于-0.3的约束条 件计入算法中，得到结构最优材料分布。
进一步的，步骤1.1)所述的相关函数包括指数模型函数和高斯模型函数。
进一步的，步骤2.1)所述的单元弹性模量幂函数插值关系表达式为其中，为的Heaviside映射函数，光滑参数从0 到9阶梯性增加，即满足收敛条件后每次增加1.5；收敛条件为相临迭代步目标 值相对差小于0.005；P为惩罚因子，一般取3；E0为所述材料的弹性模量。
进一步的，步骤2.1)所述的梯度类优化算法为准则法或MMA算法，非梯 度类优化算法为代理模型方法或遗传算法。
本发明的有益效果为：在传统密度法拓扑优化在大规模复杂结构拓扑设计 中，设计变量多，且需要耗时的灵敏度或密度过滤方法，从而严重影响优化效 率。采用本发明方法进行大规模复杂结构的拓扑优化设计，能大幅度减少设计 变量规模，高效率得到具有清晰边界的拓扑构型。该方法还继承了密度法形式 简单，便于工程化推广，易于理解和编程等优点，优化求解速度快，适用于非 梯度类求解算法，将保证复杂装备结构创新拓扑设计的研发效率。
附图说明
图1为本发明实施例提供的一种二维MBB梁结构设计域。图中：F表示 施加在结构上的荷载。
图2为二维MBB梁结构最优拓扑构型。
图3为本发明实施例提供的一种三维悬臂梁结构设计域。
图4(a)为材料体积比为7.5％时的三维悬臂梁结构最优拓扑设计图。
图4(b)为材料体积比为30％时的三维悬臂梁结构最优拓扑设计图。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施例。
一种基于材料场缩减级数展开的结构拓扑优化方法，该拓扑优化方法通过 定义一种考虑相关性的有界材料场，采用谱分解方法变换为一系列待定系数的 线性组合，并以这些待定系数作为设计变量，采用单元密度插值模型和梯度类 优化算法对拓扑优化问题进行求解，进而高效率获得带有清晰边界的拓扑构型。
第一步，对设计域材料场进行离散化和缩减级数展开
1.1)根据结构的实际情况和尺寸要求确定二维或三维设计域，在设计域中 均匀选取若干观察点，定义具有空间相关性的有界材料场函数。图1为二维MBB 梁结构的设计域，设计域长为180mm，宽为30mm，选取均匀分布观察点数N ＝2700个；图3为三维悬臂梁结构的设计域，选取均匀分布观察点数N＝6570 个；所述的材料场函数界限为[-1,1]，材料场任意两点之间的相关性采用依赖于 两点空间距离的相关函数进行定义，表达式为C(x1,x2)＝exp(-||x1-x2||2/lc 2)。
1.2)确定相关长度，计算所有观察点之间的相关性，构造对角线为1的对 称正定相关矩阵；图1的材料场相关长度lc＝2mm和8mm，图3的材料场相 关长度lc＝6mm。
1.3)对特征矩阵进行特征值分解，将特征值从大到小排序，依照截断准则 选取前几阶特征值；所述的截断准则为：选取的特征值之和占所有特征值总和 的99.9999％。
1.4)对材料场进行缩减级数展开，即其中，η为待定 的级数展开系数，Λ为1.3)中选取的特征值组成的对角矩阵，ψ为相应的特征 向量，C(x)为1.1)中相关函数得到的相关性向量。
第二步，对结构进行拓扑优化
2.1)首先，对设计域划分有限元单元网格，图1设计域划分有限元网格数 NE＝43200个，图3设计域划分有限元网格数NE＝93312个。建立有限元单元 弹性模量与材料场的幂函数插值关系为其中，为 的Heaviside映射函数，光滑参数从0到9阶梯性增加，即满足收敛条件后 每次增加1.5；收敛条件为相临迭代步目标值相对差小于0.005。其次，在设计 域中施加荷载和约束边界，进行有限元分析。最后，建立结构拓扑优化模型， 使结构的整体柔顺性最小化。
a)约束一：要求每个观察点的材料场函数值不超过1；
b)约束二：确定结构材料用量，作为材料体积约束上限；图1的材料体积 比为50％，图3的材料体积比为7.5％和30％；
c)设计变量：材料场的缩减级数展开系数η，η取值为-100和100之间。
2.2)根据步骤2.1)建立的拓扑优化模型，对目标函数和约束条件进行灵敏 度分析；利用梯度类优化算法或者非梯度类优化算法进行迭代求解，在迭代过 程中采用紧约束处理策略，仅将当前观察点材料场函数值大于-0.3的约束条件 计入优化算法中，得到结构最优拓扑构型，分别见图2和图4。
本发明的实质是引入具有空间相关性的材料场，采用谱分解方法将连续的 材料场进行变换，达到减少设计变量的目的，并能内在地避免棋盘格式和网格 依赖性。其对前述各实施例所记载的优化模型、方法、方案进行修改，或者对 其中部分或者全部方法特征进行等同替换(例如采用其他幂函数插值关系、改 变目标函数或约束具体形式等)，并不使相应方法与方案的本质脱离本发明各 实施例方法与方案的范围。
Claims (5)
1.一种基于材料场缩减级数展开的结构拓扑优化方法，其特征在于，该结构拓扑优化方法主要包括材料场缩减级数展开和结构拓扑优化建模两部分，步骤如下：
第一步，对设计域材料场进行离散化和缩减级数展开
1.1)根据结构的实际情况和尺寸要求确定二维或三维设计域，定义具有空间相关性的有界材料场函数，在设计域中均匀选取若干观察点对材料场进行离散化；所述的观察点个数控制在10,000个以内；所述的材料场函数界限为[-1,1]，材料场任意两点之间的相关性采用依赖于两点空间距离的相关函数进行定义，即C(x1,x2)＝exp(-||x1-x2||2/lc 2)；其中，x1和x2为两点空间位置，lc为相关长度，|| ||为2-norm范数；
1.2)确定相关长度，计算所有观察点之间的相关性，构造对角线为1的对称正定相关矩阵；所述的相关长度不大于设计域长边尺寸的25％；
1.3)对步骤1.2)中的对称正定相关矩阵进行特征值分解，将特征值从大到小排序，依照截断准则选取前几阶特征值；所述的截断准则为：选取的特征值之和占所有特征值总和的99.9999％；
1.4)对材料场进行缩减级数展开，即其中，η为待定的级数展开系数，Λ为1.3)中选取的特征值组成的对角矩阵，ψ为相应的特征向量，C(x)为步骤1.1)中相关函数得到的相关向量；
第二步，对结构进行拓扑优化
2.1)首先，对设计域划分有限元单元网格，建立有限元单元弹性模量与材料场的幂函数插值关系；其次，在设计域中施加荷载和约束边界，进行有限元分析；最后建立结构拓扑优化模型，优化目标为使结构的整体刚度最大化或者整体柔顺性最小化，约束条件和设计变量如下：
a)约束条件一：要求每个观察点的材料场函数值不超过1；
b)约束条件二：确定结构材料用量，作为材料体积约束上限；所述的材料用量为设计域体积的5％-50％；
c)设计变量：材料场的缩减级数展开系数η，η取值为-100和100之间；
2.2)根据步骤2.1)建立的结构拓扑优化模型，对优化目标和约束条件进行灵敏度分析；利用梯度类优化算法或者非梯度类优化算法进行迭代求解，在迭代过程中采用紧约束处理策略，仅将当前观察点材料场函数值大于-0.3的约束条件计入算法中，得到结构最优材料分布。
2.根据权利要求1所述的一种基于材料场缩减级数展开的结构拓扑优化方法，其特征在于，步骤1.1)所述的相关函数包括指数模型函数和高斯模型函数。
3.根据权利要求1或2所述的一种基于材料场缩减级数展开的结构拓扑优化方法，其特征在于，步骤2.1)所述的单元弹性模量幂函数插值关系表达式为其中，为的Heaviside映射函数，光滑参数从0到9阶梯性增加，即满足收敛条件后每次增加1.5；收敛条件为相临迭代步目标值相对差小于0.005；P为惩罚因子；E0为所述材料的弹性模量。
4.根据权利要求1或2所述的一种基于材料场缩减级数展开的结构拓扑优化方法，其特征在于，步骤2.1)所述的梯度类优化算法为准则法或MMA算法，非梯度类优化算法为代理模型方法或遗传算法。
5.根据权利要求3所述的一种基于材料场缩减级数展开的结构拓扑优化方法，其特征在于，步骤2.1)所述的梯度类优化算法为准则法或MMA算法，非梯度类优化算法为代理模型方法或遗传算法。
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910341472.6A CN110110413B (zh) | 2019-04-26 | 2019-04-26 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
| KR1020217009347A KR20210047351A (ko) | 2019-04-26 | 2019-08-12 | 재료장 축소 급수전개에 기초하는 구조 토폴로지 최적화 방법 |
| EP19870048.6A EP3754521A4 (en) | 2019-04-26 | 2019-08-12 | STRUCTURAL TOPOLOGY OPTIMIZATION PROCESS BASED ON A MATERIAL FIELD REDUCED SERIES EXPANSION |
| JP2020522826A JP6951001B2 (ja) | 2019-04-26 | 2019-08-12 | 材料場縮減級数展開に基づく構造トポロジー最適化方法 |
| US16/757,669 US20210073428A1 (en) | 2019-04-26 | 2019-08-12 | Structure topology optimization method based on material-field reduced series expansion |
| PCT/CN2019/100131 WO2020215533A1 (zh) | 2019-04-26 | 2019-08-12 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
| AU2019442319A AU2019442319B2 (en) | 2019-04-26 | 2019-08-12 | Structural topology optimization method based on material-field reduction series expansion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910341472.6A CN110110413B (zh) | 2019-04-26 | 2019-04-26 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110110413A true CN110110413A (zh) | 2019-08-09 |
| CN110110413B CN110110413B (zh) | 2022-11-18 |
Family
ID=67486709
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910341472.6A Active CN110110413B (zh) | 2019-04-26 | 2019-04-26 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US20210073428A1 (zh) |
| EP (1) | EP3754521A4 (zh) |
| JP (1) | JP6951001B2 (zh) |
| KR (1) | KR20210047351A (zh) |
| CN (1) | CN110110413B (zh) |
| AU (1) | AU2019442319B2 (zh) |
| WO (1) | WO2020215533A1 (zh) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110609975A (zh) * | 2019-09-20 | 2019-12-24 | 大连理工大学 | 一种基于Kriging代理模型的并行计算全局优化算法 |
| CN110852011A (zh) * | 2019-11-08 | 2020-02-28 | 大连理工大学 | 一种基于序列Kriging代理模型的结构非梯度拓扑优化方法 |
| CN111027110A (zh) * | 2019-11-27 | 2020-04-17 | 中国科学院光电技术研究所 | 一种连续体结构拓扑与形状尺寸综合优化方法 |
| CN111090942A (zh) * | 2019-12-18 | 2020-05-01 | 浙江大学 | 基于拓扑优化的高灵敏度压阻式单轴力传感器设计方法 |
| WO2020215533A1 (zh) * | 2019-04-26 | 2020-10-29 | 大连理工大学 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
| CN112182941A (zh) * | 2020-10-23 | 2021-01-05 | 南京大学 | 一种间隔收缩隔声结构拓扑优化方法 |
| CN112257319A (zh) * | 2020-10-26 | 2021-01-22 | 大连理工大学 | 一种基于非梯度拓扑优化的声学超材料设计方法 |
| CN112818576A (zh) * | 2021-01-28 | 2021-05-18 | 华中科技大学 | 一种曲线纤维复合结构设计多层级优化方法 |
| CN113806982A (zh) * | 2021-09-17 | 2021-12-17 | 北京航空航天大学杭州创新研究院 | 变结构可穿戴柔性热电器件基底拓扑优化方法 |
| WO2022087770A1 (zh) * | 2020-10-26 | 2022-05-05 | 大连理工大学 | 一种基于非梯度拓扑优化的声学超材料设计方法 |
| CN115630556A (zh) * | 2022-12-07 | 2023-01-20 | 东南大学 | 一种基于顶点法的电机拓扑优化方法 |
| CN117133385A (zh) * | 2023-07-27 | 2023-11-28 | 中南大学 | 一种弹性超材料非梯度多目标拓扑优化方法及系统 |
Families Citing this family (56)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112417599B (zh) * | 2020-11-20 | 2023-08-15 | 北京动力机械研究所 | 基于拓扑优化的航空发动机用传动壳体结构设计方法 |
| CN112380633B (zh) * | 2020-12-08 | 2023-09-26 | 安徽江淮汽车集团股份有限公司 | 车身空腔中的加强件设计方法、终端及存储介质 |
| CN112765865B (zh) * | 2021-02-04 | 2022-05-20 | 上海交通大学 | 一种控制金属粉床增材制造热变形的支撑结构设计方法 |
| CN112966419B (zh) * | 2021-03-16 | 2022-09-13 | 武汉大学 | 联合拓扑优化与形状优化的多场耦合作用的结构设计方法 |
| CN113094944B (zh) * | 2021-03-22 | 2024-04-19 | 华中科技大学 | 一种微流道散热器及其细观尺度拓扑优化设计方法 |
| CN113239584B (zh) * | 2021-04-26 | 2022-03-11 | 云南大学 | 一种优化增材制造方法及系统 |
| CN113204906B (zh) * | 2021-05-11 | 2022-10-14 | 华中科技大学 | 一种考虑结构稳定性的多相材料拓扑优化设计方法和系统 |
| CN113345536B (zh) * | 2021-05-31 | 2022-10-11 | 山东大学 | 一种基于极限各向异性点阵材料的结构拓扑优化方法 |
| CN113268842B (zh) * | 2021-06-18 | 2024-05-10 | 西安交通大学 | 一种像素式拓扑优化结果的高保真矢量图转换方法 |
| CN113536623B (zh) * | 2021-06-24 | 2022-03-25 | 河海大学 | 一种材料不确定性结构稳健性拓扑优化设计方法 |
| CN113673123B (zh) * | 2021-07-02 | 2023-09-01 | 华南理工大学 | 一种实体-壳耦合结构拓扑优化方法 |
| CN113673125B (zh) * | 2021-07-16 | 2023-08-22 | 华南理工大学 | 一种用于装配式防撞护栏的连接件的设计方法 |
| CN113515824B (zh) * | 2021-07-29 | 2023-08-08 | 山东大学 | 一种筋条布局与基板形状协同的拓扑优化设计方法 |
| CN113782116B (zh) * | 2021-08-04 | 2024-02-02 | 华中科技大学 | 一种自由曲面热超材料结构及其设计和制备方法 |
| CN113761717B (zh) * | 2021-08-13 | 2024-03-08 | 中南大学 | 一种用于工业优化问题的自动变量约简方法 |
| CN113742970B (zh) * | 2021-08-30 | 2023-12-26 | 华中科技大学 | 一种点阵结构优化设计的基于水平集的变密度方法 |
| CN113887095B (zh) * | 2021-09-09 | 2024-07-26 | 华中科技大学 | 一种基于等几何分析的渐进式结构拓扑优化方法 |
| CN113779802B (zh) * | 2021-09-15 | 2024-02-02 | 湘潭大学 | 基于无网格efgm和等几何分析耦合的结构拓扑优化方法 |
| CN113821887B (zh) * | 2021-09-22 | 2023-06-27 | 湘潭大学 | 基于无网格efgm和plsm的各向异性结构热力耦合拓扑优化方法 |
| CN114030201B (zh) * | 2021-09-26 | 2022-07-15 | 上海交通大学 | 拓扑优化气囊辅助成型复合材料梯形夹芯结构制作工艺 |
| CN113987860B (zh) * | 2021-10-08 | 2022-06-07 | 中山大学 | 一种基于动态高斯核卷积滤波的结构优化方法和装置 |
| CN113868928B (zh) * | 2021-10-19 | 2024-04-12 | 大连理工大学 | 一种面向结构拓扑优化的显式最小尺寸控制方法 |
| CN114048616B (zh) * | 2021-11-17 | 2024-06-14 | 大连理工大学人工智能大连研究院 | 基于径向积分边界元法辨识非线性热导率的修正lm方法 |
| CN114117877B (zh) * | 2021-11-26 | 2024-04-09 | 西安交通大学 | 一种基于等几何粒子描述的拓扑优化方法 |
| CN114357546B (zh) * | 2021-12-21 | 2024-06-25 | 武汉大学 | 梯度随机点阵/多孔结构多尺度优化设计方法 |
| CN114254507B (zh) * | 2021-12-21 | 2024-06-14 | 华中科技大学 | 一种压电执行器稳健性等几何拓扑优化方法及其应用 |
| CN114372361B (zh) * | 2021-12-23 | 2024-05-28 | 北京理工大学 | 基于bddc区域分解并行算法的粗网格选取方法 |
| CN114282372B (zh) * | 2021-12-27 | 2024-06-14 | 华中科技大学 | 一种等几何应力拓扑优化方法及其应用 |
| CN114329664B (zh) * | 2021-12-30 | 2024-06-14 | 华中科技大学 | 一种面向复合材料的等几何拓扑优化方法及其应用 |
| CN114330066B (zh) * | 2021-12-30 | 2024-07-05 | 兰州大学 | 一种同位素热光伏系统中晶元尺寸的确定方法及装置 |
| CN114266184B (zh) * | 2021-12-31 | 2024-05-07 | 西安交通大学 | 一种适配风扇特性曲线的电子设备异型风道设计方法 |
| CN114462274B (zh) * | 2022-01-21 | 2024-10-11 | 哈尔滨工业大学 | 一种功能梯度压电复合材料集成器的设计方法 |
| CN114595655B (zh) * | 2022-03-07 | 2024-03-12 | 西安交通大学 | 一种可调制狭窄空间冷却气流形态的屏蔽罩形状设计方法 |
| CN114707373B (zh) * | 2022-03-07 | 2024-04-09 | 西安交通大学 | 一种高压直流gil电场等效及其结构快速优化方法 |
| CN114757002B (zh) * | 2022-03-25 | 2022-11-08 | 中国船舶重工集团公司第七�三研究所 | 一种约束阻尼板壳的结构-材料分层拓扑优化设计方法 |
| CN114722463B (zh) * | 2022-03-28 | 2024-09-24 | 大连理工大学 | 一种基于各向异性材料场叠加的薄壁结构加筋设计方法 |
| CN114781205B (zh) * | 2022-03-28 | 2024-06-28 | 华中科技大学 | 一种结构拓扑优化模型的灵敏度分析方法及应用 |
| CN114756934B (zh) * | 2022-04-12 | 2024-04-02 | 大连交通大学 | 一种三维多尺度超材料结构优化设计方法 |
| CN114595617B (zh) * | 2022-05-10 | 2022-07-22 | 山东科技大学 | 一种微脉冲激光雷达光路参数拓扑结构的设计方法 |
| CN114912409A (zh) * | 2022-05-13 | 2022-08-16 | 南京邮电大学 | 一种基于三维拓扑优化的芯片被动散热的热沉设计方法 |
| CN114936497B (zh) * | 2022-05-24 | 2024-07-09 | 河北工业大学 | 一种基于缩减基的igbt抗热冲击设计方法及系统 |
| CN114996879B (zh) * | 2022-06-23 | 2024-04-09 | 湖南大学 | 一种针对柔性夹钳机构的区间场几何不确定性拓扑优化方法 |
| CN115081113B (zh) * | 2022-06-23 | 2024-09-13 | 华中科技大学 | 一种基于多水平集函数参数化的网格加强筋结构优化方法 |
| CN115358011B (zh) * | 2022-08-31 | 2023-03-28 | 西北工业大学 | 一种分体式飞行汽车乘客舱结构设计优化方法 |
| CN115495863B (zh) * | 2022-11-22 | 2023-03-14 | 北京理工大学深圳汽车研究院(电动车辆国家工程实验室深圳研究院) | 一种平滑边界表达的散热拓扑优化方法 |
| CN116541910B (zh) * | 2023-06-07 | 2024-02-13 | 黄理鑫 | 一种用于生物冻存的传热模块及其设计、制造方法 |
| CN116843656B (zh) * | 2023-07-06 | 2024-03-15 | 安徽正汇汽配股份有限公司 | 钢带管的涂塑控制方法及其系统 |
| CN116644535B (zh) * | 2023-07-27 | 2023-09-29 | 华东交通大学 | 一种基于基频最大化的柔性铰链拓扑优化方法 |
| CN116644620B (zh) * | 2023-07-27 | 2023-10-17 | 华东交通大学 | 一种基于等几何分析的柔顺机构应力约束拓扑优化方法 |
| CN117057038B (zh) * | 2023-08-15 | 2024-04-05 | 合肥工业大学 | 一种面向机翼的单循环可靠性拓扑优化设计方法 |
| CN117350116B (zh) * | 2023-10-13 | 2024-09-06 | 湖南大学 | 一种考虑惯性载荷的热力耦合拓扑优化方法 |
| CN117610380B (zh) * | 2024-01-19 | 2024-04-09 | 中国科学院苏州生物医学工程技术研究所 | 基于正交实验联合拓扑优化的机器人轻量化设计方法 |
| CN117807813B (zh) * | 2024-02-29 | 2024-05-17 | 浙江大学 | 一种二维拓扑收敛的针织物样式仿真实现方法 |
| CN117892414B (zh) * | 2024-03-18 | 2024-05-24 | 华东交通大学 | 简谐激励下的各向异性材料结构稳健性拓扑优化设计方法 |
| CN117972910B (zh) * | 2024-03-29 | 2024-06-21 | 湖南大学 | 一种多轴智能底盘的转向系统协同设计方法 |
| CN118378554B (zh) * | 2024-06-27 | 2024-08-30 | 华东交通大学 | 基于应力约束的多自由度柔顺机构稳健性拓扑优化方法 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1178413A2 (en) * | 2000-08-03 | 2002-02-06 | Kabushiki Kaisha Toyota Chuo Kenkyusho | Method and system for supporting user in analyzing performance of object, using generalized and specialized models on computer |
| US20060073333A1 (en) * | 1997-09-09 | 2006-04-06 | David Anderson | Coated particles, methods of making and using |
| US20080300831A1 (en) * | 2006-12-19 | 2008-12-04 | Board Of Governors For Higher Education, State Of Rhode Island And Providence | System and method for finite element based on topology optimization |
| CN101950319A (zh) * | 2010-10-09 | 2011-01-19 | 大连理工大学 | 一种高速立式加工中心大跨度横梁拓扑优化设计方法 |
| WO2013176852A1 (en) * | 2012-05-19 | 2013-11-28 | Eric Olsen | Residue number arithmetic logic unit |
| CN103528634A (zh) * | 2013-10-23 | 2014-01-22 | 山东大学 | 一种科里奥利质量流量计云传输数字信号处理装置及方法 |
| CN105681383A (zh) * | 2014-11-21 | 2016-06-15 | 腾讯科技(成都)有限公司 | 应用升级数据包下载方法、装置、终端、服务器及系统 |
| CN106709215A (zh) * | 2017-02-21 | 2017-05-24 | 北京航空航天大学 | 一种基于级数展开的连续体结构非概率可靠性拓扑优化方法 |
| CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
| CN107633106A (zh) * | 2017-07-25 | 2018-01-26 | 东南大学 | 一种基于全局差分法的非均匀温度场下热模态灵敏度分析方法 |
| CN109063357A (zh) * | 2018-08-15 | 2018-12-21 | 大连理工大学 | 基于拓扑优化的自适应多步变域的汽车构件焊点布局优化方法 |
| CN109670200A (zh) * | 2018-11-13 | 2019-04-23 | 华中科技大学 | 一种等几何材料密度场结构拓扑优化方法 |
| CN114722463A (zh) * | 2022-03-28 | 2022-07-08 | 大连理工大学 | 一种基于各向异性材料场叠加的薄壁结构加筋设计方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2763058B1 (en) * | 2013-01-30 | 2021-10-20 | Honda Research Institute Europe GmbH | Optimizing the design of physical structures/objects |
| CN104361246A (zh) * | 2014-11-24 | 2015-02-18 | 河南理工大学 | 基于距离场的功能梯度材料隐式建模方法 |
| WO2017096259A1 (en) * | 2015-12-03 | 2017-06-08 | The Regents Of The University Of Michigan | Topology optimization using reduced length boundaries on structure segments of different thicknesses |
| CN107315872A (zh) * | 2017-06-23 | 2017-11-03 | 华中科技大学 | 一种高效的结构频率响应拓扑优化方法 |
| CN110110413B (zh) * | 2019-04-26 | 2022-11-18 | 大连理工大学 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
-
2019
- 2019-04-26 CN CN201910341472.6A patent/CN110110413B/zh active Active
- 2019-08-12 EP EP19870048.6A patent/EP3754521A4/en not_active Withdrawn
- 2019-08-12 KR KR1020217009347A patent/KR20210047351A/ko not_active Application Discontinuation
- 2019-08-12 AU AU2019442319A patent/AU2019442319B2/en not_active Ceased
- 2019-08-12 JP JP2020522826A patent/JP6951001B2/ja active Active
- 2019-08-12 US US16/757,669 patent/US20210073428A1/en not_active Abandoned
- 2019-08-12 WO PCT/CN2019/100131 patent/WO2020215533A1/zh unknown
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060073333A1 (en) * | 1997-09-09 | 2006-04-06 | David Anderson | Coated particles, methods of making and using |
| EP1178413A2 (en) * | 2000-08-03 | 2002-02-06 | Kabushiki Kaisha Toyota Chuo Kenkyusho | Method and system for supporting user in analyzing performance of object, using generalized and specialized models on computer |
| US20080300831A1 (en) * | 2006-12-19 | 2008-12-04 | Board Of Governors For Higher Education, State Of Rhode Island And Providence | System and method for finite element based on topology optimization |
| CN101950319A (zh) * | 2010-10-09 | 2011-01-19 | 大连理工大学 | 一种高速立式加工中心大跨度横梁拓扑优化设计方法 |
| WO2013176852A1 (en) * | 2012-05-19 | 2013-11-28 | Eric Olsen | Residue number arithmetic logic unit |
| CN103528634A (zh) * | 2013-10-23 | 2014-01-22 | 山东大学 | 一种科里奥利质量流量计云传输数字信号处理装置及方法 |
| CN105681383A (zh) * | 2014-11-21 | 2016-06-15 | 腾讯科技(成都)有限公司 | 应用升级数据包下载方法、装置、终端、服务器及系统 |
| CN106709215A (zh) * | 2017-02-21 | 2017-05-24 | 北京航空航天大学 | 一种基于级数展开的连续体结构非概率可靠性拓扑优化方法 |
| CN107633106A (zh) * | 2017-07-25 | 2018-01-26 | 东南大学 | 一种基于全局差分法的非均匀温度场下热模态灵敏度分析方法 |
| CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
| CN109063357A (zh) * | 2018-08-15 | 2018-12-21 | 大连理工大学 | 基于拓扑优化的自适应多步变域的汽车构件焊点布局优化方法 |
| CN109670200A (zh) * | 2018-11-13 | 2019-04-23 | 华中科技大学 | 一种等几何材料密度场结构拓扑优化方法 |
| CN114722463A (zh) * | 2022-03-28 | 2022-07-08 | 大连理工大学 | 一种基于各向异性材料场叠加的薄壁结构加筋设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| 周克民: "结构拓扑优化的一些基本概念和研究方法", 《力学与实践》 * |
| 徐涛等: "基于预条件Lanczos算法的结构拓扑修改静态重分析方法", 《吉林大学学报(工学版)》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020215533A1 (zh) * | 2019-04-26 | 2020-10-29 | 大连理工大学 | 一种基于材料场缩减级数展开的结构拓扑优化方法 |
| CN110609975B (zh) * | 2019-09-20 | 2023-02-17 | 大连理工大学 | 一种基于Kriging代理模型的并行计算全局优化方法 |
| CN110609975A (zh) * | 2019-09-20 | 2019-12-24 | 大连理工大学 | 一种基于Kriging代理模型的并行计算全局优化算法 |
| CN110852011B (zh) * | 2019-11-08 | 2022-09-20 | 大连理工大学 | 一种基于序列Kriging代理模型的结构非梯度拓扑优化方法 |
| CN110852011A (zh) * | 2019-11-08 | 2020-02-28 | 大连理工大学 | 一种基于序列Kriging代理模型的结构非梯度拓扑优化方法 |
| CN111027110A (zh) * | 2019-11-27 | 2020-04-17 | 中国科学院光电技术研究所 | 一种连续体结构拓扑与形状尺寸综合优化方法 |
| CN111090942A (zh) * | 2019-12-18 | 2020-05-01 | 浙江大学 | 基于拓扑优化的高灵敏度压阻式单轴力传感器设计方法 |
| CN112182941A (zh) * | 2020-10-23 | 2021-01-05 | 南京大学 | 一种间隔收缩隔声结构拓扑优化方法 |
| CN112182941B (zh) * | 2020-10-23 | 2023-11-24 | 南京大学 | 一种间隔收缩隔声结构拓扑优化方法 |
| CN112257319A (zh) * | 2020-10-26 | 2021-01-22 | 大连理工大学 | 一种基于非梯度拓扑优化的声学超材料设计方法 |
| CN112257319B (zh) * | 2020-10-26 | 2022-09-20 | 大连理工大学 | 一种基于非梯度拓扑优化的声学超材料设计方法 |
| WO2022087770A1 (zh) * | 2020-10-26 | 2022-05-05 | 大连理工大学 | 一种基于非梯度拓扑优化的声学超材料设计方法 |
| CN112818576B (zh) * | 2021-01-28 | 2024-04-19 | 华中科技大学 | 一种曲线纤维复合结构设计多层级优化方法 |
| CN112818576A (zh) * | 2021-01-28 | 2021-05-18 | 华中科技大学 | 一种曲线纤维复合结构设计多层级优化方法 |
| CN113806982A (zh) * | 2021-09-17 | 2021-12-17 | 北京航空航天大学杭州创新研究院 | 变结构可穿戴柔性热电器件基底拓扑优化方法 |
| CN113806982B (zh) * | 2021-09-17 | 2023-08-29 | 北京航空航天大学杭州创新研究院 | 变结构可穿戴柔性热电器件基底拓扑优化方法 |
| CN115630556A (zh) * | 2022-12-07 | 2023-01-20 | 东南大学 | 一种基于顶点法的电机拓扑优化方法 |
| CN117133385A (zh) * | 2023-07-27 | 2023-11-28 | 中南大学 | 一种弹性超材料非梯度多目标拓扑优化方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20210073428A1 (en) | 2021-03-11 |
| AU2019442319B2 (en) | 2022-01-13 |
| AU2019442319A1 (en) | 2021-04-08 |
| KR20210047351A (ko) | 2021-04-29 |
| CN110110413B (zh) | 2022-11-18 |
| EP3754521A4 (en) | 2021-01-27 |
| WO2020215533A1 (zh) | 2020-10-29 |
| JP2021516801A (ja) | 2021-07-08 |
| EP3754521A1 (en) | 2020-12-23 |
| JP6951001B2 (ja) | 2021-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
| CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
| CN111737839B (zh) | 基于动态进化率和自适应网格的beso拓扑优化方法及其应用 | |
| CN113191040B (zh) | 一种考虑结构稳定性的单材料结构拓扑优化方法和系统 | |
| CN111737835B (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN109344524B (zh) | 一种薄板结构加强筋分布优化方法 | |
| CN113158527B (zh) | 一种基于隐式fvfd计算频域电磁场的方法 | |
| CN107273613A (zh) | 一种基于应力惩罚和自适应体积的结构拓扑优化设计方法 | |
| CN109190233A (zh) | 一种结构拓扑优化方法 | |
| CN107729648A (zh) | 一种基于Shepard插值的曲线纤维复合结构设计瀑布型多级优化方法 | |
| CN112613206B (zh) | 一种基于各向异性体调和场的附面层网格生成方法 | |
| CN113887095B (zh) | 一种基于等几何分析的渐进式结构拓扑优化方法 | |
| CN113204906A (zh) | 一种考虑结构稳定性的多相材料拓扑优化设计方法和系统 | |
| CN115688276A (zh) | 一种基于离散伴随方法的飞行器外形自动化优化方法、系统、设备、介质 | |
| CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
| CN108763777B (zh) | 基于泊松方程显式解的vlsi全局布局模型建立方法 | |
| CN113779802A (zh) | 基于无网格efgm和等几何分析耦合方法的结构拓扑优化技术 | |
| CN113051796B (zh) | 一种应用于增材制造的结构拓扑优化设计方法 | |
| CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
| CN117556672A (zh) | 面向智能制造中结构应力最小化设计的高效拓扑优化方法 | |
| CN113515824B (zh) | 一种筋条布局与基板形状协同的拓扑优化设计方法 | |
| CN111079279A (zh) | 一种多构型晶格结构多尺度拓扑优化设计方法 | |
| CN117454502A (zh) | 基于Wolfe准则修正牛顿法的三维网格质量优化方法 | |
| CN116842799A (zh) | 用于载具轻量化设计的连续体多材料结构拓扑优化方法 | |
| CN105608267A (zh) | 一种多变量全局优化算法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240626 Address after: Unit 5A091, Office Building 1, Gaobei Sixteen B2, No. 16 Gaobei North Sixth Road, Songpingshan Community, Xili Street, Nanshan District, Shenzhen City, Guangdong Province, 518057 Patentee after: Shenzhen Hongtuo Weishe Technology Co.,Ltd. Country or region after: China Address before: 116024 No. 2 Ling Road, Ganjingzi District, Liaoning, Dalian Patentee before: DALIAN University OF TECHNOLOGY Country or region before: China |