振 动 与 冲 击
第 26 卷 第 1 期 JO U R N A L O F V IB R A TIO N A N D SH O C K V ol.26 N o.1 2007
基于神经网络响应面法的随机结构动力可靠度分析
基金项目:中国工程物理研究院院内基金资助项目(20030653)
收稿日期: 2005 - 11 - 10 修改 稿收到 日期 :2006 - 01 - 04
第一作者 陈 颖 男,硕士生,工程师,1975 年 3 月生
陈颖,宁佐贵
(中国工程物理研究院结构力学研究所,绵阳 621900)
摘 要
在对神经网络响应面法的原理和算法进行研究的基础上,建立了基于神经网络响应面法的随机结构动力
可靠度分析方法。首先,基于首次超越破坏准则,参照静力可靠度的功能函数模式,建立了随机结构的动力可靠度功能函
数;然后引入响应面法,以三层 BP 神经网络作为拟合函数,推导了功能函数的拟合表达式;最后结合一次二阶矩方法求
解可靠指标。算例分析表明了本文方法有较好的计算精度和计算效率,在复杂结构的动力可靠度分析中有较强的实用
价值。
关键词: 神经网络,响应面法,随机结构,动力可靠度
中图分类号: TU 311,O 242 文献标识码:A
由于工程实际结构的复杂性和所用材料在统计上
的离散性以及测量、加工、制造误差的存在,结构参数
在许多方面都具有随机性,会给结构动力响应分析结
果带来不容忽略的影响,在一定条件下还可能成为主
导因素
[ 1 ,2]
。因此具有随机参数的结构(以下简称随机
结构)在随机振动载荷下的动力可靠度分析问题越来
越受到国内外学者的重视。
目前常用的随机结构动力响应和可靠度分析是基
于随机有限元方法进行
[1 ~4]
,但它需要推导复杂的计
算公式,特别是其中求解结构特性关于各随机量的偏
导往往是非常麻烦的。而且随机有限元难以与现有有
限元分析商用软件如 A N SY S、N ASTRA N 等结合,需要
搭建自己的计算分析平台,大大限制了其在复杂工程
结构中的实用性。
近年来响应面法(R SM )
[ 5,6 ]
的蓬勃兴起,为该问题
的解决提供了新的思路。本文在建立了随机结构动力
可靠度功能函数的基础上,引入响应面法来对功能函
数进行拟合求解,使整个可靠度分析过程中无需对随
机参数作数学处理,并可使用现有商用软件进行有限
元分析。在响应面函数的构建问题上,对新兴的神经
网络响应面法
[ 7,8 ]
进行了研究,引入三层 BP 神经网络
作为响应面函数,推导了相应的功能函数拟合表达式
和可靠指标求解公式。最后通过算例分析表明了本文
方法的有效性。
1 随机结构的动力可靠度功能函数的建立
在随机结构的静力可靠度分析中,常用到的结构
功能函数为:
Z=L-Y (1)
其中 L 为界限值,Y 为某种响应值,求可靠度实际
上就是求 Z > 0 的概率。
现将这种模式扩展到动力分析中。在随机振动载
荷下,结构的响应是随时间变化的随机过程,因此功能
函数也与时间有关,称为功能随机过程
[10]
:
Z(t) = L - Y(t) (2)
在某个时间区域[0,T] 内式(2) 进行功能极小化
变换,得到:
Z
min
(t)=min(L-Y(t))=L-y
max
t∈[0,T]
(t) (3)
y
max
t∈ [0,T]
(t)是结构响应 y(t)在时间区域[0,T]内的
极大值。本文讨论动载荷为平稳随机过程的情况。采
用首次超越破坏准则
[12]
(即结构的响应首次超过了一
个规定的上限值或下限值,则认为结构遭到破坏),参
照确定性结构在平稳随机载荷下响应极大值的推导过
程,将 y
max
t∈[0,T]
(t)表示为:
y
max
t∈[0,T]
(t) = pσ
y
(4)
上式 中, p 为 随 机变 量, 其均 值、方差和分布函数分
别为:
μ
p
= (2lnv
+
0
T)
1
2
+
0.5772
(2lnv
+
0
T)
1
2
σ
2
p
=
π
2
6
(2lnv
+
0
T)
-1
(5)
F(p) = exp - v
+
0
Texp
-
p
2
2
()[]
(6)
v
+
0
为结构的正斜率期望穿零率,计算方法详见文献
[12]。
σ
y
为结构在随机振动载荷下的响应总均方根值。
对于确定性结构来说,σ
y
是定值;对于随机结构来说,
由于结构的材料常数、尺寸、阻尼比等是随机变量,对
于其中的每一组样本取值,都对应了一种响应总均方
根值,因此 σ
y
应描述为随结构参数取值而变化的随机
变量。