CN114925565A - 纤维增强橡胶复合材料的纤维角度优化方法 - Google Patents
纤维增强橡胶复合材料的纤维角度优化方法 Download PDFInfo
- Publication number
- CN114925565A CN114925565A CN202210527417.8A CN202210527417A CN114925565A CN 114925565 A CN114925565 A CN 114925565A CN 202210527417 A CN202210527417 A CN 202210527417A CN 114925565 A CN114925565 A CN 114925565A
- Authority
- CN
- China
- Prior art keywords
- fiber
- optimization
- composite material
- rubber composite
- reinforced rubber
- 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
- 239000000835 fiber Substances 0.000 title claims abstract description 155
- 238000005457 optimization Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000002131 composite material Substances 0.000 title claims abstract description 58
- 238000013461 design Methods 0.000 claims abstract description 49
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 239000000463 material Substances 0.000 claims abstract description 30
- 238000006073 displacement reaction Methods 0.000 claims abstract description 17
- 230000035945 sensitivity Effects 0.000 claims abstract description 17
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000009795 derivation Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 1
- 239000013013 elastic material Substances 0.000 description 1
- 239000002657 fibrous material Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
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
- 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
- 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
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Moulding By Coating Moulds (AREA)
- Reinforced Plastic Materials (AREA)
Abstract
本发明涉及一种纤维增强橡胶复合材料的纤维角度优化方法,属于橡胶复合材料设计领域。步骤包括:一、通过将纤维层等效成纤维平面单元来构建有限元模型;二、采用改进的离散多材料优化方法建立纤维增强橡胶复合材料的纤维角度优化模型和有限元模型;三、通过非线性有限元分析得到节点位移;四、根据优化模型和节点位移采用伴随方法计算灵敏度信息;五、根据灵敏度信息采用优化算法进行优化迭代,得到最优结果。本发明可以分别考虑橡胶基体的超弹性与纤维层的几何非线性来解决纤维增强橡胶复合材料纤维角度的优化问题,同时在离散材料优化方法中添加了映射函数提高了纤维收敛的速度。为纤维增强橡胶复合材料设计提供一种有效且高效的优化方法。
Description
技术领域
本发明涉及橡胶复合材料设计领域,特别涉及一种纤维增强橡胶复合材料的纤维角度优化方法。
背景技术
由于纤维增强橡胶复合材料具有低密度和高强重比等优秀的力学性能,在航空航天、土木工程、汽车、软组织工程等工程领域的应用日益广泛,例如航空轮胎、空气弹簧、人工组织等。纤维增强橡胶复合材料通过按特定方向排列纤维,可以在纤维方向上获得更高的强度。因此,纤维增强橡胶复合材料结构的整体力学性能在很大程度上取决于纤维方向的设计。然而,纤维增强橡胶复合材料涉及非线性大变形问题,目前关于复合材料纤维角度优化的研究停留在线性材料层面,在纤维增强橡胶复合材料的优化设计中无法适用,橡胶复合材料设计只能通过试验得到可用的纤维角度,试验过程需要花费大量的时间和资源并且不容易得到较优的纤维角度设计方案。因此,急需一种可以提高橡胶复合材料的设计能力的复合材料的优化方法,同时具有重要的理论与应用价值。
发明内容
本发明的目的在于提供一种纤维增强橡胶复合材料的纤维角度优化方法,解决了现有的纤维角度优化方法无法应用于纤维增强橡胶复合材料的问题。本发明为纤维增强橡胶复合材料设计提供一种有效且高效的优化方法。
本发明的上述目的通过以下技术方案实现:
纤维增强橡胶复合材料的纤维角度优化方法,包括以下步骤:
步骤一、给定设计域和边界条件,通过将纤维层等效成纤维平面单元的纤维增强橡胶复合材料建模方法来构建有限元模型;
步骤二、采用改进的离散多材料优化方法构建纤维增强橡胶复合材料的纤维角度优化模型和有限元模型;
步骤三、通过非线性有限元分析得到节点位移;
步骤四、根据优化模型和节点位移采用伴随方法引入拉格朗日乘子得到伴随方程,进而得到每个设计变量的灵敏度信息;
步骤五、根据灵敏度信息采用优化算法对所构建的纤维增强橡胶复合材料的纤维角度优化问题进行求解并判断收敛性,若满足收敛条件,得到最优结果。
步骤二所述的采用改进的离散多材料优化方法构建纤维增强橡胶复合材料的纤维角度优化模型和有限元模型,具体是:构造离散多材料优化方法的纤维平面单元材料矩阵的插值函数形式来实现将离散的角度优化问题转化为连续设计变量优化问题,即设计组件的每个待选角度的材料矩阵乘以权系数后叠加得到优化形式的材料矩阵;并添加设计变量映射函数提高纤维角度优化的收敛效率;得到纤维增强橡胶复合材料纤维角度优化的优化模型;目标函数为结构的柔度,约束条件为满足非线性有限元的平衡方程和设计变量的取值范围。
步骤四所述的根据优化模型和节点位移采用伴随方法引入拉格朗日乘子得到伴随方程,进而得到每个设计变量的灵敏度信息,具体是:由于橡胶基体和纤维平面单元分别涉及超弹性和几何非线性问题,因此,引入拉格朗日乘子将纤维平面单元刚度凝聚到橡胶基体后的非线性有限元平衡方程与目标函数联立,通过伴随法得到拉格朗日乘子;分别对橡胶基体和纤维平面单元两部分的残余力进行求导联立后得到灵敏度信息。
本发明的有益效果在于:通过将纤维层等效成纤维平面单元来构建橡胶复合材料有限元模型,使得该方法可以分别考虑橡胶基体的超弹性与纤维层的几何非线性来解决纤维增强橡胶复合材料纤维角度的优化问题,同时在离散材料优化方法中添加了映射函数提高了纤维收敛的速度。本发明为橡胶复合材料设计提供一种有效且高效的优化方法,解决了现有的纤维角度优化方法不能优化纤维增强橡胶复合材料的问题。实用性强。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1为本发明的纤维增强橡胶复合材料的纤维角度优化方法的流程图;
图2为本发明的纤维增强橡胶复合材料等效建模的示意图;
图3为本发明实施例中模型几何尺寸与边界条件示意图;
图4为本发明的实施例的设计优化结果图。
具体实施方式
下面将结合附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
参见图1至图4所示,本发明的纤维增强橡胶复合材料的纤维角度优化方法,为纤维增强橡胶复合材料设计提供一种有效且高效的优化方法。包括如下步骤:
步骤一、给定设计域和边界条件,通过将纤维层等效成纤维平面单元的纤维增强橡胶复合材料建模方法来构建有限元模型;
步骤二、构建纤维角度优化模型和有限元模型:构造离散多材料优化方法的纤维平面单元材料矩阵的插值函数形式来实现将离散的角度优化问题转化为连续设计变量优化问题,即设计组件的每个待选角度的材料矩阵乘以权系数后叠加得到优化形式的材料矩阵。并添加设计变量映射函数提高纤维角度优化的收敛效率。得到纤维增强橡胶复合材料纤维角度优化的优化模型:目标函数为结构的柔度,约束条件为满足非线性有限元的平衡方程和设计变量的取值范围。
步骤三、进行非线性有限元分析得到节点位移:通过坐标转换矩阵将纤维平面单元的刚度叠加到橡胶基体上,来实现消除纤维平面单元的自由度。采用迭代求解方法对非线性有限元平衡方程进行线性化处理,得到纤维增强橡胶复合材料的线性化平衡方程。
步骤四、根据优化模型和节点位移采用伴随方法计算灵敏度信息:由于橡胶基体和纤维平面单元分别涉及超弹性和几何非线性问题,直接求解设计变量的灵敏度非常困难。因此,引入拉格朗日乘子将纤维平面单元刚度凝聚到橡胶基体后的非线性有限元平衡方程与目标函数联立,通过伴随法得到拉格朗日乘子。分别对橡胶基体和纤维平面单元两部分的残余力进行求导联立后得到灵敏度信息。
步骤五、根据灵敏度信息进行优化迭代和收敛性判断:根据步骤四得到的灵敏度信息采用优化算法进行纤维角度优化,通过更新设计变量来实现在待选角度中转换及调整,实现结构不断更新优化,在每次结构纤维角度更新后进行收敛性判断。若满足收敛条件,则优化过程结束,输出当前的设计结果作为最优的纤维角度;否则,返回步骤三,继续对结构进行优化。
实施例1:
参见图1所示,一种纤维增强橡胶复合材料的纤维角度优化方法,包括以下步骤:
步骤一、构建有限元模型。给定设计域和边界条件,对设计域中的橡胶基体进行网格划分,在每个橡胶基体单元内创建纤维平面单元,采用纤维平面单元来等效纤维层构建纤维增强橡胶复合材料有限元模型,图2为纤维增强橡胶复合材料等效建模的示意图,图2中的(a)部分为纤维平面单元的几何模型,(b)部分为纤维平面单元的等效模型。
步骤二、构建纤维角度优化模型和有限元模型。采用改进的双值编码参数化方法(简称BCP法)构建散多材料优化方法的纤维平面单元材料矩阵的插值函数形式为:
其中,Di为第i个待确定纤维角度组件的材料矩阵,nc为待选角度的数量,β为设计变量映射函数惩罚因子,p为权系数惩罚因子,xik为优化的设计变量,为了消除初始结构对优化结果的影响,初始设计变量为0,Dj是纤维平面单元为第j个待选角度时的材料矩阵,nv为BCP法中设计变量的数量,定义为:
sjk为BCP法中的参数,具体表达式为:
通过在BCP方法中添加设计变量映射函数——tanh函数,并对惩罚因子β进行适当的控制来提高纤维角度优化方法的收敛效果与收敛效率。
得到纤维增强橡胶复合材料纤维角度优化的优化列式:
步骤三、进行非线性有限元分析得到节点位移。通过坐标转换矩阵将纤维平面单元的刚度叠加到橡胶基体上,实现消除纤维平面单元的自由度。通过Newton-Raphson迭代方法进行线性化处理得到纤维增强橡胶复合材料的平衡方程表达式为:
[k]e{Δu}e=R(U) (5)
其中,{Δu}e为橡胶单元的位移增量,[k]e为橡胶单元叠加纤维平面单元的切线刚度矩阵,由基体项[kM]e和纤维项[kR]e两部分组成。R(U)为橡胶单元叠加纤维平面单元的残余力向量,具体表达式如下:
其中,BN为橡胶单元应变与位移的关系矩阵,BL为纤维平面单元的线性应变与位移的关系矩阵,BNL为纤维平面单元的非线性应变与位移的关系矩阵,BG和G分别表示橡胶和基体确定非线性应变与位移的关系的中间矩阵。DM和DR分别为橡胶单元和纤维平面单元的材料矩阵,ΣM和ΣR分别为橡胶单元和纤维平面单元的第二基尔霍夫应力矩阵,SM和SR分别为橡胶单元和纤维平面单元的第二基尔霍夫应力的列向量形式。CT为纤维平面单元相对于所在橡胶单元的相对位置矩阵,L为纤维平面单元整体坐标到局部坐标的变换矩阵。
步骤四、计算灵敏度信息。通过目标函数对设计变量求导得到:
其中,权系数ωik对设计变量xij的导数为:
拉格朗日乘子λ通过计算伴随方程得到:
步骤五、优化迭代和收敛性判断。根据步骤四得到的灵敏度信息采用MMA优化算法进行纤维角度优化,通过更新设计变量来实现在待选角度中转换及调整,实现结构不断更新优化,在每次结构纤维角度更新后进行收敛性判断,收敛条件为:
其中,wij为第i个待确定纤维角度组件的第j个待选角度的权系数,ε为允许的收敛误差。若满足收敛条件,则优化过程结束,输出当前的设计结果作为最优的纤维角度;否则,返回步骤三,继续对结构进行优化。
本发明纤维增强橡胶复合材料的纤维角度优化方法的有益效果在于:通过将纤维层等效成纤维平面单元来构建橡胶复合材料有限元模型,使得该方法可以分别考虑橡胶基体的超弹性与纤维层的几何非线性来解决纤维增强橡胶复合材料纤维角度的优化问题,同时在离散材料优化方法中添加了映射函数提高了纤维收敛的速度。本发明为橡胶复合材料设计提供一种有效且高效的设计方法。
实施例2:
以纤维增强橡胶复合悬臂板结构为例,参见图3所示,为本实施例的模型几何尺寸、有限元模型与边界条件,单元尺寸为5mm。其中改进的BCP方法中β的取值从0开始,在优化进行到第十代时增大到5,此后每两代增加1。优化的允许收敛误差设置为0.01%。规定每个橡胶基体单元对应的纤维平面单元的纤维角度都进行优化设计。橡胶基体采用Neo-Hookean本构模型,其中A10是拟合Neo-Hookean本构模型的材料参数,K为体积模量,纤维为线弹性材料。橡胶基体和纤维的材料参数如表1所示。本实施例最优的纤维角度优化结果如图4所示。本发明方法经过实施例26步迭代后得到优化设计结果。优化结果的目标末端柔度为4142.0,相比于所有角度选为90°结构的末端柔度4483.9降低了7.6%。
表1橡胶与纤维材料参数
以上所述仅为本发明的优选实例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡对本发明所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种纤维增强橡胶复合材料的纤维角度优化方法,其特征在于:包括以下步骤:
步骤一、给定设计域和边界条件,通过将纤维层等效成纤维平面单元的纤维增强橡胶复合材料建模方法来构建有限元模型;
步骤二、采用改进的离散多材料优化方法构建纤维增强橡胶复合材料的纤维角度优化模型和有限元模型;
步骤三、通过非线性有限元分析得到节点位移;
步骤四、根据优化模型和节点位移采用伴随方法引入拉格朗日乘子得到伴随方程,进而得到每个设计变量的灵敏度信息;
步骤五、根据灵敏度信息采用优化算法对所构建的纤维增强橡胶复合材料的纤维角度优化问题进行求解并判断收敛性,若满足收敛条件,得到最优结果。
2.根据权利要求1所述的纤维增强橡胶复合材料的纤维角度优化方法,其特征在于:步骤二所述的采用改进的离散多材料优化方法构建纤维增强橡胶复合材料的纤维角度优化模型和有限元模型,具体是:构造离散多材料优化方法的纤维平面单元材料矩阵的插值函数形式来实现将离散的角度优化问题转化为连续设计变量优化问题,即设计组件的每个待选角度的材料矩阵乘以权系数后叠加得到优化形式的材料矩阵;并添加设计变量映射函数提高纤维角度优化的收敛效率;得到纤维增强橡胶复合材料纤维角度优化的优化模型;目标函数为结构的柔度,约束条件为满足非线性有限元的平衡方程和设计变量的取值范围。
3.根据权利要求1所述的纤维增强橡胶复合材料的纤维角度优化方法,其特征在于:步骤四所述的根据优化模型和节点位移采用伴随方法引入拉格朗日乘子得到伴随方程,进而得到每个设计变量的灵敏度信息,具体是:由于橡胶基体和纤维平面单元分别涉及超弹性和几何非线性问题,因此,引入拉格朗日乘子将纤维平面单元刚度凝聚到橡胶基体后的非线性有限元平衡方程与目标函数联立,通过伴随法得到拉格朗日乘子;分别对橡胶基体和纤维平面单元两部分的残余力进行求导联立后得到灵敏度信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210527417.8A CN114925565B (zh) | 2022-05-16 | 纤维增强橡胶复合材料的纤维角度优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210527417.8A CN114925565B (zh) | 2022-05-16 | 纤维增强橡胶复合材料的纤维角度优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114925565A true CN114925565A (zh) | 2022-08-19 |
CN114925565B CN114925565B (zh) | 2024-06-28 |
Family
ID=
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
CN110175430A (zh) * | 2019-06-04 | 2019-08-27 | 成都数象科技有限公司 | 一种叶片极限状态下多目标自动优化方法 |
WO2020211012A1 (zh) * | 2019-04-17 | 2020-10-22 | 大连理工大学 | 一种面向混杂纤维复合材料板壳结构的快速协同优化方法 |
CN112364546A (zh) * | 2020-12-03 | 2021-02-12 | 华中科技大学 | 一种基于双边滤波的纤维增强复合材料结构优化设计方法 |
CN113515850A (zh) * | 2021-05-20 | 2021-10-19 | 大连理工大学 | 一种考虑纤维连续的纤维增强复合材料结构布局优化设计方法 |
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590325A (zh) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | 一种基于Shepard插值的纤维增强复合材料结构优化方法 |
WO2020211012A1 (zh) * | 2019-04-17 | 2020-10-22 | 大连理工大学 | 一种面向混杂纤维复合材料板壳结构的快速协同优化方法 |
CN110175430A (zh) * | 2019-06-04 | 2019-08-27 | 成都数象科技有限公司 | 一种叶片极限状态下多目标自动优化方法 |
CN112364546A (zh) * | 2020-12-03 | 2021-02-12 | 华中科技大学 | 一种基于双边滤波的纤维增强复合材料结构优化设计方法 |
CN113515850A (zh) * | 2021-05-20 | 2021-10-19 | 大连理工大学 | 一种考虑纤维连续的纤维增强复合材料结构布局优化设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ho-Huu et al. | Optimization of laminated composite plates for maximizing buckling load using improved differential evolution and smoothed finite element method | |
Do et al. | Material optimization of tri-directional functionally graded plates by using deep neural network and isogeometric multimesh design approach | |
CN112100774A (zh) | 一种基于变密度法的应力和应变能双约束的拓扑优化方法 | |
Gao et al. | Topology optimization of micro-structured materials featured with the specific mechanical properties | |
Haj-Ali et al. | A multi-scale constitutive formulation for the nonlinear viscoelastic analysis of laminated composite materials and structures | |
Pi et al. | Nonlinear inelastic analysis of steel beam-columns. I: Theory | |
Santos et al. | A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures | |
Haj-Ali et al. | Nonlinear micromechanical formulation of the high fidelity generalized method of cells | |
Atluri | On the hybrid stress finite element model for incremental analysis of large deflection problems | |
Kim et al. | Coupled stability analysis of thin-walled composite beams with closed cross-section | |
Kim et al. | A co-rotational 8-node assumed strain shell element for postbuckling analysis of laminated composite plates and shells | |
Verma et al. | Geometrically nonlinear flexural analysis of multilayered composite plate using polynomial and non-polynomial shear deformation theories | |
CN112131770B (zh) | 一种考虑可靠性的功能梯度连续体结构轻量化设计方法 | |
CN114925565B (zh) | 纤维增强橡胶复合材料的纤维角度优化方法 | |
CN117292777A (zh) | 一种基于锯齿效应的复合材料层合板层间应力预测方法、装置、介质及设备 | |
CN114925565A (zh) | 纤维增强橡胶复合材料的纤维角度优化方法 | |
Carrera et al. | Comparisons between 1d (beam) and 2d (plate/shell) finite elements to analyze thin walled structures | |
Komurcu et al. | Towards composite suspension control arm: Conceptual design, structural analysis, laminate optimization, manufacturing, and experimental testing | |
Topal et al. | Frequency optimization of laminated folded composite plates | |
Kim et al. | A 4-node co-rotational ANS shell element for laminated composite structures | |
CN113779825A (zh) | 一种基于单位分解及移动控制点的结构优化设计方法 | |
Kim et al. | Postbuckling analysis of composite panels with imperfection damage | |
Bellahkim et al. | The Effect of Boundary Conditions and Stacking Sequence on the Nonlinear Behavior of Laminated Composite Plates in Bending | |
Ju-Fen et al. | An improved arc-length method and application in the post-buckling analysis for composite structures | |
Kim et al. | Buckling strength prediction of CFRP cylindrical panels using finite element method |
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 |