CN111625910A - 一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 - Google Patents
一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 Download PDFInfo
- Publication number
- CN111625910A CN111625910A CN202010476305.5A CN202010476305A CN111625910A CN 111625910 A CN111625910 A CN 111625910A CN 202010476305 A CN202010476305 A CN 202010476305A CN 111625910 A CN111625910 A CN 111625910A
- Authority
- CN
- China
- Prior art keywords
- stress
- fatigue
- constraint
- formula
- static strength
- 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
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/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于静强度和疲劳约束的柔顺机构拓扑优化方法,包括以下步骤:评定柔顺机构的静强度和疲劳强度;将机构的静强度和疲劳约束转化为不同的应力约束,将单元的各类应力约束凝聚成不同的全局应力约束;以输出位移最大化为目标函数,以静强度、疲劳和体积作为约束,构建基于静强度和疲劳约束的柔顺机构拓扑优化模型;进行灵敏度分析,利用映射过滤修正灵敏度,求解拓扑优化问题,迭代至满足收敛准则,获得最优的柔顺机构拓扑构型。本发明利用带符号的von Mises和修正的Goodman疲劳准则进行基于静强度和疲劳约束的柔顺机构拓扑优化,获得的柔顺机构能够满足静强度和疲劳强度要求,并且有效地避免类铰链结构。
Description
技术领域
本发明涉及柔顺机构设计相关技术领域,特别是涉及一种基于静强度和疲劳约束的柔顺机构拓扑优化方法。
背景技术
柔顺机构是采用柔性元件的变形传递运动、力或能量的一种新型免装配的机构。与传统刚性机构相比,它具有以下优点:1)节约材料,制造成本低;2)无需装配,精度高;3)不需要润滑;4)容易集成现代驱动技术。由于柔顺机具有以上优点,使得它在微机电系统、微/纳制造、精密定位、精密操作、精密和超精密加工领域中具有广泛应用。
目前,拓扑优化方法为柔顺机构的常用设计方法,该方法是将设计域离散为有限单元,在给定的边界条件下,根据优化算法删除部分单元,寻求最优的柔顺机构构型,以达到机构的性能最优。目前,现有拓扑优化方法获得的柔顺机构往往难以满足静强度和疲劳强度要求,因此不能满足实际工程的需要。
发明内容
本发明的目的是提供一种基于静强度和疲劳约束的柔顺机构拓扑优化方法,以解决上述现有技术存在的问题,使柔顺机构满足静强度和疲劳要求,满足实际工程的需要。
为实现上述目的,本发明提供了如下方案:本发明提供一种基于静强度和疲劳约束的柔顺机构拓扑优化方法,包括以下步骤:
S1、初始化参数:设定柔顺机构的设计域、边界条件、材料参数、作用正弦载荷的参数,设置输入、输出的弹簧刚度、设计变量的过滤半径,将初始设计域离散成N个有限元单元,设置移动渐近算法(MMA)的初始变量;
S2、采用固体各向同性材料惩罚模型计算结构的总装刚度矩阵,采用准静态分析法求解最大的载荷列阵作用下的结构节点位移场;
S3、采用应力松弛方法对柔顺机构单元应力进行惩罚,然后获得带符号的单元vonMises应力幅值;采用带符号的von Mises应力和修正的Goodman准则评定机构的疲劳强度,采用带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度,将柔顺机构的静强度和疲劳约束转化为不同的应力约束,采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束,采用自适应约束缩放方法修正P范数应力;
S4、以柔顺机构的输出端位移最大化为目标函数,以机构的静强度、疲劳及体积为约束,构建基于静强度和疲劳约束的柔顺机构拓扑优化模型;
S5、由有限元分析结果,求解柔顺机构的输出位移、静强度和疲劳强度,进行目标函数和约束的灵敏度分析,采用映射过滤技术修正灵敏度;
S6、采用移动渐近算法MMA更新设计变量,判断优化算法是否满足收敛条件,若不满足,转至步骤S2;若满足,则拓扑优化过程结束,输出满足疲劳和静强度约束的柔顺机构拓扑构型。
优选的,采用准静态分析法求解最大的载荷列阵作用下的结构节点位移场
KU=Fmax
式中:K为结构总装刚度矩阵;U为结构的单元节点位移列阵;Fmax为正弦载荷的最大载荷列阵。
优选的,步骤S3中,为了避免应力约束优化问题的奇异解现象,采用应力松弛方法对单元应力进行惩罚,公式为
优选的,采用带符号的von Mises等效应力度量单元的应力均值和应力幅值表示公式为:
优选的,所述修正的Goodman疲劳准则表示为:
式中:L1,e和L2,e为应力约束,σut为拉伸强度;σNf为表示无限循环次数下的许用临界应力幅值。
所述许用临界应力幅值σNf由Basquin方程求得
σNf=σ'f(2Nf)b
式中:σ'f为材料的疲劳强度系数;b为材料的疲劳强度指数。
优选的,所述带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度表示为:
优选的,采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束Li,PN(i=1,2,3,4),计算公式为
式中,P为P范数参数;
由于L3,e和L4,e的值有可能为负值,导致不用采用P范数进行最大近似处理,因此对静强度约束进行改写
采用自适应约束缩放方法修正P范数应力,计算公式为
Li,max=γiLi,PN i=1,2,3,4
γi表示自适应约束缩放系数公式表示为
优选的,基于静强度和疲劳约束的柔顺机构拓扑优化模型为
Fmax=KU
V≤f*V0
0<ρmin≤ρe≤1,e=1,2,3...N
式中,Uout为柔顺机构的输出位移,为作用在柔顺机构输出端的虚拟载荷列阵,为虚拟载荷阵列作用产生的节点位移列阵,Fmax为作用在柔顺机构输入端的载荷列阵,U为输入端载荷F作用产生的节点位移矢量,V是优化后的结构总体积,V0是优化前的结构总体积,f*是允许的材料体积比,N为有限元的个数。
优选的,目标函数和体积约束的灵敏度公式表示为
式中,v0为实体单元的体积。
优选的,任一Li,max应力约束对设计变量的灵敏度表示为
其中,λe为单元的伴随列阵,其满足以下伴随方程
单元应力对设计变量和单元结点位移的灵敏度分别为
单元刚度矩阵对设计变量的灵敏度为
优选的,映射过滤,先采用密度过滤方法对单元密度进行过滤,公式表示为:
式中:rj为单元j与单元e质心之间的距离。
采用密度过滤,拓扑优化结构的边界仍然存在中间密度的灰度单元,进一步采用Heaviside映射修正过滤的单元密度,从而使得单元密度向0和1两端集中
本发明公开了以下技术效果:
(1)采用修正的Goodman准则评定机构的疲劳强度,采用带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度,将静强度和疲劳转化不同的应力约束。
(2)采用基于静强度和疲劳约束的柔顺机构拓扑优化设计方法,拓扑优化获得的柔顺机构构型能够同时满足静强度和疲劳要求,能够有效地抑制类铰链和应力集中现象,使材料的应力分布较为均匀。
(3)采用P范数将所有单元的静强度约束和疲劳约束凝聚化成全局静强度约束和疲劳约束,采用自适应约束缩放法使得各类的全局P应力更加接近与之对应的应力的最大值,克服了静强度和疲劳引起的巨大约束数目导致的约束的局部性,大大地减少计算量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为按照本发明基本流程图;
图2为用于示范柔性夹持器设计域、作用载荷与边界条件的示意图;
图3为基于静强度和疲劳约束的柔顺机构拓扑优化方法获得的柔顺夹持器构型;
图4为柔顺夹持器的应力分布图;
图5为柔顺夹持器的修正的Goodman图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
参考图1,本发明提供一种基于静强度和疲劳约束的柔顺机构拓扑优化设计方法,包括如下步骤:
步骤1、进行参数初始化:设定柔顺机构的设计域、边界条件、材料参数、作用正弦载荷的参数,设置输入、输出的弹簧刚度、设计变量的过滤半径,将初始设计域离散成N个有限元单元,设置移动渐近算法(MMA)的初始变量;
步骤2、采用固体各向同性材料惩罚模型计算材料结构的总装刚度矩阵,采用准静态分析法求解最大的载荷列阵作用下的结构节点位移场
KU=Fmax (1)
式中:K为结构总装刚度矩阵;U为结构的单元节点位移列阵;Fmax为正弦载荷的最大载荷列阵。
步骤3、采用应力松弛方法对单元的应力进行惩罚,然后求得带符号的单元vonMises应力;采用带符号的von Mises应力和修正的Goodman准则评定机构的疲劳强度,采用带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度,将柔顺机构的静强度和疲劳约束转化为不同的应力约束,采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束,采用自适应约束缩放方法使得各类的全局P应力更加接近与之对应的应力的最大值;
为了避免应力约束优化问题的奇异解现象,采用应力松弛方法对单元应力进行惩罚
由于von Mises等效应力的值总为正,如果采用它进行疲劳强度评估,使得负的von Mises等效应力均值符号丢失,与实际的结构响应不符。因此,采用带符号的von Mises等效应力度量单元的应力均值和应力幅值
采用修正的Goodman疲劳准则进行柔顺机构的疲劳强度评估
式中:σut为拉伸强度;σNf为表示无限循环次数下的许用临界应力幅值,当循环次数Nf>107即认为是无限次循环。
许用应力幅值σNf可由Basquin方程求得
σNf=σ'f(2Nf)b (7)
式中,σ'f为材料的疲劳强度系数,b为材料的疲劳强度指数。
采用带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度
传统的p范数方法是以单元的体积作为比例因子,具有内置缩放的作用,有利于应力约束拓扑优化问题收敛,但是采用国际单位制时,由于单元体积的数值会特别小,采用传统的p范数方法导致比例因子内置缩放作用消失。改进的p范数方法以单元密度作为比例因子,不管采用何种单位制,变化范围为0到1,都能够具有很好的内置缩放作用。采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束Li,PN(i=1,2,3,4)
式中,P为P范数参数。
由于L3,e和L4,e的值有可能为负值,从而导致不用采用P范数进行最大近似处理,因此对静强度约束进行改写:
采用自适应约束缩放法使得各类的全局P应力更加接近与之对应的应力的最大值
Li,max=γiLi,PN i=1,2,3,4 (13)
γi表示自适应约束缩放系数
步骤4、以柔顺机构的输出端位移最大化为目标函数,以机构的静强度、疲劳及体积为约束,构建基于静强度和疲劳约束的柔顺机构拓扑优化模型为
式中,Uout为柔顺机构的输出位移,为作用在柔顺机构输出端的虚拟载荷列阵,为虚拟载荷矢量作用产生的节点位移列阵,Fmax为作用在柔顺机构输入端的载荷列阵,U为输入端载荷F作用产生的节点位移矢量,V是优化后的结构总体积,V0是优化前的结构总体积,f*是允许的材料体积比,N为有限元的个数。
步骤5、由有限元分析结构,求解柔顺机构的输出端位移、静强度约束、疲劳约束和体积约束,进行目标函数和约束的灵敏度分析,采用映射过滤技术修正灵敏度;
目标函数和体积约束的灵敏度为
式中,v0为实体单元的体积。
任一Li,max应力约束对设计变量单元密度的灵敏度为
由有限元平衡方程(1)两边求导得
将式(19)带入式(18),可得
引入伴随向量:
将式(21)带入式(20),可得
其中,λe为单元的伴随列阵。
由公式(10)和公式(15)可得
由公式(10)、公式(15)则有
由公式(10)可得
根据公式(15)可得
式中:a和b为任意两个实数;ε为无穷小量。
则由公式(11)和公式(12)可得
由公式(4)可得von Mises等效应力的幅值和均值对单元应力幅值和均值的灵敏度
由公式(2)得到单元应力对设计变量和单元结点位移的灵敏度分别为
单元刚度矩阵对设计变量的灵敏度为
密度过滤方法表示为
式中:rj为单元j与单元e质心之间的距离。
采用密度过滤,拓扑优化结构的边界仍然存在中间密度的灰度单元,进一步采用Heaviside映射修正过滤的单元密度,从而使得单元密度向0和1两端集中
步骤6、采用移动渐近算法(MMA)更新设计变量,判断优化算法是否满足收敛条件,若不满足,转至步骤3;若满足,则拓扑优化过程结束,获得满足疲劳和静强度约束的柔顺机构拓扑构型。
设置两次连续迭代设计变量的绝对值小于0.0001作为收敛条件
||ρn-ρn-1||≤0.0001 (46)
式中,ρn为第n次迭代的单元密度向量,ρn-1为第n-1次迭代的单元密度向量。
如图2-图5所示,以二维柔顺夹持器为例进一步解释本发明。夹持器的设计域、边界条件及输入输出端如图2所示,设计域大小为80mm×80mm,作用载荷为最大值Fmax=120N和最小值Fmin=-48N的正弦载荷,施加在最左端中间处,输入弹簧刚度Kin和输出弹簧刚度Kout分别为1×103N/mm和1.5×103N/mm,材料选用普通1020碳素钢,弹性模量为E0=203GPa,泊松比为0.3,屈服强度σy为165MPa,拉伸强度σut为300MPa,疲劳强度系数σ'f为300MPa,疲劳强度指数b为-0.075,P范数参数P取值为8;由于设计域对称,以一半作为设计域并且离散成3200个单元,设计变量的过滤半径rmin为2.0。
参考图3—图5,由柔性夹持器的应力分布图可知,该拓扑结构很好地抑制了类铰链结构和应力集中现象,应力分布较为均匀,最大应力较小。由修正的Goodman图可知,获得机构结构的带符号的von Mises应力幅值和均值均落在静强度和疲劳约束所围成的区域中,使得柔顺夹持器同时满足了静强度和疲劳要求,能够满足实际工作需求。
本发明提供的一种基于静强度和疲劳约束的柔顺机构拓扑优化设计方法中,考虑恒定振幅的比例载荷作用下,采用修正的Goodman疲劳准则评定柔顺机构的疲劳强度,以带符号的von Mises应力幅值和均值之和的最大值表征柔顺机构的静强度;将结构静强度和疲劳约束转化为不同的应力约束,采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束,以柔顺机构的输出位移最大化为目标函数,以静强度、疲劳和体积作为约束,构建基于静强度和疲劳约束的柔顺机构拓扑优化模型。利用映射过滤方法修正灵敏度以避免数值不稳定性现象,采用移动近似线算法求解拓扑优化问题,更新迭代知道满足收敛条件,获得最优的机构构型。通过柔顺夹持器算例结果表明,一种基于静强度和疲劳约束的柔顺机构拓扑优化方法获得的柔顺机构构型能够同时满足疲劳和静强度要求,并且有效地抑制类铰链和应力集中现象。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (10)
1.一种基于静强度和疲劳约束的柔顺机构拓扑优化方法,其特征在于:包括以下步骤:
S1、初始化参数:设定柔顺机构的设计域、边界条件、材料参数和正弦载荷参数,设置输入、输出的弹簧刚度、设计变量的过滤半径,将初始设计域离散成N个有限元单元,设置移动渐近算法MMA的初始变量;
S2、采用固体各向同性材料惩罚模型计算结构的总装刚度矩阵,采用准静态分析法求解最大的载荷列阵作用下的结构节点位移场;
S3、采用应力松弛方法对柔顺机构单元应力进行惩罚,然后获得带符号的单元vonMises应力幅值;采用带符号的von Mises应力和修正的Goodman准则评定机构的疲劳强度,采用带符号的单元von Mises应力幅值和均值之和的最大值表征机构的静强度,将柔顺机构的静强度和疲劳约束转化为不同的应力约束,采用改进的P范数将所有单元的各类应力约束分别凝聚成不同的全局应力约束,采用自适应约束缩放方法修正P范数应力;
S4、以柔顺机构的输出端位移最大化为目标函数,以机构的静强度、疲劳及体积为约束,构建基于静强度和疲劳约束的柔顺机构拓扑优化模型;
S5、由有限元分析结果,求解柔顺机构的输出位移、静强度和疲劳强度,进行目标函数和约束的灵敏度分析,采用映射过滤技术修正灵敏度;
S6、采用移动渐近算法MMA更新设计变量,判断优化算法是否满足收敛条件,若不满足,转至步骤S2;若满足,则拓扑优化过程结束,输出满足疲劳和静强度约束的柔顺机构拓扑构型。
2.根据权利要求1所述的基于静强度和疲劳约束的柔顺机构拓扑优化方法,其特征在于:采用准静态分析法求解最大的载荷列阵作用下的结构节点位移场,公式为:
KU=Fmax
式中:K为结构总装刚度矩阵;U为结构的单元节点位移列阵;Fmax为正弦载荷的最大载荷列阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010476305.5A CN111625910B (zh) | 2020-05-29 | 2020-05-29 | 一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010476305.5A CN111625910B (zh) | 2020-05-29 | 2020-05-29 | 一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111625910A true CN111625910A (zh) | 2020-09-04 |
CN111625910B CN111625910B (zh) | 2022-06-07 |
Family
ID=72271140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010476305.5A Active CN111625910B (zh) | 2020-05-29 | 2020-05-29 | 一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111625910B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114722655A (zh) * | 2022-03-14 | 2022-07-08 | 哈尔滨理工大学 | 一种基于局部有限寿命疲劳约束条件的结构拓扑优化方法 |
CN115952628A (zh) * | 2023-03-09 | 2023-04-11 | 中国航发四川燃气涡轮研究院 | 一种计算航空发动机涡轮盘盘心等效应力的方法 |
CN116644535A (zh) * | 2023-07-27 | 2023-08-25 | 华东交通大学 | 一种基于基频最大化的柔性铰链拓扑优化方法 |
CN116644620A (zh) * | 2023-07-27 | 2023-08-25 | 华东交通大学 | 一种基于等几何分析的柔顺机构应力约束拓扑优化方法 |
CN117892414A (zh) * | 2024-03-18 | 2024-04-16 | 华东交通大学 | 简谐激励下的各向异性材料结构稳健性拓扑优化设计方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629761A (zh) * | 2012-03-14 | 2012-08-08 | 广东电网公司珠海供电局 | 一种avc系统电网模型维护方法及其系统 |
WO2013091085A1 (en) * | 2011-12-23 | 2013-06-27 | The Royal Institution For The Advancement Of Learning/Mcgill University | Bone replacement implants with mechanically biocompatible cellular material |
CN105653741A (zh) * | 2014-11-09 | 2016-06-08 | 许亚夫 | 一种高速系统结构载荷优化方法 |
CN106446387A (zh) * | 2016-09-18 | 2017-02-22 | 同济大学 | 一种高速磁悬浮列车制动磁铁连接装置的疲劳寿命预测方法 |
CN107563097A (zh) * | 2017-09-26 | 2018-01-09 | 华中科技大学 | 一种基于应力的启发式结构离散拓扑优化设计方法 |
CN108009323A (zh) * | 2017-11-22 | 2018-05-08 | 北京航空航天大学 | 一种基于多工况失效载荷的复合材料结构不确定优化方法 |
CN109508495A (zh) * | 2018-11-12 | 2019-03-22 | 华东交通大学 | 一种基于k-s函数的柔顺机构全局应力约束拓扑优化方法 |
CN109583091A (zh) * | 2018-11-30 | 2019-04-05 | 长沙理工大学 | 基于自适应约束的柔性机构拓扑优化设计方法 |
CN109614640A (zh) * | 2018-10-30 | 2019-04-12 | 中国电力科学研究院有限公司 | 一种大型风电机组轮毂疲劳寿命预测方法及系统 |
CN110414165A (zh) * | 2019-08-01 | 2019-11-05 | 华东交通大学 | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 |
CN110705167A (zh) * | 2019-10-09 | 2020-01-17 | 华东交通大学 | 具有规定的真实输出位移的柔顺机构拓扑优化设计方法 |
-
2020
- 2020-05-29 CN CN202010476305.5A patent/CN111625910B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013091085A1 (en) * | 2011-12-23 | 2013-06-27 | The Royal Institution For The Advancement Of Learning/Mcgill University | Bone replacement implants with mechanically biocompatible cellular material |
CN102629761A (zh) * | 2012-03-14 | 2012-08-08 | 广东电网公司珠海供电局 | 一种avc系统电网模型维护方法及其系统 |
CN105653741A (zh) * | 2014-11-09 | 2016-06-08 | 许亚夫 | 一种高速系统结构载荷优化方法 |
CN106446387A (zh) * | 2016-09-18 | 2017-02-22 | 同济大学 | 一种高速磁悬浮列车制动磁铁连接装置的疲劳寿命预测方法 |
CN107563097A (zh) * | 2017-09-26 | 2018-01-09 | 华中科技大学 | 一种基于应力的启发式结构离散拓扑优化设计方法 |
CN108009323A (zh) * | 2017-11-22 | 2018-05-08 | 北京航空航天大学 | 一种基于多工况失效载荷的复合材料结构不确定优化方法 |
CN109614640A (zh) * | 2018-10-30 | 2019-04-12 | 中国电力科学研究院有限公司 | 一种大型风电机组轮毂疲劳寿命预测方法及系统 |
CN109508495A (zh) * | 2018-11-12 | 2019-03-22 | 华东交通大学 | 一种基于k-s函数的柔顺机构全局应力约束拓扑优化方法 |
CN109583091A (zh) * | 2018-11-30 | 2019-04-05 | 长沙理工大学 | 基于自适应约束的柔性机构拓扑优化设计方法 |
CN110414165A (zh) * | 2019-08-01 | 2019-11-05 | 华东交通大学 | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 |
CN110705167A (zh) * | 2019-10-09 | 2020-01-17 | 华东交通大学 | 具有规定的真实输出位移的柔顺机构拓扑优化设计方法 |
Non-Patent Citations (7)
Title |
---|
侯杰 等: "考虑机械连接载荷和疲劳性能的装配结构拓扑优化设计方法", 《科学通报》, vol. 64, no. 1, 31 January 2019 (2019-01-31), pages 79 - 86 * |
刘沐宇 等: "考虑焊接残余应力的构造细节疲劳性能评估", 《武汉理工大学学报(交通科学与工程版)》, vol. 41, no. 5, 31 October 2017 (2017-10-31), pages 713 - 718 * |
占金青 等: "热弹性结构全局应力约束拓扑优化设计", 《机械科学与技术》, vol. 38, no. 9, 30 September 2019 (2019-09-30), pages 1386 - 1392 * |
吉国明 等: "承受轴压载荷的加筋板的准静态分析", 《机械强度》, vol. 35, no. 3, 30 June 2013 (2013-06-30), pages 308 - 311 * |
梅帅 等: "开放式结构拓扑优化软件设计与研发", 《机械工程学报》, vol. 51, no. 9, 31 May 2015 (2015-05-31), pages 144 - 152 * |
王选 等: "基于改进的双向渐进结构优化法的应力约束拓扑优化", 《力学学报》, vol. 50, no. 2, 31 March 2018 (2018-03-31), pages 385 - 394 * |
蔡守宇 等: "基于光滑变形隐式曲线的模型重构、应力分析与优化设计一体化方法", 《计算机辅助设计与图形学学报》, no. 9, 30 September 2018 (2018-09-30), pages 1765 - 1772 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114722655A (zh) * | 2022-03-14 | 2022-07-08 | 哈尔滨理工大学 | 一种基于局部有限寿命疲劳约束条件的结构拓扑优化方法 |
CN115952628A (zh) * | 2023-03-09 | 2023-04-11 | 中国航发四川燃气涡轮研究院 | 一种计算航空发动机涡轮盘盘心等效应力的方法 |
CN116644535A (zh) * | 2023-07-27 | 2023-08-25 | 华东交通大学 | 一种基于基频最大化的柔性铰链拓扑优化方法 |
CN116644620A (zh) * | 2023-07-27 | 2023-08-25 | 华东交通大学 | 一种基于等几何分析的柔顺机构应力约束拓扑优化方法 |
CN116644535B (zh) * | 2023-07-27 | 2023-09-29 | 华东交通大学 | 一种基于基频最大化的柔性铰链拓扑优化方法 |
CN116644620B (zh) * | 2023-07-27 | 2023-10-17 | 华东交通大学 | 一种基于等几何分析的柔顺机构应力约束拓扑优化方法 |
CN117892414A (zh) * | 2024-03-18 | 2024-04-16 | 华东交通大学 | 简谐激励下的各向异性材料结构稳健性拓扑优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111625910B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111625910B (zh) | 一种基于静强度和疲劳约束的柔顺机构拓扑优化方法 | |
CN110414165B (zh) | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 | |
CN109508495B (zh) | 一种基于k-s函数的柔顺机构全局应力约束拓扑优化方法 | |
CN102720295B (zh) | 一种基于索穹顶张拉和承载全过程分析的预应力确定方法 | |
CN111695188B (zh) | 一种地锚式悬索桥动力特性的快速精细分析方法 | |
CN110008512A (zh) | 一种考虑承载特性的负泊松比点阵结构拓扑优化方法 | |
CN109614580A (zh) | 基于在线Xgboost算法的抗震混合试验模型更新方法 | |
CN107766670B (zh) | 周期性手征蜂窝结构材料等效弹性模量预测方法 | |
CN112417692B (zh) | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 | |
Toi et al. | Finite element analysis of two-dimensional electrochemical–mechanical response of ionic conducting polymer–metal composite beams | |
Chan | Buckling analysis of structures composed of tapered members | |
CN113076677B (zh) | 基于五次埃尔米特形函数的柔性体结构高阶非线性有限元数值模拟方法 | |
CN114297779A (zh) | 一种飞行器气动力系数快速预测方法、系统及设备 | |
CN113505405B (zh) | 等效荷载获取方法、基于等效荷载的拓扑优化方法及系统 | |
CN115659843A (zh) | 一种利用人工智能技术构建显示动力学模型的方法 | |
CN112818583B (zh) | 等效静荷载获取方法、拓扑优化方法及系统 | |
CN115906333A (zh) | 一种桁架结构的几何非线性等效板动力学建模与响应分析方法 | |
CN112989275B (zh) | 一种用于网络大规模控制系统的多方向方法 | |
CN114756934A (zh) | 一种三维多尺度超材料结构优化设计方法 | |
CN111159946B (zh) | 一种基于最小势能原理的非连续性问题分区求解方法 | |
Rahmad et al. | Single-and multi-objective optimization of a low-speed airfoil using genetic algorithm | |
Chiandussi | On the solution of a minimum compliance topology optimisation problem by optimality criteria without a priori volume constraint specification | |
Wang et al. | Optimization and sensitivity analysis of space frames allowing for large deflection | |
CN116911124A (zh) | 一种拉压结构拓扑优化方法 | |
CN113051678B (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 |