CN107220461A - 一种变刚度复合材料板壳结构高效优化方法 - Google Patents

一种变刚度复合材料板壳结构高效优化方法 Download PDF

Info

Publication number
CN107220461A
CN107220461A CN201710492492.4A CN201710492492A CN107220461A CN 107220461 A CN107220461 A CN 107220461A CN 201710492492 A CN201710492492 A CN 201710492492A CN 107220461 A CN107220461 A CN 107220461A
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
variation rigidity
composite panel
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
Application number
CN201710492492.4A
Other languages
English (en)
Other versions
CN107220461B (zh
Inventor
郝鹏
袁枭桀
王雨田
刘晨
王博
牛飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201710492492.4A priority Critical patent/CN107220461B/zh
Publication of CN107220461A publication Critical patent/CN107220461A/zh
Application granted granted Critical
Publication of CN107220461B publication Critical patent/CN107220461B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/24—Sheet material

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

一种变刚度复合材料板壳结构高效优化方法
技术领域
本发明属于航空航天结构中承力结构设计技术领域，涉及一种变刚度复合材料板壳结构高效优化方法。
背景技术
航天器承力结构轻量化可以减少制造结构所需的资源和能耗，这无疑能减少巨大的成本。由于较高的比刚度和比强度，复合材料板壳结构在现代航天工业复杂载荷承力结构的应用中愈发普遍。目前纤维增强复合材料层合板多采用平行顺直纤维铺放成型，且同层的纤维角度固定。考虑到制造工艺性约束，工程中常采用0°、90°和±45°的铺层方向，铺层纤维缺乏更为丰富的空间分布形式，也使结构的设计空间受限，不能完全发挥复合材料板壳的承载优势。
随着复合材料制造技术的发展，纤维自动铺放技术使得纤维路径为曲线的变刚度复合材料板壳的制备成为可能。国外的纤维自动铺放技术发展很快，有较多的大型纤维铺放机投入使用。相比之下，国内的纤维铺放成型技术和装备尚处于原理性研究和工程样机研制阶段。国内外相关技术和设备水平的提升，为纤维铺放技术创造了基本条件。
相比变刚度复合材料板壳结构的制造，其纤维路径设计也充满挑战。对于复合材料薄壁板壳结构，其抗屈曲性能直接决定了结构服役安全。但在传统的基于有限元方法的变刚度复合材料板壳结构屈曲分析中，由于单元离散导致变刚度设计中的纤维铺层角度不再连续光滑。该处理方式在分析大型结构的力学行为时将可能导致严重的预测误差。为保证仿真计算的结果精度，只能采用增加单元数和节点数的方式，导致模型自由度数超高、分析效率低下。更重要的是，由于传统有限元方法无法提供解析灵敏度信息，只能采用差分方法来近似计算，导致计算量极大且容易陷入局部最优解，常导致优化失败，难以挖掘变刚度复合材料板壳结构的承载潜力。因此，现有的变刚度复合材料板壳结构优化方法无法实现曲线纤维路径的快速设计，即使应用代理模型技术，直接导致优化设计能力不足、产品研发周期极长。相比之下，等几何方法利用NURBS基函数和控制点代替传统有限元方法中的形函数和节点，使结构的设计、分析和优化模型具有相同的几何描述，使得CAD和CAE纳入到统一的框架下。等几何方法容易构造高阶协调单元，特别适合研究对近似函数有高阶连续性要求的曲线铺层问题，可推导获得全解析的灵敏度信息，为后续优化设计提供高效手段。
发明内容
本发明主要解决现有技术中变刚度复合材料板壳结构优化设计效率低的技术问题，提出一种针对变刚度复合材料板壳结构的高效优化方法，达到显著提高变刚度复合材料板壳结构的承载效率、大幅降低产品研发周期的目的。
为了达到上述目的，本发明技术方案为：
一种变刚度复合材料板壳结构高效优化方法，包括以下步骤：
步骤100，建立变刚度复合材料板壳结构的等几何屈曲设计模型，包括以下子步骤：
步骤101，对特定的曲线纤维铺层，基于等几何方法实现纤维路径的参数化；所述的曲线纤维铺层包括：线性梯度函数、高阶曲面等高线函数、高阶Bézier曲线函数、拉格朗日多项式函数等。
步骤102，建立控制点角度与全场纤维路径的映射关系，并进行收敛性分析，较有限元方法以较少的计算时间保证屈曲分析的精确性。
步骤103，选定设计变量，建立变刚度复合材料板壳结构的等几何屈曲设计模型，得到等几何分析所需的输入数据；所述的设计变量为描述纤维路径的函数以及参数变量，所述的参数变量为描述纤维路径的角度。
步骤200，根据得到的等几何屈曲设计模型，进行变刚度复合材料板壳结构的等几何屈曲分析，得到屈曲载荷与屈曲模态，包括以下子步骤：
步骤201，根据有限元基本理论，将变刚度复合材料板壳结构的单元刚度矩阵组装成等几何屈曲设计模型的刚度矩阵K。
步骤202，对等几何屈曲设计模型进行静力分析，得到应力分布，并以此求得等几何屈曲设计的几何刚度矩阵KG。
步骤203，求解屈曲分析的控制方程(K-λKG)ai＝0，得到等几何屈曲设计模型的屈曲载荷与屈曲模态，其中，λ表示屈曲载荷，ai表示屈曲模态。
步骤300，根据屈曲分析控制方程(K-λKG)ai＝0，推导变刚度复合材料板壳结构屈曲载荷的全解析灵敏度，包括以下子步骤：
步骤301，计算刚度矩阵K对设计变量的导数，根据求导链式法则，转换为材料的偏轴刚度矩阵对设计变量的导数。
步骤302，计算几何刚度矩阵KG对设计变量的导数，同样根据求导链式法则，转换为变刚度复合材料板壳结构的应力分布矩阵σ对设计变量的导数。
其中，运用伴随法求解变刚度复合材料板壳结构的平衡方程。
步骤303，将步骤300所述的屈曲分析控制方程对设计变量求导，将步骤301中刚度矩阵K对设计变量的导数与步骤302中几何刚度阵KG对设计变量的导数代入，得到变刚度复合材料板壳结构屈曲载荷的全解析灵敏度，并校验其正确性。
步骤304，根据曲率函数约束纤维变化，防止纤维路径出现大幅度折角，推导并计算曲率函数对设计变量的全解析灵敏度。
步骤400，采用梯度类局部优化算法，利用步骤100和200提供的等几何屈曲设计模型以及步骤300所提供的全解析灵敏度信息进行曲线纤维路径优化，以屈曲载荷值作为优化目标，考虑实际工艺制造约束条件，得到具有更高承载效率的变刚度复合材料板壳结构的最优设计。所述的梯度类局部优化算法包括：移动渐近线法、最速下降法、可行方向法、单纯形法、序列线性规划法、序列二次规划法等。
本发明的有益效果为：本发明提供的变刚度复合材料板壳结构高效优化方法，针对现有的变刚度板壳结构传统的纤维路径优化方法效率低下的缺点，利用了等几何方法快速精确的分析优势，以少量的控制点实现屈曲性能的精确分析，更重要的是，全解析灵敏度可以在保证灵敏度信息精度的同时，避免了差分灵敏度计算的高耗时问题，从而大幅提高变刚度复合材料板壳结构的纤维铺放路径优化效率、降低产品研发周期。本发明有望成为将来我国运载火箭设计等航空航天领域中变刚度复合材料板壳结构的主要优化方法之一。
附图说明
图1为本发明实施例提供的变刚度复合材料板壳结构高效优化方法的实现流程图；
图2为本发明实施例提供的变刚度复合材料方形板结构示意图；
图3为本发明实施例提供的变刚度复合材料板壳结构纤维路径铺放示意图；
图4为本发明实施例提供的等几何方法与有限元方法计算时间对比图；
图5为本发明实施例提供的变刚度复合材料板壳结构前5阶屈曲载荷与屈曲模态；
图6为本发明实施例提供的有限元差分方法、等几何差分方法以及等几何解析灵敏度方法的优化迭代曲线示意图。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚，下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。
图1为本发明实施例提供的变刚度复合材料板壳结构高效优化方法的实现流程图。如图1所示，本发明实施例提供的变刚度复合材料板壳结构高效优化方法包括：
步骤100，变刚度复合材料方形板边长a＝254mm，如图2所示，共铺设20层，单层厚度0.15mm，弹性常数E1＝181GPa,E2＝10.270GPa,G12＝G13＝7.170GPa,G23＝3.780GPa,υ12＝0.28，承受大小为1的均布轴压载荷和剪切载荷。优化初始铺层信息为[<60°|15°>/<-60°|-15°>]5s。建立相应的变刚度复合材料板壳结构等几何屈曲设计模型，包括以下子步骤：
步骤101，以线性梯度函数铺层形式为例，给定纤维路径之后沿y轴平移得到的纤维铺放形式，其它层则将第一层关于x轴对称得到，如图3所示。
其中，θ(x)表示变刚度复合材料板壳结构中任意位置的纤维路径角度，x位纤维路径函数的自变量即结构位置信息，a表示结构的尺寸参数，T0和T1分别表示纤维路径在板中心与边缘的角度，如图1中所示。
步骤102，建立控制点角度与全场纤维路径的映射关系，进行收敛性分析，得到如图4所示的对比曲线。相比有限元方法，等几何方法能以较少的计算时间即可保证屈曲分析的精确性，分析效率。
步骤103，选定设计变量为每层的T0和T1，建立变刚度复合材料板壳结构的等几何屈曲设计模型，得到等几何分析所需的输入数据。；
步骤200，根据得到的等几何屈曲设计模型，进行变刚度复合材料板壳结构的等几何屈曲分析，得到屈曲载荷与屈曲模态。
步骤201，计算等几何单元积分点上纤维角度对刚度的影响，并累加成单元刚度矩阵，根据有限元基本理论，将变刚度复合材料板壳结构的单元刚度矩阵组装成等几何屈曲设计模型的刚度矩阵K；所述的刚度矩阵K由n个单元刚度矩阵kn组成。
步骤202，对等几何屈曲设计模型进行静力分析，得到应力分布，并以此求得模型的几何刚度矩阵KG；所述的几何刚度矩阵KG由n个单元的几何刚度矩阵组成。
步骤203，根据屈曲分析的控制方程(K-λKG)ai＝0，得到等几何屈曲设计模型的屈曲载荷与屈曲模态，如图5所示；其中，λ表示屈曲载荷，K表示结构的总体刚度矩阵，KG表示结构的几何刚度矩阵，ai表示屈曲模态。
步骤300，推导变刚度复合材料板壳结构屈曲响应的解析灵敏度，具体包括以下子步骤：
步骤301，计算刚度矩阵K对设计变量的导数，根据求导链式法则，问题归结于材料的偏轴刚度矩阵对设计变量的导数。
其中，kn表示单元刚度矩阵，Ti表示i个设计变量，numgas表示单元数，B表示应变矩阵，表示材料的偏轴刚度矩阵，w1表示等几何控制点垂直方向权函数，w2表示等几何控制点水平方向权函数，J表示雅可比矩阵。
步骤302，计算几何刚度矩阵KG对设计变量的导数，同样根据求导链式法则，问题归结于变刚度复合材料板壳结构的应力分布矩阵σ对设计变量的导数。
其中kgn表示单元的几何刚度矩阵，G表示几何刚度应变矩阵，[σ]表示应力分布矩阵，h表示变刚度复合材料层合板壳的厚度。
其中，需要采用伴随法求解结构的平衡方程。
其中，u表示结构的受力作用下产生的位移。
步骤303，推导并计算结构屈曲响应对设计变量的解析灵敏度，根据下式可得到全解析灵敏度信息：
其中，d表示结构位移向量。
步骤304，设置纤维的曲率函数约束纤维变化，曲率函数κ不超过0.1，由下述公式推导并计算曲率函数对设计变量的全解析灵敏度；
其中，κ表示曲线纤维路径的曲率函数，θ表示结构任意位置的曲线纤维路径的角度；θx表示θ对x的一阶导数；
步骤400，采用梯度类算法，利用步骤100和200提供的等几何屈曲设计模型以及步骤300所提供的灵敏度信息进行曲线纤维优化，采用移动渐进算法，以一阶屈曲载荷值作为优化目标，考虑曲率约束条件，优化迭代曲线如图6所示。相比差分方法计算灵敏度，全解析灵敏度方法不仅可以显著提高计算效率(只需14步即可获得最终优化结果，而差分法将近50步曲线仍然振荡)，单次分析时长全解析灵敏法也只是差分的1/10～1/5。最终优化结果解析灵敏度法提高到了38.66kN，而差分法只能达到34.54kN。相比传统有限元分析，等几何分析计算效率和计算精度都有极大提升，不仅为后续优化节约了大量的计算时间，优化后结构的承载效率也得到大幅提高。本发明有望成为未来我国运载火箭、导弹设计等航空航天领域中变刚度复合材料结构优化设计的方法之一。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (5)

1.一种变刚度复合材料板壳结构高效优化方法，其特征在于以下步骤：
步骤100，建立变刚度复合材料板壳结构的等几何屈曲设计模型，包括以下子步骤：
步骤101，对曲线纤维铺层，基于等几何方法实现纤维路径的参数化；
步骤102，建立控制点角度与全场纤维路径的映射关系，并进行收敛性分析；
步骤103，选定设计变量，建立变刚度复合材料板壳结构的等几何屈曲设计模型，得到等几何分析所需的输入数据；所述的设计变量为描述纤维路径的函数及描述纤维路径的角度；
步骤200，根据得到的等几何屈曲设计模型，进行变刚度复合材料板壳结构的等几何屈曲分析，得到屈曲载荷与屈曲模态，包括以下子步骤：
步骤201，根据有限元基本理论，将变刚度复合材料板壳结构的单元刚度矩阵组装成等几何屈曲设计模型的刚度矩阵K，其中，刚度矩阵K由n个单元刚度矩阵kn组成；
步骤202，对等几何屈曲设计模型进行静力分析，得到应力分布，并以此求得等几何屈曲设计的几何刚度矩阵KG，其中，几何刚度矩阵KG由n个单元的几何刚度矩阵kgn组成；
步骤203，求解屈曲分析的控制方程(K-λKG)ai＝0，得到等几何屈曲设计模型的屈曲载荷与屈曲模态，其中，λ表示屈曲载荷，ai表示屈曲模态；
步骤300，根据屈曲分析控制方程(K-λKG)ai＝0，推导变刚度复合材料板壳结构屈曲载荷的全解析灵敏度，包括以下子步骤：
步骤301，计算刚度矩阵K对设计变量的导数，根据求导链式法则，转换为材料的偏轴刚度矩阵对设计变量的导数；
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>k</mi> <mi>n</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mi>u</mi> <mi>m</mi> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </munderover> <msup> <mi>B</mi> <mi>T</mi> </msup> <mfrac> <mrow> <mo>&amp;part;</mo> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>Bw</mi> <mn>1</mn> </msub> <msub> <mi>w</mi> <mn>2</mn> </msub> <mo>|</mo> <mi>J</mi> <mo>|</mo> </mrow>
其中，kn表示单元刚度矩阵，Ti表示i个设计变量，numgas表示单元数，B表示应变矩阵，表示材料的偏轴刚度矩阵，w1表示等几何控制点垂直方向权函数，w2表示等几何控制点水平方向权函数，J表示雅可比矩阵；
步骤302，计算几何刚度矩阵KG对设计变量的导数，同样根据求导链式法则，转换为变刚度复合材料板壳结构的应力分布矩阵对设计变量的导数；
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>k</mi> <mrow> <mi>g</mi> <mi>n</mi> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mi>u</mi> <mi>m</mi> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </munderover> <msup> <mi>G</mi> <mi>T</mi> </msup> <mfrac> <mrow> <mo>&amp;part;</mo> <mo>&amp;lsqb;</mo> <mi>&amp;sigma;</mi> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>Ghw</mi> <mn>1</mn> </msub> <msub> <mi>w</mi> <mn>2</mn> </msub> <mo>|</mo> <mi>J</mi> <mo>|</mo> </mrow>
其中，kgn表示单元的几何刚度矩阵，G表示几何刚度应变矩阵，[σ]表示应力分布矩阵，h表示变刚度复合材料层合板壳的厚度；
运用伴随法求解变刚度复合材料板壳结构的平衡方程；
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>K</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mi>u</mi> <mo>+</mo> <mi>K</mi> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>u</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mrow>
其中，u表示结构的受力作用下产生的位移；
步骤303，将步骤300所述的屈曲分析控制方程对设计变量求导，将步骤301中刚度矩阵K对设计变量的导数与步骤302中几何刚度阵KG对设计变量的导数代入，根据下式得到变刚度复合材料板壳结构屈曲载荷的全解析灵敏度信息，并校验其正确性；
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;lambda;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>d</mi> <mi>T</mi> </msup> <msub> <mi>K</mi> <mi>G</mi> </msub> <mi>d</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msup> <mi>d</mi> <mi>T</mi> </msup> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>K</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mi>d</mi> <mo>-</mo> <msup> <mi>&amp;lambda;d</mi> <mi>T</mi> </msup> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>K</mi> <mi>G</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mi>d</mi> <mo>&amp;rsqb;</mo> </mrow>
其中，d表示结构位移向量；
步骤304，根据曲率函数约束纤维变化，由下述公式推导并计算曲率函数对设计变量的全解析灵敏度；
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;kappa;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>tan</mi> <mn>2</mn> </msup> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>x</mi> <mo>,</mo> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;theta;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>x</mi> <mo>,</mo> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow>
其中，κ表示曲线纤维路径的曲率函数，θ表示结构任意位置的曲线纤维路径的角度；θx表示θ对x的一阶导数；
步骤400，采用梯度类局部优化算法，根据等几何屈曲设计模型以及全解析灵敏度信息进行曲线纤维路径优化，以屈曲载荷值作为优化目标，考虑实际工艺制造约束条件，得到变刚度复合材料板壳结构的最优设计。
2.根据权利要求1所述的一种变刚度复合材料板壳结构高效优化方法，其特征在于，步骤101中所述的曲线纤维铺层形式包括：线性梯度函数、高阶曲面等高线函数、高阶Bézier曲线函数、拉格朗日多项式函数。
3.根据权利要求1或2所述的一种变刚度复合材料板壳结构高效优化方法，其特征在于，步骤304中所述的曲率函数不超过0.1。
4.根据权利要求1或2所述的一种变刚度复合材料板壳结构高效优化方法，其特征在于，步骤400中所述的梯度类局部优化算法包括：移动渐近线法、最速下降法、可行方向法、单纯形法、序列线性规划法、序列二次规划法。
5.根据权利要求3所述的一种变刚度复合材料板壳结构高效优化方法，其特征在于，步骤400中所述的梯度类局部优化算法包括：移动渐近线法、最速下降法、可行方向法、单纯形法、序列线性规划法、序列二次规划法。
CN201710492492.4A 2017-06-26 2017-06-26 一种变刚度复合材料板壳结构高效优化方法 Active CN107220461B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710492492.4A CN107220461B (zh) 2017-06-26 2017-06-26 一种变刚度复合材料板壳结构高效优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710492492.4A CN107220461B (zh) 2017-06-26 2017-06-26 一种变刚度复合材料板壳结构高效优化方法

Publications (2)

Publication Number Publication Date
CN107220461A true CN107220461A (zh) 2017-09-29
CN107220461B CN107220461B (zh) 2019-10-11

Family

ID=59951230

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710492492.4A Active CN107220461B (zh) 2017-06-26 2017-06-26 一种变刚度复合材料板壳结构高效优化方法

Country Status (1)

Country Link
CN (1) CN107220461B (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107729659A (zh) * 2017-10-19 2018-02-23 吴尘瑾 层合板铺层处理方法及装置
CN108009357A (zh) * 2017-12-01 2018-05-08 湖南大学 一种复合材料板形变量的分析方法
CN108363828A (zh) * 2017-12-28 2018-08-03 中国空间技术研究院 一种变刚度复合材料的建模方法
CN109766613A (zh) * 2018-12-29 2019-05-17 西安交通大学 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法
CN110083900A (zh) * 2019-04-17 2019-08-02 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN110955941A (zh) * 2019-11-29 2020-04-03 华中科技大学 基于向量场的复合材料结构优化设计方法及设备
CN111027151A (zh) * 2019-12-02 2020-04-17 大连理工大学 一种复合材料异形壳的纤维路径与几何形状一体化设计方法
CN111247080A (zh) * 2017-10-04 2020-06-05 蒂森克虏伯座椅电梯有限公司 规划平台升降机的方法
CN111259549A (zh) * 2020-01-16 2020-06-09 中南大学 一种任意边界条件下板组合结构屈曲分析的半解析方法
CN111444579A (zh) * 2020-03-11 2020-07-24 华中科技大学 一种考虑可制造性的复合材料结构优化设计方法
WO2020211012A1 (zh) * 2019-04-17 2020-10-22 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN111898295A (zh) * 2020-07-10 2020-11-06 西北工业大学 一种变刚度复合材料层合板的有限元建模方法
CN112182923A (zh) * 2020-09-08 2021-01-05 武汉大学 含孔复合材料结构件整体化成形方法
CN112818470A (zh) * 2021-02-09 2021-05-18 大连理工大学 一种基结构的优化方法、装置、计算机设备和存储介质
CN112861304A (zh) * 2020-12-04 2021-05-28 中国航空工业集团公司成都飞机设计研究所 一种面向压缩稳定性的变刚度复合材料结构设计方法
CN112989648A (zh) * 2021-02-04 2021-06-18 西安理工大学 一种协同拓扑构型与纤维路径的柔性机构优化设计方法
CN113434921A (zh) * 2021-07-05 2021-09-24 西安交通大学 一种考虑介纳观尺度效应的结构等几何拓扑优化方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930079A (zh) * 2012-10-08 2013-02-13 西北工业大学 一种分析复合材料层合板层间损伤的方法
CN104484531A (zh) * 2014-12-18 2015-04-01 大连理工大学 考虑多源不确定性的加筋板壳结构可靠度优化方法
CN104866673A (zh) * 2015-05-28 2015-08-26 大连理工大学 一种轴压加筋柱壳的开口补强方法
CN105224750A (zh) * 2015-10-10 2016-01-06 北京工业大学 一种基于响应面的新型空间可展铰链中单簧片结构优化方法
CN106156449A (zh) * 2016-08-31 2016-11-23 中航沈飞民用飞机有限责任公司 一种复合材料机翼壁板优化设计方法
CN106384384A (zh) * 2016-09-18 2017-02-08 上海理工大学 一种三维产品模型的形状优化算法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930079A (zh) * 2012-10-08 2013-02-13 西北工业大学 一种分析复合材料层合板层间损伤的方法
CN104484531A (zh) * 2014-12-18 2015-04-01 大连理工大学 考虑多源不确定性的加筋板壳结构可靠度优化方法
CN104866673A (zh) * 2015-05-28 2015-08-26 大连理工大学 一种轴压加筋柱壳的开口补强方法
CN105224750A (zh) * 2015-10-10 2016-01-06 北京工业大学 一种基于响应面的新型空间可展铰链中单簧片结构优化方法
CN106156449A (zh) * 2016-08-31 2016-11-23 中航沈飞民用飞机有限责任公司 一种复合材料机翼壁板优化设计方法
CN106384384A (zh) * 2016-09-18 2017-02-08 上海理工大学 一种三维产品模型的形状优化算法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PENG HAO等: "Isogeometric buckling analysis of composite variable-stiffness panels", 《COMPOSITE STRUCTURES》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111247080A (zh) * 2017-10-04 2020-06-05 蒂森克虏伯座椅电梯有限公司 规划平台升降机的方法
US11597634B2 (en) 2017-10-04 2023-03-07 T.K. Home Solutions B.V. Method of planning platform lift
CN107729659A (zh) * 2017-10-19 2018-02-23 吴尘瑾 层合板铺层处理方法及装置
CN108009357B (zh) * 2017-12-01 2021-01-29 湖南大学 一种复合材料板形变量的分析方法
CN108009357A (zh) * 2017-12-01 2018-05-08 湖南大学 一种复合材料板形变量的分析方法
CN108363828A (zh) * 2017-12-28 2018-08-03 中国空间技术研究院 一种变刚度复合材料的建模方法
CN108363828B (zh) * 2017-12-28 2021-10-08 中国空间技术研究院 一种变刚度复合材料的建模方法
CN109766613B (zh) * 2018-12-29 2020-10-27 西安交通大学 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法
CN109766613A (zh) * 2018-12-29 2019-05-17 西安交通大学 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法
WO2020211012A1 (zh) * 2019-04-17 2020-10-22 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN110083900A (zh) * 2019-04-17 2019-08-02 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN110955941A (zh) * 2019-11-29 2020-04-03 华中科技大学 基于向量场的复合材料结构优化设计方法及设备
CN110955941B (zh) * 2019-11-29 2022-03-18 华中科技大学 基于向量场的复合材料结构优化设计方法及设备
CN111027151A (zh) * 2019-12-02 2020-04-17 大连理工大学 一种复合材料异形壳的纤维路径与几何形状一体化设计方法
CN111027151B (zh) * 2019-12-02 2022-09-20 大连理工大学 一种复合材料异形壳的纤维路径与几何形状一体化设计方法
CN111259549A (zh) * 2020-01-16 2020-06-09 中南大学 一种任意边界条件下板组合结构屈曲分析的半解析方法
CN111259549B (zh) * 2020-01-16 2022-02-15 中南大学 一种任意边界条件下板组合结构屈曲分析的半解析方法
CN111444579B (zh) * 2020-03-11 2022-04-12 华中科技大学 一种考虑可制造性的复合材料结构优化设计方法
CN111444579A (zh) * 2020-03-11 2020-07-24 华中科技大学 一种考虑可制造性的复合材料结构优化设计方法
CN111898295A (zh) * 2020-07-10 2020-11-06 西北工业大学 一种变刚度复合材料层合板的有限元建模方法
CN112182923A (zh) * 2020-09-08 2021-01-05 武汉大学 含孔复合材料结构件整体化成形方法
CN112182923B (zh) * 2020-09-08 2023-06-02 武汉大学 含孔复合材料结构件整体化成形方法
CN112861304A (zh) * 2020-12-04 2021-05-28 中国航空工业集团公司成都飞机设计研究所 一种面向压缩稳定性的变刚度复合材料结构设计方法
CN112861304B (zh) * 2020-12-04 2022-06-28 中国航空工业集团公司成都飞机设计研究所 一种面向压缩稳定性的变刚度复合材料结构设计方法
CN112989648A (zh) * 2021-02-04 2021-06-18 西安理工大学 一种协同拓扑构型与纤维路径的柔性机构优化设计方法
CN112989648B (zh) * 2021-02-04 2022-11-25 西安理工大学 一种协同拓扑构型与纤维路径的柔性机构优化设计方法
CN112818470A (zh) * 2021-02-09 2021-05-18 大连理工大学 一种基结构的优化方法、装置、计算机设备和存储介质
CN112818470B (zh) * 2021-02-09 2023-02-03 大连理工大学 一种基结构的优化方法、装置、计算机设备和存储介质
CN113434921A (zh) * 2021-07-05 2021-09-24 西安交通大学 一种考虑介纳观尺度效应的结构等几何拓扑优化方法

Also Published As

Publication number Publication date
CN107220461B (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN107220461B (zh) 一种变刚度复合材料板壳结构高效优化方法
Hao et al. Design of manufacturable fiber path for variable-stiffness panels based on lamination parameters
Caliri Jr et al. A review on plate and shell theories for laminated and sandwich structures highlighting the Finite Element Method
JP7058902B2 (ja) ハイブリッド繊維複合材料の板巻きシェル構造に対する高速協調最適化方法
CN104866673B (zh) 一种轴压加筋柱壳的开口补强方法
Saka et al. Mathematical and metaheuristic applications in design optimization of steel frame structures: an extensive review
Xu et al. Robust topology optimization for multiple fiber-reinforced plastic (FRP) composites under loading uncertainties
CN105512385B (zh) 采用序列幂函数插值方法实现结构多材料拓扑优化的方法
CN107590325B (zh) 一种基于Shepard插值的纤维增强复合材料结构优化方法
Teng et al. An adaptively refined XFEM with virtual node polygonal elements for dynamic crack problems
CN106709215B (zh) 一种基于级数展开的连续体结构非概率可靠性拓扑优化方法
Zhou et al. Two-scale buckling topology optimization for grid-stiffened cylindrical shells
Raju et al. Optimal postbuckling design of variable angle tow composite plates
CN110955941B (zh) 基于向量场的复合材料结构优化设计方法及设备
Muc et al. Design of plates with curved fibre format
Liu et al. Second-order direct analysis of domelike structures consisting of tapered members with I-sections
CN109241562A (zh) 基于多尺度有限元方法的微结构材料弹性性能测定方法
Ma et al. Generative design of stiffened plates based on homogenization method
CN108763778A (zh) 一种实体材料和类桁架微结构材料一体化非概率可靠性拓扑优化方法
Duan et al. Discrete material selection and structural topology optimization of composite frames for maximum fundamental frequency with manufacturing constraints
Qin et al. Isogeometric vibration and buckling analyses of curvilinearly stiffened composite laminates
Chu et al. Explicit topology optimization of novel polyline-based core sandwich structures using surrogate-assisted evolutionary algorithm
Kussmaul et al. A novel computational framework for structural optimization with patched laminates
CN108491591A (zh) 一种高温环境下曲线加筋板有限元分析方法
CN103065015B (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