CN112417692B - 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 - Google Patents
基于载荷不确定性的材料结构多尺度拓扑优化设计方法 Download PDFInfo
- Publication number
- CN112417692B CN112417692B CN202011329988.8A CN202011329988A CN112417692B CN 112417692 B CN112417692 B CN 112417692B CN 202011329988 A CN202011329988 A CN 202011329988A CN 112417692 B CN112417692 B CN 112417692B
- Authority
- CN
- China
- Prior art keywords
- microstructure
- load
- adopting
- macrostructure
- design
- 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/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供基于载荷不确定性的材料结构多尺度拓扑优化设计方法,包括如下步骤:定义设计域、初始化设计参数,依次进行材料微观结构和宏观结构有限元分析,采用摄动法计算结构柔顺度的期望值和方差,构建基于载荷不确定的材料结构多尺度拓扑优化模型,灵敏度分析,利用映射过滤方法修正目标函数和约束灵敏度,采用OC算法更新设计变量,判断是否满足收敛准则,最后输出宏观拓扑型构和材料微观结构拓扑型构;本发明采用一阶摄动法对不确定性载荷进行量化,有效地降低统计性指标计算的复杂程度,提高计算效率,实现了宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计,获得的周期性多空复合材料结构具有更好的稳健性。
Description
技术领域
本发明涉及材料结构多尺度优化设计相关技术领域,特别是涉及基于载荷不确定性的材料结构多尺度拓扑优化设计方法。
背景技术
为了满足现代工业产品对“结构超轻量化”、“功能特殊化”与“性能集成化”的需求,同时考虑材料微观结构优化和宏观结构优化进行材料结构多尺度拓扑优化设计,获取最优的材料微观结构和宏观结构多尺度设计,获得周期性多空复合材料结构;它具有高比刚度、高比强度、吸能等优异的力学性能,广泛应用于航空航天、汽车工业、建筑行业等相关领域。
现有的材料结构多尺度拓扑优化设计方法是考虑确定性载荷工况,即载荷的大小与方向均为确定的,但在实际工程中不可避免存在载荷不确定性因素。因此,目前急需一种基于载荷不确定的材料结构多尺度拓扑优化设计方法,以解决上述问题。
发明内容
本发明的目的是提供基于载荷不确定的材料结构一体化拓扑优化方法,采用一阶摄动法对不确定性载荷进行量化,有效地减低统计性指标计算的复杂程度,提高计算效率,实现了宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计,获得的周期性多空复合材料结构具有更好的稳健性。
本发明提供基于载荷不确定性的材料结构多尺度拓扑优化设计方法,包括如下步骤:
S1、定义设计域,给定初始测试应变条件,设定宏观结构与材料微观结构的尺寸、边界、材料参数和不确定性载荷参数,设置设计变量的过滤半径,设置迭代准则,初始化宏观结构与材料微观结构的设计变量;
S2、对微观结构边界添加周期性边界条件,在给定的初始测试应变条件下进行微观结构有限元分析,计算材料微观结构的宏观等效弹性张量;
S3、采用摄动方法对不确定性载荷的参数进行量化,采用随机有限元法得到宏观结构的单元节点位移矢量,计算宏观结构柔顺度的期望和方差;
S4、将宏观结构柔顺度的期望值与方差进行加权得到加权和,以加权和最小化作为优化目标,以宏观结构和材料微观结构的体积份数作为约束,采用变密度法与能量均匀化方法构建考虑载荷不确定性的材料结构多尺度拓扑优化模型;
S5、根据所述材料结构多尺度拓扑优化模型,计算宏观结构与材料微观结构优化设计的目标函数和约束灵敏度,采用Heaviside映射过滤方法对目标函数和约束灵敏度进行修正,获得宏观结构和材料微观结构拓扑构型;
S6、采用优化准则法OC更新宏观结构与材料微观结构设计变量,判断是否满足收敛准则;若不满足,转至步骤S2,若满足,输出最优的宏观结构与材料微观结构拓扑型构;
优选地,所述S2包括以下步骤:
S2-1、在微观结构边界添加周期性边界条件,将宏观位移场和微观材料内的周期性波动场相加得到材料微观结构的位移场;
S2-2、在给定的初始测试应变条件下,将位移场分为四个部分,建立材料微观结构的有限元平衡方程进行有限元分析,并采用能量均匀化方法求得材料微观结构的宏观等效弹性张量;
优选地,所述S3包括以下步骤:
S3-1、采用摄动法,以一阶泰勒展开式表示具有不确定性大小向量和方向向量的作用载荷向量,对不确定性载荷的参数进行量化,得到不确定性量化参数;
S3-2、采用固体各向同性材料惩罚模型SIMP表示宏观结构的材料分布,由材料微观结构的宏观等效弹性张量建立宏观结构有限元平衡方程,采用随机有限元法,通过一阶泰勒展开式求得宏观结构的单元节点位移矢量;
S3-3、将不确定性量化参数和高阶函数的平均值代入结构柔顺度的不确定性表达式中,求得柔顺度的期望E(C)和方差Var(C);
优选地,所述S4构建基于载荷不确定的材料结构多尺度拓扑优化模型为如下方法:
find{x y}
s.t.K(x,y)ui=fi(i=0,1,2)
0≤yj≤1 j=1,...,n
为宏观结构的柔顺度,x和y分别为宏观设计变量和微观设计变量,统一表示为单元密度,ui为宏观结构的位移场,f0为载荷列阵的均值矩阵,f1为不确定载荷大小对随机变量的偏导数,f2为不确定载荷方向对随机变量的偏导数;α为加权系数,K(x,y)为宏观结构的总刚度矩阵,Ω0和Y0分别为宏观和微观结构的初始体积,Ωe和Ye分布为宏观和微观结构的单元体积,和分别为宏观和微观结构允许的材料体积比,N和n分别为宏观和微观结构的有限元的数目;
优选地,所述S5包括以下步骤:
S5-1、采用拉格朗日乘子法计算宏观结构柔顺度相对宏观设计变量的灵敏度,计算柔顺度对微观结构设计变量的灵敏度,并且求解宏观结构和材料微观结构体积约束灵敏度;
S5-2、采用基于Heaviside函数的映射过滤方法对目标函数和宏观结构材料微观结构约束灵敏度进行修正,首先采用密度过滤方法对单元密度进行过滤,通过相邻单元密度的加权平均方法,获得过滤后的单元密度;
S5-3、采用Heaviside映射过滤法对过滤后的单元密度进行修正,使单元密度值向0和1两端集中,获得宏观结构和材料微观结构拓扑构型;
优选地,所述S6采用优化准则法OC算法更新宏观结构与材料微观结构设计变量,求解拓扑优化问题,判断是否满足收敛准则,若不满足,转至步骤S2,若满足,则材料结构多尺度拓扑优化过程结束,输出最优的宏观结构和材料微观结构拓扑构型。
本发明公开了以下技术效果:
(1)本发明基于载荷不确定性的材料结构多尺度拓扑优化设计,采用一阶摄动法对不确定性载荷进行量化;采用变密度法与能量均匀化方法构建基于载荷不确定性的材料结构多尺度拓扑优化模型;采用伴随矩阵技术进行目标函数和约束灵敏度分析;获得最优的宏观结构拓扑和材料微观结构拓扑;
(2)采用一阶摄动法对不确定性载荷进行量化,计算考虑载荷大小和方向不确定条件下的柔顺度期望值和方差表达式,有效地减低统计性指标计算的复杂程度,提高计算效率;
(3)考虑载荷的方向和大小的不确定,进行宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计,获得的周期性多空复合材料结构,具有更好的稳健性。
附图说明
图1为本发明基于载荷不确定性的材料结构多尺度拓扑优化设计方法流程示意图;
图2为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的二维立柱结构设计域图;
图3为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型图;
图4为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的材料微观结构拓扑图;
图5为载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的3×3微观结构图;
图6为载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型;
图7为载荷不确定性的材料结构多尺度拓扑优化设计获得的材料微观结构拓扑;
图8为载荷不确定性的材料结构多尺度拓扑优化设计获得的3×3微观结构图。
具体实施方式
下面将结合本实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1所示,本发明提供基于载荷不确定性的材料结构多尺度拓扑优化设计方法,包括如下步骤:
S1、定义设计域,给定初始测试应变条件,设定宏观结构与材料微观结构的尺寸、边界条件、材料参数和不确定载荷的大小和方向参数,设置设计变量的过滤半径,设置迭代准则,初始化宏观结构与材料微观结构的设计变量;
S2、对微观结构边界添加周期性边界,给定的初始测试应变条件下进行微观结构有限元分析,计算材料微观结构的宏观等效弹性张量;
S2-1、在微结构边界施加周期性边界,将宏观位移场和微观材料内的周期性波动场相加得到微观结构的位移场:
求得材料微观结构所有法线方向上边界的位移场为:
k+与k-表示微观结构内一对相反的两个边界的法线方向,由上式相减消除材料微观结构内部的周期性波动位移场,可得:
S2-2、将位移场全局位移分为四个部分,建立材料微观结构的有限元平衡方程,进行有限元分析:
Kkl(k,l=1,2,3,4)表示材料微观结构刚度矩阵分块矩阵,表示给定的位移,χ2表示微观结构的内部节点的未知位移,χ3和χ4表示位于相反边界的节点的未知位移,为材料微观结构相反边界的位移差值,且满足其中F1为给定位移节点出的反力,F2=0,由周期假设有F3+F4=0,且Kij=Kji,则上式得:
采用能量均匀化方法求得材料微观结构的宏观等效弹性张量:
为材料微观结构的宏观等效弹性张量,Y0表示材料微观结构的面积,Yj表示微观材料结构有限元单元的面积,dY表示材料结构中的微元体积,T表示矩阵转置的符号,χ0(ij)表示由施加在微观结构上的初始测试应变引起的位移,χij表示为未知的位移,由初始测试应变引起的微观结构内部产生的位移;b为应变位移矩阵,C(y)为基材料的弹性矩阵,n为材料微观结构的有限元的数目;
S3、考虑载荷大小和方向不确定性,采用摄动方法对不确定性载荷的参数进行量化,采用随机有限元法求解宏观结构的单元节点位移矢量,计算宏观结构柔顺度的期望和方差;
S3-1、采用摄动法,将一阶泰勒展开式表示具有不确定性大小向量μ和方向向量θ的作用载荷向量F:
F(μ,θ)≈f0+f1·Δμ+f2·Δθ
Δμ和Δθ分别为载荷大小和方向随机变化量称为载荷不确定参数,f0为载荷列阵的均值矩阵,f1为不确定载荷大小对随机变量的偏导数,f2为不确定载荷方向对随机变量的偏导数;
S3-2、采用SIMP(Solid Isotropic Material with Penalization,固体各向同性材料惩罚模型)表示宏观结构的材料分布,并建立宏观结构有限元平衡方程:
KU(μ,θ)=F(μ,θ)
K为宏观结构的总装刚度矩阵,U为宏观结构的单元节点位移矢量,宏观结构的节点位移矢量通过一阶泰勒展开式求得为:
U(μ,θ)=u0+u1·Δμ+u2·Δθ
u0为位移的均值矢量,u1和u2分别位移矢量对载荷大小和方向随机变量的偏导数,由下式求得:
ui=K-1fi(i=0,1,2);
S3-3、将载荷不确定参数Δμ和Δθ高阶函数的平均值代入结构柔顺度C的不确定性表达式中,求得柔顺度的期望E(C)和方差Var(C):
E(Δμ)2为载荷大小的高阶函数的均值,E(Δθ)2为载荷方向的高阶函数的均值,u2为位移矢量对随机变量的偏导数:
E(C2)为柔顺度平方的均值,E(C)2柔顺度均值的平方,Var(Δμ)为不确定载荷大小的方差,Var(Δθ)为不确定载荷方向的方差,Var(ΔμΔθ)为不确定载荷大小与方向乘积的方差;
S4、以结构柔顺度的期望值与标准差进行加权得到加权和,以加权和最小化作为优化目标,以宏观结构和材料微观结构的体积份数作为约束,采用变密度法与能量均匀化方法构建基于载荷不确定的材料结构多尺度拓扑优化模型,构建基于载荷不确定的材料结构多尺度拓扑优化模型方法如下:
find{x y}
s.t.K(x,y)ui=fi(i=0,1,2)
0≤yj≤1 j=1,...,n,
为宏观结构的柔顺度,x和y分别为宏观设计变量和微观设计变量,ui为宏观结构的位移场,α为加权系数,K(x,y)为宏观结构的总刚度矩阵,Ω0和Y0分别为宏观和微观结构的初始体积,Ωe和Ye分别为宏观和微观结构的单元体积,和分别为宏观和微观结构允许的材料体积比,N和n分别为宏观和微观结构的有限元的数目;
S5、根据所述材料结构多尺度拓扑优化模型,计算宏观结构与材料微观结构优化设计的目标函数和约束灵敏度,采用Heaviside映射过滤方法对目标函数和约束灵敏度进行修正,获得清晰的宏观结构和材料微观结构拓扑构型;
S5-1、引入拉格朗日乘子计算宏观结构柔顺度约束灵敏度:
α为加权系数,λi(i=0,1,2)表示为拉格朗日乘子,结构柔顺度对宏观设计变量的约束灵敏度求得为:
其中A0=f0 T+fT 1·E(Δμ)+f2 T·E(Δθ),
A1=f0 T·E(Δμ)+fT 1·E(Δμ)2+f2 T·E(Δθ)E(Δμ),
A2=f0 T·E(Δθ)+fT 1·E(Δμ)E(Δθ)+f2 T·E(Δθ)2,
B0=f1·Var(Δμ)·f1 T+f2·Var(Δθ)·f2 T,
B1=f0·Var(Δμ)·f0 T+f1·Var(Δμ)2·f1 T+f2·Var(ΔμΔθ)·f2 T,
B2=f0·Var(Δθ)·f0 T+f1·Var(ΔμΔθ)·f1 T+f2·Var(Δθ)2·f2 T,
结构柔顺度的约束灵敏度公式为:
将拉格朗日乘子λi代入结构柔顺度的约束灵敏度公式得:
其中
宏观体积约束对宏观设计变量的约束灵敏度为:
结构柔顺度对微观结构设计变量的约束灵敏度为:
其中
微观体积约束对微观设计变量的约束灵敏度为:
5-2、采用基于Heaviside函数的映射过滤方法对目标函数和宏观结构、材料微观结构的约束灵敏度进行修正,首先采用密度过滤方法对单元密度进行过滤,通过相邻单元密度xj的加权平均,获得过滤后的单元密度x%e:
Ωe表示与单元e质心距离小于过滤半径rmin的所有单元的设计变量集合,wj为权重系数,其表达式为
w(xj)=R-|xj-xi|
R为最小过滤半径,xj和xi分别表示单元j和单元i的坐标;
5-3、采用密度过滤方法的拓扑优化结果中仍存在灰度单元,采用Heaviside映射过滤法对过滤后的单元密度进行修正,使单元的密度值向0和1两端集中,得到修正后的单元密度,获取清晰的宏观结构和材料微观结构拓扑构型;
S6、采用优化准则法OC算法更新宏观结构与材料微观结构设计变量,求解拓扑优化问题,迭代至满足收敛为止,获得最优的宏观结构拓扑和材料微观结构拓扑,获取宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计,获得的周期性多空复合材料结构具有更好的稳健性,判断是否满足收敛准则,若不满足,转至步骤S2,若满足,则材料结构多尺度拓扑优化过程结束,输出最优的宏观结构和材料微观结构拓扑构型。
为了进一步验证本发明基于载荷不确定性的材料结构多尺度拓扑优化设计方法的有效性,本实施例以二维立柱结构为例对本发明进行解释;
立柱结构的设计域、边界条件如图2所示,宏观结构的设计域尺寸长×高为500mm×250mm,设计域下端为固定边界,载荷作用F在宏观结构上边界的中心处,所用材料的弹性模量E为1Gpa,泊松比v为0.3,宏观设计域离散为50×50个平面四边形单元,材料微观结构的设计域长×高为0.01mm×0.01mm,材料微观结构的设计域50×50个平面四边形单元,载荷F的大小为1kN,作用方向与水平方向的夹角为-π/2(垂直向下),载荷大小标准差为0,载荷方向标准差为0.3。宏观结构和材料微观结构的允许的体积份数为25%,最小过滤半径rmin为2倍单元尺寸,惩罚系数P为3。
图3、图4和图5分别为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型、材料微观结构拓扑和3×3微观结构图,立柱宏观结构的材料主要分布在竖直载荷线附近,形成一种柱状结构,以承受宏观结构的竖直方向的承载,材料微观结构的材料主要分布在两根竖直柱状结构。图6、图7和图8分别为基于载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型、材料微观结构拓扑和3×3微观结构图。与确定性载荷条件下的材料结构多尺度拓扑优化结果相比,基于载荷不确定性的材料结构多尺度拓扑优化获得的立柱宏观结构底部出现的一定角度的开口,以此来提高立柱结构的水平方向的承载能力;基于载荷不确定性的材料结构多尺度拓扑优化获得的材料微观结构拓扑构型形成更多网格结构,材料分布会更加均匀,并非集中在竖直的两根柱状结构附近,这种材料微观结构使得立柱宏观结构具有更好水平方向的承载能力。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (5)
1.基于载荷不确定性的材料结构多尺度拓扑优化设计方法,包括如下步骤:
S1、定义设计域,给定初始测试应变条件,设定宏观结构与材料微观结构的尺寸、边界、材料参数和不确定性载荷参数,设置设计变量的过滤半径,设置迭代准则,初始化宏观结构与材料微观结构的设计变量;
S2、对材料微观结构边界添加周期性边界条件,在给定的初始测试应变条件下进行所述材料微观结构有限元分析,计算材料微观结构的宏观等效弹性张量;
S3、采用摄动方法对所述不确定性载荷参数进行量化,采用随机有限元法得到所述宏观结构的单元节点位移矢量,计算宏观结构柔顺度的期望值和方差;
S4、将所述宏观结构柔顺度的期望值与方差进行加权得到加权和,以所述加权和最小化作为优化目标,以宏观结构和材料微观结构的体积份数作为约束,采用变密度法与能量均匀化方法构建基于载荷不确定性的材料结构多尺度拓扑优化模型;
S5、根据所述材料结构多尺度拓扑优化模型,计算宏观结构与材料微观结构优化设计的目标函数和约束灵敏度,采用Heaviside映射过滤方法对所述目标函数和约束灵敏度进行修正,获得宏观结构和材料微观结构拓扑构型;
S6、采用优化准则法OC更新宏观结构与材料微观结构设计变量,判断是否满足收敛准则;若不满足,转至步骤S2,若满足,输出最优的宏观结构与材料微观结构拓扑构型;
所述S3包括以下步骤:
S3-1、采用摄动法,以一阶泰勒展开式表示具有不确定性大小向量和方向向量的作用载荷向量,对所述不确定性载荷参数进行量化,得到不确定性量化参数;
S3-2、采用固体各向同性材料惩罚模型SIMP表示宏观结构的材料分布,由材料微观结构的宏观等效弹性张量建立宏观结构有限元平衡方程,采用随机有限元法,通过一阶泰勒展开式求得宏观结构的单元节点位移矢量;
2.根据权利要求1所述基于载荷不确定性的材料结构多尺度拓扑优化设计方法,其特征在于,所述S2包括以下步骤:
S2-1、在微观结构边界添加周期性边界条件,将宏观位移场和微观材料内的周期性波动场相加得到材料微观结构的位移场;
S2-2、在给定的初始测试应变条件下,将位移场分为四个部分,建立材料微观结构的有限元平衡方程进行有限元分析,并采用能量均匀化方法求得材料微观结构的宏观等效弹性张量。
4.根据权利要求1所述基于载荷不确定性的材料结构多尺度拓扑优化设计方法,其特征在于,所述S5包括以下步骤:
S5-1、采用拉格朗日乘子法计算宏观结构柔顺度相对宏观设计变量的灵敏度,计算柔顺度对微观结构设计变量的灵敏度,并且求解宏观结构和材料微观结构体积约束灵敏度;
S5-2、采用基于Heaviside函数的映射过滤方法对目标函数和宏观结构材料微观结构约束灵敏度进行修正,首先采用密度过滤方法对单元密度进行过滤,通过相邻单元密度的加权平均方法,获得过滤后的单元密度;
S5-3、采用Heaviside映射过滤法对过滤后的单元密度进行修正,使单元密度值向0和1两端集中,获得宏观结构和材料微观结构拓扑构型。
5.根据权利要求1所述基于载荷不确定性的材料结构多尺度拓扑优化设计方法,其特征在于,所述S6采用优化准则法OC算法更新宏观结构与材料微观结构设计变量,求解拓扑优化问题,判断是否满足收敛准则,若不满足,转至步骤S2,若满足,则材料结构多尺度拓扑优化过程结束,输出最优的宏观结构和材料微观结构拓扑构型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011329988.8A CN112417692B (zh) | 2020-11-24 | 2020-11-24 | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011329988.8A CN112417692B (zh) | 2020-11-24 | 2020-11-24 | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112417692A CN112417692A (zh) | 2021-02-26 |
CN112417692B true CN112417692B (zh) | 2022-08-12 |
Family
ID=74778182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011329988.8A Active CN112417692B (zh) | 2020-11-24 | 2020-11-24 | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112417692B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113094944B (zh) * | 2021-03-22 | 2024-04-19 | 华中科技大学 | 一种微流道散热器及其细观尺度拓扑优化设计方法 |
CN113343512B (zh) * | 2021-05-11 | 2024-05-14 | 华中科技大学 | 一种基于Mobile-U-Net的多尺度拓扑优化设计方法 |
CN113536623B (zh) * | 2021-06-24 | 2022-03-25 | 河海大学 | 一种材料不确定性结构稳健性拓扑优化设计方法 |
CN114595617B (zh) * | 2022-05-10 | 2022-07-22 | 山东科技大学 | 一种微脉冲激光雷达光路参数拓扑结构的设计方法 |
CN117454705A (zh) * | 2023-11-07 | 2024-01-26 | 北京航空航天大学 | 机翼结构/材料多尺度气动弹性优化方法、装置及介质 |
CN117373579B (zh) * | 2023-11-09 | 2024-05-17 | 哈尔滨理工大学 | 一种时域动载荷作用下多微结构多尺度并行拓扑优化方法及系统 |
CN117892414B (zh) * | 2024-03-18 | 2024-05-24 | 华东交通大学 | 简谐激励下的各向异性材料结构稳健性拓扑优化设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844917A (zh) * | 2017-01-10 | 2017-06-13 | 西北工业大学 | 一种基于支反力方差约束的机床底座拓扑优化设计方法 |
CN107341316A (zh) * | 2017-07-13 | 2017-11-10 | 西北工业大学 | 设计相关压力载荷作用下的结构形状‑拓扑联合优化方法 |
CN111310377A (zh) * | 2020-02-21 | 2020-06-19 | 北京航空航天大学 | 一种基频和频率间隔混合约束下的连续体结构非概率可靠性拓扑优化设计方法 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412997B (zh) * | 2013-08-14 | 2016-04-13 | 北京航空航天大学 | 一种不确定载荷下的桁架结构稳健设计方法 |
US10391530B2 (en) * | 2013-09-27 | 2019-08-27 | The Regents Of The University Of California | Liquid-repellent surfaces made of any materials |
CN104821596A (zh) * | 2015-05-21 | 2015-08-05 | 国家电网公司 | 一种基于内模控制的混合直流输电系统拓扑及控制方法 |
EP3379434B1 (en) * | 2017-03-22 | 2022-09-28 | Tata Consultancy Services Limited | A system and method for design of additively manufactured products |
CN107577845B (zh) * | 2017-08-02 | 2018-09-14 | 北京航空航天大学 | 一种基于非概率可靠性的确定骨架结构最优拓扑形式的方法 |
CN107563010B (zh) * | 2017-08-08 | 2020-09-25 | 西北工业大学 | 基于形状特征的多尺度结构材料一体化设计方法 |
CN107609320B (zh) * | 2017-10-30 | 2019-02-15 | 西安科技大学 | 一种桁架非概率可靠性形状优化设计方法 |
CN107844676B (zh) * | 2017-12-18 | 2021-07-09 | 北京工业大学 | 一种基于多性能约束的结构拓扑优化设计方法 |
US11132478B2 (en) * | 2018-02-01 | 2021-09-28 | Toyota Motor Engineering & Manufacturing North America, Inc. | Methods for combinatorial constraint in topology optimization using shape transformation |
US10984593B2 (en) * | 2018-10-31 | 2021-04-20 | The Hong Kong University Of Science And Technology | Methods of high-definition cellular level set in B-splines for modeling and topology optimization of three-dimensional cellular structures |
CN109670207B (zh) * | 2018-11-22 | 2021-04-06 | 华中科技大学 | 一种面向多种多孔材料结构的动力学一体化设计方法 |
CN109583091B (zh) * | 2018-11-30 | 2022-11-29 | 长沙理工大学 | 基于自适应约束的柔性机构拓扑优化设计方法 |
US11288417B2 (en) * | 2019-04-23 | 2022-03-29 | Autodesk, Inc. | Topology optimization of structure with multiple targets |
CN110414165B (zh) * | 2019-08-01 | 2022-06-10 | 华东交通大学 | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 |
CN110705167B (zh) * | 2019-10-09 | 2022-12-27 | 华东交通大学 | 具有规定的真实输出位移的柔顺机构拓扑优化设计方法 |
CN111475976B (zh) * | 2020-03-30 | 2022-07-26 | 浙江大学 | 考虑混合不确定性的颗粒增强材料构件稳健拓扑优化方法 |
CN111523264B (zh) * | 2020-04-02 | 2023-06-02 | 三峡大学 | 一种具有极限弹性性能的多相材料微结构拓扑优化方法 |
-
2020
- 2020-11-24 CN CN202011329988.8A patent/CN112417692B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844917A (zh) * | 2017-01-10 | 2017-06-13 | 西北工业大学 | 一种基于支反力方差约束的机床底座拓扑优化设计方法 |
CN107341316A (zh) * | 2017-07-13 | 2017-11-10 | 西北工业大学 | 设计相关压力载荷作用下的结构形状‑拓扑联合优化方法 |
CN111310377A (zh) * | 2020-02-21 | 2020-06-19 | 北京航空航天大学 | 一种基频和频率间隔混合约束下的连续体结构非概率可靠性拓扑优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112417692A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112417692B (zh) | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 | |
CN110414165B (zh) | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 | |
CN110795873B (zh) | 一种考虑尺寸控制的跨尺度拓扑优化方法 | |
CN113434921B (zh) | 一种考虑介纳观尺度效应的结构等几何拓扑优化方法 | |
CN110008512B (zh) | 一种考虑承载特性的负泊松比点阵结构拓扑优化方法 | |
CN102720295B (zh) | 一种基于索穹顶张拉和承载全过程分析的预应力确定方法 | |
EP2866019A1 (en) | Air-water-soil skeleton coupled calculation device, coupled calculation method, and coupled calculation program | |
CN112765856B (zh) | 一种功能梯度多孔结构拓扑优化的混合水平集方法 | |
Cao et al. | Form-finding analysis of suspension bridges using an explicit iterative approach | |
CN113191040A (zh) | 一种考虑结构稳定性的单材料结构拓扑优化方法和系统 | |
CN107766670B (zh) | 周期性手征蜂窝结构材料等效弹性模量预测方法 | |
CN109858071B (zh) | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 | |
Ye et al. | The modified force-density method for form-finding of membrane structures | |
CN113158510B (zh) | 一种索网结构的优化找形方法 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Borgart | An approximate calculation method for air inflated cushion structures for design purposes | |
CN115600383B (zh) | 一种不确定性数据驱动计算力学方法、存储介质及产品 | |
CN116226986A (zh) | 一种拉索无应力长度计算方法 | |
CN108052729B (zh) | 超疏水表面微结构的逆向设计方法 | |
CN115544623A (zh) | 一种基于逆吊法的空间拱桥找形方法 | |
CN113139313B (zh) | 基于有限元理论的等效轻量化结构分析方法 | |
CN110110387B (zh) | 一种索结构有限元找形分析方法 | |
CN108345729B (zh) | 一种基于群论的对称索杆结构找形方法 | |
CN117892414B (zh) | 简谐激励下的各向异性材料结构稳健性拓扑优化设计方法 | |
Ye et al. | The modified dynamic relaxation method for the form-finding of membrane structures |
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 |