一种FRP结构件的结构与工艺一体化设计方法技术

本发明专利技术涉及一种FRP结构件的结构与工艺一体化设计方法，包括：1、建立FRP结构件的有限元模型，参数化铺层角度和铺层厚度设计变量，建立材料插值模型；2、对结构件进行有限元分析，提取灵敏度计算所需参数，建立铺层角度和总厚度设计变量的目标函数及约束条件；3、计算目标函数和约束条件的灵敏度；4、根据灵敏度，采用数学规划方法迭代求出目标函数最优解，更新铺层角度和总厚度设计变量的参数化模型；5、重复步骤2～4，直至结果收敛或达到最大迭代次数，得到FRP结构件最优铺层结构。与现有技术相比，本发明专利技术实现了结构与工艺的一体化设计，在满足性能要求和制造约束的前提下，有效提高FRP结构件的材料利用率和优化设计效率。

【技术实现步骤摘要】
一种FRP结构件的结构与工艺一体化设计方法
本专利技术涉及复合材料结构件设计方法，尤其是涉及一种FRP结构件的结构与工艺一体化设计方法。
技术介绍
纤维复合材料具有比模量高，比强度高，耐腐蚀，抗疲劳性好，减震性好，密度低等优点。此外，纤维复合材料成型性好，可设计性强，可以通过改变纤维、基体种类、纤维体积分数、铺层角度、铺层厚度、铺层顺序以及结构拓扑等参数，达到结构件在实际使用过程中的工程性能要求，如刚度、强度、模态等性能指标。纤维复合材料可设计性强的特点使得其在轻量化方面具有巨大潜力，在诸如航空航天、汽车工业、船舶工程等学科领域得到了广泛的应用。然而，不同于传统各向同性的金属材料，纤维复合材料具有各向异性，因此传统的结构优化方法不能够充分发挥其轻量化潜力。以往的研究多基于工程经验进行简单替代，或仅对复合材料进行厚度优化，或运用代理模型的方法进行参数优化，或采用分级式优化，轻量化效果有限。有必要针对这些问题，深入研究纤维复合材料设计变量的参数化方法，采用合理的优化算法，对FRP结构件进行结构优化。对于纤维增强复合材料(FRP)，现有商业软件如Optistruct需经过自由尺寸优化、尺寸优化等过程进行分级式优化，该过程的缺点是铺层角度和铺层厚度的优化不能有效的统一起来进行优化，使得优化结构的优化效果并不显著，此外，该方法不能考虑拓扑结构会产生相同角度铺层过多连续、铺层结构不对称及中间空心等缺陷，导致优化结果与实际产品性能差别大。因此，针对FRP复合材料，应同时将铺层角度和铺层厚度作为优化设计变量，并在此基础上，引入制造工艺约束，实现复合FRP结构件的结构与工艺一体化设计，充分发挥FRP复合材料的性能。经过对现有技术的检索发现，中国专利文献号CN106202597A(公开日为2016年12月07日)公开了一种复合材料加筋壁板结构优化分析方法，该技术对结构进行分级优化，第一级进行拓扑优化，第二级在第一级优化的基础上进行铺层优化，确定铺层角度、厚度，结构减重明显。但该技术将拓扑优化与铺层优化分开进行，使得优化结果具有局限性，不能充分发挥复合材料的优异性能。中国专利文献CN105711106A(公开日为2016年6月29日)公开了一种轻量化的复合材料汽车尾门，采用了模块化的设计理念，充分发挥了纤维增强复合材料可设计性强、易成型复杂形状的特点，在保证汽车安全性的前提下，减少了汽车尾门的零部件数量，降低了汽车车身的重量。但是该技术对于结构设计是基于工程经验积累基础上进行的，缺乏理论方法的支持，不能保证设计结构为最优结构。
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种FRP结构件的结构与工艺一体化设计方法。本专利技术的目的可以通过以下技术方案来实现：一种FRP结构件的结构与工艺一体化设计方法，包括以下步骤：S1、建立FRP结构件的有限元模型，参数化铺层角度和铺层厚度设计变量，得到参数化模型，并赋予初始值，通过罚函数法建立材料插值模型；S2、对FRP结构件的各工况进行有限元分析，从有限元分析结果中提取灵敏度计算所需参数，建立铺层角度设计变量和总厚度设计变量的目标函数及约束条件；S3、计算得到目标函数、约束条件对于铺层角度设计变量及总厚度设计变量的灵敏度；S4、根据步骤S3得到的灵敏度，采用数学规划方法迭代求出步骤S2中目标函数的最优解，更新铺层角度和铺层厚度设计变量的参数化模型；S5、重复步骤S2～S4，直至结果收敛或达到最大迭代次数，得到FRP结构件最优的铺层结构。优选的，所述参数化铺层角度和铺层厚度设计变量具体包括：如果备选角度c存在于设计子域j铺层l当中，则xjlc＝1，否则xjlc＝0，以实现铺层角度设计变量的参数化；如果单元e的l层有铺层材料，则ρel＝1，否则ρel＝0，以实现铺层厚度设计变量的参数化；铺层角度及铺层厚度设计变量的参数化模型为：其中，Eel为属于设计子域j的单元e在第l层铺层的本构矩阵，xjlc为设计子域j在l层备选角度c下的铺层角度，ρel为单元e在l层的铺层厚度，E0是一种假想的弹性模量很低的人工材料，目的是避免在优化过程中材料本构矩阵发生奇异，Ec为备选角度c的材料本构矩阵，nc为备选角度的个数。优选的，所述通过罚函数法建立材料插值模型具体为：采用RAMP插值模型建立材料插值模型。优选的，所述RAMP插值模型为：其中，q为第一惩罚因子，p为第二惩罚因子。优选的，所述步骤S2中以结构件的各工况的加权柔度值最小为目标函数，所述约束条件包括：刚度约束，体积分数约束，相同铺层角度连续铺设层数约束，每层的铺层厚度约束，铺层结构对称约束。优选的，所述相同铺层角度连续铺设层数约束包括：其中，t为第t层铺层，xjlc为设计子域j在第l层铺层备选角度c下的铺层角度，ns是相同铺层角度连续铺设的最大层数。优选的，所述每层的铺层厚度约束是用于防止铺层结构产生中间空心，包括：其中，ρel为单元e在l层的铺层厚度，ρe为单元e的总厚度设计变量，β表示第三惩罚因子，s(l)是l层正交后的铺层坐标系，具体为：其中，nl为铺层总数。优选的，所述铺层结构对称约束包括：对称平面的对称单元具有相同的材料本构矩阵。优选的，所述目标函数对铺层角度设计变量的灵敏度为：其中，C为结构件整体柔度，xjlc为设计子域j在l层备选角度c下的铺层角度，q为第一惩罚因子，e为第e个单元，i为铺层界面上的积分点编号，Pj为第j个设计子域，nl为总铺层数，ni为积分点总个数，T为矩阵的转置，Veli和εeli分别为第i点的积分点权重系数及应变矩阵，Ec为备选角度c的材料本构矩阵；所述约束条件对铺层角度设计变量的灵敏度为：其中，dr为某一点的位移，χeli为在第i点的应变向量。所述目标函数对于总厚度设计变量的灵敏度为：其中，式中，ρe表示单元e的总厚度设计变量，Ne为属于单元e的过滤半径范围内的单元集合，ω(Xe)和ω(Xi)是权重因子项，Xe、Xi分别表示单元e和单元i的中心坐标，为第i个单元进行密度过滤后新的总厚度设计变量，ρil为根据该新的总厚度设计变量而设计的每一层的厚度变量，ve为中心单元的体积，vi为过滤半径范围内第i个单元的体积；约束条件对于总厚度设计变量的灵敏度为：优选的，所述步骤S2中通过密度过滤算法对总厚度设计变量进行处理，用于避免拓扑优化时数值不稳定现象的产生。与现有技术相比，本专利技术具有以下优点：1、对FRP材料进行了铺层角度及铺层厚度的参数化，该参数化方法可在单次优化过程中同时实现铺层角度及铺层厚度的优化，通过材料插值模型等技术实现了FRP结构件铺层角度和厚度的一体化设计，从而实现了结构与工艺的一体化设计，对复合材料结构件优化设计提供理论层次的指导，克服了以往仅能进行分级式优化或者基于经验的非最优FRP结构件问题，能够有效的提升FRP结构件的优化设计效率，产生更高的经济效益。2、由于增大了铺层角度、总厚度设计变量的结合空间，故较现有商业软件相比具有更高的材料利用率。3、通过约束条件增加对制造工艺的约束，避免了通过现有商业软件进行优化而产生的如中间空心、相同铺层角度连续铺设过多、铺层结构不对称等缺陷，从而更加符合工程实际要求。附图说明图1为本专利技术的流程示意图；图2为实施例中发动机罩的有限元模型示意图本文档来自技高网...

