一种复合材料层合结构的非概率动力可靠性评估方法
【专利摘要】本发明公开了一种面向全寿命周期的复合材料层合结构非概率动力可靠性评估方法。该方法综合考虑复材层合结构载荷的动态波动效应以及材料性能的累积退化，基于有限元思想和区间数学方法，构建了该结构动特性的非概率区间过程模型；进而结合首次穿越理论和轻质层合结构的失效判定准则，定义了结构非概率动力可靠性指标，并探索了高效稳健的求解策略。本发明在进行可靠度计算过程中合理表征了不确定性对结构全生命周期内动力安全的综合影响，为确保对其开展精细化设计提供必要的数据参考。
【专利说明】
-种复合材料层合结构的非概率动力可靠性评估方法
技术领域
[0001] 本发明设及轻质复材层合结构的安全性评估技术领域，特别设及考虑不确定性、 载荷交变及材料退化共同作用下结构安全性能的有效认知与定量表征，为合理制定复合材 料层合结构的选材及铺层方案设计提供重要的理论支持。
【背景技术】
[0002] 由于质量轻、强度高、多功能，使复合材料层合板成为具有优良特性的工程结构 件。此外，鉴于其自身具有结构设计裕度大、性能匹配度高等天然优势，复材层合结构不仅 被大量应用于航空航天、船舶、兵器等军工领域，还作为最基本和最主要的构件频繁出现在 民用客机、汽车、建筑等民用结构系统中。因此，针对复合材料层合结构的力学特性分析与 安全态势评估技术研究具有重要的理论意义与工程实用价值。
[0003] 然而，复合材料由于成型工艺W及加工批次的差异，材料缺陷导致的分散性不可 避免。此外，结构的服役环境复杂多样，外部激励的未确知性同样客观存在。加之随着服役 时间的累积，材料性能的退化与载荷的交变效应严重影响着复合材料层合结构的实际使用 性能，如若无法有效准确地预估结构的真实服役状态，将导致重大的安全事故发生。综合上 述情况，针对轻质复材层合板结构开展不确定性传播分析与动力可靠性评估方法研究已受 到学术界和工程界的高度重视。
[0004] 当前，国内外学者与工程技术人员对复合材料结构的可靠性分析方法研究主要集 中在两个方面：（1)基于概率统计理论的结构不确定性传播影响预测技术；（2)基于准静态 假设的可靠度计算方法研究。上述工作一定程度上丰富了复合材料结构的分析与强度理 论，但是忽略了随机方法对样本信息的依赖性W及时间累积效应下结构失效事件的相关 性，大大限制了其理论的工程实用化进程。
[0005] 由于实际工程中贫信息、少数据的情况时有发生，建立W非概率理论框架为基础 的不确定性表征技术、时变可靠度建模与求解技术具有显著的现实意义。目前，相关研究工 作尚不成熟，针对复合材料层合结构的方案设计经常无法严格满足所需的应用要求，亦或 是安全冗余度过大，造成严重的资源浪费与时间成本损耗。鉴于此，本发明将重点探究静动 力不确定性作用下轻质层合结构安全性能的退化历程，为其轻量化、多功能设计提供理论 保障。
【发明内容】
[0006] 本发明解决的技术问题是：克服现有技术的不足，提供一种针对轻质复合材料层 合板结构的安全性评价方法，充分考虑实际工程问题中普遍存在的静动力不确定性因素， W提出的非概率动力可靠性度量评判结构安全与否的量化指标，所得到的强度校核结果更 加符合真实情况，工程适用性更强。
[0007] 本发明采用的技术方案实现步骤如下：
[000引第一步:根据复材层合结构的本构特征及载荷边界条件，基于最小势能原理构建 有限元列式，推导出如下动力方程：
庚中， 矣多，S分别表示节点加速度、速度和位移，M，C，K分别为总体质量矩阵、总体阻尼矩阵和总 体刚度矩阵，F为载荷列向量，N为节点个数。运里，复合材料层合板结构总体刚度矩阵的计 算如下巧A击夫巧，
[0009]
[0010]其中，V代表全域体积，Vk表示第k层板的体积，n为层板总数，A为层板面积，B为应 变矩阵。
[0011] 第二步：引入强度比R，结合Tsai-Wu强度失效判定准则，代入到步骤一中的动力方 程中，可得到临界失效载荷的显式表达式：
[0012] F 化 max+Fu〇imax〇 jmax =li，j = l，2，...6
[0013]
[0014] 其中，表示强度特征量，O= [01，02，Ti2]%应力向量，Oimax对应最危险点处 的应力值，角标i和j分别表示单元对应的自由度计数指标。运里所述的强度比R其物理意义 为极限载荷与真实载荷的比值，具体公式为：
[0015]
[0016] 其中，分别定义中间参量Coe_A = Fij〇i〇j和Coe_B = Fi〇i，〇i和Oj分别对应计数指标i 和j处的单元应力。强度特征量Fi和Fu取决于材料工程强度参数{扣，杞，￥*，￥。，52}，分别表示 材料绑向掠伸端底、纵向店缩端底、描向掠伸端底、横向压缩强度W及剪切强度，即：
[0017;
[0018；
[0019] 第=步:综合考虑存在于第一步和第二步中结构参数中的不确定性效应，不确定 性效应包括材料分散性引起的M，C，K的变动、载荷F的未确知性引起响应足矣S的变化、强 度判定准则的模糊性导致强度特征量Fi和Fu的不一致。引入区间向量XExI = Ui,X2,...， Xm似及区间过程模型X(t)GX(t)i=化(〇龙(〇，...而(0)，得到有限样本信息条件下 结构静动力不确定性参数的数学表达：
[0020]
[0021] X(t)U/L=(Xi(t)U/L，X2(t)U/L, . . .，Xn(t)U/L) = (Xi(t)c±Xl(t)r，X2(t)c±X2(t )T, . . .，Xm(t)e±Xm(t)T)其中，上标U代表参量的取值上界，上标L代表参量的取值下界，上标 C代表中屯、值，上标r代表半径。静动力不确定性参数包括:载荷边界参数、材料特征参数、测 量精度参数W及设计准则(设计许用值)参数等，为便于计算，通常表示为标准化形式，具体 如下：
[0022] X= [xL'X。] = [xC-xT'xC+xr]
[0023] =义。+义叮-1,1]
[0024] X Cl
[0025] 和
[0026] X(t) = [XHt),X'^(t)] = [X"(t)-X"(t),X"(t)+X"(t)]
[0027] =X"(t)+X"(t)[-l,l]
[002引=xc(t)+XT(t)X 号2
[0029] 其中，X。= . . .，Xm。)和 X'= (XiT,X2T，. . .，?/)表示区间向量 X 的均值和半 径，XC(t) = (XlC(t)，X2C(t)，. . .，XnC(t))和 XC(t) = (Xir(t)，X2T(t)，. . .，XnT(t))表示区间过 程X(t)的均值和半径点点G曰4,曰4定义为所有元素包含在[-1，U内的4维向量集合，符号 "X"定义为两个向量各对应元素相乘的算子，乘积仍为维数为4的向量。
[0030] 第四步:将第=步中表征的静动力不确定参数代入到第一步和第二步的方程中， 形成考虑时间效应的区间集格式，即给定任意时刻有：
[0031]
。
[0032] 其中，矿(.V，AXO)，各Yv, 乂(〇)/; I (X，X( t))分别为考虑不确定性影响下的加速度、速 度和位移历程，11^，乂(*))，(：1^，乂(*))，1(1^，乂(*))分别表示不确定性影响下的总体质量、 总体阻尼和总体刚度矩阵历程。
[0033] 对上述连续时刻下的平衡方程进行时间离散化处理，结合化wmark迭代求解算法 和区间差分格式，实现任意离散点处临界许用载荷的上下界计算，即：
[0034]
[0035]
[0036] 进而，获得完整寿命期内复材层合结构临界许用载荷巧的区间历程，即 上界端胃.托X似)和下界蠕f片，X(W)。上式中，Si(x，X(t))的计算采用节点位移的差分迭 代格式，针对确定性问题，Si(x，X(t))采用退化的S(tk)表示，计算如下：
[0037]
[003引其中，tk表示第k个离散时刻，即tk = k A t，S(tk)表示tk时刻下的节点位移向量，S (0)，如0)，<j(0)表示初始位移向量、速度向量和加速度向量，j和k分别表示迭代离散过程的 计数指标。并且：
[0039]
[0040]
[0041]
[0042] 第五步：根据实际加载历程Pi(t)，构建与步骤四获得的临界许用载荷历程 间的应力-强度区间过程干设模型，并定义时变极限状态函数如下：
[0043] …
[0044] 引入非概率区间过程理论，实现任意离散时刻极限状态中屯、值G。、半径值护W及任 意微小时间增量[i A t，（i+l) A t]内协方差CovcQ A t，（i+l) A t)和相关系数PcQ A t，（i+ 1) At)的显式表达。运里，时变极限状态函数G(t，x，X(t))的定义及相关性特征的求解过程 如下：
[0045] 极限状态函数G(t，x，X(t))同样被定义为一个区间过程，即G(t，x，X(t))GGi(t, x，X(t))。对于任意给定瞬时时刻ti = iAt，极限状态将退化为一个区间量G(ti，x，X(ti))， 有限个离散区间量的可行范围被约束在一个超立方体域QD内。引入非概率区间过程理论， 实现任意离散时刻极限状态中屯、值护、半径值护W及任意微小时间增量[i A t，（i+1) A t]内 协方差CovgQ A t, (i+1) A t)和相关系数f>G(i A t, (i+1) A t)的显式表达。其中，CoVG(i A t, (i+1) A t)和PG(i A t，（i+1) A t)的定义需借助标准化处理手段，并转换工作坐标系至化1， 化），即：
[0046] G(ti,x,X(ti))G[GHti,x,X(ti)),G^ti,x,X(ti))]
[0047] =GC(ti，x，X(ti))+Gr(ti，x，X(ti))XUi
[004引 G(ti+i，x，X(ti+i))G[GL(ti+ii，x，X(ti+i))，GU(ti+i，x，X(ti+i))]
[0049] =GC(ti+i，x，X(ti+i))+护（ti+i，x，X(ti+i))XU2
[0050] 其中，GC(ti，x，X(ti))和护(ti，x，X(ti))表示ti时刻极限状态函数的均值和半径， (ti+i，x，X(ti+i))和护(ti+i，x，X(ti+i))表示心拥刻极限状态函数的均值和半径。
[0051] 第六步:将首次穿越理论与步骤五中建立的复合材料层合结构极限状态的应力- 强度区间过程干设模型相结合，可W得到任意时间区间内的穿越可能度：
其中，Posl ?}表示事件发生 的可能性度量，EiAt表示穿越事件，即事件A: Q A t)时刻结构安全G(i A t)〉0与事件B: (a+ I) A t)时刻结构失效G((i+1) A t)<0的交事件，A t表示微小时间增量。
[0052] 第屯步:遍历所有时间段内的穿越可能度化s{EiAt}，定义复合材料层板结构的非 概率动力前靠底计當指标，
[0化3]
[0054]其中，Rs(T)表示整个生命周期T内的动力可靠度，GQ A t)为时刻i A t的极限状态 函数，巧,,,U，X(/A〇)和pi(i A t)分别对应时刻i A t的临界许用载荷和真实加载，求解上式 即实现复材层合结构动力安全态势的有效评估。运里，Rs(T)的求解计算需借助时间离散化 方法进行简化近似，通过遍历每一个微小时间增量内结构发生穿越破坏的可能性指标，经 叠加运算后得到可靠度的等效表达式如下：
[0化5]
[0056]其中，Pf (T)表示整个生命周期内的失效度，
表 示结构在初始时刻即发生失效的可能度，A t = tw-ti表示微小时间增量，其取值设定为整 个服役寿命的1/1000。
[0057] 本发明与现有技术相比的优点在于:本发明提供了含静动力不确定性复合材料层 合结构可靠性分析的新思路，弥补和完善了传统基于概率理论的静态可靠性分析方法的局 限性。所构建的动力可靠性度量模型，一方面可大幅减小对样本信息的依赖性，另一方面可 有效计及并量化材料退化与载荷交变作用下结构安全性能的折减效应。在对结构极限载荷 进行定界计算时，高效稳健的迭代算法与区间差分格式相结合，可确保载荷信息输入条件 的可信性;基于首次穿越的可靠度计算方法更为真实地反映了时间相关性对结构安全的综 合影响，为复合材料结构精细化强度分析理论的完善作出了积极贡献。
【附图说明】
[0058] 图1是本发明针对复合材料层合结构的非概率动力可靠性评估流程图；
[0059] 图2是本发明针对复合材料层合结构有限元建模过程中坐标变换示意图；
[0060] 图3是本发明定义的极限状态相关性函数所对应的几何可行域示意图；
[0061 ]图4是本发明采用的首次穿越方法示意图；
[0062] 图5是本发明提出的微小时间段内穿越失效可能度计算方法示意图；
[0063] 图6是本发明实施例中复合材料层合结构的几何模型示意图；
[0064] 图7是本发明实施例中临界许用载荷与真实载荷的区间过程干设情况示意图。
【具体实施方式】
[0065] 如图1所示，本发明提出了一种复合材料层合结构的非概率动力可靠性评估方法， 包括W下步骤：
[0066] (1)根据复材层合结构的本构特征及载荷边界条件，基于最小势能原理构建有限 元列式，推导出如下动力方程：
其 中，表式S分别表示节点加速度、速度和位移，M，C，K分别为总体质量矩阵、总体阻尼矩阵和 总体刚度矩阵，F为载荷列向量，N为节点个数。运里，复合材料层合板结构总体刚度矩阵的 计算如下巧A式亲巧：
[0067]
[0068] 其中，V代表全域体积，Vk表示第k层板的体积，n为层板总数，A为层板面积，B为应 变矩阵。上式中，应变矩阵B= {Bi}如下表达式：
[0069] e =BSe= [Bi B2 B3 B4]Se
[0070] 和
[0071]
[0072] 其中，e为应变向量，SB为单元位移场，Ni表示第i个节点自由度对应的形函数。弹性 矩阵D可通过下式获得：
[0073] 〇wz = [T]-i〇=[T]-i[C]e = [T]-i[C][T]-Tewz=W]exyz
[0074] 其中，Oxyz和Exyz表示单元在全局坐标系下的应力场和应变场，T为坐标转置矩阵 (如图2所示）。
[0075] (2)引入强度比R，结合Tsai-Wu强度失效判定准则，代入到步骤一中的动力方程 中，可得到临界失效载荷的显式表达式：
[0076] Fi〇imax+Fij〇imax〇jmax=l(i , j = l ,2, ---6)
[0077]
[007引其中，尸1和枯表示强度特征量，O= [01，02, Ti2]T为应力向量，Oimax对应最危险点处 的应力值，角标i和j分别表示单元对应的自由度计数指标。运里所述的强度比R其物理意义 为极限载荷与真实载荷的比值，具体公式为：
[0079]
[0080] 其中，分别定义中间参量Coe_A = Fij〇i〇j和Coe_B = Fi〇i，〇i和Oj分别对应计数指标i 和j处的单元应力。强度特征量Fi和Fu取决于材料工程强度参数{沿，也，￥*，片，52}，分别表示 材料纵向扮化端睛.幼向圧縮端睛喊向扮化端睛横向压缩强度W及剪切强度，即：
[0081]
[0082]
[0083] (3)综合考虑存在于第一步和第二步中结构参数中的不确定性效应，不确定性效 应包括材料分散性引起的M，C，K的变动、载荷F的未确知性引起响应泉矣S的变化、强度判 定准则的模糊性导致强度特征量Fi和Fu的不一致。引入区间向量xGxI=(X1，X2, . . .，Xm)W 及区间过程模型X(t)GX(t)I = (Xl(t)，X2(t)，...，Xn(t))，得到有限样本信息条件下结构 静动力不确定性参数的数学表达：
[0084]
[00化]X(t)U/L=(Xi(t)UA，X2(t)U/L, . . .，Xn(t)U/L) = (Xi(t)c±Xl(t)r，X2(t)c±X2(t )T，. . .，Xm(t)e±Xm(t)T)其中，上标U代表参量的取值上界，上标L代表参量的取值下界，上标 C代表中屯、值，上标r代表半径。静动力不确定性参数包括:载荷边界参数、材料特征参数、测 量精度参数W及设计准则(设计许用值)参数等，为便于计算，通常表示为标准化形式，具体 如下：
[0086] X= [yL'xU] = [xC-xT'xC+xr]
[0087] =x""+x^[-l, 1]
[008引=xe+xTX 号I
[0089] 和
[0090] X(t) = [XHt),X'^(t)] = [X"(t)-X"(t),X"(t)+X"(t)]
[0091] =X"(t)+X"(t)[-l,l]
[0092] =X"(t)+X"(t)X|2
[0093] 其中，X^= . . . ,?!^)和 X'= (XiT,X2T，. . .，?/)表示区间向量 X 的均值和半 径，XC(t) = (XlC(t)，X2C(t)，. . .，XnC(t))和 XC(t) = (Xir(t)，X2T(t)，. . .，XnT(t))表示区间过 程X(t)的均值和半径点点G曰4,曰4定义为所有元素包含在[-1，U内的4维向量集合，符号 "X"定义为两个向量各对应元素相乘的算子，乘积仍为维数为4的向量。
[0094] (4)将第=步中表征的静动力不确定参数代入到第一步和第二步的方程中，形成 考虑时间泌脉的反间隹格击.邮给宙件意时刻有，
[0095]
[0096] 其中，沪O，义脚、沪U义卿卢(X，X( t))分别为考虑不确定性影响下的加速度、速 度和位移历程，11^，乂(*))，(：1^，乂(*))，1(1^，乂(*))分别表示不确定性影响下的总体质量、 总体阻尼和总体刚度矩阵历程。
[0097] 对上述连续时刻下的平衡方程进行时间离散化处理，结合化wmark迭代求解算法 和区间差分格式，实现任意离散点处临界许用载荷的上下界计算，即：
[009引
[0099]
[0100] 进而，获得完整寿命期内复材层合结构临界许用载荷巧4^批，巧〇)的区间历程，即 上界端似巧日下界瑞帛从义把)）。上式中，Si(x，X(t))的计算采用节点位移的差分迭 代格式，针对确定性问题，Si(x，X(t))采用退化的S(tk)表示，计算如下：
[0101]
[0102]其中，tk表示第k个离散时刻，即tk = kAt，S(tk)表示tk时刻下的节点位移向量，S (0)，j(0)，新0)表示初始位移向量、速度向量和加速度向量，j和k分别表示迭代离散过程的 计数指标。并且：
[0103]
[0104]
[0105]
[0106] 根据节点位移的动态边界信息，可知应力场界限为：[0107] 〇L(tk)=min{D(x，X(t))B(x，X(t))SL(tk)WP〇u(tk)=max{D(x，X(t))B(x，X(t))SU (tk)}[0108] 因此，任意离散点处临界许用载荷的上下界可进一步表示为：[0109]
[
[
[0112] (5)根据实际加载历程Pi(t)，构建与步骤四获得的临界许用载荷历程 的)间的应力-强度区间过程干设模型，并定义时变极限状态函数如下：
[0113]
[0114] 引入非概率区间过程理论，实现任意离散时刻极限状态中屯、值G。、半径值护W及任 意微小时间增量[i A t, (i+1) A t]内协方差CovgQ A t, (i+1) A t)和相关系数f>G(i A t, Q+ 1) At)的显式表达。运里，时变极限状态函数G(t，x，X(t))的定义及相关性特征的求解过程 如下：
[0115] 极限状态函数G(t，x，X(t))同样被定义为一个区间过程，即G(t，x，X(t))GGi(t, x，X(t))。对于任意给定瞬时时刻ti = iAt，极限状态将退化为一个区间量G(ti，x，X(ti))， 有限个离散区间量的可行范围被约束在一个超立方体域Q n内。
[0116] 基于区间数学理论，分别定义出任意给定时刻ti下中屯、值护(*1，义，乂(*1))和半径护 (ti,x,X(ti))如下：
[0117]
[011 引
[0119]
[0120]
[0121] 利用正则化手段，首先得到：
[0122] G(ti，x，X(ti))G[GL(ti，x，X(ti))，GU(ti，x，X(ti))]
[0123] =GC(ti,x,X(ti))+Gr(ti,x,X(ti)) XUi
[0124] G(ti+i,x,X(ti+i)) G [G^(ti+ii,x,X(ti+i)) ,G^(ti+i,x,X(ti+i))]
[0125] =GC(ti+i,x,X(ti+i))+Gr(ti+i,x,X(ti+i)) XU2
[0126] 其中，护(ti，x，X(ti))和护(ti，x，X(ti))表示ti时刻极限状态函数的均值和半径，护 (ti+i，X，X( ti+i))和护(ti+i，X，X( ti+i))表示ti+拥刻极限状态函数的均值和半径。
[0127] 从几何角度不难发现，存在无数多个不同形状的偏转矩形域包含于标准方形域 内，而运些矩形域形状的改变与其对应的极限状态相关性具有映射关系（如图3所示）。于 是，定义协力'差巧1猶"[0乂「,（1 A t. (i +1 ) A t)末时目卖^系猶日G(i A t. (i +1 ) A t)々n下；
[012 引
[0129]
[0130] 其中，d表示如图3所示矩形域边长的一半，馬和化,分别是标准区间变量化和化的 方差（化，二A ^ = I )，PG(i A t，（i + 1) A t)是一个无量纲量，其大小代表了G(ti)和G(ti+1) 的线性相关度。
[0131] 综上，本发明实现了对复合材料层合结构强度失效的定量表征，为后续动力可靠 性建模及求解提供了必要的理论依据。
[0132] (6)将首次穿越理论(如图4所示）与与步骤五中建立的复合材料层合结构极限状 态的应力-强度区间过程干设模型相结合，可W得到任意时间区间内的穿越可能度：
，其中，Posl ?}表示事件发生
的可能性度量，EiAt表示穿越事件，即事件A: Q A t)时刻结构安全G(i A t)〉0与事件B: (a+ I) A t)时刻结构失效G((i+1) A t)<0的交事件，A t表示微小时间增量。运里,引入面积比思 想(如图5所示），化s{EiAt}可定义为穿越几何条件与极限状态可行域的干设面积与过程中 总可行域(偏掉巧形）而巧，比，邮：
[0133]
[0134] 斗Sf的计算通常是一个分段函数，需要结合几何边界与可行域的相交条件分类讨论。
[0135] (7)遍历所有时间段内的穿越可能度Pos{EiAt}，定义复合材料层板结构的非概率
时变可告曲A咎化枯-
[0136]
[0137] 其中，Rs(T)表示整个生命周期T内的动力可靠度，GQ A t)为时刻i A t的极限状态 函数，巧;lf.(x，X(z'A/))和pi(i A t)分别对应时亥Iji A t的临界许用载荷和真实加载，求解上式 即实现复材层合结构动力安全态势的有效评估。运里，Rs(T)的求解计算需借助时间离散化 方法进行简化近似，通过遍历每一个微小时间增量内结构发生穿越破坏的可能性指标，经
叠加运算 巨婚至Il前靖:睹斷锋就車化才加下.
[013 引
[0139] 其中，Pf(T)表示整个生命周期内的失效度
表 示结构在初始时刻即发生失效的可能度，A t = tw-ti表示微小时间增量，其取值设定为整 个服役寿命的1/1000。
[0140] 实施例；
[0141] 为了更充分地了解该发明的特点及其对工程实际的适用性，本发明针对如图6所 示24层复合材料层合板结构进行非概率动力可靠性分析。其采用对称铺层方案[0/0/0/0/ 0/0/-0/-0/-0/-0/-0/-0]对称。该结构受到作用于几何中屯、的集中载荷，并采用四边固支加 W约束。板的长度巧日宽度b均为100mm，每层板厚度为t = 0.147mm。层板材料具有横向各向 同性，其密度为P=1.38X 103kg/m3。该层合板工程强度参数信息如表1所示，表2列出了层板 模量的非概率静动力不确定性特征。
[0142] 表1
[0143]
[0146] 根据本发明提出的方法可W快速获得临界许用载荷的区间过程表达式，即： 句,,,.(0，￡;(化￡2(化知(〇，成.馬;）。分别讨论目=15°和目=45°两种铺层形式，对应的载荷工 况分别为：
[0149]于是，可获得上述两种动态载荷工况复材层合结构强度约束下的极限状态函数G
[0147]
[014 引 (t) Ie=I日。和G(t)|e=4日。。图7显示了临界许用载荷与真实载荷的区间过程干设情况。基于前述 的动力可靠性评价方法，两种工况对应的可靠度结果分别为：（1)9 = 15° ,Rs(T) = O.53; (2) 0 = 45° ,Rs(T)=O.85。其中，全寿命周期为T = 20年。从结果可W看出，不同铺层形式下复合 材料层合结构的强度性能差异明显，也为后续最优设计提供了较大的控制裕度。此外，动力 可靠度可W有效量化结构的安全性，可靠度越高，结构越安全，反之越危险。
[0150] 综上所述，本发明提出了一种复合材料层合结构的非概率动力可靠性评估方法。 首先，基于有限元思想和最小势能原理，构建并解析结构的平衡方程;其次，引入非概率静 动力不确定性参量，结合数值迭代算法与有限差分思想，实现复材层合结构临界许用载荷 的区间过程表述，进而构建极限状态函数的显式表达式;基于首次穿越理论和离散化策略， 完成任意微小时间段内穿越失效可能度的几何定义与求解;最后，通过遍历所有时间区间 的穿越可能度函数，叠加运算后计算动力可靠度指标，构造出复合材料层合结构动力安全 的量化判据。
[0151] W上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制；其可扩展应 用于含缺陷结构的优化设计领域，凡采用等同变换或者等效替换而形成的技术方案，均落 在本发明权利保护范围之内。
[0152] 本发明未详细阐述部分属于本领域技术人员的公知技术。
【主权项】
1. 一种复合材料层合结构的非概率动力可靠性评估方法，其特征在于实现步骤如下： 第一步:根据复材层合结构的本构特征及载荷边界条件，基于最小势能原理构建有限 元列式，建立如下动力方程v，其 中，#，武s分别表示节点加速度、速度和位移，M，C，K分别为总体质量矩阵、总体阻尼矩阵和 总体刚度矩阵，F为载荷列向量，N为节点个数； 第二步：引入强度比R，结合Tsai-Wu强度失效判定准则，代入到第一步中建立的动力方 程中，得到临界失效载荷的显式表达式： Fi〇imax+Fij〇imax〇jmax- 1 ? , j - 1，2，···6其中，Fi和Fij表示强度特征量，σ = ^^^^^^为应力向量…^对应最危险点处的应 力值，角标i和j分别表示单元对应的自由度计数指标； 第三步:综合考虑存在于第一步和第二步中结构参数中的不确定性效应，不确定性效 应包括材料分散性引起的Μ，C，K的变动、载荷F的未确知性引起响应忒次δ的变化、强度判 定准则的模糊性导致强度特征量Fi和Fij的不一致；引入区间向量...，x m)以 及区间过程模型= 得到有限样本信息条件下结构 静动力不确定性参数的数学表达：=UV,<厂，尤:(/)'，.….义⑴，±冬(，广土… 其中，上标U代表参量的取值上界，上标L代表参量的取值下界，上标c代表中心值，上标 r代表半径;静动力不确定性参数包括:载荷边界参数、材料特征参数、测量精度参数以及设 计准则即设计许用值参数，为便于计算，通常表示为标准化形式； 第四步:将第三步中表征的静动力不确定参数代入到第一步和第二步的方程中，形成 考虑时间效应的区间集格式，即给定任意时刻有：其中，^(^^^^(.^'(/:^^^⑴彡分别为考虑不确定性影响下的加速度~速度和 位移历程，11(1，(〇)，(：1(1，(〇)，1(1(13(〇)分别表示不确定性影响下的总体质量、总体 阻尼和总体刚度矩阵历程； 对上述连续时刻下的平衡方程进行时间离散化处理，结合Newmark迭代求解算法和区 间差分格式，实现任意离散点处临界许用载荷的上下界计算，即：进而，获得完整寿命期内复材层合结构临界许用载荷⑴)的区间历程，即上界 校界(X，))和下界界〇, Xft)); 第五步:根据实际加载历程PHt)，构建与第四步获得的临界许用载荷区间 历程间的应力-强度区间过程干涉模型，并建立时变极限状态函数如下：引入非概率区间过程理论，实现任意离散时刻极限状态中心值Ge、半径值(7以及任意微 小时间增量[i Δ t，（i+Ι) Δ t]内协方差Covc(i Δ t，（i+Ι) Δ t)和相关系数PG(i Δ t，（i+Ι) Δ t)的显式表达； 第六步:将首次穿越理论与步骤五中建立的复合材料层合结构时变极限状态函数相结 合，得到任意时间内的穿越可能度：办4匙} W1心丨印~卜< 〇}，其 中，P〇s{ · }表示事件发生的可能性度量，EiAt表示穿越事件，即事件A:(iAt)时刻结构安全 G(i △ t)>0与事件B: ((i+1) △ t)时刻结构失效G((i+1) △ t)〈0的交事件，△ t表示微小时间 增量； 第七步:遍历所有时间段内的穿越可能度P〇s{ElAt}，计算复合材料层板结构的非概率 动力可靠度计算指标：其中，RS(T)表示整个生命周期T内的动力可靠度，G(i Δ t)为时刻i Δ t的极限状态函数， (谢))和Ρ? Δ t)分别对应时亥Iji Δ t的临界许用载荷和真实加载，求解上式即实 现复材层合结构动力安全态势的有效评估。2. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第一步中总体刚度矩阵的计算如下式：其中，V代表全域体积，Vk表示第k层板的体积，η为层板总数，A为层板面积，B为应变矩 阵。3. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第二步中强度比R的公式为：其中，分别定义中间参量(：〇θ_Α = Ρυσ?σ」和Coe_B = Fi〇i，(^P〇j分别对应计数指标i和j 处的单元应力，强度特征量Fi和Fij取决于材料工程强度参数{Xt，X。，Yt，Y。，S 2}，分别表示材 料纵向拉伸强度、纵向压缩强度、横向拉伸强度、横向压缩强度以及剪切强度，即：4. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第三步中区间向量X和区间过程X( t)的标准化形式可表示为： x= [xL,xu] = [xc-xr,xc+xr] = xc+xr[-l, 1] = χ°+χΓΧξι 和 x(t) = [xL(t)，xu(t)] = [xc(t)-xr(t)，xc(t)+xr(t)] = xc(t)+xr(t)[-l，l] =xc(t)+xr(tm2 其中，Z 和χ^^χ/,χ/，· · ·，χ/)表示区间向量X的均值和半径，Xlt) = (Xic(t)，X2c(t)，· · ·，xnc(t))和Xc(t) = (Xir(t)，X2r(t)，· · ·，xnr(t))表示区间过程X(t)的 均值和半径山，ξ2εΞ4，Ξ4定义为所有元素包含在[-1，1]内的4维向量集合，符号"X"定义 为两个向量各对应元素相乘的算子，乘积仍为维数为4的向量。5. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第四步中δΗχ,Χα))的计算采用节点位移的差分迭代格式，针对确定性问题， δ^,χα))采用退化的S( tk)表示，计算如下：其中，tk表示第k个离散时刻，即tk = kA t，5(tk)表示tk时刻下的节点位移向量，δ(〇)， 3(0)，1(0)表示初始位移向量、速度向量和加速度向量，j和k分别表示迭代离散过程的计数 指标。6. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第五步中时变极限状态函数G(t，x，X(t))的定义及相关性特征的求解过程如 下： 所述第五步中时变极限状态函数G(t，x，X(t))同样被定义为一个区间过程，即G(t，x，X (〇)^1(〖^，(〇)，因此，对于任意给定瞬时时刻〖1 = 1八〖，极限状态时，时变极限状态函 数G(t，x，X(t))将退化为一个区间量6(〖1^3(〖 1))，有限个离散区间量的可行范围被约束 在一个超立方体域Ωη内；此外，要想获得协方差函数Covc(i A t，（i+l) △ t)和相关系数函数 pG(i Δ t，（i+1) Δ t)的显示表达，需借助标准化处理手段，并转换工作坐标系至(l^Us)，即： G(ti,x,X(ti)) e [GL(ti,x,X(ti)) ,Gu(ti,x,X(ti))] = Gc(ti,x,X(ti) )+Gr(ti,x,X(ti)) XUi G(ti+i,x,X(ti+i)) e [GL(ti+ii,x,X(ti+i)) ,Gu(ti+i,x,X(ti+i))] = Gc(ti+i,x,X(ti+i))+Gr(ti+i,x,X(ti+i)) XU2 其中，和表示U时刻极限状态函数的均值和半径， (ti+1，X，X( ti+1))和ti+1，X，X( ti+1))表示ti+1时刻极限状态函数的均值和半径。7. 根据权利要求1所述的一种复合材料层合结构的非概率动力可靠性评估方法，其特 征在于:所述第七步中非概率动力可靠度RS(T)的求解计算需借助时间离散化方法进行简 化近似，通过遍历每一个微小时间增量内结构发生穿越破坏的可能性指标，经叠加运算后 得到可靠度的等效表达式如下：其中，MT)表示整个生命周期内的失效度，/5m(o) = Pav(f:.L(.v, A(0)^/^(())1表示结 构在初始时刻即发生失效的可能度；△ t = ti+i-ti表示微小时间增量，其取值设定为整个服 役寿命的1/1000。
【文档编号】G06F19/00GK105956368SQ201610251078
【公开日】2016年9月21日
【申请日】2016年4月21日
【发明人】王磊, 王晓军, 陈潇, 王睿星, 陈贤佳, 郑宇宁, 邱志平
【申请人】北京航空航天大学