CN108629140B - 一种基于测地线距离的带孔复合材料结构设计优化方法 - Google Patents
一种基于测地线距离的带孔复合材料结构设计优化方法 Download PDFInfo
- Publication number
- CN108629140B CN108629140B CN201810456168.1A CN201810456168A CN108629140B CN 108629140 B CN108629140 B CN 108629140B CN 201810456168 A CN201810456168 A CN 201810456168A CN 108629140 B CN108629140 B CN 108629140B
- Authority
- CN
- China
- Prior art keywords
- design
- composite material
- fiber
- geodesic distance
- material structure
- 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.)
- Active
Links
Images
Classifications
-
- 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
- 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)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Moulding By Coating Moulds (AREA)
Abstract
本发明属于复合材料结构设计优化领域,并具体公开了一种基于测地线距离的带孔复合材料结构设计优化方法,包括如下步骤:在带孔复合材料结构的实体区域内均匀定义离散设计点,设定各设计点处纤维角度的初始值;对结构进行网格划分,确定各单元中心点的坐标;计算各设计点到各单元中心点的测地线距离,基于测地线距离求得各单元中心点处的纤维角度;建立刚度矩阵,并计算得到整体位移向量;根据整体位移向量得到敏度,然后更新纤维角度θi直至满足优化终止条件。本发明通过计算单元中心点到设计点之间的测地线距离,并基于测地线距离得到带孔复合材料结构的连续变角度纤维布局,改善了带孔复合材料结构的连续变角度纤维布局,提升了结构性能。
Description
技术领域
本发明属于复合材料结构设计优化领域,更具体地,涉及一种基于测地线距离的带孔复合材料结构设计优化方法。
背景技术
与常规的定刚度复合材料结构相比,曲线纤维复合结构具有极高的可设计性。除了常规的结构外形设计外,复合材料结构的设计还包括纤维布局的设计。例如,通过设计纤维的铺放角度使结构具有变刚度的特性,进而获得更优的结构性能。随着自动纤维铺放和增材制造技术的发展,曲线纤维复合结构逐渐应用于航空航天、汽车、高端装备等领域。
对于复合材料结构中纤维角度布局的设计,由于要保证纤维铺设角的空间连续性以便于加工制造,所以连续变角度纤维布局描述模型是必不可少的。目前已有的优化设计方法中,基于Shepard插值的纤维增强复合材料结构优化方法能保证纤维角度的空间连续变化布局,例如CN107590325A公开的一种基于Shepard插值的纤维增强复合材料结构优化方法,CN107729648A公开的一种基于Shepard插值的曲线纤维复合结构设计瀑布型多级优化方法。但上述方法涉及的Shepard插值基于欧氏距离,只能解决不含孔洞的复合材料结构设计,而未考虑带孔的复合材料结构设计问题。在实际的工程结构中,大部分复合材料结构带有孔洞,或称为非凸形状。为解决该类复合材料结构的纤维角度优化问题,需研究设计一种新的优化方法,以实现带孔洞复合材料结构的纤维角度的优化。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于测地线距离的带孔复合材料结构设计优化方法,其通过计算单元中心点到设计点之间的测地线距离,并基于测地线距离得到带孔复合材料结构的连续变角度纤维布局,解决了现有的基于欧氏距离的Shepard插值方法的局限性,改善了带孔复合材料结构的连续变角度纤维布局,提升了结构性能。
为实现上述目的,本发明提出了一种基于测地线距离的带孔复合材料结构设计优化方法,其包括如下步骤:
S1在带孔复合材料结构的实体区域内均匀定义一系列离散设计点pi,并设定各设计点pi处纤维角度θi的初始值,其中,i=1,2,...,n,n为正整数;
S2对带孔复合材料结构进行网格划分,提取节点数据和单元数据,确定各单元中心点xc的坐标;
S4建立依赖于纤维角度θe的刚度矩阵,并计算得到整体位移向量;
S5根据整体位移向量得到敏度,然后更新纤维角度θi直至满足优化终止条件,以此完成带孔复合材料结构的设计优化。
作为进一步优选的,测地线距离||xc-pi||g采用如下方式获得:
2)利用二维线性插值根据步骤1)计算得到的各单元中心点xc和对应背景网格节点x间的测地线距离d计算单元中心点xc和设计点pi之间的测地线距离||xc-pi||g。
作为进一步优选的,速度场F(x)具体为:
作为进一步优选的,测地线距离||xc-pi||g具体采用如下公式计算:
其中,分别表示从单元中心点xc到四个背景网格节点A1、A2、A3、A4的测地线距离,这四个背景网格节点分别分布在设计点pi左上、左下、右上、右下的位置,(a1,b1)为A1的坐标,(a1,b2)为A2的坐标,(a2,b1)为A3的坐标,(a2,b2)为A4的坐标,(a,b)为pi的坐标。
作为进一步优选的,更新纤维角度θi直至满足优化终止条件具体为:
(1)更新纤维角度θi:
(2)根据更新后的纤维角度计算目标函数值,判断目标函数值是否满足收敛条件,若否,返回步骤(1)再次更新纤维角度,若是,则结束,此时对应的纤维角度即为优化结果,以此完成带孔复合材料结构的设计优化。
作为进一步优选的,所述的收敛条件为:
其中,ck为第k次更新后的目标函数值,ck-1为第k-1次更新后的目标函数值,k为纤维角度更新的次数。
其中,θmin和θmin分别为θi的上下界,θmin=0,θmax=θmin+π-ε,ε=1×10-8。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:
1.本发明改变现有优化方法,通过计算单元中心点到设计点之间的测地线距离,并基于测地线距离得到带孔复合材料结构的连续变角度纤维布局,解决了现有的基于欧氏距离的Shepard插值方法的局限性,改善了带孔复合材料结构的连续变角度纤维布局,提升了结构性能。
2.由于带孔复合材料结构孔洞的存在,使得现有的基于欧氏距离的优化方法不在适用于该结构,因为现有的方法忽视了孔洞两侧单元中心点处纤维角度和设计点处纤维角度的不相关性,导致孔洞两侧的纤维分布不合理而使设计的结构刚度偏小,本发明通过研究获得了一种新的优化方法,可适用于带孔复合材料结构的设计优化,该方法相比于现有方法,可将设计的结构刚度大大提升,提高约30%,并且孔周围的纤维角度分布更合理。
3.本发明方法通过划分背景网格并定义速度场,以区分实体区域与孔洞区域,然后基于上述区分对应计算出各区域中设计点到单元中心点的测地线距离,可以较好地解决带孔复合材料结构的纤维角度设计问题,改善孔周围的纤维分布,提升结构的力学性能,可获得最优的连续变角度纤维布局。
附图说明
图1是本发明的基于测地线距离的带孔复合材料结构设计优化方法的流程图;
图2是本发明较佳实施方式提供的平面悬臂梁结构优化设计例子示意图;
图3是图2中的优化设计例子关于设计点处的纤维角度的优化结果;
图4是图2中的优化设计例子关于有限单元中心点处的纤维角度的优化结果;
图5是基于欧氏距离优化方法得到的有限单元中心点处纤维角度的优化结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明实施例提供的一种基于测地线距离的带孔复合材料结构设计优化方法,包括如下步骤:
S1在带孔复合材料结构的实体区域内均匀定义一系列离散设计点pi,并设定各设计点pi处纤维角度θi的初始值,其中,i=1,2,...,n,n为正整数;
S2对带孔复合材料结构进行网格划分,提取节点数据和单元数据,确定各单元中心点xc的坐标;
S4建立依赖于纤维角度θe的刚度矩阵,并计算得到整体位移向量;
S5根据整体位移向量得到敏度,然后更新纤维角度θi直至满足优化终止条件,以此完成带孔复合材料结构的设计优化。
以下通过具体实施例对本发明进行详细说明,本实施例中,以带有均布载荷的悬臂梁柔度最小化的优化问题为例来解释本发明。如图2所示,在给定的1m×2m设计域内设置纤维角度初始值为90°,区域左侧边界施加位移约束,垂直于梁方向的面内载荷f=1均匀地施加在区域上边界右侧0.5m范围内,对带孔的纤维增强悬臂梁结构进行纤维角度布局优化,使其柔度最小化,包括如下步骤:
S1在带孔复合材料结构的实体区域内均匀定义一系列离散设计点pi,并设定各设计点pi处纤维角度θi的初始值,其中,i=1,2,...,200,设计点均匀排布,如图3所示,设计为10×20,给定pi处纤维角度θi的初始值为90°,即纤维与水平线的夹角为90°;
S2对带孔复合材料结构进行网格划分,提取节点数据和单元数据,确定各单元中心点xc的坐标,节点数为807,单元数Ne为700;
具体的,如图4所示,可在ANSYS有限元分析软件中选用四边形单元对带孔结构划分有限元网格,提取节点数据和单元数据,其中节点数据包括节点编号、节点坐标,单元数据包括单元编号和每个单元对应的四个节点的编号,根据每个单元对应节点的横坐标平均值和纵坐标平均值,即可得到各单元中心点xc的横纵坐标,各单元中心点xc的横坐标等于该单元对应的所有节点的横坐标的平均值,各单元中心点xc的纵坐标等于该单元对应的所有节点的纵坐标的平均值,其中可以以纤维的铺设平面为二维坐标平面,以铺设平面上任意一点为坐标原点,以从左往右的水平线为X轴,以从前往后的与X轴垂直的水平线为Y轴,具体可根据实际需要进行设计。
S3计算各单元中心点xc到各设计点pi的测地线距离||xc-pi||g,基于测地线距离||xc-pi||g求得各单元中心点xc处的纤维角度其中,Ixc为单元中心点xc的影响域内所有设计点索引集合,影响域半径Rc的选取要保证至少有一个设计点pi在单元中心点xc的影响域内,优选的Rc=0.3以满足要求。由于单元中心点xc和设计点pi的位置坐标固定,在优化迭代前计算得到的测地线距离||x-pi||g和权函数wi(x)均被保存,以节省优化时间。
具体的,测地线距离||xc-pi||g采用如下方式计算获得:
在带孔复合材料结构的实体区域内划分正方形的背景网格,正方形网格的各顶点称为节点,x为背景网格节点的坐标,x∈Ω表示节点在实体区域内,表示节点在孔洞内,其中Ω代表结构的实体区域;在该网格背景下,定义速度场F(x),表示如下:
F(x)与d之间可采用Eikonal非线性偏微分方程进行表示:
其中,d表示从单元中心点xc到背景网格节点x之间的测地线距离,表示梯度,|·|表示欧几里德范数,当单元中心点xc与背景网格节点x重合时d=0,通过求解非线性偏微分方程即可求解获得测地线距离d,而对于非线性偏微分方程的求解方法而言,其可以采用现有的常规求解方法进行求解,例如可利用快速推进法求解,其为现有技术,在此不赘述,具体求解方法本发明不做限定,只需能求解出非线性偏微分方程中待求解的参数即可。例如,单元中心点xc为(1,0.65),背景网格节点x为(1,0.35),分布在实体区域内,则F(x)=1,利用快速推进法求解方程得d=1.175;再如,单元中心点xc为(1,0.65),背景网格节点x为(1,0.5),分布在孔洞区域内,则F(x)=1×10-8,利用快速推进法求解方程得d=2.684×106。
利用二维线性插值计算单元中心点xc和设计点pi之间的测地线距离||xc-pi||g,具体计算如下:
其中,d1、d2、d3、d4分别表示从单元中心点xc到四个背景网格节点A1、A2、A3、A4的测地线距离,其通过求解非线性偏微分方程获得,这四个背景网格节点分别分布在设计点pi左上、左下、右上、右下位置,例如设计点pi分布在背景网格中,那么该背景网格的四个节点即为这四个节点,再如设计点pi分布在背景网格上,那么设计点pi左上、左下、右上、右下位置对应的与设计点pi相邻的四个节点即为这四个节点,其中,A1的坐标为(a1,b1)、A2的坐标为(a1,b2)、A3的坐标为(a2,b1)、A4的坐标为(a2,b2),pi的坐标为(a,b)。
S4建立依赖于纤维角度θe的刚度矩阵,并计算得到整体位移向量;
S5根据整体位移向量得到敏度,然后更新纤维角度θi直至满足优化终止条件。
其中,步骤S4和S5采用的计算方法与背景技术中介绍的现有专利的计算方法一致,但输入的原始数据存在本质区别,本发明的原始数据是测地线距离,而现有专利的原始数据是欧氏距离,因此数据处理细节存在区别。
具体的,步骤S4中,利用步骤S2已划分的有限元网格,单元数Ne为700,在每个单元e(e=1,2,3,......,700)上建立依赖于单元中心处纤维角度θe的刚度矩阵Ke(θe),其中:
式中,B为位移应变矩阵,D(θe)是依赖于单元中心处纤维角度的单元弹性矩阵,A为结构的面积,Ae为单元e的面积;
单元弹性矩阵D(θe)计算如下:
D(θe)=T(θe)D0T(θe)T
式中,D0为纤维未旋转时原始弹性矩阵,T(θe)为旋转矩阵:
其中,Ex和Ey是杨氏模量,Gxy是剪切模量,νxy和νyx是泊松比,满足条件νxyEy=νyxEx;
再通过组装单元刚度矩阵Ke(θe)得到整体刚度矩阵K,其具体操作就是把所有单元刚度矩阵Ke(θe)扩大为整体刚度矩阵同阶的方块阵将单元刚度矩阵的元素存放至单元扩大矩阵中对应节点整体编号的行和列,其余元素为零,扩大后的矩阵也称为单元贡献矩阵,将所有单元贡献矩阵相加,即得到整体刚度矩阵,计算公式如下:
根据Ku=f计算得到整体位移向量u,其中,f为外力向量,可根据实际需要进行选择与设定。
式中,ue为单元e对应的位移向量,单元e的编号为ne,则ue为整体位移向量u第ne行的元素。
利用基于敏度共轭映射的优化算法更新设计变量θi直到收敛,其更新格式如下:
其中,e为自然常数;
其中,θmin和θmin分别为θi的上下界,θmin=0,θmax=θmin+π-ε,ε=1×10-8,其为一个极小正值,用来消除插值时由于三角函数计算造成的数值错误;
获得更新后的设计变量(即纤维角度),例如第k+1次更新后的θi值为然后通过有限元分析得到设计变量第k+1次更新后的整体位移向量uk+1,利用公式ck+1=fTuk+1得到ck+1,直至满足如下收敛准则,完成整个优化过程:
即当满足上述三个条件中的任意一个条件即完成优化,其中,ck为第k次更新后的目标函数值,ck-1为第k-1次更新后的目标函数值,k为纤维角度更新的次数。
综上,带孔复合材料结构的优化设计模型可以概括为:
findθi(i=1,2,…,200)
min c=fTu
s.t.Ku=f
θmin≤θi≤θmax
其中,离散设计点处的纤维角度值θi为设计变量,目标函数为柔度c,设计目标使结构的柔度c最小化,设计约束为平衡方程Ku=f和θi的上下界θmax、θmin。
本发明的优化结果如下:基于测地线距离的带孔复合材料结构设计点处的细微角度优化结果如图3所示,有限单元中心点处的纤维角度优化结果如图4所示,对应的结构柔度值为105.29;而采用现有的基于欧式距离的优化方法对带孔复合材料结构进行优化,其在有限单元中心点处的纤维角度优化结果如图5所示,对应的结构柔度值为151.64。比较本发明与现有的优化方法,可以看出本发明方法得到的结构刚度更大(柔度更小),约提高30%,并且孔周围的纤维角度分布更合理,因此本发明的方法更适用于带孔复合材料结构的设计优化。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于测地线距离的带孔复合材料结构设计优化方法,其特征在于,包括如下步骤:
S1在带孔复合材料结构的实体区域内均匀定义一系列离散设计点pi,并设定各设计点pi处纤维角度θi的初始值,其中,i=1,2,...,n,n为正整数;
S2对带孔复合材料结构进行网格划分,提取节点数据和单元数据,确定各单元中心点xc的坐标;
其中,测地线距离||xc-pi||g采用如下方式获得:
1)在带孔复合材料结构的实体区域内划分正方形的背景网格,定义速度场F(x),根据非线性偏微分方程|▽d|F(x)=1求解各单元中心点xc和对应背景网格节点x间的测地线距离d;
2)利用二维线性插值根据步骤1)计算得到的各单元中心点xc和对应背景网格节点x间的测地线距离d计算单元中心点xc和设计点pi之间的测地线距离||xc-pi||g;
S4建立依赖于纤维角度θe的刚度矩阵,并计算得到整体位移向量;
S5根据整体位移向量得到敏度,然后更新纤维角度θi直至满足优化终止条件,以此完成带孔复合材料结构的设计优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810456168.1A CN108629140B (zh) | 2018-05-14 | 2018-05-14 | 一种基于测地线距离的带孔复合材料结构设计优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810456168.1A CN108629140B (zh) | 2018-05-14 | 2018-05-14 | 一种基于测地线距离的带孔复合材料结构设计优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108629140A CN108629140A (zh) | 2018-10-09 |
CN108629140B true CN108629140B (zh) | 2020-07-10 |
Family
ID=63693018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810456168.1A Active CN108629140B (zh) | 2018-05-14 | 2018-05-14 | 一种基于测地线距离的带孔复合材料结构设计优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108629140B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109460577B (zh) * | 2018-10-10 | 2020-11-24 | 华中科技大学 | 一种保证制造质量的变刚度复合材料结构优化设计方法 |
CN110955941B (zh) * | 2019-11-29 | 2022-03-18 | 华中科技大学 | 基于向量场的复合材料结构优化设计方法及设备 |
CN111444579B (zh) * | 2020-03-11 | 2022-04-12 | 华中科技大学 | 一种考虑可制造性的复合材料结构优化设计方法 |
CN111723457B (zh) * | 2020-05-08 | 2022-06-17 | 华中科技大学 | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 |
CN112818576B (zh) * | 2021-01-28 | 2024-04-19 | 华中科技大学 | 一种曲线纤维复合结构设计多层级优化方法 |
CN114407350B (zh) * | 2022-01-11 | 2023-07-21 | 西北工业大学 | 连续纤维增强复合材料3d打印填充路径规划方法和装置 |
CN115876883B (zh) * | 2022-12-29 | 2024-03-29 | 南京航空航天大学 | 一种复合材料层合板分层损伤位置检测方法及检测系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101811580A (zh) * | 2010-01-18 | 2010-08-25 | 中国航空工业集团公司北京航空制造工程研究所 | 一种复合材料结构中曲面铺层的下料方法 |
CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
-
2018
- 2018-05-14 CN CN201810456168.1A patent/CN108629140B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101811580A (zh) * | 2010-01-18 | 2010-08-25 | 中国航空工业集团公司北京航空制造工程研究所 | 一种复合材料结构中曲面铺层的下料方法 |
CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
Non-Patent Citations (3)
Title |
---|
"A cascadic multilevel optimization algorithm for the design of composite structures with curvilinear fiber based on Shepard interpolation";Qi Xia等;《Composite Structures》;20180104;第209-219页 * |
"Finding shortest paths on surfaces using level sets propagation";Ron Kimmel等;《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》;19950630;第17卷(第6期);第635-640页 * |
"两种不同铝内衬碳纤维缠绕复合材料气瓶受力状况的对比分析";周丹 等;《玻璃钢/复合材料》;20160630(第6期);第76-80页第1-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN108629140A (zh) | 2018-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108629140B (zh) | 一种基于测地线距离的带孔复合材料结构设计优化方法 | |
CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
CN110941924B (zh) | 一种多组件系统集成一体化的多尺度拓扑优化设计方法 | |
CN110955941B (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
CN111723457B (zh) | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 | |
CN112765856B (zh) | 一种功能梯度多孔结构拓扑优化的混合水平集方法 | |
CN108228970B (zh) | 结构动力学分析显式异步长并行计算方法 | |
CN107729648A (zh) | 一种基于Shepard插值的曲线纤维复合结构设计瀑布型多级优化方法 | |
CN112446163B (zh) | 基于参数化水平集的能量有限元拓扑优化方法 | |
CN109670207B (zh) | 一种面向多种多孔材料结构的动力学一体化设计方法 | |
CN109460577B (zh) | 一种保证制造质量的变刚度复合材料结构优化设计方法 | |
CN111125963A (zh) | 基于拉格朗日积分点有限元的数值仿真系统及方法 | |
CN113204906A (zh) | 一种考虑结构稳定性的多相材料拓扑优化设计方法和系统 | |
CN111079327A (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
CN113505405B (zh) | 等效荷载获取方法、基于等效荷载的拓扑优化方法及系统 | |
CN110717208B (zh) | 一种基于连续梯度微结构的多尺度频率响应拓扑优化方法 | |
CN116362079B (zh) | 一种基于新型插值模型的多材料结构拓扑优化方法 | |
CN111127491B (zh) | 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备 | |
CN115630542B (zh) | 一种薄壁加筋结构的加筋布局优化方法 | |
CN111444579A (zh) | 一种考虑可制造性的复合材料结构优化设计方法 | |
CN116187074A (zh) | 基于等几何的各向异性周期结构材料多尺度拓扑优化方法 | |
CN113779825B (zh) | 一种基于单位分解及移动控制点的结构优化设计方法 | |
CN112364546B (zh) | 一种基于双边滤波的纤维增强复合材料结构优化设计方法 | |
CN115081113A (zh) | 一种基于多水平集函数参数化的网格加强筋结构优化方法 | |
CN113705060A (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 |