一种基于自适应配点的结构模糊不确定性分析方法技术

本发明专利技术公开了一种基于自适应配点的结构模糊不确定性分析方法，属于结构不确定性分析领域。首先，得到隶属度轴上的区间；其次，通过模糊截集策略得到模糊不确定变量的截集区间；然后，通过区间不确定性分析方法得到截集区间下的结构响应上下界；最后，通过两点梯形公式和三点梯形公式计算结构响应的上下界与隶属度轴所围成的面积并判断是否收敛，如果收敛，则输出结构响应的模糊分布，否则，进行下一次循环直至收敛。数值算例表明，基于自适应配点的结构模糊不确定性分析方法能够得到精确的结构响应的模糊分布，为结构模糊不确定性分析提供了一种新的方法。

A Fuzzy Uncertainty Analysis Method for Structures Based on Adaptive Collocation

The invention discloses a structural fuzzy uncertainty analysis method based on adaptive collocation points, which belongs to the field of structural uncertainty analysis. Firstly, the interval on the axis of membership degree is obtained; secondly, the interval of cut sets of fuzzy uncertain variables is obtained by using the strategy of fuzzy cut sets; secondly, the upper and lower bounds of structural response under the interval uncertainty analysis method are obtained; finally, the upper and lower bounds of structural response are calculated by two-point trapezoidal formula and three-point trapezoidal formula. The area surrounded by the membership axis is judged to be convergent. If convergence occurs, the fuzzy distribution of the output structure response will be obtained. Otherwise, the next cycle will be carried out until convergence. The numerical examples show that the method of structural fuzzy uncertainty analysis based on adaptive collocation can obtain the accurate fuzzy distribution of structural response, which provides a new method for structural fuzzy uncertainty analysis.

【技术实现步骤摘要】
一种基于自适应配点的结构模糊不确定性分析方法
本专利技术涉及结构模糊不确定性分析的
，特别涉及基于自适应配点的结构模糊不确定性分析方法。
技术介绍
工程中存在多种多样的不确定性，为了避免由于这些不确定性因素导致的结构失效，许多结构不确定性分析方法被提出。总体上，三种不确定性模型被用来处理这些不确定性，即概率模型、区间模型和模糊模型。在概率模型中，不确定变量被看做随机变量，完善的概率理论给随机不确定性分析提供了良好的数学基础，然而概率模型需要大量的样本信息来构造不确定变量精确的概率密度分布，然而工程中往往得不到充足的样本。因此，发展了非概率模型。区间模型是一种常用的非概率模型，在区间模型中，只需要确定不确定变量的上下界。常用的区间不确定分析方法有区间摄动法、顶点法、蒙特卡洛法、优化法。模糊模型是近年来发展的又一种非概率模型。随着对不确定性因素的进一步认识，人们发现由于系统本身的复杂性以人们现有知识或技能的局限性所产生的一些不确定性因素中，有些只能根据已有实验和专家经估计出这些参数在某值附近这样模糊的信息。这类不确定性广泛地存在于工程对象的几何特征、材料特性、荷载及边界条件等方面，将其按照模糊因素处理更为合适。通过截集策略，模糊模型可以转化为一系列的区间。常见的模糊不确定分析方法基于标准模糊算法，其将隶属度函数等分为m段。在工程实际中，m的取值通常很小，以保证计算效率，然后这样往往会导致计算所得的结果不精确，特别是隶属度函数为非线性时。基于此，本专利技术提出一种基于自适应配点的结构模糊不确定性分析方法以改善传统方法的精度。
技术实现思路
本专利技术要解决的技术问题是：克服现有方法的不足，提供一种基于自适应配点的结构模糊不确定性分析方法。基于自适应配点的结构模糊不确定性分析方法计算精度较高，是现有方法很好的一个补充。本专利技术采用的技术方案为：基于自适应配点的结构模糊不确定性分析方法，可用于工程结构的模糊不确定性分析，其实现步骤如下：步骤一：根据所作的现实的抽象、假设情况，知识的缺乏情况，几何尺寸和加载条件，材料特性确定输入参数的模糊分布，利用合理表征贫信息、少数据条件下的模糊不确定性参数的集合；步骤二：设置隶属度轴上的区间[a,b]，a,b的初始值分别为0和1；步骤三：通过截集策略得到在隶属度a,b,(a+b)/2下的截集区间步骤四：进行区间不确定性分析，得到截集区间下结构响应u的上下界；步骤五：通过两点梯形公式和三点梯形公式计算结构响应上界与隶属度轴围成的面积，计算公式如下：其中为两点梯形公式所得面积，为三点梯形公式所得面积，为区间[a,b]的中点，分别是隶属度为a,b,c时求得的结构响应u的上界；步骤六：比较与判断是否收敛，如果收敛，则进行步骤七，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二到六直至每一个子区间都能满足收敛条件。步骤七：分别用两点和三点梯形公式计算结构响应下界与隶属度轴围成的面积，计算公式如下：步骤八：比较SI与SII，判断是否收敛，如果收敛，则输出结构响应结果，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二、步骤三、步骤四、步骤七、步骤八直至每一个子区间都能满足收敛条件。进一步的，所述步骤二中设置隶属度轴上的区间[a,b]，a,b的初始值分别为0和1。进一步的，所述步骤三中通过截集策略得到在隶属度a,b,(a+b)/2下的截集区间，以截集a为例，其区间上下界计算方法如下：Xa＝min{X|X∈R,μ(X)≥a}进一步的，所述步骤四中采用顶点法、泰勒展开法、模特卡洛法、优化法进行截集区间下的区间不确定性分析。进一步的，所述步骤五中通过两点梯形公式和三点梯形公式计算结构响应上界与隶属度轴围成的面积，计算公式如下：进一步的，所述步骤六中判断是否收敛的计算公式为：其中，ξ是设定的收敛因子。进一步的，所述步骤七中分别用两点和三点梯形公式计算结构响应下界与隶属度轴围成的面积，计算公式如下：进一步的，所述步骤八中判断是否收敛的计算公式为：|SI-SII|≤x本专利技术与现有技术相比的优点在于：本专利技术提供了结构模糊不确定性分析的一种新思路，弥补和完善了传统基于标准模糊算法的结构模糊不确定性分析方法，所采用的基于自适应配点的结构模糊不确定性分析方法，能够在隶属度函数变化剧烈的地方配置更多的点，在隶属度函数变化平缓的地方配置较少的点，得到较为准确的结构响应的模糊分布。在提高精度的同时，不降低效率。附图说明图1是本专利技术针对基于自适应配点的结构模糊不确定性分析方法的总体流程图；图2是本专利技术中十杆桁架结构模型示意图；图3是本专利技术中模糊不确定变量的隶属函数图，其中，图3(a)为弹性模量的隶属度函数，图3(b)为载荷的隶属度函数；图4是本专利技术中基于自适应配点及标准模糊算法所得结构响应的隶属度函数图。具体实施方式下面结合附图以及具体实施例进一步说明本专利技术。本专利技术提出了一种基于自适应配点的结构模糊不确定性分析方法，具体步骤如下：步骤一：根据所作的现实的抽象、假设情况，知识的缺乏情况，几何尺寸和加载条件，材料特性确定输入参数的模糊分布，利用合理表征贫信息、少数据条件下的模糊不确定性参数的集合。步骤二：设置隶属度轴上的区间[a,b]，a,b的初始值分别为0和1；。步骤三：通过截集策略得到在隶属度a,b,(a+b)/2下的截集区间以截集a为例，其区间上下界计算方法如下：Xa＝min{X|X∈R,μ(X)≥a}步骤四：进行区间不确定性分析，得到截集区间下结构响应u的上下界。步骤五：通过两点梯形公式和三点梯形公式计算结构响应上界与隶属度轴围成的面积，计算公式如下：其中为两点梯形公式所得面积，为三点梯形公式所得面积，为区间[a,b]的中点，分别是隶属度为a,b,c时求得的结构响应u的上界；步骤六：比较与判断是否收敛，如果收敛，则进行步骤七，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二到六直至每一个子区间都能满足收敛条件,收敛准则为；其中ξ是设定的收敛因子，是一个常数。步骤七：分别用两点和三点梯形公式计算结构响应下界与隶属度轴围成的面积，计算公式如下：步骤八：比较SI与SII，判断是否收敛，如果收敛，则输出结构响应结果，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二、三、四、七、八直至每一个子区间都能满足收敛条件,收敛准则为：|SI-SII|≤x。实施例1：为了更充分地了解该专利技术的特点及其对工程实际的适用性，本专利技术针对图2所示十杆桁架结构进行了结构模糊不确定性分析。实施例中载荷P和材料的弹性模量E被选作模糊不确定变量，其隶属度函数如图3所示。各杆的长度均为1米，截面积均为0.0001平方米。结构响应中，关注结点2的位移。ξ设为10-6。表1列出了标准模糊算法与自适应配点法的计算结果，其中，λ是截集水平，uλ是在λ截集下所得的结构响应的下界，是在λ截集下所得的结构响应的上界。标准模糊算法将隶属度[0,1]区间等分为10段。而自适应配点法在区间[0,1]自适应撒点，先求采用两点梯形公式和三点梯形公式计算区间[0,1]内的面积，即：显然不满足因此将本文档来自技高网...

