CN113191048A - 一种等距曲线纤维增强结构的优化设计方法 - Google Patents
一种等距曲线纤维增强结构的优化设计方法 Download PDFInfo
- Publication number
- CN113191048A CN113191048A CN202110459680.3A CN202110459680A CN113191048A CN 113191048 A CN113191048 A CN 113191048A CN 202110459680 A CN202110459680 A CN 202110459680A CN 113191048 A CN113191048 A CN 113191048A
- Authority
- CN
- China
- Prior art keywords
- function
- fiber reinforced
- reinforced structure
- fiber
- 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.)
- Pending
Links
- 239000000835 fiber Substances 0.000 title claims abstract description 111
- 238000013461 design Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005457 optimization Methods 0.000 claims abstract description 32
- 230000035945 sensitivity Effects 0.000 claims abstract description 24
- 238000006073 displacement reaction Methods 0.000 claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000009833 condensation Methods 0.000 claims description 3
- 230000005494 condensation Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 11
- 239000002131 composite material Substances 0.000 abstract description 3
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 abstract description 2
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 3
- 239000003733 fiber-reinforced composite Substances 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000009727 automated fiber placement Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 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]
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)
Abstract
本发明公开了一种等距曲线纤维增强结构的优化设计方法,属于结构优化设计领域,所述设计方法包括:确定纤维增强结构的设计域,将权重系数作为设计变量,利用水平集函数的等值线来描述纤维铺放路径;进行结构有限元分析得到整体位移向量;利用整体位移向量计算目标函数及其关于设计变量的灵敏度,利用水平集函数梯度计算约束函数及其关于设计变量的灵敏度;基于灵敏度的优化算法更新设计变量获得新的纤维铺放路径,反复迭代,直至满足优化终止条件,获得最优的等距曲线纤维增强结构。本发明极大改善了复合结构的可制造性,避免制造过程中纤维丝束之间出现间隙和发生重叠,同时也有较高的优化效率。
Description
技术领域
本发明属于结构优化设计技术领域,更具体地,涉及一种等距曲线纤维增强结构的优化设计方法。
背景技术
相比于一般金属材料,纤维增强复合材料具有更大的比模量、更高的比强度和更佳的抗疲劳性,因此其被广泛应用于航空航天、车辆工程、土木建筑等领域。随着先进制造技术的发展,特别是自动纤维铺放技术的出现,这类结构中的纤维可沿曲线铺放,由此获得了变刚度的特性。曲线纤维增强结构的设计自由度更大,工程技术人员通过优化纤维的铺设路径,可以充分发挥材料潜能来满足结构力学性能的要求。需要注意的是,为了避免这类结构的纤维在自动铺放的过程中出现重叠和间隙的制造缺陷,优化得到的纤维铺设路径应尽量平行且等距。
现有技术中,专利CN201710758619.2公开了一种基于Shepard插值的纤维增强复合材料结构优化设计方法,该方法利用Shepard插值构建一个全局连续的函数来计算纤维角度,保证了整个结构设计域内纤维角度的空间连续性。但该技术不易从纤维角度布局提取出连续纤维的铺设路径,更无法获得等距的曲线纤维铺放路径。专利CN202010166397.7公开了一种考虑可制造性的曲线纤维增强结构优化设计方法,该方法利用滤波算法对设计域内的纤维角度进行过滤处理,从而得到满足制造约束的纤维角度布局。同样地,该技术虽改善了这类结构的可制造性,但不能严格保证曲线纤维铺放路径等距。专利CN201811179771.6公开了一种保证制造质量的变刚度复合材料结构优化设计方法,该方法利用数学公式来表达两类制造约束(即曲线纤维铺放路径等距以及曲率合适),并在优化模型中引入约束条件,保证了优化后的纤维角度布局满足制造工艺的要求。但该技术中,大量的约束条件会降低计算效率,使设计过程的时间成本大大增加。以上技术的共同点在于优化结果均为纤维角度布局,需进行后处理来提取纤维铺放路径。
为方便获取曲线纤维铺放路径,专利CN201911206379.0公开了一种基于向量场的复合材料结构优化设计方法,该方法通过径向基函数来构建一个无散场,再利用无散场来描述纤维角度布局,而优化后的无散场的场线即可作为纤维铺放路径。该技术的优化效率较高,但仍无法获得等距的曲线纤维铺放路径。专利CN202010383646.8公开了一种纤维曲线铺放变刚度结构优化设计的水平集方法,该方法利用参数化水平集函数的切线方向来描述整体纤维角度布局,基本可避免最优结构的纤维出现重叠或存在间隙的情况。但该技术不能直接通过水平集函数的等值线来获得等距的纤维铺设路径,仍存在后处理较为麻烦的问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种等距曲线纤维增强结构的优化设计方法,其目的在于,以水平集函数各处的梯度值为约束条件,获取等距的曲线纤维铺设路径,从而提高纤维增强结构的制作效率。
为实现上述目的,按照本发明的一个方面,提供了一种等距曲线纤维增强结构的优化设计方法,包括:
S1:将纤维增强结构的设计域中各个单元对应的径向基函数对应的权重系数作为设计变量;利用水平集函数的等值线描述所述纤维增强结构的纤维铺放路径;
S2:对所述纤维增强结构进行有限元分析得到整体位移向量;
S3:利用所述整体位移向量计算目标函数及其对应的目标灵敏度;利用各单元中心点处的水平集函数梯度计算约束函数及其对应的约束灵敏度;
S4:基于所述目标函数、所述目标灵敏度、所述约束函数和所述约束灵敏更新所述设计变量;利用所述更新后的设计变量更新所述纤维铺放路径;
S5:当更新后的纤维铺放路径对应的结构柔度误差不满足优化终止条件时,针对更新后的纤维铺放路径重复S2-S4;当更新后的纤维铺放路径对应的结构柔度误差满足所述优化终止条件时,将对应的纤维增强结构作为所述等距曲线纤维增强结构。
在其中一个实施例中,所述S1包括:
将各个所述单元的节点作为径向基函数的中心点Pi,并将各个Pi处的权重系数αi作为所述设计变量,其中,i=1,2,...,N,N为所述单元上节点的数量;
结合所述设计变量αi和径向基函数得到各个所述单元中心点处的水平集函数ψ(xe);利用所述水平集函数ψ(xe)的等值线描述所述纤维铺放路径。
在其中一个实施例中,所述ψ(xe)的表达式为:
φ(r)=(max(0,1-r))4(4r+1)
其中,pi=(xi,yi)为所述径向基函数的中心点坐标;xe=(xe,ye)为所述单元的中心点坐标,dmin为自定义参数。
在其中一个实施例中,所述S2包括:
利用公式θe=actan(vey/vex)计算得到各所述单元中心点处的纤维角度θe,其中,vex和vey分别为所述单元中心点处的向量ve的两个分量,e=1,2,...,n,n为单元的数量;
利用所述纤维角度θe建立各所述单元的刚度矩阵Ke,将各个所述单元对应的刚度矩阵Ke组装为整体刚度矩阵K;
利用公式Ku=f求解出所述整体位移向量u,f为外力向量,从而完成所述纤维增强结构的有限元分析。
在其中一个实施例中,所述S3包括:
在其中一个实施例中,所述S3中:
在其中一个实施例中,所述S4包括:
将所述目标函数、所述约束函数、所述目标灵敏度、所述约束灵敏度结合移动渐进线方法更新所述设计变量αi得到更新后的设计变量αi′;
利用所述更新后的设计变量αi′求得所述单元中心点处的更新向量ve和更新水平集函数ψ(xe)′,进而获得更新后的纤维铺放路径;
其中,各所述单元中心点处的更新向量ve′的计算式为:
在其中一个实施例中,所述优化终止条件设置为:
其中,cxrr为柔度误差,k为当前迭代次数,ε为柔度误差的下限值,kmax为迭代次数的上限值。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
本发明提供一种等距曲线纤维增强结构的优化设计方法,利用水平集函数的等值线描述纤维铺放路径;以结构刚度最大化作为设计目标,引入与水平集函数梯度相关的约束条件,使用移动渐进线方法更新设计变量,直至满足优化终止条件,得到最优的曲线纤维增强结构。本发明能够保证得到的纤维增强结构具有等距的曲线纤维铺放路径,极大改善了这类结构的可制造性,避免制造过程中纤维丝束之间出现间隙和发生重叠,同时也有较高的优化效率。
附图说明
图1是本发明一实施例提供的一个简支梁纤维增强结构优化设计的结构示意图;
图2是本发明一实施例提供的等距曲线纤维增强结构的优化设计方法流程图;
图3是图2中关于纤维铺放路径的初始设计图;
图4是图2中采用本发明提供的方法时关于纤维铺放路径的优化结果;
图5是图2中未考虑水平集函数梯度约束时关于纤维铺放路径的优化结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1,本实施方式以带有集中载荷的平面简支梁结构柔度最小化的优化问题为例来解释本发明。在给定的1m×3m矩形设计域D内设置水平的初始纤维铺放路径,固定区域左下角,限制区域右下角竖直方向自由度,并在区域上边界中心处施加集中力F=1N。对纤维增强简支梁结构进行纤维铺放路径优化,使其刚度最大化。
请参阅图2的流程图,本发明提供一种等距曲线纤维增强结构的优化设计方法,包括:
S1:将纤维增强结构的设计域中各个单元的权重系数作为设计变量;利用水平集函数的等值线描述纤维增强结构的纤维铺放路径;
具体的,将给定的设计域D划分为一定数量的正方形单元,单元边长为0.05m,单元数n为20×60,单元节点数N为21×61,在每个单元中心点处定义向量ve(e=1,2,...,n),以单元节点作为径向基函数的中心点Pi(i=1,2,...,N),这里选用紧支径向基函数,其计算式为:
φ(r)=(max(0,1-r))4(4r+1),
其中参数r可根据径向基函数中心点坐标pi=(xi,yi)和单元中心点坐标xe=(xe,ye)计算得到,具体计算式为:
这里取值为1,并在每个Pi处定义权重系数αi,作为优化中的设计变量,根据单元中心点处的水平集函数计算式:
可推导出单元中心处的水平集函数梯度计算式为:
进而向量ve的计算式为:
其中两个偏导数的计算式分别为:
其中r关于x和y的偏导数分别为
S2:对纤维增强结构进行有限元分析得到整体位移向量;
具体的,通过向量ve的两个分量vex和vey计算得到各单元中心点处的纤维角度θe,其计算式为θe=actan(vey/vex),再利用纤维角度θe建立各单元的刚度矩阵Ke,其计算式为:其中Ωe为单元面积,B为位移应变矩阵,D(θe)为单元弹性矩阵,通过组装Ke得到整体刚度矩阵K,然后利用公式Ku=f求解出所述整体位移向量u,f为外力向量,从而完成所述纤维增强结构的有限元分析。S3:利用整体位移向量计算目标函数及其对应的目标灵敏度;利用各单元中心点处的水平集函数梯度计算约束函数及其对应的约束灵敏度。
其中Ke关于θe的偏导数计算式为:
以及θe关于αi的偏导数计算式为:
其中de关于αi的偏导数计算式为其中de关于αi的偏导数计算式为:
S4:基于目标函数、目标灵敏度、约束函数和约束灵敏更新设计变量;利用更新后的设计变量更新纤维铺放路径;
具体的,根据目标函数、约束函数以及它们关于设计变量的灵敏度,结合移动渐进线方法(Method of Moving Asymptotes,通常简称为MMA)更新设计变量αi,MMA为现有的优化算法,在此不赘述,利用更新后的αi求得单元中心点处新的向量ve以及水平集函数ψ(xe),并绘制水平集函数的等值线,即得到新的纤维铺放路径。
S5:当更新后的纤维铺放路径不满足优化终止条件时,针对更新后的纤维铺放路径重复S2-S4;当更新后的纤维铺放路径满足优化终止条件时,将对应的纤维增强结构作为等距曲线纤维增强结构。
具体的,重复S2-S4,每一次重复称为一次迭代,直至满足优化终止条件,得到等距曲线纤维增强结构的最优设计,优化终止条件设定为
其中cerr为柔度误差,k为当前迭代次数,ε为柔度误差的下限值,这里取0.1%,kmax为迭代次数的上限值,这里取1000。
本发明较佳实施例的优化结果如下:优化后的纤维铺放路径如图4所示,柔度值为95.58。作为对比,未考虑水平集函数梯度约束时优化得到的纤维铺放路径如图5所示,柔度值为80.74,此时得到的结构虽刚度较大(柔度较小),但纤维铺放路径之间不等距。因此,从避免制造缺陷的角度,本发明提供的基于水平集的等距曲线纤维增强结构优化设计方法更优,可保证最优结构的纤维丝束在制造过程中不会产生重叠和间隙。
本领域的技术人员容易理解,以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种等距曲线纤维增强结构的优化设计方法,其特征在于,包括:
S1:将纤维增强结构的设计域中各个单元对应的径向基函数对应的权重系数作为设计变量;利用水平集函数的等值线描述所述纤维增强结构的纤维铺放路径;
S2:对所述纤维增强结构进行有限元分析得到整体位移向量;
S3:利用所述整体位移向量计算目标函数及其对应的目标灵敏度;利用各单元中心点处的水平集函数梯度计算约束函数及其对应的约束灵敏度;
S4:基于所述目标函数、所述目标灵敏度、所述约束函数和所述约束灵敏更新所述设计变量;利用所述更新后的设计变量更新所述纤维铺放路径;
S5:当更新后的纤维铺放路径对应的结构柔度误差不满足优化终止条件时,针对更新后的纤维铺放路径重复S2-S4;当更新后的纤维铺放路径对应的结构柔度误差满足所述优化终止条件时,将对应的纤维增强结构作为所述等距曲线纤维增强结构。
2.如权利要求1所述的等距曲线纤维增强结构的优化设计方法,其特征在于,所述S1包括:
将各个所述单元的节点作为径向基函数的中心点Pi,并将各个Pi处的权重系数αi作为所述设计变量,其中,i=1,2,...,N,N为所述单元上节点的数量;
结合所述设计变量αi和径向基函数得到各个所述单元中心点处的水平集函数ψ(xe);利用所述水平集函数ψ(xe)的等值线描述所述纤维铺放路径。
4.如权利要求1所述的等距曲线纤维增强结构的优化设计方法,其特征在于,所述S2包括:
利用公式θe=actan(vey/vex)计算得到各所述单元中心点处的纤维角度θe,其中,vex和vey分别为所述单元中心点处的向量ve的两个分量,e=1,2,...,n,n为单元的数量;
利用所述纤维角度θe建立各所述单元的刚度矩阵Ke,将各个所述单元对应的刚度矩阵Ke组装为整体刚度矩阵K;
利用公式Ku=f求解出所述整体位移向量u,f为外力向量,从而完成所述纤维增强结构的有限元分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459680.3A CN113191048A (zh) | 2021-04-27 | 2021-04-27 | 一种等距曲线纤维增强结构的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459680.3A CN113191048A (zh) | 2021-04-27 | 2021-04-27 | 一种等距曲线纤维增强结构的优化设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113191048A true CN113191048A (zh) | 2021-07-30 |
Family
ID=76979490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110459680.3A Pending CN113191048A (zh) | 2021-04-27 | 2021-04-27 | 一种等距曲线纤维增强结构的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113191048A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779825A (zh) * | 2021-08-24 | 2021-12-10 | 华中科技大学 | 一种基于单位分解及移动控制点的结构优化设计方法 |
CN114692455A (zh) * | 2022-03-28 | 2022-07-01 | 江苏理工学院 | 一种基于分割等距重组曲线的复合材料层合板的建模方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111723457A (zh) * | 2020-05-08 | 2020-09-29 | 华中科技大学 | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 |
-
2021
- 2021-04-27 CN CN202110459680.3A patent/CN113191048A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111723457A (zh) * | 2020-05-08 | 2020-09-29 | 华中科技大学 | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 |
Non-Patent Citations (2)
Title |
---|
HONDA, SHINYA; NARITA, YOSHIHIRO: "Vibration design of laminated fibrous composite plates with local anisotropy induced by short fibers and curvilinearfibers", 《COMPOSITE STRUCTURES》, vol. 93, no. 2, 31 December 2011 (2011-12-31), XP027477730, DOI: 10.1016/j.compstruct.2010.07.003 * |
王琥;李启迪;李光耀;: "变刚度复合材料结构设计方法及不确定性分析研究进展", 《机械工程学报》, no. 08, 1 April 2019 (2019-04-01) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779825A (zh) * | 2021-08-24 | 2021-12-10 | 华中科技大学 | 一种基于单位分解及移动控制点的结构优化设计方法 |
CN113779825B (zh) * | 2021-08-24 | 2024-04-19 | 华中科技大学 | 一种基于单位分解及移动控制点的结构优化设计方法 |
CN114692455A (zh) * | 2022-03-28 | 2022-07-01 | 江苏理工学院 | 一种基于分割等距重组曲线的复合材料层合板的建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111723457B (zh) | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 | |
CN107220461B (zh) | 一种变刚度复合材料板壳结构高效优化方法 | |
CN113191048A (zh) | 一种等距曲线纤维增强结构的优化设计方法 | |
CN108763658B (zh) | 基于等几何方法的组合薄壁结构固有频率设计方法 | |
CN110955941B (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
CN111950149A (zh) | 基于参数化水平集法的连续体结构非概率拓扑优化方法 | |
CN107729648A (zh) | 一种基于Shepard插值的曲线纤维复合结构设计瀑布型多级优化方法 | |
CN108629140B (zh) | 一种基于测地线距离的带孔复合材料结构设计优化方法 | |
CN108829914B (zh) | 一种frp结构件的结构与工艺一体化设计方法 | |
CN107590325A (zh) | 一种基于Shepard插值的纤维增强复合材料结构优化方法 | |
JP6736108B2 (ja) | しわを抑制する空間薄膜構造治具形状を最適化する設計方法 | |
CN109460577B (zh) | 一种保证制造质量的变刚度复合材料结构优化设计方法 | |
CN111079327B (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
Kussmaul et al. | A novel computational framework for structural optimization with patched laminates | |
CN108491591A (zh) | 一种高温环境下曲线加筋板有限元分析方法 | |
CN112364546B (zh) | 一种基于双边滤波的纤维增强复合材料结构优化设计方法 | |
CN113505405B (zh) | 等效荷载获取方法、基于等效荷载的拓扑优化方法及系统 | |
CN111444579B (zh) | 一种考虑可制造性的复合材料结构优化设计方法 | |
CN117497108A (zh) | 基于主应力方向改进的连续纤维复合材料多尺度优化方法 | |
CN112131770A (zh) | 一种考虑可靠性的功能梯度连续体结构轻量化设计方法 | |
CN113779825A (zh) | 一种基于单位分解及移动控制点的结构优化设计方法 | |
Opgenoord et al. | Aeroelastic tailoring using additively manufactured lattice structures | |
Zhang et al. | Combining density-based approach and optimization refinement in the design of morphing airfoil structures | |
CN115238396B (zh) | 三维气动外形反设计方法及装置 | |
CN112818576A (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 |