CN113442441B - 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 - Google Patents
纤维取向与结构并行优化的连续纤维3d打印路径规划方法 Download PDFInfo
- Publication number
- CN113442441B CN113442441B CN202110875031.1A CN202110875031A CN113442441B CN 113442441 B CN113442441 B CN 113442441B CN 202110875031 A CN202110875031 A CN 202110875031A CN 113442441 B CN113442441 B CN 113442441B
- Authority
- CN
- China
- Prior art keywords
- fiber
- printing
- sub
- path
- optimized
- 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
- 239000000835 fiber Substances 0.000 title claims abstract description 146
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000010146 3D printing Methods 0.000 title claims abstract description 53
- 238000005457 optimization Methods 0.000 title claims abstract description 48
- 239000000463 material Substances 0.000 claims abstract description 73
- 238000007639 printing Methods 0.000 claims abstract description 66
- 239000003733 fiber-reinforced composite Substances 0.000 claims abstract description 27
- 239000002131 composite material Substances 0.000 claims abstract description 25
- 238000013461 design Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000004519 manufacturing process Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 3
- 108020001568 subdomains Proteins 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920001169 thermoplastic Polymers 0.000 description 2
- 239000004416 thermosoftening plastic Substances 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000007847 structural defect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
- B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
- B29C64/00—Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
- B29C64/30—Auxiliary operations or equipment
- B29C64/386—Data acquisition or data processing for additive manufacturing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y50/00—Data acquisition or data processing for additive manufacturing
-
- 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/10—Additive manufacturing, e.g. 3D printing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Materials Engineering (AREA)
- Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
Abstract
一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,构建纤维取向与复合材料结构并行优化模型,以材料密度和纤维角度作为设计变量,对其过滤得到纤维取向与复合材料结构并行优化结构;将复杂的优化结构离散成简单几何形状的结构,采用拓扑学思想将每个子区域抽象成点,根据子区域所属优化结构的位置关系将点与点连接,组成含有优化结构特征信息的连通图,将路径规划归为寻找连通图中的哈密顿路径;将材料密度值作为其纤维角度的权重因子,得到优化结构每个子区域内纤维轨迹方向;将子区域打印路径根据哈密顿路径顺序连接,生成打印代码;本发明发挥连续纤维增强复合材料的各向异性力学性能,满足3D打印工艺的要求。
Description
技术领域
本发明属于结构优化、复合材料和增材制造交叉技术领域,具体涉及一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法。
技术背景
连续纤维增强复合材料作为一种各向异性材料,是用于先进结构制造的重要材料。连续纤维增强复合材料3D打印工艺打破了传统模压、铺放技术在纤维角度方向上的约束,以路径规划信息打印模型,可以实现纤维取向的精细化控制和自由设计。采用熔融沉积成形(FDM)技术类似的原理,将纤维干丝与热塑性丝材在线浸渍并通过打印喷头进行打印,其力学性能根据纤维路径的打印角度产生相应变化。由于连续纤维增强复合材料沿纤维方向的力学性能远优于垂直纤维方向的力学性能,因此3D打印工艺中纤维路径打印方向对构件的整体性能有很大影响。
连续纤维增强复合材料3D打印工艺使材料、结构并行优化和一体化成形成为可能。虽然国内外对基于复合材料各向异性的优化设计有一定的理论研究,但是并未和3D打印工艺技术紧密结合,优化设计中包含的微观纤维分布、宏观拓扑结构信息与制造技术的融合仍面临巨大挑战,无法发挥出先进成形工艺与纤维增强性能的巨大潜力,缺乏相应的路径规划方法。
目前针对连续纤维增强复合材料3D打印工艺技术尚不完善,其路径规划方法多采用栅格轮廓填充、轮廓偏置路径填充、混合路径填充等传统FDM工艺手段,未能充分考虑连续纤维增强复合材料的各向异性力学性能。对于纤维取向与复合材料结构的并行优化设计,其模型中包含了有限元离散网格下各单元材料密度和纤维角度,充分考虑了纤维方向对结构性能的影响。然而,这种结构优化方法支持任意形状输出作为优化结果,其结果并非直接可行,通过传统的3D打印路径规划方法不能同时实现优化结构的宏观拓扑几何特征和微观纤维取向,还会出现打印过程中转角过小、路径跳转、路径重叠等问题,从而严重影响优化结构的力学性能,限制了连续纤维增强复合材料3D打印工艺的发展。
发明内容
为了克服上述现有技术的缺陷,本发明的目的在于提供一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,得到既能实现宏观拓扑几何结构特征,又充分考虑微观纤维增强方向特性的无重叠、无跳转的连续纤维3D打印路径,充分发挥连续纤维增强复合材料的各向异性力学性能,满足连续纤维增强复合材料3D打印工艺的要求。
为了达到上述目的,本发明采取的技术方案为:
一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,包括以下步骤:
1)构建纤维取向与复合材料结构并行优化模型,以材料密度x和纤维角度θ作为设计变量,其数学模型如下:
式中,目标函数c代表最小柔度值;U和F分别表示整体位移向量和整体载荷向量;K表示整体刚度矩阵;ue和ke分别表示单元位移向量和单元刚度矩阵;xmin表示相对最低密度;p表示惩戒因子;N表示有限元划分网格的数量;V(x)和V0分别表示材料体积和设计域初始体积;f表示设置的复合材料容许体积比;vf表示连续纤维在复合材料内的体积分数;通过求解得到有限元单元网格下的材料密度xe和纤维角度θe;
2)对步骤1)得到的有限元单元网格下的材料密度xe和纤维角度θe进行过滤处理,设置材料过滤密度xset,当xe<xset时,令xe等于0,此时材料密度xe不为0的单元网格组成宏观拓扑几何结构;根据连续纤维增强复合材料的正交各向异性,将步骤1)中得到的纤维角度θe的取值范围为从[-2π,2π]调整至[-π/2,π/2]或[0,π],得到纤维取向与结构并行优化的结果,即优化结构;
3)将步骤2)得到的优化结构根据其宏观几何特征划分子区域,由此将复杂的优化结构离散成有限个具有简单几何形状特征的结构;采用拓扑学思想将每个子区域抽象成点,并根据子区域所属优化结构的位置关系将点与点之间相连接,相互连接的点图组成了含有优化结构特征信息的连通图,对优化结构的路径规划问题归为图论中寻找连通图中的哈密顿路径问题,即在连通图中寻找不重复遍历所有点的路径;
4)对步骤3)划分的子区域内纤维轨迹方向进行计算,根据子区域的几何特征将其分割成n个区间,每个区间内包含ne个单元网格;然后将步骤2)处理后的并行优化结果代入,区间内单元网格的材料密度xe值作为其纤维角度θe的权重因子,由此得到每个区间的纤维轨迹方向θi,如式(3)所示;
式中,θj表示相邻区间的轨迹方向,ε表示区间之间的角度偏差,根据子区域的几何特征ε值相应变化;当步骤2)中纤维角度θe取值范围的不同时,得到的区间纤维轨迹角度θi也会有所不同,应以取值范围边界和区间内θe值相隔较远原则选择特定取值范围,使计算得到的纤维轨迹角度θi值尽可能与所在区间内并行优化计算结果θe相匹配;另外,在子区域的局部网格内可能出现并行优化计算得到的纤维角度和子区域的几何特征不符,此时计算的区间纤维轨迹角度应不包含这些局部单元网格,得到优化结构每个子区域内的纤维轨迹方向;
5)根据步骤4)得到的子区域纤维轨迹方向对子区域内材料进行铺放,其铺放间距即为该子区域的不同区间内3D打印扫描间距hi,i=1,2,…,n;根据连续纤维增强复合材料3D打印工艺,扫描间距h和纤维体积分数vf存在比例关系,而纤维体积分数vf作为并行优化的参数在步骤1)中已经确定,因此将3D打印工艺最大纤维含量下扫描间距h(vfmax)和并行优化中设置的纤维含量下扫描间距h(vf)作为3D打印工艺的制造约束,即h(vfmax)≤hi≤h(vf),以满足连续纤维增强复合材料结构满足目标性能需求和3D打印工艺要求,而区间内扫描间距hi应尽可能接近并行优化中设置纤维含量下扫描间距h(vf);
6)根据步骤5)的制造约束对扫描间距hi进行调节,在满足子区域结构特征的基础上结合步骤3)得到的含有优化结构特征信息的连通图确定铺放在子区域内的铺放纤维的条数,使连通图中基于哈密顿路径相连接的子区域内所包含的路径条数相等,由此得到优化结构的每个子区域连续纤维3D打印路径;
7)将步骤6)得到的子区域打印路径根据步骤3)得到的哈密顿路径顺序连接,连接时根据具体材料和设备特性确定允许最小转角半径rmin作为另一制造约束,使子区域连接时的打印半径rprint≥rmin;
8)经过以上的路径规划后,根据基体材料进给量和打印距离的比例将路径输出,生成打印G-code代码用于3D打印机进行打印。
所述的步骤3)中如果含有优化结构特征信息的连通图中不存在哈密顿路径,或所含有的哈密顿路径无法规划打印路径,则将连通图中已有的部分点即子区域分割为多个点,即多个子区域;或在连通图中增添点,即增加新的子区域,具体指在优化结构内增加新的结构,直至含有优化结构特征信息的连通图中至少存在一种哈密度路径作为连续纤维3D打印路径的规划依据。
所述的步骤5)中根据每个区间纤维轨迹方向θi和扫描间距hi初步得到的子区域的打印路径,采用不同的子区域分割区间范围,得到的打印路径也会有所不同;如果打印路径不能展现子区域的形状特征,则回到步骤4)的调整分割区间范围,直至每个子区域内都至少存在一种打印路径可以实现优化结构的宏观几何特征和3D工艺打印要求。
本发明的有益效果为:
本发明将连续纤维增强复合材料、结构优化与3D打印工艺相结合,完成对纤维取向与复合材料结构并行优化的连续纤维3D打印路径规划。与现有技术相比,本发明将哈密顿路径思想引入3D打印路径规划,将复杂结构的打印路径规划问题看作连通图哈密度路径问题,基于哈密顿路径生成的打印路径在打印过程中没有跳转点,极大地提高了对复杂构件的打印效率;同时充分考虑纤维取向与复合材料结构并行优化中纤维角度和材料密度,基于宏观优化结构的拓扑几何特征、微观纤维分布方向与3D打印工艺制造约束规划连续纤维打印路径,得到的打印路径可以充分发挥连续纤维增强复合材料各向异性的力学性能优势,并满足3D打印工艺的要求。
本发明具有良好的适用性,将各向异性复合材料的方向特性引入结构优化和3D打印工艺,在打印过程中实现复合材料纤维取向的精确调控,有效解决目前连续纤维增强复合材料轻量化、高承载复杂结构的制造问题,从而推进连续纤维增强复合材料结构优化及3D打印技术的发展。
附图说明
图1是本发明的流程图。
图2是本发明纤维取向与复合材料结构并行优化的示意图。
图3是本发明对优化结构进行子区域划分处理的示意图。
图4是本发明含有优化结构特征信息的哈密顿路径示意图。
图5是本发明对子区域进行材料铺放的示意图。
图6是本发明路径规划的示意图。
图7是本发明基于纤维取向与复合材料结构并行优化的3D打印构件。
具体实施方式
以下结合附图和实施例对本发明作进一步的详细说明
参照图1,一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,包括以下步骤:
1)构建纤维取向与复合材料结构并行优化模型,以材料密度x和纤维角度θ作为设计变量,其数学模型如下:
式中,目标函数c代表最小柔度值;U和F分别表示整体位移向量和整体载荷向量;K表示整体刚度矩阵;ue和ke分别表示单元位移向量和单元刚度矩阵;xmin表示相对最低密度;p表示惩戒因子;N表示有限元划分网格的数量;V(x)和V0分别表示材料体积和设计域初始体积;f表示设置的复合材料容许体积比;vf表示连续纤维在复合材料内的体积分数;通过求解得到有限元单元网格下的材料密度xe和纤维角度θe;
2)对步骤1)得到的有限元单元网格下的材料密度xe和纤维角度θe进行过滤处理,设置材料过滤密度xset,当xe<xset时,令xe等于0,此时材料密度xe不为0的单元网格组成宏观拓扑几何结构;根据连续纤维增强复合材料的正交各向异性,将步骤1)中得到的纤维角度θe的取值范围为从[-2π,2π]调整至[-π/2,π/2]或[0,π],两种取值范围在步骤4)中进行选取,得到纤维取向与结构并行优化的结果,即优化结构;
本实施例通过对复合材料的结构设计得到纤维取向与复合材料结构并行优化的结构模型,如图2所示,以经典模型MBB梁为例进行说明,通过有限元分析和优化求解得到每个单元网格内的材料密度和纤维角度,通过过滤处理将多余单元网格删除,留下图2中的阴影部分即纤维增强复合材料MBB梁的优化结构;由于对纤维取向与复合材料结构的优化为并行处理,其结构包含的每个单元网格内都含有其相对应的纤维角度值;
3)将步骤2)得到的优化结构根据其宏观几何特征划分子区域,由此将复杂的优化结构离散成有限个具有简单几何形状特征的结构;采用拓扑学思想将每个子区域抽象成点,并根据子区域所属优化结构的位置关系将点与点之间相连接,相互连接的点图组成了含有优化结构特征信息的连通图,对优化结构的路径规划问题归为图论中寻找连通图中的哈密顿路径问题,即在连通图中寻找不重复遍历所有点的路径;
考虑优化结构的复杂性,如果含有优化结构特征信息的连通图中不存在哈密顿路径,或所含有的哈密顿路径无法规划打印路径,则将连通图中已有的部分点即子区域分割为多个点,即多个子区域;或在连通图中增添点,即增加新的子区域,具体指在优化结构模型内增加新的结构,直至含有优化结构特征信息的连通图中至少存在一种哈密度路径作为连续纤维3D打印路径的规划依据;
本实施例根据优化结构的几何形状特征划分子区域,如图3所示,对优化结构的几何形状进行判断和解析,根据形状特征选取图2中的细杆结构作为区域划分的基础,由此将复杂的优化结构拆解为16个简单的杆状结构,如图3中框选出的16个子区域;
本实施例采用拓扑学思想将每个子区域等价成点,点与点的连接关系代表优化结构的交点信息,由此将复杂的拓扑几何结构转化为连通图,如图4所示,连通图内的每个节点都代表其数字相对应的子区域,由此将路径规划问题转化为建立连通图中的哈密顿路径问题,哈密顿路径为不重复遍历所有点的路径,这里指通过所有子区域的打印路径,如图4所示,构建新的子区域(*1、*2、*3)作为哈密顿路径中新的连接点,子区域*1为子区域3的衍生区域,子区域*2为子区域2的衍生区域,子区域*3为子区域16的衍生区域;以连通图中的节点3作为起始点,依次连接子区域节点3,16,1,2,*1,11,14,15,*3,13,12,7,6,*2,4,5,8,9,10,得到含有优化结构特征信息的哈密顿路径,定义该哈密顿路径为优化结构的路径规划依据;
4)由于纤维角度θe在并行优化的计算结果中为离散值,因此不能直接得到打印路径上的纤维轨迹方向,对步骤3)划分的子区域内纤维轨迹方向进行计算,根据子区域的几何特征将其分割成n个区间,每个区间内包含ne个单元网格;然后将步骤2)处理后的并行优化结果代入,区间内单元网格的材料密度xe值作为其纤维角度θe的权重因子,由此得到每个区间的纤维轨迹方向θi,如式(3)所示;
式中,θj表示相邻区间的轨迹方向,ε表示区间之间的角度偏差,根据子区域的几何特征ε值相应变化;当步骤2)中纤维角度θe取值范围的不同时,得到的区间纤维轨迹角度θi也会有所不同,应以取值范围边界和区间内θe值相隔较远原则选择特定取值范围,使计算得到的纤维轨迹角度θi值尽可能与所在区间内并行优化计算结果θe相匹配,以最大程度上发挥连续纤维增强复合材料的各向异性性能优势;另外,在子区域的局部网格内可能出现并行优化计算得到的纤维角度和子区域的几何特征不符,此时计算的区间纤维轨迹角度应不包含这些局部单元网格,以提高优化结构的可制造性;得到优化结构每个子区域内的纤维轨迹方向,为子区域内连续纤维路径的规划提供依据;
本实施例对步骤3)划分得到的子区域进行几何形状判断和解析,根据形状特征将子区域分割为有限个区间,如图5所示1/2MBB梁子区域材料铺放,对每个区间所包含的单元纤维优化角度θe进行修正,以材料密度xe值作为其纤维角度θe的权重因子代入计算得到每个区间的纤维轨迹方向θi;如图5所示,部分子区域3和子区域12根据其形状特征以不同区间形式分割,以此保证铺放材料的整体方向和子区域的结构几何特征相匹配,不同区间内纤维铺放方向根据计算θi值相应变化;
5)根据步骤4)得到的子区域纤维轨迹方向对子区域内材料进行铺放,其铺放间距即为该子区域的不同区间内3D打印扫描间距hi,i=1,2,…,n;根据连续纤维增强复合材料3D打印工艺,扫描间距h和纤维体积分数vf存在比例关系,而纤维体积分数vf作为并行优化的参数在步骤1)中已经确定,因此将3D打印工艺最大纤维含量下扫描间距h(vfmax)和并行优化中设置的纤维含量下扫描间距h(vf)作为3D打印工艺的制造约束,即h(vfmax)≤hi≤h(vf),以满足连续纤维增强复合材料结构满足目标性能需求和3D打印工艺要求,其中,hi小于h(vfmax)保证了打印路径在子区域内不会出现路径叠加,提高制造结构的力学性能,而区间内扫描间距hi应尽可能接近并行优化中设置纤维含量下扫描间距h(vf),以减少材料损耗;
根据每个区间纤维轨迹方向θi和扫描间距hi初步得到的子区域的打印路径,采用不同的子区域分割区间范围,得到的打印路径也会有所不同;如果打印路径不能展现子区域的形状特征,则回到步骤4)的调整分割区间范围,直至每个子区域内都至少存在一种打印路径可以实现优化结构的宏观几何特征和3D工艺打印要求;
本实施例对于3D打印连续纤维增强复合材料,纤维体积分数通过3D打印工艺参数的材料系数、打印间距和层厚相对应,打印间距与纤维体积分数呈负相关关系,和如图5所示的局部结构放大图,轨迹之间的距离为打印间距,打印间距hi需要小于并尽可能接近步骤1)中代入计算的纤维体积分数所对应的打印间距值h(vf),以满足路径规划的优化结构的力学性能;且需要小于3D打印工艺最大体积分数所对应的打印间距值h(vfmax),以满足3D打印工艺的制造要求;
6)根据步骤5)的制造约束对扫描间距hi进行调节,在满足子区域结构特征的基础上结合步骤3)得到的含有优化结构特征信息的连通图确定铺放在子区域内的铺放纤维的条数,使连通图中基于哈密顿路径相连接的子区域内所包含的路径条数相等,由此得到优化结构的每个子区域连续纤维3D打印路径;
7)将步骤6)得到的子区域打印路径根据步骤3)得到的哈密顿路径顺序连接,连接时根据具体材料和设备特性确定允许最小转角半径rmin作为另一制造约束,使子区域连接时的打印半径rprint≥rmin,以此避免打印过程中转角过小造成的结构缺陷;根据哈密顿路径的特性,根据以上步骤得到的打印路径在打印复杂几何结构时不存在喷头跳转,极大地提高了零件打印效率和材料成形效果;
本实施例根据步骤3)的哈密顿路径信息按顺序依次连接各子区域,相连子区域之间应具有相同纤维根数,可根据步骤5)的打印间距约束调整轨迹间距,使相邻子区域内的轨迹一一对应;同时,设置另一制造约束最小打印半径,使子区域连接处路径回转半径不低于最小打印半径,保证优化结构的打印质量,得到纤维取向与复合材料结构并行优化的打印路径,如图6所示;
8)经过以上的路径规划后,根据基体材料进给量和打印距离的比例将路径输出,生成打印G-code代码用于3D打印机进行打印;
本实施例根据基体材料进给量和打印距离的比例将路径输出,生成打印G-code代码用于3D打印机进行打印,基于纤维取向与复合材料结构并行优化的3D打印构件如图7所示,最终得到既能实现宏观拓扑几何结构特征,又充分考虑微观连续纤维分布方向的无重叠、无跳转的连续纤维3D打印路径,满足3D打印工艺的要求。
Claims (3)
1.一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,其特征在于,包括以下步骤:
1)构建纤维取向与复合材料结构并行优化模型,以材料密度x和纤维角度θ作为设计变量,其数学模型如下:
式中,目标函数c代表最小柔度值;U和F分别表示整体位移向量和整体载荷向量;K表示整体刚度矩阵;ue和ke分别表示单元位移向量和单元刚度矩阵;xmin表示相对最低密度;p表示惩戒因子;N表示有限元划分网格的数量;V(x)和V0分别表示材料体积和设计域初始体积;f表示设置的复合材料容许体积比;vf表示连续纤维在复合材料内的体积分数;通过求解得到有限元单元网格下的材料密度xe和纤维角度θe;
2)对步骤1)得到的有限元单元网格下的材料密度xe和纤维角度θe进行过滤处理,设置材料过滤密度xset,当xe<xset时,令xe等于0,此时材料密度xe不为0的单元网格组成宏观拓扑几何结构;根据连续纤维增强复合材料的正交各向异性,将步骤1)中得到的纤维角度θe的取值范围为从[-2π,2π]调整至[-π/2,π/2]或[0,π],得到纤维取向与结构并行优化的结果,即优化结构;
3)将步骤2)得到的优化结构根据其宏观几何特征划分子区域,由此将复杂的优化结构离散成有限个具有简单几何形状特征的结构;采用拓扑学思想将每个子区域抽象成点,并根据子区域所属优化结构的位置关系将点与点之间相连接,相互连接的点图组成了含有优化结构特征信息的连通图,对优化结构的路径规划问题归为图论中寻找连通图中的哈密顿路径问题,即在连通图中寻找不重复遍历所有点的路径;
4)对步骤3)划分的子区域内纤维轨迹方向进行计算,根据子区域的几何特征将其分割成n个区间,每个区间内包含ne个单元网格;然后将步骤2)处理后的并行优化结果代入,区间内单元网格的材料密度xe值作为其纤维角度θe的权重因子,由此得到每个区间的纤维轨迹方向θi,如式(3)所示;
式中,θj表示相邻区间的轨迹方向,ε表示区间之间的角度偏差,根据子区域的几何特征ε值相应变化;当步骤2)中纤维角度θe取值范围的不同时,得到的区间纤维轨迹角度θi也会有所不同,应以取值范围边界和区间内θe值相隔较远原则选择特定取值范围,使计算得到的纤维轨迹角度θi值尽可能与所在区间内并行优化计算结果θe相匹配;另外,在子区域的局部网格内可能出现并行优化计算得到的纤维角度和子区域的几何特征不符,此时计算的区间纤维轨迹角度应不包含这些局部单元网格,得到优化结构每个子区域内的纤维轨迹方向;
5)根据步骤4)得到的子区域纤维轨迹方向对子区域内材料进行铺放,其铺放间距即为该子区域的不同区间内3D打印扫描间距hi,i=1,2,…,n;根据连续纤维增强复合材料3D打印工艺,扫描间距h和纤维体积分数vf存在比例关系,而纤维体积分数vf作为并行优化的参数在步骤1)中已经确定,因此将3D打印工艺最大纤维含量下扫描间距h(vfmax)和并行优化中设置的纤维含量下扫描间距h(vf)作为3D打印工艺的制造约束,即h(vfmax)≤hi≤h(vf),以满足连续纤维增强复合材料结构满足目标性能需求和3D打印工艺要求,而区间内扫描间距hi应尽可能接近并行优化中设置纤维含量下扫描间距h(vf);
6)根据步骤5)的制造约束对扫描间距hi进行调节,在满足子区域结构特征的基础上结合步骤3)得到的含有优化结构特征信息的连通图确定铺放在子区域内的铺放纤维的条数,使连通图中基于哈密顿路径相连接的子区域内所包含的路径条数相等,由此得到优化结构的每个子区域连续纤维3D打印路径;
7)将步骤6)得到的子区域打印路径根据步骤3)得到的哈密顿路径顺序连接,连接时根据具体材料和设备特性确定允许最小转角半径rmin作为另一制造约束,使子区域连接时的打印半径rprint≥rmin;
8)经过以上的路径规划后,根据基体材料进给量和打印距离的比例将路径输出,生成打印G-code代码用于3D打印机进行打印。
2.根据权利要求1所述的一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,其特征在于:所述的步骤3)中如果含有优化结构特征信息的连通图中不存在哈密顿路径,或所含有的哈密顿路径无法规划打印路径,则将连通图中已有的部分点即子区域分割为多个点,即多个子区域;或在连通图中增添点,即增加新的子区域,具体指在优化结构内增加新的结构,直至含有优化结构特征信息的连通图中至少存在一种哈密度路径作为连续纤维3D打印路径的规划依据。
3.根据权利要求1所述的一种纤维取向与结构并行优化的连续纤维3D打印路径规划方法,其特征在于:所述的步骤5)中根据每个区间纤维轨迹方向θi和扫描间距hi初步得到的子区域的打印路径,采用不同的子区域分割区间范围,得到的打印路径也会有所不同;如果打印路径不能展现子区域的形状特征,则回到步骤4)的调整分割区间范围,直至每个子区域内都至少存在一种打印路径可以实现优化结构的宏观几何特征和3D工艺打印要求。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110875031.1A CN113442441B (zh) | 2021-07-30 | 2021-07-30 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
PCT/CN2021/129415 WO2023005052A1 (zh) | 2021-07-30 | 2021-11-08 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110875031.1A CN113442441B (zh) | 2021-07-30 | 2021-07-30 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113442441A CN113442441A (zh) | 2021-09-28 |
CN113442441B true CN113442441B (zh) | 2022-05-06 |
Family
ID=77817820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110875031.1A Active CN113442441B (zh) | 2021-07-30 | 2021-07-30 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113442441B (zh) |
WO (1) | WO2023005052A1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113442441B (zh) * | 2021-07-30 | 2022-05-06 | 西安交通大学 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
CN114274500B (zh) * | 2021-12-23 | 2022-09-30 | 西安交通大学 | 一种基于绝对零刚度结构的隔振鞋中底的3d打印制作方法 |
CN115194931B (zh) * | 2022-09-14 | 2022-12-30 | 中电建冀交高速公路投资发展有限公司 | 混凝土3d打印路径的规划方法、装置、设备及存储介质 |
CN117565404B (zh) * | 2023-12-04 | 2024-05-28 | 重庆大学溧阳智慧城市研究院 | 基于Voronoi多边形骨架的3D打印多孔结构路径规划方法 |
CN117690532A (zh) * | 2023-12-22 | 2024-03-12 | 华中科技大学 | 变厚度蒙皮的纤维增强复合结构拓扑优化设计方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766656A (zh) * | 2019-01-25 | 2019-05-17 | 北京航空航天大学 | 一种基于拓扑优化的梯度点阵结构设计方法 |
CN111319268A (zh) * | 2020-02-20 | 2020-06-23 | 西北工业大学 | 一种考虑增材制造打印方向的自支撑结构优化设计方法 |
CN111950149A (zh) * | 2020-08-13 | 2020-11-17 | 北京航空航天大学 | 基于参数化水平集法的连续体结构非概率拓扑优化方法 |
CN112883616A (zh) * | 2021-02-26 | 2021-06-01 | 山东大学 | 一种面向纤维增强结构的3d打印喷头路径优化方法 |
CN112966410A (zh) * | 2021-02-03 | 2021-06-15 | 西安交通大学 | 一种适应于临界角可变的增材制造自支撑结构拓扑优化方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9688028B2 (en) * | 2013-03-22 | 2017-06-27 | Markforged, Inc. | Multilayer fiber reinforcement design for 3D printing |
US11084223B2 (en) * | 2018-08-20 | 2021-08-10 | Lawrence Livermore National Security, Llc | Optimal toolpath generation system and method for additively manufactured composite materials |
US20200156323A1 (en) * | 2018-11-20 | 2020-05-21 | Arevo, Inc. | Systems and methods for optimization of design and tool paths for additive manufacturing |
CN110001067B (zh) * | 2019-03-27 | 2022-01-18 | 北京机科国创轻量化科学研究院有限公司 | 一种连续纤维增强复合材料3d打印路径规划方法 |
US20200356638A1 (en) * | 2019-05-07 | 2020-11-12 | Toyota Motor Engineering & Manufacturing North America, Inc. | Orientation optimization in components fabricated with anisotropic material properies |
US11009853B2 (en) * | 2019-07-29 | 2021-05-18 | Toyota Motor Engineering & Manufacturing North America, Inc. | Method of tool path generation for additive manufacturing with vector distribution |
CN110955941B (zh) * | 2019-11-29 | 2022-03-18 | 华中科技大学 | 基于向量场的复合材料结构优化设计方法及设备 |
CN111444579B (zh) * | 2020-03-11 | 2022-04-12 | 华中科技大学 | 一种考虑可制造性的复合材料结构优化设计方法 |
CN112989648B (zh) * | 2021-02-04 | 2022-11-25 | 西安理工大学 | 一种协同拓扑构型与纤维路径的柔性机构优化设计方法 |
CN113442441B (zh) * | 2021-07-30 | 2022-05-06 | 西安交通大学 | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 |
-
2021
- 2021-07-30 CN CN202110875031.1A patent/CN113442441B/zh active Active
- 2021-11-08 WO PCT/CN2021/129415 patent/WO2023005052A1/zh unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766656A (zh) * | 2019-01-25 | 2019-05-17 | 北京航空航天大学 | 一种基于拓扑优化的梯度点阵结构设计方法 |
CN111319268A (zh) * | 2020-02-20 | 2020-06-23 | 西北工业大学 | 一种考虑增材制造打印方向的自支撑结构优化设计方法 |
CN111950149A (zh) * | 2020-08-13 | 2020-11-17 | 北京航空航天大学 | 基于参数化水平集法的连续体结构非概率拓扑优化方法 |
CN112966410A (zh) * | 2021-02-03 | 2021-06-15 | 西安交通大学 | 一种适应于临界角可变的增材制造自支撑结构拓扑优化方法 |
CN112883616A (zh) * | 2021-02-26 | 2021-06-01 | 山东大学 | 一种面向纤维增强结构的3d打印喷头路径优化方法 |
Non-Patent Citations (3)
Title |
---|
Mei Yulin ; Wang Xiaoming.A level set method for structural topology optimization and its application.《advances in Engineering Software》.2004, * |
增材制造中的结构设计与路径规划问题;翟晓雅;《博士电子期刊》;20210520;1-107 * |
高性能树脂基复合材料轻质结构3D打印与性能研究;田小永;候章浩;张俊康;李涤尘;《航空制造技术》;20170515;34-39 * |
Also Published As
Publication number | Publication date |
---|---|
CN113442441A (zh) | 2021-09-28 |
WO2023005052A1 (zh) | 2023-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113442441B (zh) | 纤维取向与结构并行优化的连续纤维3d打印路径规划方法 | |
Padhi et al. | Optimization of fused deposition modeling process parameters using a fuzzy inference system coupled with Taguchi philosophy | |
Fernandez et al. | Optimal design of fiber reinforced composite structures and their direct ink write fabrication | |
Wong et al. | Additive manufacturing of fiber-reinforced polymer composites: A technical review and status of design methodologies | |
Maute et al. | Adaptive topology optimization of shell structures | |
Li et al. | Design for additive manufacturing from a force-flow perspective | |
Sales et al. | Function-aware slicing using principal stress line for toolpath planning in additive manufacturing | |
CN110083900A (zh) | 一种面向混杂纤维复合材料板壳结构的快速协同优化方法 | |
CN112818470B (zh) | 一种基结构的优化方法、装置、计算机设备和存储介质 | |
CN113191077B (zh) | 一种基于连续纤维复材3d打印的变纤维含量拓扑优化方法 | |
US20200159186A1 (en) | Additive manufacturing system and method having toolpath analysis | |
CN117150858B (zh) | 含裂纹连续纤维增强树脂基复合材料三维有限元建模方法 | |
CN113192576B (zh) | 一种2.5d编织复合材料单胞模型的建模方法 | |
US20160121558A1 (en) | Method for defining fiber trajectories from a vector field | |
US9481135B2 (en) | Method for defining fiber trajectories from a transfer surface | |
CN110126266A (zh) | 一种三维物体制造方法 | |
Zhang et al. | A nodal-based optimization method for the design of continuous fiber-reinforced structures | |
CN111027151B (zh) | 一种复合材料异形壳的纤维路径与几何形状一体化设计方法 | |
Zhang et al. | A 3D printing tool-path generation strategy based on the partition of principal stress field for fused filament fabrication | |
CN105138750A (zh) | 一种可扩展多线交叉机身线框模型构建方法 | |
Naresh et al. | Design and development of alternate layer printing method to reduce the porosity in FDM printing process | |
CN115688462A (zh) | 一种正高斯曲面正交标架铺丝成型轨迹规划设计方法 | |
CN114741919A (zh) | 一种基于力流管载荷路径的3d打印填充设计方法 | |
CN116330706A (zh) | 一种复合材料自动铺放的工艺路径生成系统 | |
CN112364546B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231205 Address after: 518000 301-a10, building 3, Huahan science and Technology Industrial Park, No. 19, LanJin 4th Road, Heping community, Pingshan street, Pingshan District, Shenzhen, Guangdong Patentee after: Shenzhen Yunjiang Zhizao Technology Co.,Ltd. Address before: Beilin District Xianning West Road 710049, Shaanxi city of Xi'an province No. 28 Patentee before: XI'AN JIAOTONG University |
|
TR01 | Transfer of patent right |