【技术保护点】
1.基于自适应配点的结构模糊不确定性分析方法，可用于工程结构的模糊不确定性分析，其特征在于，实现步骤如下：步骤一：根据所作的现实的抽象、假设情况，知识的缺乏情况，几何尺寸和加载条件，材料特性确定输入参数的模糊分布，利用

【技术特征摘要】
1.基于自适应配点的结构模糊不确定性分析方法，可用于工程结构的模糊不确定性分析，其特征在于，实现步骤如下：步骤一：根据所作的现实的抽象、假设情况，知识的缺乏情况，几何尺寸和加载条件，材料特性确定输入参数的模糊分布，利用合理表征贫信息、少数据条件下的模糊不确定性参数的集合；步骤二：设置隶属度轴上的区间[a,b]，a,b的初始值分别为0和1；步骤三：通过截集策略得到在隶属度a,b,(a+b)/2下的截集区间步骤四：进行区间不确定性分析，得到截集区间下结构响应u的上下界；步骤五：通过两点梯形公式和三点梯形公式计算结构响应上界与隶属度轴围成的面积，计算公式如下：其中为两点梯形公式所得面积，为三点梯形公式所得面积，为区间[a,b]的中点，分别是隶属度为a,b,c时求得的结构响应u的上界；步骤六：比较与判断是否收敛，如果收敛，则进行步骤七，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二到六直至每一个子区间都能满足收敛条件；步骤七：分别用两点和三点梯形公式计算结构响应下界与隶属度轴围成的面积，计算公式如下：步骤八：比较SI与SII，判断是否收敛，如果收敛，则输出结构响应结果，否则，二分区间[a,b]为[a,c]和[c,b]，返回步骤二并更新隶属度轴上的区间，重复步骤二、步骤三、步骤四、步骤七、步骤八直至...

【专利技术属性】
技术研发人员：王磊，熊闯，石庆贺，王晓军，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11