CN105893716A - 一种基于分形理论的结构断裂非概率可靠性分析方法 - Google Patents
一种基于分形理论的结构断裂非概率可靠性分析方法 Download PDFInfo
- Publication number
- CN105893716A CN105893716A CN201610390402.6A CN201610390402A CN105893716A CN 105893716 A CN105893716 A CN 105893716A CN 201610390402 A CN201610390402 A CN 201610390402A CN 105893716 A CN105893716 A CN 105893716A
- Authority
- CN
- China
- Prior art keywords
- alpha
- bounded
- unknown uncertainty
- vector
- interval
- 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
Links
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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于分形理论的结构断裂非概率可靠性分析方法。该方法首先考虑有限样本条件下载荷、材料特性、几何尺寸等参数的不确定性效应,将不确定性参数区间量化,并在不确定性参数的区间内配点;采用自相似网格和常规网格离散含裂纹结构,建立含裂纹结构的配点型区间应力强度因子分析模型,求解得到应力强度因子的区间范围;考虑断裂韧性的不确定性效应,建立非概率应力强度因子干涉模型,得到含裂纹结构的可靠度。本发明精确、高效的获得结构的可靠度,为结构的设计提供客观有效的数据。
Description
技术领域
本发明适用于结构断裂的可靠性分析,具体涉及一种基于分形理论和非概率集合理论的结构断裂可靠性分析方法。
背景技术
在工程实际中,机械设备和金属结构的构件中往往存在因制造、使用或材料本身缺陷所致的宏观裂纹。这时要确定构件能否继续安全使用,最重要的就是判断裂纹是否会失稳扩展从而导致结构和设备的破坏。应力强度因子反映了裂纹尖端附近区域的应力场和位移场,是裂纹扩展趋势和裂纹扩展推动力的度量。按断裂力学的观点:裂纹尖端的应力强度因子若是小于材料的断裂韧性,则构件是安全的,否则,构件是危险的。
大量的实际工程表明:由于作用在结构上的外载荷不确定地波动和组成结构的材质、结构工艺的内在不均匀性,使得虽然是同一批制造的同一类型结构在同一工况下体现出不同的效能,结果导致结构疲劳寿命可相差数倍之多,存在着相当大的分散性。同样就一种结构群的每个结构而言,其承受的载荷—时间历程、全寿命期内所受的最大载荷、决定临界裂纹尺寸的材料断裂韧度、描述裂纹扩展速率的曲线(表达式的参数)、乃至结构产生的裂纹形态,均是不确定的。由此可见,必须要利用可靠性理论对含不确定性参数结构的安全性进行研究。
当前,国内外学者与工程技术人员对含不确定性问题的方法大致分为三种:概率理论、模糊理论及非概率凸集合理论。同样,用于处理结构可靠性分析的模型也可分为:概率可靠性模型、模糊可靠性模型及非概率可靠性模型。目前,概率可靠性模型是最成功且应用最为普遍的可靠性模型;模糊可靠性模型也使得结构的设计及可靠性分析得到了很大的提高;非概率可靠性模型近年来也得到了很大的发展,对另外两种可靠性模型起到了一定的弥补作用。从应用条件来看,前两种可靠性模型以大量的统计数据为基础,以便得到不确定参数的概率密度函数或隶属度函数。事实上,由于试件精确数据的缺乏,概率密度函数或隶属度函数不能被精确的得到。人们仅能通过经验或较少的数据对它们进行一定的数据处理及假设。这样一些人为因素将带来模型的误差,从而影响结构可靠性的确定。因此,这两种模型对结构可靠性问题的分析有一定的局限性。非概率可靠性模型所需信息较少,因此对其进行深入研究有很大的科研及实际意义。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供一种基于分形理论和非概率集合理论的结构断裂可靠性分析方法,充分考虑实际工程问题中普遍存在的不确定性因素,结合分形理论与非概率集合理论,所得到的设计结果更加符合真实情况,工程适用性更强。
本发明采用的技术方案为:一种基于分形理论的结构断裂非概率可靠性分析方法,具体实现步骤如下:
第一步:区间不确定性参数向量是以区间形式表示的不确定性参数向量,利用区间不确定性参数向量表征贫信息、少数据条件下的结构参数和载荷的不确定性,表示为:
i=1,2,…,m
其中,和分别为区间不确定性参数向量α的上、下界,和i=1,2,…,m分别为第i个区间不确定性参数的上、下界,m为区间不确定性参数的个数,为区间不确定性参数向量α的中心值,为区间不确定性参数向量α的半径,和△αi,i=1,2,…,m分别为第i个区间不确定性参数的中心值和半径;
区间不确定性参数向量还可表示为:
其中,e∈Ξm,Ξm定义为所有元素包含在[-1,1]内的m维向量集合,符号“×”定义为两个向量各对应元素相乘的算子,乘积仍为维数为m的向量。
第二步:将第一步中的区间不确定性参数向量处理为一元区间不确定性参数向量,一个m维的区间不确定性参数向量变为m个一元区间不确定性参数向量,m维的区间不确定性参数向量中每一维都是区间不确定性参数,一元区间不确定性参数向量中只有其中一维是区间不确定性参数,其它m-1维是确定性参数;
一元区间不确定性参数向量表示为:
αi=αc+△α×Xi
其中,Xi=(0,…,x,…,0)T,x处于第i行;αi为一元区间不确定性参数向量,角标i表示αi中第i个分量为区间不确定性参数。由此可见,一个m维的区间不确定性参数向量通过处理变为m个一元区间不确定性参数向量。
第三步:在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集。配点原则是采用高斯积分点在区间内配点,区间内的Gauss积分点记为xk,表示为:
其中,xk为区间内配置的第k个高斯积分点,q为区间内配点的个数。
第四步:根据含裂纹结构的几何模型及裂纹位置,用人工边界Γ,将几何模型划分为常规区域Ω和靠近裂纹尖端的分形区域D,其中人工边界Γ是圆形边界,圆心在裂纹尖端端点,半径是r,0≤r≤a,a为裂纹长度;
根据分形理论的自相似性,在分形区域D内构造比例系数为ξ的自相似单元。自相似单元的层数为k,k为大于等于1的正整数,比例系数0<ξ<1;
含裂纹结构的配点型区间应力强度因子分析模型表示为:
K(α)u=f(α)
其中,K(α)为含裂纹结构的区间不确定性结构刚度矩阵,f(α)为含裂纹结构的区间不确定性结构节点载荷向量,u为含裂纹结构的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数;
K(α)、u、f(α)分别表示为:
其中,ur分别为区域Ω内节点的区间不确定性结构刚度矩阵、节点载荷向量和节点位移向量,um分别为边界Γ上主节点的区间不确定性结构刚度矩阵、节点载荷向量和节点位移向量,和为区域Ω和边界Γ上的区间不确定性结构耦合刚度矩阵。分别为区域D第1层单元主节点的区间不确定性结构刚度矩阵、节点载荷向量,分别为区域D第1层单元从节点的区间不确定性结构刚度矩阵、节点载荷向量,a为分形区域D内的区间不确定性广义坐标向量,和为区域D第1层单元的区间不确定性结构耦合刚度矩阵,为区域D第1层单元从节点的转换矩阵,和分别为分形区域D内的第2至k层单元的有界不确定性结构刚度矩阵、节点载荷向量。通过求解u的上下界可以直接得到应力强度因子KΙ,ΙΙ的上界和下界,其中和分别为Ι型平面裂纹的应力强度因子的上界和下界,和分别为ΙΙ型平面裂纹的应力强度因子的上界和下界。
第五步:根据第三步得到的区间配点方案和第四步建立的配点型区间应力强度因子分析模型,求解应力强度因子的上界和下界;
具体求解时,采用最佳平方逼近多项式逼近结构响应函数,表示为:
其中,Tj(x)=cos(jarccosx),-1≤x≤1,0≤j≤r,为正交多项式系;Tj(xk)为正交多项式系Tj(x)在第k个高斯积分点处对应的函数值;为第k个高斯积分点对应的结构响应;为r阶第一类Chebyshev多项式,角标i表示针对第i个一元区间不确定性参数向量,采用最佳平方逼近多项式逼近结构响应函数;q为配点个数;
求解x∈[-1,1]的最小值点和最大值点,分别记为和重复上述过程,直到i遍历完1~m时,就能得到具有m个元素的最值点向量,记为和将Xmin和Xmax分别带入结构响应函数中,得到结构响应的近似区间估计为按下式计算:
其中,为响应的下界,为响应的上界;和分别为第i个一元区间不确定性参数向量在区间[-1,1]内的最小值和最大值点,Xmin和Xmax分别是由最小值点和最大值点构成的最值点向量;αc为区间不确定性参数向量的中值向量;△α为区间不确定性参数向量的区间半径向量;响应的中值由下式给出
第六步:根据断裂韧性的不确定性效应和结构断裂可靠性的功能函数,建立非概率应力强度因子干涉模型,基于断裂准则和非概率应力强度因子干涉模型,对结构断裂非概率可靠性进行度量,得到结构断裂的可靠度;
结构断裂可靠性的功能函数表示为:
M(KΙ,ΙΙ,Kc)=Kc-KΙ,ΙΙ
其中,应力强度因子KΙ,ΙΙ和断裂韧性Kc均为区间变量,即和和分别为应力强度因子的下界和上界,和分别为断裂韧性的下界和上界。当M(KΙ,ΙΙ,Kc)>0时,结构安全,裂纹稳定不扩展;当M(KΙ,ΙΙ,Kc)<0时,结构失效,裂纹将发生不稳定扩展;当M(KΙ,ΙΙ,Kc)=0时,为临界状态。
本发明与现有技术相比的优点在于:
(1)本发明提出了一种基于分形理论的区间应力强度因子求解方法,该方法根据分形理论的自相似性,在裂纹尖端形成无限细化的网格,可无限逼近裂纹尖端,提高了计算的精确性,克服现有方法精度不高的缺点。
(2)本发明可以处理贫数据、少信息条件下的结构可靠性问题,无需知道不确定性参数的概率分布,只要知道不确定性参数的上下界就可以预测结构的可靠度,有更强的工程适用性。
附图说明
图1是本发明方法实现流程图;
图2是本发明中含裂纹结构的分区示意图;
图3是本发明中含裂纹结构常规区域的离散示意图;
图4是本发明中含裂纹结构分形区域的离散示意图;
图5是本发明中非概率应力强度因子干涉模型示意图;
图6是本发明二维区间变量的断裂韧性和应力强度因子干涉模型示意图;
图7是本发明实施例中带一条单边裂纹的弹性板几何模型示意图;
图8是本发明实施例中不同变异系数下弹性板的非概率可靠度。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,本发明提出了一种基于分形理论的结构断裂非概率可靠性分析方法,其具体实现步骤是:
(1)区间不确定性参数向量是以区间形式表示的不确定性参数向量,利用区间不确定性参数向量表征贫信息、少数据条件下的结构参数和载荷的不确定性,表示为:
i=1,2,…,m
其中,和分别为区间不确定性参数向量α的上、下界,和,i=1,2,…,m分别为第i个区间不确定性参数的上、下界,m为区间不确定性参数的个数,为区间不确定性参数向量α的中心值,为区间不确定性参数向量α的半径,和△αi,i=1,2,…,m分别为第i个区间不确定性参数的中心值和半径;
区间不确定性参数向量还可表示为:
其中,e∈Ξm,Ξm定义为所有元素包含在[-1,1]内的m维向量集合,符号“×”定义为两个向量各对应元素相乘的算子,乘积仍为维数为m的向量。
(2)将(1)中的区间不确定性参数向量处理为一元区间不确定性参数向量。取(2)式中e的第i(1≤i≤m)个元素为x,其他元素为0,记为:
Xi=(0,…,x,…,0)T (3)
1 i m
其中,x∈[-1,1]。由式(2)和(3)得到一元区间不确定性参数向量为:
αi=αc+△α×Xi (4)
其中,αi为一元区间不确定性参数向量,角标i表示αi中第i个分量为区间不确定性参数。由此可见,一个m维的区间不确定性参数向量通过处理变成m个一元区间不确定性参数向量。m维的区间不确定性参数向量中每一维都是区间不确定性参数,一元区间不确定性参数向量中只有其中一维是区间不确定性参数,其它m-1维是确定性参数。
(3)在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集,配点原则是采用高斯积分点在区间内配点,在[-1,1]上配置q个Gauss积分点记为xk,表示为:
其中,xk为区间内配置的第k个高斯积分点,q为区间内配点的个数。
(4)根据含裂纹结构的几何模型及裂纹位置,用人工边界Γ,将几何模型划分为常规区域Ω和靠近裂纹尖端的分形区域D,见图2。对常规区域Ω的几何模型进行离散,见图3(仅表示了一半模型)。建立常规区域Ω内的区间不确定性结构响应求解模型,表示为:
KR(α)uR=fR(α) (5)
其中,KR(α)为常规区域Ω内的区间不确定性结构刚度矩阵,表示为:
其中,为区域Ω内节点的区间不确定性结构刚度矩阵,为边界Γ上主节点的区间不确定性结构刚度矩阵,和为区域Ω和边界Γ上的区间不确定性结构耦合刚度矩阵;
fR(α)为常规区域Ω内的区间不确定性结构节点载荷向量,表示为:
其中,为区域Ω内节点的区间不确定性结构节点载荷向量,为边界Γ上主节点的区间不确定性结构节点载荷向量;
uR为区域Ω内的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数。表示为:
其中,ur为区域Ω内节点的区间不确定性结构节点位移向量,um为边界Γ上主节点的区间不确定性结构节点位移向量;
在分形区域D内,采用比例系数为ξ的自相似网格进行离散,建立k层自相似单元。其中,0<ξ<1,k为大于等于1的正整数。如图4所示;
以裂纹尖端端点为原点构造极坐标系,裂纹尖端位移场的William’s一般解具体表示为:
其中,u,v分别为直角坐标系下裂纹尖端沿x和y方向的位移分量,G为剪切模量,r为极坐标系下节点的极径,n为William’s级数的项数。和n=1,2,…,为广义坐标,fn,ij(n,θ),i,j=1,2,具体表达式为:
其中,θ为极坐标系下节点的极角,κ为常数,对于平面应变问题,κ=3-4ν,平面应力问题,κ=(3-ν)/(1+ν),其中ν为泊松比;
利用裂纹尖端位移场的William’s一般解作为整体插值函数,将分形区域D内的区间不确定性结构节点位移向量表示为:
uS=TSa (15)
其中,uS为分形区域D内的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数,TS为分形区域D内的转换矩阵,a为分形区域D内的区间不确定性广义坐标向量,也是区间不确定性参数向量α=(αi)的函数,均为分形区域D内的区间不确定性广义坐标;
根据(9)、(10)式,区间不确定性平面裂纹应力强度因子KΙ,ΙΙ与区间不确定性广义坐标有关,表示为:
其中,KI为Ι型平面裂纹的区间不确定性应力强度因子,KΙΙ为ΙΙ型平面裂纹的区间不确定性应力强度因子,二者均是区间不确定性参数向量α=(αi)的函数。由此可见,区间不确定性应力强度因子可由区间不确定性广义坐标得到;
在区间不确定性参数约束的条件下,建立分形区域D内第1层单元的区间不确定性结构响应求解模型,表示为:
其中,分别为区域D第1层单元主节点的区间不确定性结构刚度矩阵、节点载荷向量,um为边界Γ上主节点的区间不确定性节点位移向量,分别为区域D第1层单元从节点的区间不确定性结构刚度矩阵、节点载荷向量和节点位移向量,和为区域D第1层单元的区间不确定性结构耦合刚度矩阵;
利用分形区域D内的转换矩阵TS和区间不确定性广义坐标向量a,将表示为:
其中,为区域D第1层单元从节点的转换矩阵。进而,将分形区域D内第1层单元的区间不确定性结构响应求解模型表示为:
或
其中,I为单位矩阵;
分形区域D内的第k1层单元,其中2≤k1≤k,k为分形区域D内自相似单元的总层数,其区间不确定性结构响应求解模型表示为:
其中,为分形区域D内第k1层单元的区间不确定性结构刚度矩阵,为分形区域D内第k1层单元的区间不确定性结构节点载荷向量,为分形区域D内第k1层单元的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数;
将上式表示成广义坐标的形式:
其中,为分形区域D内第k1层单元的转换矩阵。根据分形区域D内单元的自相似性,每一层单元的刚度矩阵相等,即:
其中,为分形区域D内第2层单元的区间不确定性结构刚度矩阵;
根据自相似性,将(23)式中的表示为:
其中,为分形区域D内第2层单元的转换矩阵。Diag[ηi]为对角元素为ηi的对角矩阵,ηi具体表示为:
其中,ξ为比例系数,ni具体表示为:
其中,1≤i≤2n,n为William’s级数的项数。结合(22)、(23)、(24)、(25)四式,将分形区域D内的第2至k层的区间不确定性结构刚度矩阵叠加,表示为:
其中,1≤j≤2n,n为William’s级数的项数,为分形区域D内的第2至k层的区间不确定性结构刚度矩阵,αij和[kij(α)]具体表示为:
同样,分形区域D内第2至k层的区间不确定性结构节点载荷向量表示为:
其中,为分形区域D内第2至k层的区间不确定性结构节点载荷向量,[fij(α)]具体表示为:
其中,为分形区域D内第2层单元的区间不确定性结构节点载荷向量;
叠加式(21)、(28)、(31)建立整个分形区域D内的区间不确定性结构响应求解模型,表示为:
KS(α)uS=fS(α) (33)
其中,KS(α)为分形区域D内的区间不确定性结构刚度矩阵,fS(α)为分形区域D内的区间不确定性结构节点载荷向量,uS为分形区域D内的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数;
KS(α)、uS、fS(α)分别表示为:
组合(5)和(33)的区间不确定性结构响应求解模型,建立含裂纹结构的区间不确定性结构响应求解模型,表示为:
K(α)u=f(α) (37)
其中,K(α)为含裂纹结构的区间不确定性结构刚度矩阵,f(α)为含裂纹结构的区间不确定性结构节点载荷向量,u为含裂纹结构的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数;
K(α)、u、f(α)分别表示为:
(5)根据第(3)步得到的区间配点方案采用最佳平方逼近多项式逼近结构响应函数的方法求解(37)式,得到区间不确定性节点位移向量u的上界和下界;
引入r阶第一类Chebyshev多项式,其正交多项式系{Tj(x)}和最佳平方逼近函数Pr(x)为:
Tj(x)=cos(jarccosx),-1≤x≤1,0≤j≤r (41)
其中,j为非负整数,aj为逼近函数展开式系数。Pr(x)为r阶第一类Chebyshev多项式;
由Gauss积分点求出多项式系数并代入(42),可进一步得到:
其中,Tj(xk)为正交多项式系Tj(x)在第k个高斯积分点处对应的函数值;为第k个高斯积分点对应的结构响应;为r阶第一类Chebyshev多项式,角标i表示针对第i个一元区间不确定性参数向量,采用最佳平方逼近多项式逼近结构响应函数;q为配点个数;
简记为:
其中:
T(x)=[1 T1(x) T2(x) … Tr(x)]T (47)
先考虑如何求解的最值,对式(44)关于x求导并令导数为零,得:
求解式(48)的根,并联合和根据连续函数在闭区间上的最值定理,可得一元逼近函数的最小值点和最大值点,分别记为和
重复上述过程,直到i遍历完1~m时,就能得到具有m个元素的最值点向量,记为和将Xmin和Xmax分别带入结构响应函数中,得到结构响应的近似区间估计为按下式计算
其中,和分别为第i个一元区间不确定性参数向量在区间[-1,1]内的最小值和最大值点,Xmin和Xmax分别是由最小值点和最大值点构成的最值点向量;αc为区间不确定性参数向量的中值向量;△α为区间不确定性参数向量的区间半径向量;响应的中值由下式给出根据(49)和(50)式得到的区间不确定性节点位移向量的上界和下界,提取相应的区间不确定性广义坐标的上界和下界根据(17)式得到区间不确定性应力强度因子KΙ,ΙΙ的上界和下界。其中KΙ为Ι型平面裂纹的区间不确定性应力强度因子,KΙΙ为ΙΙ型平面裂纹的区间不确定性应力强度因子。应力强度因子的中值由给出。
(6)考虑断裂韧性的不确定性效应,将其表示为:
其中和分别表示断裂韧性的下界和上界。其中值表示为
将结构断裂非概率可靠性的功能函数表示为:
M(KΙ,ΙΙ,Kc)=Kc-KΙ,ΙΙ (52)
当M(KΙ,ΙΙ,Kc)>0时,结构安全,裂纹稳定不扩展;当M(KΙ,ΙΙ,Kc)<0时,结构失效,裂纹将发生不稳定扩展;当M(KΙ,ΙΙ,Kc)=0时,为临界状态;
由于断裂韧性和应力强度因子都是区间变量,根据(51)和(52)式,二者可能发生干涉的情况,如图5所示,即为非概率应力强度因子干涉模型;
基于断裂准则和非概率应力强度因子干涉模型,对结构断裂非概率可靠性进行度量,将图5所示的干涉关系转换为二维区间变量的断裂韧性和应力强度因子干涉的关系,如图6所示。则结构断裂的非概率可靠度即为安全区域面积与变量总区域面积之比:
其中,R为结构断裂非概率可靠度,Ssafe为安全区域的面积,Ssum为变量区域的总面积。
实施例:
为了更充分地了解该发明的特点及其对工程实际的适用性,本发明针对如图7所示的带一条单边裂纹的弹性板进行非概率可靠性分析。图7中弹性板的宽为w,高h=200cm,裂纹长度为a,弹性模量E=2×105MPa,泊松比ν=0.167,受到均布拉力F的作用。由于制造和测量误差,裂纹长度a、均布拉力F和弹性板宽w均为区间不确定性参数,裂纹长度a的中心值为ac=5cm,均布拉力F的中心值为Fc=0.3kN/cm,弹性板宽w的中心值为wc=40cm,并且有a=[ac-βac,ac+βac],F=[Fc-βFc,Fc+βFc],w=[wc-βwc,wc+βwc],β为可以变化的变异系数,分别取0.05,0.10,0.15,0.20,0.25,0.30。断裂韧性的区间范围为:Kc=[1.5,2.7]。本例中需要预测弹性板断裂的非概率可靠度。
本例中裂纹类型为I型,因此,KΙΙ=0。用一个圆心在裂纹尖端端点,半径为r=3cm的圆将弹性板分为常规区域Ω和分形区域D。由于弹性板是对称结构,故取一半模型进行分析。对常规区域Ω采用四节点的四边形等参单元离散,共有32个单元,47个节点。对分形区域D采用比例系数为ξ=0.5的自相似单元离散,建立k=10层自相似单元。裂纹尖端位移场William’s一般解取10项,即n=10。通过编程得到分形区域D内的有界不确定性结构静力响应求解模型。最佳平方逼近多项式的项数r=5,配点个数取为q=5。表1给出了不同变异系数下,根据本发明所提出的方法得到的弹性板的非概率可靠度。
表1
β | R(%) |
0.05 | 98.62 |
0.10 | 95.36 |
0.15 | 91.38 |
0.20 | 86.58 |
0.25 | 80.76 |
0.30 | 73.57 |
图8显示了本文方法得到的结构断裂非概率可靠度。从图中可以看出,随着变异系数β的增大,结构的可靠度逐渐降低,也就是说,随着不确定性参数区间变大,结构的可靠度逐渐降低。这与工程实际情况吻合,由此证明本发明方法所得的结果是可信的。此外,工程实际中很难获知载荷、结构参数等区间不确定性参数的真实概率分布情况,只能得到其分布的上下界,此时,本发明方法的易用性和有效性优势更加得以凸显。此外,本发明提出的一种基于分形理论求解区间应力强度因子的半解析方法,计算精度高,弥补了现有方法不精确的缺点,且通过对实施例的计算发现,本发明方法与传统有限元方法相比计算效率提高了90%,存储容量减少了60%。以上实例验证了本发明方法针对结构断裂可靠度分析的可行性与精确性。
本发明可以精确、高效的获得结构断裂的可靠度,为结构后续设计提供客观有效的数据。
以上仅是本发明的具体步骤,对本发明的保护范围不构成任何限制。
本发明未详细阐述部分属于本领域技术人员的公知技术。
Claims (7)
1.一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于实现步骤如下:
第一步:区间不确定性参数向量是以区间形式表示的不确定性参数向量,利用区间不确定性参数向量表征贫信息、少数据条件下的结构参数和载荷的不确定性,表示为:
其中,和分别为区间不确定性参数向量α的上、下界,和α i,i=1,2,…,m分别为第i个区间不确定性参数的上、下界,m为区间不确定性参数的个数,为区间不确定性参数向量α的中心值,为区间不确定性参数向量α的半径,和Δαi,i=1,2,…,m分别为第i个区间不确定性参数的中心值和半径;
第二步:将第一步中的区间不确定性参数向量处理为一元区间不确定性参数向量,一个m维的区间不确定性参数向量变为m个一元区间不确定性参数向量,m维的区间不确定性参数向量中每一维都是区间不确定性参数,一元区间不确定性参数向量中只有其中一维是区间不确定性参数,其它m-1维是确定性参数;
第三步:在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集;
第四步:根据含裂纹结构的几何模型及裂纹位置,用人工边界Γ,将几何模型划分为常规区域Ω和靠近裂纹尖端的分形区域D,基于分形理论的自相似性,在分形区域D内构造比例系数为ξ的自相似单元,建立含裂纹结构的配点型区间应力强度因子分析模型;
第五步:根据第三步得到的区间配点方案和第四步建立的配点型区间应力强度因子分析模型,求解应力强度因子的上界和下界K I,II,其中和K I分别为I型平面裂纹的应力强度因子的上界和下界,和K II分别为II型平面裂纹的应力强度因子的上界和下界;
第六步:根据断裂韧性的不确定性效应和结构断裂可靠性的功能函数,建立非概率应力强度因子干涉模型,基于断裂准则和非概率应力强度因子干涉模型,对结构断裂非概率可靠性进行度量,得到结构断裂的可靠度。
2.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第一步中区间不确定性参数向量还可表示为:
其中,e∈Ξm,Ξm定义为所有元素包含在[-1,1]内的m维向量集合,符号“×”定义为两个向量各对应元素相乘的算子,乘积仍为维数为m的向量。
3.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第二步中一元区间不确定性参数向量表示为:
αi=αc+Δα×Xi
其中,Xi=(0,…,x,…,0)T,x处于第i行;αi为一元区间不确定性参数向量,角标i表示α中第i个分量为区间不确定性参数,由此可见,一个m维的区间不确定性参数向量通过处理变为m个一元区间不确定性参数向量。
4.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第三步中在一元区间不确定性参数向量的区间内配点的原则是采用高斯积分点在区间内配点,区间内的Gauss积分点记为xk,表示为:
其中,xk为区间内配置的第k个高斯积分点,q为区间内配点的个数。
5.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第四步中人工边界Γ是圆形边界,圆心在裂纹尖端端点,半径是r,其中0≤r≤a,a为裂纹长度,分形区域D内自相似单元的层数为k,k为大于等于1的正整数,比例系数0<ξ<1;
含裂纹结构的配点型区间应力强度因子分析模型表示为:
K(α)u=f(α)
其中,K(α)为含裂纹结构的区间不确定性结构刚度矩阵,f(α)为含裂纹结构的区间不确定性结构节点载荷向量,u为含裂纹结构的区间不确定性结构节点位移向量,也是区间不确定性参数向量α=(αi)的函数;
K(α)、u、f(α)分别表示为:
其中,ur分别为区域Ω内节点的区间不确定性结构刚度矩阵、节点载荷向量和节点位移向量,um分别为边界Γ上主节点的区间不确定性结构刚度矩阵、节点载荷向量和节点位移向量,和为区域Ω和边界Γ上的区间不确定性结构耦合刚度矩阵,分别为区域D第1层单元主节点的区间不确定性结构刚度矩阵、节点载荷向量,分别为区域D第1层单元从节点的区间不确定性结构刚度矩阵、节点载荷向量,a为分形区域D内的区间不确定性广义坐标向量,和为区域D第1层单元的区间不确定性结构耦合刚度矩阵,为区域D第1层单元从节点的转换矩阵,和分别为分形区域D内的第2至k层单元的有界不确定性结构刚度矩阵、节点载荷向量,通过求解u的上下界可以直接得到应力强度因子KI,II的上界和下界K I,II,其中和K I分别为I型平面裂纹的应力强度因子的上界和下界,和K II分别为II型平面裂纹的应力强度因子的上界和下界。
6.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第五步中对配点型区间应力强度因子分析模型求解时,采用最佳平方逼近多项式逼近结构响应函数,表示为:
其中,Tj(x)=cos(j arccos x),-1≤x≤1,0≤j≤r,为正交多项式系;Tj(xk)为正交多项式系Tj(x)在第k个高斯积分点处对应的函数值;为第k个高斯积分点对应的结构响应;为r阶第一类Chebyshev多项式,角标i表示针对第i个一元区间不确定性参数向量,采用最佳平方逼近多项式逼近结构响应函数;q为配点个数;
求解x∈[-1,1]的最小值点和最大值点,分别记为和重复上述过程,直到i遍历完1~m时,就能得到具有m个元素的最值点向量,记为和将Xmin和Xmax分别带入结构响应函数中,得到结构响应的近似区间估计为按下式计算:
u=umin=u(αc+Δα×Xmin)
其中,u为响应的下界,为响应的上界;和分别为第i个一元区间不确定性参数向量在区间[-1,1]内的最小值和最大值点,Xmin和Xmax分别是由最小值点和最大值点构成的最值点向量;αc为区间不确定性参数向量的中值向量;Δα为区间不确定性参数向量的区间半径向量;响应的中值由下式给出
7.根据权利要求1所述的一种基于分形理论的结构断裂非概率可靠性分析方法,其特征在于:所述第六步中结构断裂非概率可靠性的功能函数表示为:
M(KI,II,Kc)=Kc-KI,II
其中,应力强度因子KI,II和断裂韧性Kc均为区间变量,即和 K I,II和分别为应力强度因子的下界和上界,K c和分别为断裂韧性的下界和上界,当M(KI,II,Kc)>0时,结构安全,裂纹稳定不扩展;当M(KI,II,Kc)<0时,结构失效,裂纹将发生不稳定扩展;当M(KI,II,Kc)=0时,为临界状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390402.6A CN105893716B (zh) | 2016-06-02 | 2016-06-02 | 一种基于分形理论的结构断裂非概率可靠性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390402.6A CN105893716B (zh) | 2016-06-02 | 2016-06-02 | 一种基于分形理论的结构断裂非概率可靠性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105893716A true CN105893716A (zh) | 2016-08-24 |
CN105893716B CN105893716B (zh) | 2018-06-29 |
Family
ID=56710925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610390402.6A Expired - Fee Related CN105893716B (zh) | 2016-06-02 | 2016-06-02 | 一种基于分形理论的结构断裂非概率可靠性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105893716B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844945A (zh) * | 2017-01-19 | 2017-06-13 | 电子科技大学 | 一种考虑认知不确定性的多状态系统概率重要度分析方法 |
CN108846149A (zh) * | 2018-04-20 | 2018-11-20 | 北京航空航天大学 | 一种基于多源不确定性的结构分布式动态载荷识别的方法 |
CN108920786A (zh) * | 2018-06-20 | 2018-11-30 | 北京航空航天大学 | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 |
CN109374449A (zh) * | 2018-09-25 | 2019-02-22 | 南京航空航天大学 | 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法 |
CN110795820A (zh) * | 2019-09-24 | 2020-02-14 | 武汉大学 | 工程结构裂纹问题求解方法以及装置 |
CN112257197A (zh) * | 2020-10-19 | 2021-01-22 | 郑州轻工业大学 | 一种大型铸锻件微缺陷工作应力评估方法 |
CN115146482A (zh) * | 2022-08-01 | 2022-10-04 | 暨南大学 | 基于子区间分析与Chebyshev稀疏模型的结构可靠性分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110159248A1 (en) * | 2009-12-30 | 2011-06-30 | Airbus Operations, S.L. | Aircraft fuselage frame of composite material with stabilising ribs |
CN103324782A (zh) * | 2013-05-29 | 2013-09-25 | 北京航空航天大学 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
CN103366065A (zh) * | 2013-07-17 | 2013-10-23 | 北京航空航天大学 | 一种基于区间可靠性的飞行器热防护系统尺寸优化设计方法 |
CN104459776A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油大学(华东) | 一种断裂分形特征优化计算方法 |
-
2016
- 2016-06-02 CN CN201610390402.6A patent/CN105893716B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110159248A1 (en) * | 2009-12-30 | 2011-06-30 | Airbus Operations, S.L. | Aircraft fuselage frame of composite material with stabilising ribs |
CN103324782A (zh) * | 2013-05-29 | 2013-09-25 | 北京航空航天大学 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
CN103366065A (zh) * | 2013-07-17 | 2013-10-23 | 北京航空航天大学 | 一种基于区间可靠性的飞行器热防护系统尺寸优化设计方法 |
CN104459776A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油大学(华东) | 一种断裂分形特征优化计算方法 |
Non-Patent Citations (3)
Title |
---|
孙文彩: "结构非概率可靠性模型在断裂力学分析中的应用研究", 《应用基础与工程科学学报》 * |
张衡: "分形裂纹的断裂参数", 《华南理工大学学报(自然科学版)》 * |
胡春明: "断裂分形理论的应用", 《大庆石油地质与开发》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844945A (zh) * | 2017-01-19 | 2017-06-13 | 电子科技大学 | 一种考虑认知不确定性的多状态系统概率重要度分析方法 |
CN108846149A (zh) * | 2018-04-20 | 2018-11-20 | 北京航空航天大学 | 一种基于多源不确定性的结构分布式动态载荷识别的方法 |
CN108846149B (zh) * | 2018-04-20 | 2020-02-07 | 北京航空航天大学 | 一种基于多源不确定性的结构分布式动态载荷识别的方法 |
CN108920786A (zh) * | 2018-06-20 | 2018-11-30 | 北京航空航天大学 | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 |
CN108920786B (zh) * | 2018-06-20 | 2022-10-25 | 北京航空航天大学 | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 |
CN109374449A (zh) * | 2018-09-25 | 2019-02-22 | 南京航空航天大学 | 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法 |
CN110795820A (zh) * | 2019-09-24 | 2020-02-14 | 武汉大学 | 工程结构裂纹问题求解方法以及装置 |
CN112257197A (zh) * | 2020-10-19 | 2021-01-22 | 郑州轻工业大学 | 一种大型铸锻件微缺陷工作应力评估方法 |
CN112257197B (zh) * | 2020-10-19 | 2022-06-28 | 郑州轻工业大学 | 一种大型铸锻件微缺陷工作应力评估方法 |
CN115146482A (zh) * | 2022-08-01 | 2022-10-04 | 暨南大学 | 基于子区间分析与Chebyshev稀疏模型的结构可靠性分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105893716B (zh) | 2018-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105893716A (zh) | 一种基于分形理论的结构断裂非概率可靠性分析方法 | |
Liu et al. | A systematic approach for dynamic security assessment and the corresponding preventive control scheme based on decision trees | |
CN107133406B (zh) | 一种电力系统静态电压稳定域边界的快速搜索方法 | |
CN101635457B (zh) | 一种基于状态估计残差参数灵敏度的电网参数估计方法 | |
CN105808884A (zh) | 一种基于分形理论的有界不确定性平面裂纹应力强度因子上下界的预测方法 | |
CN105976257A (zh) | 基于隶属度函数的模糊综合评价法的电网脆弱性评估方法 | |
CN105471637A (zh) | 一种复杂网络节点重要性评估方法及系统 | |
CN108051035A (zh) | 基于门控循环单元的神经网络模型的管网漏损识别方法 | |
CN107632590B (zh) | 一种基于优先级的底事件排序方法 | |
CN107092984A (zh) | 一种基于级联失效的网络功能端节点传播预测方法 | |
CN105956336A (zh) | 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法 | |
CN105023071A (zh) | 一种基于高斯云变换和模糊时间序列的水质预测方法 | |
CN105843733A (zh) | 一种大数据平台的性能检测方法及装置 | |
CN106651660A (zh) | 基于g1‑熵权法寻找电网静态薄弱点的综合评价方法 | |
CN104898089A (zh) | 一种基于空间迁移压缩感知的被动式定位方法 | |
Cui et al. | Fuel savings model after aero-engine washing based on convolutional neural network prediction | |
Bahmanikashkooli et al. | Application of particle swarm optimization algorithm for computing critical depth of horseshoe cross section tunnel | |
CN104462414A (zh) | 一种基于拓扑结构的流程图相似性方法 | |
CN104281978A (zh) | 一种基于概率潮流的可用输电能力计算方法 | |
CN107276093B (zh) | 基于场景削减的电力系统概率潮流计算方法 | |
CN110571788A (zh) | 基于动态等值电路的静态电压稳定域边界系数的计算方法 | |
CN107798353A (zh) | 一种间歇过程监控数据处理方法 | |
CN106874611B (zh) | 一种基于超体积迭代策略的含区间参数结构响应区间的分析方法 | |
Li et al. | On the matching number of an uncertain graph | |
CN106874561A (zh) | 基于牛顿迭代的多学科不确定性传播分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180629 |
|
CF01 | Termination of patent right due to non-payment of annual fee |