CN108763604B - 一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 - Google Patents
一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 Download PDFInfo
- Publication number
- CN108763604B CN108763604B CN201810120219.3A CN201810120219A CN108763604B CN 108763604 B CN108763604 B CN 108763604B CN 201810120219 A CN201810120219 A CN 201810120219A CN 108763604 B CN108763604 B CN 108763604B
- Authority
- CN
- China
- Prior art keywords
- radial basis
- basis function
- neural network
- composite material
- material structure
- 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
Images
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]
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)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法。该方法首先在复合材料结构参数区间域内选择一系列样本点,并使用有限元的方法求解这些样本点的结构静力学响应值。将这些数据作为径向基神经网络的训练样本。根据样本点的数目选择神经网络隐藏层中径向基函数的数量,再用K均值算法得到各个径向基函数的中心值,然后再用递归最小二乘方法得神经网络隐藏层和输出层之间的权值。将训练好的径向基神经网络作为原结构响应函数的近似响应函数。再利用遗传算法求解径向基神经网络的最大和最小值作为含区间参数的复合材料结构静力学响应的上界和下界。
Description
技术领域
本发明涉及含区间不确定参数的复合材料结构静力学响应分析的技术领域,具体涉及一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点法。
背景技术
结构分析在机械工程、土木工程、车辆工程和航空航天工程中占有很重要的位置。为了保证结构能够在其经济寿命内安全、可靠地工作,在其服役之前,其应变、应力和位移等各项指标首先需要仔细地校核。结构分析在结构设计的整个环节中有着举足轻重的地位,而静力学响应分析又是结构分析中最为基本和重要的环节之一。传统的结构静力学分析往往把结构的各项参数视为确定的值,一整套分析过程之后会得到一个确定的值。
复合材料作为新型材料的典型代表,具有高的比强度和比刚度、优越的疲劳性能等优良的力学特性,其合理应用可实现结构有效减重,具有广阔的应用前景。但是由于复合材料的分散性较大,加工过程中不可控因素较多,因而复合材料结构的各项参数均存在不容忽视的不确定性的存在。这些不确定性存在可能会对分析结果造成比较大的影响。
一般描述不确定参数的方法有概率密度函数、模糊集合以及证据理论等等,但是为了得到概率密度函数曲线需要有大量的数据,一般来说很难通过实验得到如此多的数据;而利用模糊集合时需要这个领域的专家分配隶属度函数;运用证据理论时,同样也需要专家分配概率密度给各个焦元。当数据的信息量和对数据的认知较少时,区间方法在对不确定性进行表征时有较大的优势。区间方法表示不确定性时只需要知道不确定参数的上下界。
一般的求解区间不确定性传播的方法有蒙特卡洛法,摄动法,顶点法以及配点法等。运用蒙特卡洛方法进行求解时需要在不确定域内进行大量的采样,因而需要大量的有限元计算,这在实际工程中是无法实施的;摄动法为一种一阶近似方法,其只能够处理若非线性和不确定区间较小的问题;顶点法只能够处理单调问题;而配点法能够处理非线性较强的问题,并且能够得到较为精确的解,但是其处理过程较为复杂。
由于在复合材料结构静力学响应分析过程中,复合材料的参数不确定性广泛存在且不容忽视,因而建立一种精度较高,并且计算量小的含区间不确定性参数的复合材料结构的静力学分析方法有着很强的现实意义。
发明内容
本发明要解决的技术问题是:本发明提供了一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法。该方法首先在复合材料结构参数区间域内选择一系列样本点,并使用有限元的方法求解这些样本点的结构静力学响应值。将这些数据作为径向基神经网络的训练样本。根据样本点的数目选择神经网络隐藏层中径向基函数的数量,再用K均值算法得到各个径向基函数的中心值,然后再用递归最小二乘方法得神经网络隐藏层和输出层之间的权值。将训练好的径向基神经网络作为原结构响应函数的近似响应函数。再利用遗传算法求解径向基神经网络的最大和最小值作为含区间参数的复合材料结构静力学响应的上界和下界。
本发明采用的技术方案为:一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法,该方法的主要步骤如下:
第一步:确定区间不确定性变量αI以及其区间域Θ;确定各个径向基函数的方差θ。
第二步:在不确定域Θ内选择样本点α(1),...α(N)(简写为α(1):(N));利用有限元方法计算复合材料结构的静力学响应Ku(α(i))=F,i=1,...N;再根据选择的样本点的个数确定隐藏层中径向基函数的数目M。
其中||·||表示向量的二范数。输出层为隐藏层的带权求和,设输出层的数据为y(i),则输出层可以表示为:
其中wj为连接第j个径向基函数和输出层的权值,记w=[w1,...,wM]。该径向基神经网络中的w和c(j),j=1,...,M需要经过训练得到。
第四步:确定利用K均值算法计算各个径向基函数的中心c(j),j=1,...,M。首先随机初始化径向基函数的中心c(j),j=1,...,M;再计算每个样本点距离每个径向基函数中心的欧氏距离Dij:
找到离样本点α(i)最近的径向基函数中心c(j),并且令指标函数如下:
zip=1,p=j,zip=0,p≠j
如果在当第q个迭代步中有:
那么表示K均值算法已经收敛,此时输出径向基函数的数据中心c(j),j=1,...,M;否则利用以下的公式对径向基函数的中心进行更新:
直到K均值算法收敛。
第五步:利用递归最小二乘法计算隐藏层中各个径向基函数到输出层的权值w。首先定义第i个样本点映射到第j个径向基函数表示为:
为了简化符号再定义如下的向量:
φ(i)=[φi1,φi2,...,φiM]T
和定义如下的矩阵:
输出层的权值训练问题可以转换为最小二乘问题,则由该问题的法方程可以得到:
R(s+1)w(s+1)=r(s+1)
其中w(s+1)表示用N个样本点中的(s+1)得到的权值向量。然后初始化w(0)=0和R-1(0)=γ-1I,γ为一个正的小量,I∈RM×M表示单位矩阵。计算:
h(s+1)=u(α(s+1))-φT(s+1)w(s)
并利用谢尔曼-莫里森-伍德伯里公式计算:
然后利用下式更新隐藏层到输出层的权值:
w(s+1)=w(s)+R-1(s+1)φ(s+1)h(s+1)
直到用完N个样本点数据。
本发明与现有技术相比优点在于:本发明通过利用径向基神经函数网络求解含区间不确定性参数的复合材料结构的静力学响应分析问题。与蒙特卡洛方法相比,该发明精度相近,但是计算量远远小于前者;与摄动法和顶点法相比,该方法精度高,并且能处理更为复杂的问题;与配点法相比,该方法实施简单。总而言之,本发明提供了一种高精度、高效的分析含区间不确定性的复合材料结构静力学响应的方法
附图说明
图1是本发明求解含区间参数的复合材料结构的静力学响应求解的流程图。
图2是本发明构造的径向基函数神经网络。
图3是本发明应用的对象——复合材料中央翼盒结构,其中图3(a)是该中央翼盒的外观示意图,左边的矩形方块表示壁面,中央翼盒固定连接在壁面上;图3(b)为中央翼盒的俯视图,D1~D11表示响应的输出位置。
图4表示本发明应用对象——复合材料中央翼盒结构的有限元网格划分形式。
图5表示用本发明和蒙特卡洛方法得打的各个测点的纵向位移的区间图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法,包括以下的步骤:
(1)确定区间不确定性变量以及其区间;确定各个径向基函数的方差。
(2)在不确定域内选择样本点;利用有限元方法计算复合材料结构的静力学响应;再根据选择的样本点的个数确定隐藏层中径向基函数的数目。
(3)构建径向基函数神经网络。
(4)确定利用K均值算法计算各个径向基函数的中心。
(5)利用递归最小二乘法计算隐藏层中各个径向基函数到输出层的权值。
(6)利用遗传算法搜索训练好的神经网络的最大和最小值,并将其作为此含区间参数的复合材料结构的静力学响应的上界和下界。
该方法的具体步骤如下:
第一步:确定区间不确定性变量αI以及其区间域Θ;确定各个径向基函数的方差θ。
第二步:在不确定域Θ内选择样本点α(1),...α(N)(简写为α(1):(N));利用有限元方法计算复合材料结构的静力学响应Ku(α(i))=F,i=1,...N;再根据选择的样本点的个数确定隐藏层中径向基函数的数目M。
其中||·||表示向量的二范数。输出层为隐藏层的带权求和,设输出层的数据为y(i),则输出层可以表示为:
其中wj为连接第j个径向基函数和输出层的权值,记w=[w1,...,wM]。该径向基神经网络中的w和c(j),j=1,...,M需要经过训练得到。
第四步:确定利用K均值算法计算各个径向基函数的中心c(j),j=1,...,M。首先随机初始化径向基函数的中心c(j),j=1,...,M;再计算每个样本点距离每个径向基函数中心的欧氏距离Dij:
找到离样本点α(i)最近的径向基函数中心c(j),并且令指标函数如下:
zip=1,p=j,zip=0,p≠j
如果在当第q个迭代步中有:
那么表示K均值算法已经收敛,此时输出径向基函数的数据中心c(j),j=1,...,M;否则利用以下的公式对径向基函数的中心进行更新:
直到K均值算法收敛。
第五步:利用递归最小二乘法计算隐藏层中各个径向基函数到输出层的权值w。首先定义第i个样本点映射到第j个径向基函数表示为:
为了简化符号再定义如下的向量:
φ(i)=[φi1,φi2,...,φiM]T
和定义如下的矩阵:
输出层的权值训练问题可以转换为最小二乘问题,则由该问题的法方程可以得到:
R(s+1)w(s+1)=r(s+1)
其中w(s+1)表示用N个样本点中的(s+1)得到的权值向量。然后初始化w(0)=0和R-1(0)=γ-1I,γ为一个正的小量,I∈RM×M表示单位矩阵。计算:
h(s+1)=u(α(s+1))-φT(s+1)w(s)
并利用谢尔曼-莫里森-伍德伯里公式计算:
然后利用下式更新隐藏层到输出层的权值:
w(s+1)=w(s)+R-1(s+1)φ(s+1)h(s+1)
直到用完N个样本点数据。
实施例:
为了更充分地了解该发明的特点及其对工程实际的适用性,本发明针对图3所示的空天飞机的碳纤维复合材料中央翼盒结构的静力学响应分析进行求解。中央翼盒的尺寸为900mm×710mm×250mm,这个中央翼盒包含5个部分:上下翼面,翼肋,翼梁,机翼前缘和机翼后缘。整个中央翼盒都由MT300碳纤维复合材料构成,各个部分的铺层如表1所示。
表1
该复合材料中央翼盒的有限元网格如图4所示,其包含12361个单元和8723个节点。集中载荷施加在中央翼盒尖部,其大小为P=7200N。复合材料的杨氏模量和泊松比为区间变量,其范围如下:
图3(b)中的测点D1~D11的纵向位移为响应输出量。纵向位移的精确区间可由蒙特卡洛方法得到。利用基于径向基函数的方法和蒙特卡洛方法得到的位移区间如图5所示。蒙特卡洛方法共采用了100000个样本点,本发明在每一个维度上取了5个样本点。图5表面,在使用较少样本点的条件下,本发明能够得到很精确的位移响应区间。
综上所述,本发明提出了一种基于径向基神经网络的求解含区间参数复合材料结构静力学响应的方法。该方法首先在复合材料结构参数区间域内选择一系列样本点,并使用有限元的方法求解这些样本点的结构静力学响应值。将这些数据作为径向基神经网络的训练样本。根据样本点的数目选择神经网络隐藏层中径向基函数的数量,再用K均值算法得到各个径向基函数的中心值,然后再用递归最小二乘方法得神经网络隐藏层和输出层之间的权值。将训练好的径向基神经网络作为原结构响应函数的近似响应函数。再利用遗传算法求解径向基神经网络的最大和最小值作为含区间参数的复合材料结构静力学响应的上界和下界。
以上仅是本发明的具体步骤,对本发明的保护范围不构成任何限制;其可扩展应用于含区间参数的其它结构形式的复合材料结构静力学响应求解领,凡采用等同变换或者等效替换而形成的技术方案,均落在本发明权利保护范围之内。
本发明未详细阐述部分属于本领域技术人员的公知技术。
Claims (1)
1.一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法,其特征在于:首先在复合材料结构参数区间域内选择一系列样本点,并使用有限元的方法求解这些样本点的结构静力学响应值;将这些数据作为径向基神经网络的训练样本;根据样本点的数目选择神经网络隐藏层中径向基函数的数量,再用K均值算法得到各个径向基函数的中心值,然后再用递归最小二乘方法得神经网络隐藏层和输出层之间的权值;再利用遗传算法求解径向基神经网络的最大和最小值作为含区间参数的复合材料结构静力学响应的上界和下界;该方法的实现步骤如下:
第一步:确定区间不确定性变量αI以及其区间域Θ;确定各个径向基函数的方差θ;
第二步:在区间域Θ内选择样本点α(1),...α(N)(简写为α(1):(N));利用有限元方法计算复合材料结构的静力学响应Ku(α(i))=F,i=1,...N;再根据选择的样本点的个数确定隐藏层中径向基函数的数目M;
其中,||.||表示向量的二范数,输出层为隐藏层的带权求和,设输出层的数据为y(i),则输出层表示为:
其中,wj为连接第j个径向基函数和输出层的权值,记w=[w1,...,wM],该径向基神经网络中的w和c(j)需要经过训练得到,其中,j=1,...,M;
第四步:确定利用K均值算法计算各个径向基函数的数据中心c(j),j=1,...,M,首先随机初始化径向基函数的数据中心c(j),j=1,...,M;再计算每个样本点距离每个径向基函数的数据中心的欧氏距离Dij:
找到离样本点α(i)最近的径向基函数的数据中心c(j),并且令指标函数如下:
zip=1,p=j,zip=0,p≠j
如果在当第q个迭代步中有:
那么表示K均值算法已经收敛,此时输出径向基函数的数据中心c(j),j=1,...,M;否则利用以下的公式对径向基函数的中心进行更新:
直到K均值算法收敛;
第五步:利用递归最小二乘法计算隐藏层中各个径向基函数到输出层的权值w,首先定义第i个样本点映射到第j个径向基函数表示为:
为了简化符号再定义如下的向量:
φ(i)=[φi1,φi2,...,φiM]T
和定义如下的矩阵:
输出层的权值训练问题转换为最小二乘问题,则由该问题的法方程得到:
R(s+1)w(s+1)=r(s+1)
其中w(s+1)表示用N个样本点中的(s+1)得到的权值向量,然后初始化w(0)=0和R-1(0)=γ-1I,γ为一个正的小量,I∈RM×M表示单位矩阵,计算:
h(s+1)=u(α(s+1))-φT(s+1)w(s)
并利用谢尔曼-莫里森-伍德伯里公式计算:
然后利用下式更新隐藏层到输出层的权值:
w(s+1)=w(s)+R-1(s+1)φ(s+1)h(s+1)
直到用完N个样本点数据;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810120219.3A CN108763604B (zh) | 2018-02-07 | 2018-02-07 | 一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810120219.3A CN108763604B (zh) | 2018-02-07 | 2018-02-07 | 一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763604A CN108763604A (zh) | 2018-11-06 |
CN108763604B true CN108763604B (zh) | 2022-03-15 |
Family
ID=63980066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810120219.3A Active CN108763604B (zh) | 2018-02-07 | 2018-02-07 | 一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763604B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046601A (zh) * | 2019-10-25 | 2020-04-21 | 武汉数字化设计与制造创新中心有限公司 | 基于工件变形代理模型的薄壁曲面加工进给速度优化方法 |
CN111177966B (zh) * | 2019-12-30 | 2021-10-01 | 北京航空航天大学 | 一种基于贝叶斯理论的导弹结构不确定载荷区间重构方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942375A (zh) * | 2014-04-11 | 2014-07-23 | 浙江大学 | 基于区间的高速压力机滑块尺寸稳健设计方法 |
CN104537134A (zh) * | 2014-06-16 | 2015-04-22 | 北京空间飞行器总体设计部 | 一种动态级联耦合机电系统的可靠性建模与评估方法 |
CN104679956A (zh) * | 2015-02-14 | 2015-06-03 | 浙江大学 | 考虑动态特性的高速压力机底座可靠性稳健设计方法 |
CN105260532A (zh) * | 2015-09-30 | 2016-01-20 | 浙江大学 | 基于序列近似优化的薄板拉伸变压边力不确定性设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120317058A1 (en) * | 2011-06-13 | 2012-12-13 | Abhulimen Kingsley E | Design of computer based risk and safety management system of complex production and multifunctional process facilities-application to fpso's |
-
2018
- 2018-02-07 CN CN201810120219.3A patent/CN108763604B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942375A (zh) * | 2014-04-11 | 2014-07-23 | 浙江大学 | 基于区间的高速压力机滑块尺寸稳健设计方法 |
CN104537134A (zh) * | 2014-06-16 | 2015-04-22 | 北京空间飞行器总体设计部 | 一种动态级联耦合机电系统的可靠性建模与评估方法 |
CN104679956A (zh) * | 2015-02-14 | 2015-06-03 | 浙江大学 | 考虑动态特性的高速压力机底座可靠性稳健设计方法 |
CN105260532A (zh) * | 2015-09-30 | 2016-01-20 | 浙江大学 | 基于序列近似优化的薄板拉伸变压边力不确定性设计方法 |
Non-Patent Citations (6)
Title |
---|
Estimation of contour parameter uncertainties in permittivity imaging using MCMC sampling;C. Schwarzl et al;《 2008 5th IEEE Sensor Array and Multichannel Signal Processing Workshop》;20080822;全文 * |
HGMs_EP复合泡沫塑料材料设计与性能研究;何长军;《中国博士学位论文全文数据库电子期刊 信息科学辑》;20111115;第2011年卷(第11期);全文 * |
RBF神经网络在极限承载力预测中的应用;王磊 等;《长沙交通学院学报》;20041231;第20卷(第4期);全文 * |
Research on dynamic modeling and simulation of axial-flow pumping system based on RBF neural network;Qinghui Wu et al;《Neurocomputing》;20160419;第186卷;全文 * |
机翼颤振控制与不确定性优化问题研究;张军红;《中国博士学位论文全文数据库电子期刊 工程科技II辑》;20120715;第2012年卷(第7期);全文 * |
飞行器复合材料结构优化设计研究进展;王晓军 等;《中国科学》;20180131;第48卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108763604A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104866692B (zh) | 一种基于自适应代理模型的飞行器多目标优化方法 | |
CN106874573B (zh) | 一种分区变厚度复合材料层合板的设计方法 | |
CN108763604B (zh) | 一种求解含区间参数复合材料结构静力学响应的径向基神经网络配点方法 | |
CN105893716B (zh) | 一种基于分形理论的结构断裂非概率可靠性分析方法 | |
CN113935086B (zh) | 一种基于计算机视觉和深度学习的结构智能设计方法 | |
CN111105005B (zh) | 一种风电功率预测方法 | |
CN113343559B (zh) | 迭代重加权最小二乘法极限学习机响应面可靠性分析方法 | |
CN114722655A (zh) | 一种基于局部有限寿命疲劳约束条件的结构拓扑优化方法 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
CN108830407B (zh) | 多工况条件下结构健康监测中的传感器分布优化方法 | |
CN114564787A (zh) | 用于目标相关翼型设计的贝叶斯优化方法、装置及存储介质 | |
Salmalian et al. | Multi-objective evolutionary optimization of polynomial neural networks for fatigue life modelling and prediction of unidirectional carbon-fibre-reinforced plastics composites | |
Ambjorn et al. | Topology induced first-order phase transitions in lattice quantum gravity | |
CN111079279A (zh) | 一种多构型晶格结构多尺度拓扑优化设计方法 | |
CN113688465B (zh) | 一种基于载荷与状态结合的飞行器结构强度数字孪生方法 | |
CN114611310A (zh) | 一种小子样条件下基于模糊理论的参数不确定性建模方法 | |
CN114282310A (zh) | 一种基于自适应加点代理模型的气弹结构耦合寻优方法 | |
CN115600383A (zh) | 一种不确定性数据驱动计算力学方法、存储介质及产品 | |
CN114329805A (zh) | 基于自适应代理模型的连接件多学科协同设计优化方法 | |
CN108920787B (zh) | 一种基于自适应配点的结构模糊不确定性分析方法 | |
CN110826182A (zh) | 一种基于顶点法和序贯优化策略的飞行器结构气动弹性鲁棒设计方法 | |
CN115146366B (zh) | 基于Direct算法和少量样本更新的结构混合可靠性分析方法 | |
Gong et al. | A multiparameter group decision making method based on the interval-valued intuitionistic fuzzy soft set | |
Athira et al. | Incomplete pythagorean fuzzy soft sets | |
CN118051996A (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 |