【技术保护点】
1.一种FRP结构件的结构与工艺一体化设计方法，其特征在于，包括以下步骤：S1、建立FRP结构件的有限元模型，参数化铺层角度和铺层厚度设计变量，得到参数化模型，并赋予初始值，通过罚函数法建立材料插值模型；S2、对FRP结构件的各工况进行有限元分析，从有限元分析结果中提取灵敏度计算所需参数，建立铺层角度设计变量和总厚度设计变量的目标函数及约束条件；S3、计算得到目标函数、约束条件对于铺层角度设计变量及总厚度设计变量的灵敏度；S4、根据步骤S3得到的灵敏度，采用数学规划方法迭代求出步骤S2中目标函数的最优解，更新铺层角度和铺层厚度设计变量的参数化模型；S5、重复步骤S2～S4，直至结果收敛或达到最大迭代次数，得到FRP结构件最优的铺层结构。

【技术特征摘要】
1.一种FRP结构件的结构与工艺一体化设计方法，其特征在于，包括以下步骤：S1、建立FRP结构件的有限元模型，参数化铺层角度和铺层厚度设计变量，得到参数化模型，并赋予初始值，通过罚函数法建立材料插值模型；S2、对FRP结构件的各工况进行有限元分析，从有限元分析结果中提取灵敏度计算所需参数，建立铺层角度设计变量和总厚度设计变量的目标函数及约束条件；S3、计算得到目标函数、约束条件对于铺层角度设计变量及总厚度设计变量的灵敏度；S4、根据步骤S3得到的灵敏度，采用数学规划方法迭代求出步骤S2中目标函数的最优解，更新铺层角度和铺层厚度设计变量的参数化模型；S5、重复步骤S2～S4，直至结果收敛或达到最大迭代次数，得到FRP结构件最优的铺层结构。2.根据权利要求1所述的一种FRP结构件的结构与工艺一体化设计方法，其特征在于，所述参数化铺层角度和铺层厚度设计变量具体包括：如果备选角度c存在于设计子域j铺层l当中，则xjlc＝1，否则xjlc＝0，以实现铺层角度设计变量的参数化；如果单元e的l层有铺层材料，则ρel＝1，否则ρel＝0，以实现铺层厚度设计变量的参数化；铺层角度及铺层厚度设计变量的参数化模型为：其中，Eel为属于设计子域j的单元e在第l层铺层的本构矩阵，xjlc为设计子域j在l层备选角度c下的铺层角度，ρel为单元e在l层的铺层厚度，E0是一种假想的弹性模量很低的人工材料，目的是避免在优化过程中材料本构矩阵发生奇异，Ec为备选角度c的材料本构矩阵，nc为备选角度的个数。3.根据权利要求2所述的一种FRP结构件的结构与工艺一体化设计方法，其特征在于，所述通过罚函数法建立材料插值模型具体为：采用RAMP插值模型建立材料插值模型。4.根据权利要求3所述的一种FRP结构件的结构与工艺一体化设计方法，其特征在于，所述RAMP插值模型为：其中，q为第一惩罚因子，p为第二惩罚因子。5.根据权利要求1所述的一种FRP结构件的结构与工艺一体化设计方法，其特征在于，所述步骤S2中以结构件的各工况的加权柔度值最小为目标函数，所述约束条件包括：刚度约束，体积分数约束，相同铺层角度连续铺设层数约束，每层的铺层...

【专利技术属性】
技术研发人员：高云凯，徐亚男，吴驰，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海,31