CN107526898B - 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 - Google Patents
变刚度复合材料板壳结构建模分析与可靠度优化设计方法 Download PDFInfo
- Publication number
- CN107526898B CN107526898B CN201710817832.6A CN201710817832A CN107526898B CN 107526898 B CN107526898 B CN 107526898B CN 201710817832 A CN201710817832 A CN 201710817832A CN 107526898 B CN107526898 B CN 107526898B
- Authority
- CN
- China
- Prior art keywords
- reliability
- variable
- shell structure
- analysis
- constraint
- 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
- 238000000034 method Methods 0.000 title claims abstract description 133
- 238000004458 analytical method Methods 0.000 title claims abstract description 96
- 239000002131 composite material Substances 0.000 title claims abstract description 74
- 238000013461 design Methods 0.000 title claims abstract description 72
- 238000005457 optimization Methods 0.000 title claims abstract description 61
- 239000000835 fiber Substances 0.000 claims abstract description 28
- 230000035945 sensitivity Effects 0.000 claims abstract description 24
- 238000004519 manufacturing process Methods 0.000 claims description 19
- 239000010410 layer Substances 0.000 claims description 17
- 238000009826 distribution Methods 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 10
- 230000000739 chaotic effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 5
- 241000375392 Tana Species 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims 1
- 230000004044 response Effects 0.000 abstract description 7
- 210000001503 joint Anatomy 0.000 abstract description 3
- 238000012827 research and development Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 7
- 235000015842 Hesperis Nutrition 0.000 description 3
- 235000012633 Iberis amara Nutrition 0.000 description 3
- 238000005452 bending Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 239000011157 advanced composite material Substances 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003733 fiber-reinforced composite Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000009864 tensile test Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
-
- 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
- B29C70/00—Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts
- B29C70/04—Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts comprising reinforcements only, e.g. self-reinforcing plastics
- B29C70/28—Shaping operations therefor
- B29C70/30—Shaping by lay-up, i.e. applying fibres, tape or broadsheet on a mould, former or core; Shaping by spray-up, i.e. spraying of fibres on a mould, former or core
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Chemical & Material Sciences (AREA)
- Composite Materials (AREA)
- Mechanical Engineering (AREA)
Abstract
本发明涉及结构可靠度优化领域,提供一种变刚度复合材料板壳结构建模分析与可靠度优化设计方法,包括:利用一次可靠性近似方法、非线性近似函数以及二次可靠性近似方法对变刚度复合材料板壳结构进行高效可靠度优化。利用非均匀有理B样条函数对变刚度复合材料板壳纤维铺设路径进行精确建模;对变刚度板壳结构进行等几何分析,包括:基于等几何方法对变刚度板壳结构进行线性屈曲分析,推导设计变量以及随机变量对结构响应的全解析灵敏度。本发明能够实现变刚度复合材料板壳结构的建模、分析与可靠度优化的无缝对接,显著提高其可靠度优化效率及准确性,大幅缩短研发周期。
Description
技术领域
本发明属于结构力学分析方法与结构可靠度优化设计技术领域,涉及一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法。
背景技术
航空航天结构追求性能第一的特点,使其成为先进复合材料技术的率先实验和转化的战场,提高复合材料用量对促进航空航天装备的轻量化和高性能化起到了至关重要的作用。受到传统制造工艺约束,目前纤维增强复合材料层合板多采用平行顺直纤维铺放成型,且同层的纤维角度固定,工程中常采用0°、90°和±45°的铺层方向,使得设计空间受限。相比之下,变刚度设计中同层的纤维铺放路径是连续曲线所示,其纤维路径可根据实际承载需求,提高或降低局部区域的强度和刚度,改善铺层平面内的应力分布,从而有效提高整体结构性能(如屈曲荷载和破坏强度),已应用于国外飞行器机身结构、异型截面机翼结构等设计中。NASA认为结构弹性剪裁是实现未来亚音速飞行器低成本、轻量化的关键技术之一,即设计者可利用具有高强度和稳定性约束的材料,使承力结构内的刚度特性按照需求分布。伴随我国重型运载火箭和大飞机的大直径、薄壁、超轻等越发苛刻的设计需求,以及承力结构多功能性融合的趋势所导致的开口增多,面内传力路径和应力分布状态变得极为复杂,且暴露的不确定性风险显著增大。显然,变刚度复合材料板壳这类极具承载潜力的结构对于具有复杂传力路径的舱段壳体具有极大的吸引力。
复合材料结构的承载性能受材料、制造公差、载荷等多种不确定性因素的影响,特别是变刚度结构大幅延拓了设计空间,制造工艺过程引起的纤维铺设角度空间分布位置等偏差成为新的不确定性因素,导致结构不确定性风险显著增大,不确定性优化设计需求越发迫切。为提高优化效率,等几何分析方法已应用于形状优化设计,这是由于设计和分析模型采用相同的几何描述,可以很容易地直接对灵敏度进行分析,摒弃了传统有限元灵敏度计算中网格近似引起的精度问题。可以预见,相对传统有限元方法,基于等几何分析的变刚度板壳结构优化设计可在满足更高精度要求的同时大幅缩短计算时间。
现有的变刚度复合材料板壳结构传统可靠度优化方法计算量巨大且过度依赖有限元分析,优化过程中需要多次调用精细有限元分析,即使采用采用等几何分析技术,其优化效率也十分低下。而且,受限于计算资源,工程中往往采用精度较差的一次可靠度近似方法进行可靠度优化,但对于可靠度要求非常苛刻的航空航天领域,高效率、高精度的可靠性分析是十分重要且需求迫切的。该领域虽然已开展了大量工作,但仍未给出一种效率高、精度好的针对变刚度复合材料板壳结构的精确建模与可靠度优化一体化方法。
发明内容
本发明主要解决现有技术的变刚度复合材料板壳结构建模分析精度差,可靠度优化设计的计算成本巨大、精度较差、优化效率低下的技术问题,提出一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法,实现建模、分析与可靠度优化过程的无缝对接,以达到提高变刚度复合材料板壳结构建模精确性、提升可靠度优化设计的计算效率、计算精度、缩短产品设计周期的目的。
为了达到上述目的,本发明的技术方案为:
一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法,包括以下步骤:
步骤100:利用高效梯度类方法进行变刚度复合材料板壳结构优化,包括以下子步骤:
步骤101:利用等几何精确建模分析方法获得变刚度复合材料板壳结构的屈曲载荷,质量以及全解析一阶导数、二阶导数,计算当前设计点处的制造约束与可靠度约束。基于可靠度分析方法,将变刚度复合材料板壳结构屈曲载荷的确定性约束P≥Pt转换为可靠度约束其中,P为结构的屈曲载荷,Pt为预设的结构屈曲载荷,pf为P<Pt的失效概率,为预设的失效概率。
步骤102:根据步骤101的制造约束、可靠度约束以及等几何精确建模分析方法,考虑不确定因素,建立变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计框架。将单层纤维铺设角度、层厚和铺层层数之一或其组合作为设计变量,将材料性能偏差和铺层角偏差作为随机变量,考虑制造约束与屈曲载荷的可靠度约束,采用高效梯度类算法对结构重量和/或结构制造成本进行优化。所述的高效梯度类算法为内点法、有效集方法、移动渐近线法、BFGS法、近似序列规划法。
步骤200:利用一次可靠度近似方法对变刚度复合材料板壳结构进行一次可靠度分析,包括以下子步骤:
步骤201:初始化内层迭代步数k=0,将步骤102中所述的随机变量转换到标准正态空间,基于等几何精确建模分析方法,计算迭代点Uk的结构屈曲载荷P和全解析灵敏度;其中,标准正态空间为U空间。
步骤202:利用一次可靠度近似方法中的增强混沌控制法(ECC)求解下一步迭代点Uk+1,并基于等几何精确建模分析方法,计算该迭代点处的结构屈曲载荷P和全解析灵敏度。
步骤203:根据Uk、Uk+1提供的结构屈曲载荷值和全解析灵敏度构建两点非线性近似函数,通过ECC求解近两点非线性近似函数的近似设计验算点UAMPP,并判断是否满足收敛条件:若满足收敛条件,计算UAMPP处的真实灵敏度;若不满足收敛条件,k=k+1且返回步骤202。所述两点非线性近似函数采用如下方法构造,包括TANA方法、TANA1方法或TANA2方法。所述的收敛条件为Uk+1与UAMPP之间差值的欧几里得范数|Uk+1-UAMPP|≤0.001。
步骤204:基于等几何精确建模分析方法,计算UAMPP处的约束函数值g(UAMPP),并判断是否足够靠近约束:若足够靠近约束,则进行步骤300中的二次可靠度分析;若不够靠近约束,返回步骤101。所述的足够靠近约束为其中,B(Xinitial)为初始设计Xinitial处的结构屈曲载荷。
步骤300:利用步骤204中一次可靠度信息进行二次可靠度分析,一次可靠度信息包括UAMPP及该点处的约束函数值,包括以下子步骤:
步骤301:基于等几何精确建模分析方法,求解UAMPP处的结构真实灵敏度以及二阶导数;
步骤302:利用二次可靠度近似方法对变刚度复合材料板壳结构进行可靠度分析,并求解失效概率pf:若迭代停止,得到考虑可靠度约束以及制造约束的最优设计;若返回步骤101。所述的二次可靠度近似方法包括:Breitung方法、Tvedt方法和基于广义卡方分布的二次可靠度方法。
针对现有的变刚度复合材料板壳结构,基于有限元建模分析存在建模不精确,求解效率低且不能提供高阶灵敏度的缺点,传统可靠度优化方法存在效率低下、精度差的问题,本发明综合利用等几何精确建模分析方法、一次可靠度近似方法、二次可靠度近似方法和函数近似方法,在不同的优化阶段引入不同的分析模型。由于近似函数与一次可靠度近似方法的计算效率高以及等几何精确建模分析能够消除设计模型与计算模型之间的误差并提供高精度结构响应和全解析灵敏度,第一阶段优化中利用基于等几何精确建模分析的一次可靠度近似方法。利用等几何精确建模分析方法提供的真实结构响应以及灵敏度,构建两点非线性近似函数,得到近似设计验算点;并通过等几何分析校验该近似验算点响应值,从而大幅提高变刚度复合材料板壳结构的可靠度优化效率、降低计算成本。当该点的约束函数值足够靠近约束函数则开展第二阶段的二次可靠度求解工作。第二阶段优化中利用等几何的全解析一、二阶导数,无需过多的计算量即可计算更加精确的结构失效概率,保证可靠度分析的准确性。本发明有望成为我国运载火箭设计等航空航天领域中变刚度复合材料板壳结构的可靠度主要优化方法之一。
本发明能够实现变刚度复合材料板壳结构的建模、分析与可靠度优化的无缝对接,显著提高其可靠度优化效率及准确性,大幅缩短研发周期。
附图说明
图1为本发明实施例提供的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法的实现流程图。
图2为变刚度复合材料板壳结构纤维铺放路径示意图。
图3为本发明实施例提供的为变刚度复合材料板壳结构铺层示意图。
图4为本发明实施例提供的变刚度复合材料板壳结构的示意图。
图5(a)为本发明实施例提供的变刚度复合材料板壳初始设计的屈曲模态图。
图5(b)为本发明实施例提供的变刚度复合材料板壳最优化设计的屈曲模态图。
图6为本发明实施例提供的变刚度复合材料板壳结构优化迭代曲线。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚,下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。
图1为本发明实施例提供的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法的实现流程图。如图1所示,本发明实施例提供的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法包括:
步骤100:利用高效梯度类方法进行变刚度复合材料板壳结构优化,包括以下子步骤:
步骤101:利用等几何精确建模分析方法获得变刚度复合材料板壳结构的屈曲载荷,质量以及全解析一阶导数、二阶导数,计算当前设计点处的制造约束与可靠度约束。基于可靠度分析方法,将变刚度复合材料板壳结构屈曲载荷的确定性约束P≥Pt转换为可靠度约束所述的P为结构的屈曲载荷,Pt为预设的结构屈曲载荷,pf为P<Pt的失效概率,为预设的失效概率。
其中,所述的等几何精确建模分析方法包括:通过B样条曲线描述纤维铺设路径,利用控制点对纤维路径进行参数化;根据得到的参数化模型,通过屈曲分析的控制方程(K-λKG)ai=0可以得到等几何屈曲设计模型的屈曲载荷λ与屈曲模态ai。其中Ke为单元刚度阵,K为整体刚度阵,KG为几何刚度矩阵。根据屈曲分析的控制方程推导各变量对变刚度复合材料板壳结构屈曲载荷的全解析一阶导数以及二阶导数。
本发明中所考虑的制造约束为变刚度复合材料板壳结构中所铺设纤维的最大曲率小于预设值。以线性梯度函数铺层形式为例,给定纤维铺设路径之后沿y轴平移得到的纤维铺放形式,如图2所示。
其中,θ(x')表示变刚度复合材料板壳结构中任意位置的纤维路径角度,α为曲线铺放的纤维参考坐标x'与x轴之间的夹角;T0为曲线铺放的纤维在原点(x0,y0)处和x'方向的夹角;T1为曲线铺放的纤维在距离参考坐标系x',y'原点距离为h处和x'方向的夹角,单层铺层角度表示方法为α<T0,T1>。变刚度复合材料板壳结构中纤维路径的曲率可以通过如下公式获得:
FIB(x,y)=fib(x,y)-fib(x0,y0)=0
其中,cur(x,y)为(x,y)处纤维的曲率,FIB(x,y)为纤维路径场函数的等高线,fib(x,y)为纤维路径的场函数,下标x,y表示对x,y的偏导数,y′为一个自定义参数,其表示纤维路径场函数的等高线在(x,y)处对x,y的偏导数比值的相反数,y”为x对y′的导数。
失效概率可以通过可靠度分析得到,每次优化迭代都在内层嵌套了可靠度分析。预设的失效概率可以根据结构的实际服役安全性要求来选择,例如航空航天结构的可靠度要求高,其失效概率一般选为0.001,可靠度=1.0-失效概率,或者可靠度+失效概率=1.0。在变刚度复合材料板壳结构可靠度优化领域,铺层维铺设角度和材料属性是影响加变刚度复合材料板壳结构轴压承载力的主要因素,通常通过评价变刚度复合材料板壳的轴压承载力来评价变刚度复合材料板壳结构的可靠度。所以本发明将变刚度复合材料板壳结构轴压下的屈曲载荷的可靠度作为可靠度约束条件(例如,轴压下屈曲载荷的可靠度小于某个值的概率不能低于一定百分比)。
步骤102:根据步骤101的曲率约束与可靠度约束以及等几何精确建模分析方法,考虑不确定因素,建立变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计框架。将单层纤维铺设角度,层厚和铺层层数之一或其组合作为设计变量,将材料性能偏差和铺层角偏差作为随机变量,考虑制造约束与屈曲载荷的可靠度约束,采用有效集方法对结构重量和/或结构制造成本进行优化。
具体的,本发明实施例提供的方法的优化列式可以表示为:
设计变量:d=[α,T0,T1,t,n]
随机变量:X=[T0,T1,E1,E2,G12,G13,G23,υ12]
目标函数:W
约束条件:
其中,d为设计变量,α为纤维场整体旋转角度,T0和T1分别表示纤维路径角度,t为单层层厚,n为铺层数量。X为随机变量,E1、E2、G12、G23、G23为弹性常数,υ12为泊松比。W为结构重量,P为结构屈曲载荷或压溃载荷,Pt为结构设计载荷,pf为失效概率,为预设的失效概率,cur(x,y)为点(x,y)处纤维的曲率,cm为制造约束所允许的最大曲率。
所述不确定性因素包括:材料性能偏差和铺层角偏差等。具体的,材料性能偏差通过拉伸试验并统计获得,铺层角偏差通过量角器测量并统计获得。不确定性因素可以根据统计数据给定不确定性因素的随机分布类型,如正态分布、对数正态分布、均匀分布、威布尔分布等,不同分布类型的不确定性因素将对约束条件以及优化目标产生不同影响。
其中,图3为本发明实施例提供的变刚度复合材料板壳结构的示意图。本发明提供的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法可以针对曲面变刚度复合材料板壳、加筋板壳、开口板壳等复合材料板壳结构开展可靠度优化。
步骤200:利用一次可靠度近似方法对变刚度复合材料板壳结构进行一次可靠度分析,包括以下子步骤:
步骤201:初始化内层迭代步数k=0,将步骤102中所述的随机变量转换到标准正态空间,基于等几何精确建模分析方法,计算Uk的结构屈曲载荷P和全解析灵敏度;其中,标准正态空间为U空间;
所述的将随机变量转换到标准正态空间的方法指当随机变量为正态分布时,通过其均值与方差直接转换到标准正态空间;当随机变量为非正态分布时,则采用等效正态分布近似法进行转换。
步骤202:利用一次可靠度近似方法中的增强混沌控制法(ECC)求解下一步迭代点Uk+1,并基于等几何精确建模分析方法,计算当该迭代点处的结构屈曲载荷P和全解析灵敏度;
所述一次可靠度近似方法是指对实际的约束函数g(U)在U*处进行一阶泰勒展开用来近似原来的复杂约束函数,其形式为其中并将所有随机变量转换到标准正态空间,用设计验算点到原点的距离表征可靠度指标β的可靠度分析方法。其中,可靠度指标对应标准正态变量概率密度函数φ的下分位点,即安全概率Ps=Φ(β)=1-Pf,其中Φ为标准正态变量累积分布函数。可靠度指标β可由如下公式计算。
其中,U*为设计验算点,为U*的第i个元素。为约束函数在设计验算点处的导数,为随机变量xi的标准差。求解设计验算点的方法包括改进均值法(AMV),稳定转换法(STM),增强混沌控制方法(ECC),高效变循环方法(EAL)。采用稳健高效的增强混沌控制方法公式如下:
ρ=0.2,其中,为Uk+1的方向向量,Gk,u为函数g(U)在Uk处的导数,f()为merit函数,为merit函数关于U的导数,F(Uk)为一过渡变量。c为merit函数中的一个系数,ρ,为两个常数。C为对合矩阵,γ为控制因子。ωk为Uk-1与Uk迭代方向的夹角,为(F(Uk)-Uk)与(Uk-Uk-1)之间的夹角。arccos()为反余弦函数。上标符号T表示转置。
步骤203:根据Uk,Uk+1提供的结构屈曲载荷值和全解析灵敏度构建两点非线性近似函数,通过ECC求解近两点非线性近似函数的近似设计验算点(UAMPP),并判断是否满足收敛条件。若满足收敛条件,计算UAMPP处的真实灵敏度;若不满足收敛条件,k=k+1且返回步骤202。所述的两点非线性近似函数可采用如下方法构造,包括TANA方法、TANA1方法或TANA2方法。所述的收敛条件为Uk+1与UAMPP之间差值的欧几里得范数|Uk+1-UAMPP|≤0.001。
所述的构建两点非线性近似函数可通过拟合两点的函数值以及梯度值,建立约束函数的局部近似函数。采用TANA方法构造两点非线性近似函数形式如下:
其中xi,k为设计变量Xk第k次迭代是的第i个元素。上式中的非线性指数r可以通过下式求得:
步骤204:基于等几何精确建模分析方法,计算UAMPP处的约束函数值g(UAMPP),并判断是否足够靠近约束。若足够靠近约束,则进行步骤300中的二次可靠度分析。若不够靠近约束,返回步骤101。其中所述的足够靠近约束为B(Xinitial)为初始设计Xinitial处结构的屈曲载荷。
步骤300:利用步骤204中一次可靠度信息包括UAMPP及该点处的约束函数值进行二次可靠度分析,包括以下子步骤
步骤301:基于等几何精确建模分析方法,求解UAMPP处的结构真实灵敏度以及二阶导数;
步骤302:利用二次可靠度近似方法对变刚度复合材料板壳结构进行可靠度分析,并求解失效概率pf。若迭代停止,得到考虑可靠度约束以及制造约束的最优设计;若返回步骤101。所述的二次可靠度近似方法包括:Breitung方法、Tvedt方法或基于广义卡方分布的SORM。
所述的二次可靠度分析方法是指根据约束函数在设计验算点处的函数值g(U*),一阶导数以及二阶导数对其进行二阶泰勒展开,用来近似原来的复杂约束函数,其形式为基于该二阶展开形式求解更加精确的当前设计点处的失效概率,其方法包括Breitung方法、Tvedt方法或基于广义卡方分布的SORM。采用Tvedt方法求得该点的失效概率:
Pf=A1+A2+A3
其中,A1,A2,A3为构成失效概率Pf的三项。Φ()为标准正态分布下的累积分布函数,φ()标准正态分布下的概率密度函数,ki为失效曲面概率验算点的主曲率。
图4为变刚度复合材料平板结构的示意图。参照图4,考虑一个几何尺寸长为609.6mm,高为711.2mm。初始铺层数为12层,材料采用碳纤维/环氧树脂基体复合材料AS4/9773。铺层的单层厚度为0.1905mm。初始设计铺层信息为{25<20,16>,-25<-20,-16>,-12<40,70>,12<-40,-70>,-12<40,70>,12<-40,-70>}s,其中下标s表示铺层形式为对称。结构重量为1.60kg,屈曲载荷为4.26kN,目标屈曲载荷为3.0kN。预设的失效概率弹性常数和泊松比被考虑为正态分布随机变量,变异系数为0.03。纤维路径角度T0和T1也被考虑为正态分布随机变量,标准差为0.05。
优化迭代曲线如图6所示。相比基于差分法的有限元方法及等几何分析方法,基于本发明所提出的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法不仅可以快速收敛到最优解,而且迭代过程更加平稳。此外,基于二次可靠度近似方法的可靠度优化能够得到更加准确的失效概率。就单次分析的效率而言,基于本发明提出的方法调用的等几何分析次数最少,仅为传统基于一次可靠度近似方法优化策略的1/5。本发明所提出的方法在保证失效概率计算准确的情况下将结构质量降低了17.5%,相比之下,传统方法的有限元+差分法只能降低15.6%且失效概率计算不够准确。初始设计和最优设计的屈曲模态图如图5(a)和图5(b)所示。相比传统基于有限元分析以及差分法的可靠度优化方法,本发明所提出的变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法在计算效率和计算精度都有极大提升,不仅为后续变刚度复合材料板壳结构的设计以及研发节约了大量的时间,优化后结构的承载效率也得到大幅提高。本发明有望成为未来我国运载火箭、导弹设计等航空航天领域中变刚度复合材料结构优化设计的方法之一。
本发明提供的一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法。针对现有的变刚度复合材料板壳结构,基于有限元建模分析存在建模不精确,求解效率低且不能提供高阶灵敏度的缺点,传统可靠度优化方法存在效率低下、精度差的问题。本发明综合利用等几何精确建模分析方法、一次可靠度近似方法、二次可靠度近似方法和函数近似方法,在不同的优化阶段引入不同的分析模型。由于近似函数与一次可靠度近似方法的计算效率高以及等几何精确建模分析能够消除设计模型与计算模型之间的误差并提供高精度结构响应和全解析灵敏度,第一阶段优化中利用一次可靠度近似方法。利用等几何精确建模分析方法提供的真实结构响应以及灵敏度,构建两点非线性近似函数,得到近似设计验算点;并通过等几何分析校验该近似验算点响应值,从而大幅提高变刚度复合材料板壳结构的可靠度优化效率、降低计算成本。当该点足够靠近约束函数则开展第二阶段的二次可靠度求解工作。第二阶段优化中利用等几何的全解析一、二阶导数,无需过多的计算量即可计算更加精确的结构失效概率,保证可靠度分析的准确性。本发明有望成为我国运载火箭设计等航空航天领域中变刚度复合材料板壳结构的可靠度主要优化方法之一。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (4)
1.变刚度复合材料板壳结构建模分析与可靠度优化设计方法,其特征在于以下步骤:
步骤100:利用高效梯度类方法进行变刚度复合材料板壳结构优化,包括以下子步骤:
步骤101:利用等几何精确建模分析方法获得变刚度复合材料板壳结构的屈曲载荷,质量以及全解析一阶导数、二阶导数,计算当前设计点处的制造约束与可靠度约束;基于可靠度分析方法,将变刚度复合材料板壳结构屈曲载荷的确定性约束P≥Pt转换为可靠度约束其中,P为结构的屈曲载荷,Pt为预设的结构屈曲载荷,pf为P<Pt的失效概率,为预设的失效概率;
步骤102:根据步骤101的制造约束、可靠度约束以及等几何精确建模分析方法,考虑不确定因素,建立变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计框架;将单层纤维铺设角度、层厚和铺层层数之一或其组合作为设计变量,将材料性能偏差和铺层角偏差作为随机变量,考虑制造约束与屈曲载荷的可靠度约束,采用高效梯度类算法对结构重量和/或结构制造成本进行优化;
步骤200:利用一次可靠度近似方法对变刚度复合材料板壳结构进行一次可靠度分析,包括以下子步骤:
步骤201:初始化内层迭代步数k=0,将步骤102中所述的随机变量转换到标准正态空间,基于等几何精确建模分析方法,计算迭代点Uk的结构屈曲载荷P和全解析灵敏度;其中,标准正态空间为U空间;
步骤202:利用一次可靠度近似方法中的增强混沌控制法(ECC)求解下一步迭代点Uk+1,并基于等几何精确建模分析方法,计算该迭代点处的结构屈曲载荷P和全解析灵敏度;所述的增强混沌控制法ECC的公式如下:
ρ=0.2,
其中,β为可靠度指标,为Uk+1的方向向量,Gk,u为函数g(U)在Uk处的导数,f()为merit函数,为merit函数关于U的导数,F(Uk)为一过渡变量;c为merit函数中的一个系数,ρ,为两个常数;C为对合矩阵,γk和γk-1分别为第k和k-1迭代步的控制因子;ωk为Uk-1与Uk迭代方向的夹角,为(F(Uk)-Uk)与(Uk-Uk-1)之间的夹角;arccos()为反余弦函数;上标符号T表示转置;
步骤203:根据Uk、Uk+1提供的结构屈曲载荷值和全解析灵敏度构建两点非线性近似函数,通过ECC求解近两点非线性近似函数的近似设计验算点UAMPP,并判断是否满足收敛条件:若满足收敛条件,计算UAMPP处的真实灵敏度;若不满足收敛条件,k=k+1且返回步骤202;所述的收敛条件为Uk+1与UAMPP之间差值的欧几里得范数|Uk+1-UAMPP|≤0.001;
步骤204:基于等几何精确建模分析方法,计算UAMPP处的约束函数值g(UAMPP),并判断是否足够靠近约束:若足够靠近约束,则进行步骤300中的二次可靠度分析;若不够靠近约束,返回步骤101;所述的足够靠近约束为其中,B(Xinitial)为初始设计Xinitial处结构的屈曲载荷;
步骤300:利用步骤204中一次可靠度信息进行二次可靠度分析,一次可靠度信息包括UAMPP及该点处的约束函数值,包括以下子步骤:
步骤301:基于等几何精确建模分析方法,求解UAMPP处的结构真实灵敏度以及二阶导数;
步骤302:利用二次可靠度近似方法对变刚度复合材料板壳结构进行可靠度分析,并求解失效概率pf:若迭代停止,得到考虑可靠度约束以及制造约束的最优设计;若返回步骤101。
2.根据权利要求1所述的变刚度复合材料板壳结构建模分析与可靠度优化设计方法,其特征在于,步骤102中所述的高效梯度类算法为内点法、有效集方法、移动渐近线法、BFGS法、近似序列规划法。
3.根据权利要求1所述的变刚度复合材料板壳结构建模分析与可靠度优化设计方法,其特征在于,步骤203中所述两点非线性近似函数采用如下方法构造,包括TANA方法、TANA1方法或TANA2方法。
4.根据权利要求1所述的变刚度复合材料板壳结构建模分析与可靠度优化设计方法,其特征在于,步骤302中所述的二次可靠度近似方法包括:Breitung方法、Tvedt方法和基于广义卡方分布的二次可靠度方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817832.6A CN107526898B (zh) | 2017-09-13 | 2017-09-13 | 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 |
US15/850,909 US10726168B2 (en) | 2017-09-13 | 2017-12-21 | Integration method for accurate modeling and analysis and reliability-based design optimization of variable stiffness composite plate and shell structures |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817832.6A CN107526898B (zh) | 2017-09-13 | 2017-09-13 | 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107526898A CN107526898A (zh) | 2017-12-29 |
CN107526898B true CN107526898B (zh) | 2019-12-27 |
Family
ID=60736587
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710817832.6A Active CN107526898B (zh) | 2017-09-13 | 2017-09-13 | 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10726168B2 (zh) |
CN (1) | CN107526898B (zh) |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763608B (zh) * | 2018-03-23 | 2022-05-13 | 江苏理工学院 | 一种基于发生函数法的复合材料层合板可靠性评估方法 |
CN109766604B (zh) * | 2018-12-27 | 2020-10-16 | 浙江大学 | 一种基于随机等几何分析的叶片高刚度设计方法 |
CN109829261A (zh) * | 2019-04-03 | 2019-05-31 | 大连理工大学 | 一种由变截面圆管组成的轻质点阵结构及其优化设计方法 |
CN110083889B (zh) * | 2019-04-10 | 2023-05-30 | 中国民航大学 | 考虑焊接空心球节点约束的圆钢管稳定承载力的诊断方法 |
CN110222355A (zh) * | 2019-04-11 | 2019-09-10 | 北京宇航系统工程研究所 | 一种考虑热力耦合的低温推进剂贮箱防晃板分析方法 |
CN110222356B (zh) * | 2019-04-12 | 2024-01-02 | 北京工业大学 | 综合考虑稳定性与振动特性的板/壳结构轻量化拓扑优化设计方法 |
CN110135051B (zh) * | 2019-05-11 | 2022-12-06 | 西安电子科技大学 | 基于解析机电耦合模型的网状天线结构稳健性设计方法 |
CN110162895B (zh) * | 2019-05-27 | 2023-04-07 | 大连海事大学 | 一种两阶段的高能效船型优化设计方法 |
CN110210134B (zh) * | 2019-06-04 | 2022-09-13 | 合肥工业大学 | 一种白车身结构优化设计方法及结构优化设计装置 |
CN110532584B (zh) * | 2019-06-04 | 2022-12-27 | 南京依维柯汽车有限公司 | 一种基于相对灵敏度的非承载式白车身轻量化方法 |
CN110362902B (zh) * | 2019-07-02 | 2023-01-03 | 宁波大学 | 一种基于区间逐维分析的单源动载荷识别方法 |
CN110457768B (zh) * | 2019-07-18 | 2022-12-13 | 东南大学 | 考虑工艺误差下基于可靠性的mems器件参数的配置方法 |
CN110633498A (zh) * | 2019-08-16 | 2019-12-31 | 江苏理工学院 | 一种应用于cfrp汽车顶盖结构材料功能一体化设计方法 |
CN110705100B (zh) * | 2019-09-30 | 2020-10-23 | 西北工业大学 | 对称复合材料层合板弯曲刚度铺层序优化设计方法 |
CN110990939B (zh) * | 2019-10-10 | 2024-04-19 | 西北工业大学 | 防冰腔结构可靠性设计方法 |
CN110781606B (zh) * | 2019-11-07 | 2023-03-24 | 西安科技大学 | 一种针对梁结构的多设计点非概率可靠性分析方法 |
CN111027151B (zh) * | 2019-12-02 | 2022-09-20 | 大连理工大学 | 一种复合材料异形壳的纤维路径与几何形状一体化设计方法 |
CN111144048B (zh) * | 2019-12-23 | 2023-04-07 | 北京机电工程研究所 | 一种透波、隐身、防热天线罩的设计优化方法 |
CN111159895B (zh) * | 2019-12-30 | 2023-03-24 | 西安科技大学 | 一种基于模糊理论的机械元动作装配精度可靠性分析方法 |
CN111274670B (zh) * | 2019-12-31 | 2023-08-08 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种钛合金型材组合壁板轴压承载能力试验设计方法 |
CN111241672A (zh) * | 2020-01-07 | 2020-06-05 | 浙江大学 | 一种混合基函数的天线模型载荷位移测量方法 |
CN111241732B (zh) * | 2020-01-07 | 2022-03-11 | 浙江大学 | 基于子结构自由度凝聚的天线模型位移快速测量方法 |
EP4107673A4 (en) * | 2020-02-17 | 2023-08-09 | Petroliam Nasional Berhad (Petronas) | EQUIPMENT FAILURE PROBABILITY CALCULATION AND LIFE ESTIMATING METHODS AND SYSTEMS |
CN111523275B (zh) * | 2020-03-27 | 2023-04-11 | 北京航空航天大学 | 基于模糊随机参数的裕度与不确定性量化结构可靠性评估方法 |
CN111382539B (zh) * | 2020-03-30 | 2023-04-07 | 合肥市太泽透平技术有限公司 | 一种基于通流计算的透平机械叶型优化方法 |
CN111539138B (zh) * | 2020-03-31 | 2024-03-26 | 中国空间技术研究院 | 基于阶跃函数的结构动力学峰值时域响应灵敏度求解方法 |
CN111723457B (zh) * | 2020-05-08 | 2022-06-17 | 华中科技大学 | 一种纤维曲线铺放变刚度结构优化设计的水平集方法 |
CN111597654B (zh) * | 2020-05-13 | 2023-04-28 | 北京航空航天大学 | 一种复杂截面形状的铝合金型材车载水箱护栏优化设计方法 |
CN111931394B (zh) * | 2020-05-20 | 2022-07-08 | 吉林大学 | 非概率混合可靠度指数的去嵌套分析方法 |
CN111625908B (zh) * | 2020-05-25 | 2023-10-20 | 西安石油大学 | 多约束下深海连接器密封结构时变可靠性优化设计方法 |
CN111783327B (zh) * | 2020-05-28 | 2022-09-09 | 西北工业大学 | 基于支持向量机的涡轮叶片多场载荷下的可靠性分析方法 |
CN111859482B (zh) * | 2020-06-09 | 2024-05-17 | 大连理工大学 | 一种肌肉式薄壁结构轻量化设计方法 |
CN111832102B (zh) * | 2020-06-19 | 2024-02-06 | 浙江大学 | 一种高维随机场条件下的新型复合材料结构优化设计方法 |
CN111967198B (zh) * | 2020-09-11 | 2023-12-12 | 浙江大学 | 一种新型复合材料结构高承载力轻量化设计方法 |
CN112182933A (zh) * | 2020-09-28 | 2021-01-05 | 西北工业大学 | 一种变刚度复合材料圆柱壳体的建模及优化设计方法 |
CN112307654B (zh) * | 2020-10-10 | 2022-05-10 | 南昌航空大学 | 基于混合可靠性模型与神经网络响应面的pma可靠性优化方法 |
CN112257180B (zh) * | 2020-10-23 | 2024-04-26 | 中国科学院微小卫星创新研究院 | 航天复合材料层合板可靠性分析系统及方法 |
CN112199880A (zh) * | 2020-10-27 | 2021-01-08 | 北京机电工程总体设计部 | 一种天线罩参数化及强度快速分析仿真系统 |
CN112948959B (zh) * | 2020-11-24 | 2021-12-21 | 合肥工业大学 | 一种固体火箭发动机壳体应力计算方法 |
CN112926161B (zh) * | 2021-03-31 | 2022-12-06 | 华侨大学 | 一种空间曲壳结构形状优化方法 |
CN113111437B (zh) * | 2021-04-16 | 2023-07-21 | 中国航空综合技术研究所 | 数据融合的飞机运动机构可靠性评估方法 |
CN113515810B (zh) * | 2021-05-17 | 2022-08-26 | 中车长春轨道客车股份有限公司 | 基于可靠性和安全性分析的动车组转向架设计开发方法 |
CN113312717B (zh) * | 2021-05-21 | 2023-03-21 | 西北工业大学 | 一种应用于低损失压气机中介机匣的设计方法 |
CN113343352B (zh) * | 2021-05-28 | 2022-06-03 | 重庆长安汽车股份有限公司 | 一种基于鲁棒优化的底盘零部件目标分解方法 |
CN113656900B (zh) * | 2021-07-06 | 2023-06-20 | 扬州大学 | 基于灵敏度修正的收割机结构强度影响因素分析方法 |
CN113779825B (zh) * | 2021-08-24 | 2024-04-19 | 华中科技大学 | 一种基于单位分解及移动控制点的结构优化设计方法 |
CN113722779A (zh) * | 2021-09-08 | 2021-11-30 | 大连理工大学 | 一种基于薄壳结构的参数化雕刻设计方法 |
CN114048547B (zh) * | 2021-11-18 | 2022-07-05 | 汉思科特(盐城)减震技术有限公司 | 一种基于自适应代理模型的车用空气弹簧工程优化设计方法 |
CN115017653A (zh) * | 2022-07-08 | 2022-09-06 | 北京科技大学 | 基于等几何设计/分析/制造框架的设计方法 |
CN115169008B (zh) * | 2022-07-27 | 2023-01-03 | 中车成型科技(青岛)有限公司 | 混合材料车体工程化轻量化方法及系统 |
CN115270585A (zh) * | 2022-09-28 | 2022-11-01 | 西北工业大学 | 一种机床床身与地脚的协同优化设计方法及相关产品 |
CN115688311B (zh) * | 2022-10-27 | 2023-06-23 | 苏州科技大学 | 一种行星滚柱丝杠副的不确定性分析方法及系统 |
CN116011301B (zh) * | 2023-02-13 | 2024-01-23 | 东华大学 | B样条等几何状态空间有限元方法 |
CN116663362B (zh) * | 2023-05-31 | 2023-12-01 | 上海交通大学 | 一种航空航天复合材料薄壁结构装配力学建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096736A (zh) * | 2011-01-28 | 2011-06-15 | 重庆大学 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
CN103399992A (zh) * | 2013-07-22 | 2013-11-20 | 中国兵器科学研究院 | 一种基于可靠寿命的结构耐久性优化设计方法 |
CN103870681A (zh) * | 2014-02-26 | 2014-06-18 | 浙江工业大学 | 带锯床混合型可靠性分析计算方法 |
CN105404732A (zh) * | 2015-11-04 | 2016-03-16 | 中国直升机设计研究所 | 一种基于随机正态分布的复合材料层压板铺层优化方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8292226B2 (en) * | 2005-11-15 | 2012-10-23 | The Boeing Company | Weight-optimizing internally pressurized composite-body aircraft fuselages having near-elliptical cross sections |
US8182234B2 (en) * | 2008-07-09 | 2012-05-22 | General Electric Company | Spaceframe wind turbine energy converter structure |
US8756037B2 (en) * | 2011-06-20 | 2014-06-17 | The Boeing Company | Curved fiber paths for composite laminates |
US9274036B2 (en) * | 2013-12-13 | 2016-03-01 | King Fahd University Of Petroleum And Minerals | Method and apparatus for characterizing composite materials using an artificial neural network |
CN105757440B (zh) * | 2014-12-15 | 2019-04-05 | 南通中集翌科新材料开发有限公司 | 底板单元 |
US10183449B2 (en) * | 2015-09-25 | 2019-01-22 | The Boeing Company | Lamination parameter-based method for optimal design and manufacturing options |
US10814564B2 (en) * | 2017-10-11 | 2020-10-27 | Divergent Technologies, Inc. | Composite material inlay in additively manufactured structures |
-
2017
- 2017-09-13 CN CN201710817832.6A patent/CN107526898B/zh active Active
- 2017-12-21 US US15/850,909 patent/US10726168B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096736A (zh) * | 2011-01-28 | 2011-06-15 | 重庆大学 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
CN103399992A (zh) * | 2013-07-22 | 2013-11-20 | 中国兵器科学研究院 | 一种基于可靠寿命的结构耐久性优化设计方法 |
CN103870681A (zh) * | 2014-02-26 | 2014-06-18 | 浙江工业大学 | 带锯床混合型可靠性分析计算方法 |
CN105404732A (zh) * | 2015-11-04 | 2016-03-16 | 中国直升机设计研究所 | 一种基于随机正态分布的复合材料层压板铺层优化方法 |
Non-Patent Citations (2)
Title |
---|
功能梯度复合材料板壳结构的弯曲、屈曲和振动;沈惠申;《力学进展》;20040624;第53-60页 * |
复合材料薄壳稳定性问题研究;陈鲲鹏;《万方数据库》;20161114;全文 * |
Also Published As
Publication number | Publication date |
---|---|
US20190080040A1 (en) | 2019-03-14 |
CN107526898A (zh) | 2017-12-29 |
US10726168B2 (en) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107526898B (zh) | 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 | |
Kassapoglou | Design and analysis of composite structures: with applications to aerospace structures | |
Vescovini et al. | A variable-kinematic model for variable stiffness plates: Vibration and buckling analysis | |
Dey et al. | A response surface modelling approach for resonance driven reliability based optimization of composite shells | |
CN110222356B (zh) | 综合考虑稳定性与振动特性的板/壳结构轻量化拓扑优化设计方法 | |
CN112417603A (zh) | 一种飞机壁板定位夹紧装配变形预测方法 | |
Wang et al. | An efficient single-loop strategy for reliability-based multidisciplinary design optimization under non-probabilistic set theory | |
López et al. | Reliability-based design optimization of composite stiffened panels in post-buckling regime | |
Coroneos et al. | Structural analysis and optimization of a composite fan blade for future aircraft engine | |
Hao et al. | Adaptive approximation-based optimization of composite advanced grid-stiffened cylinder | |
Patni et al. | On the accuracy of localised 3D stress fields in tow-steered laminated composite structures | |
KR20120082112A (ko) | 복합재 구조 분석 시스템과 분석방법 | |
CN111950091A (zh) | 一种基于区间不确定性的复合材料层合板结构非概率可靠性拓扑优化方法 | |
CN116011301B (zh) | B样条等几何状态空间有限元方法 | |
Wang et al. | Efficient computational method for the non-probabilistic reliability of linear structural systems | |
Noevere et al. | Development of a design for manufacturing tool for automated fiber placement structures | |
Kumar et al. | Investigations on the multistability of series-connected unsymmetric laminates | |
Murugesan et al. | Multi-objective design optimization of composite stiffened panel using response surface methodology | |
Franz et al. | Variation analysis of design parameters of fibre-reinforced plastic parts | |
Andraju et al. | Three-dimensional stress analysis of variable angle tow composite laminate using hybrid brick elements | |
Kothari | An Integrated Design Tool for Tow-steered Composite Laminate in Abaqus | |
Maksimović et al. | Structural Analysis and Optimization of Layered Composite Structures: Numerical and Experimental Investigations | |
Yun et al. | An innovative adaptive Kriging-based parallel system reliability method under error stopping criterion for efficiently analyzing the global reliability sensitivity index | |
Garofano et al. | On the use of double‐double design philosophy in the redesign of composite fuselage barrel frame components | |
CN115344959A (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 |