处理有界不确定问题的区间方法在模型确认中应用.pdfVIP

处理有界不确定问题的区间方法在模型确认中的应用1 李亚利* 邱志平* * （北京航空航天大学固体力学研究所，北京，100191） 模型确认与验证活动中的误差源包括物理建模误差、离散误差、舍入误差和程序误差，而不确定性源 包括建模不确定性和参数不确定性。模型确认的衡量标准是寻找实验数据和计算结果的符合度，而衡量标 准中应该包含对计算中的不确定性和误差的估计。对于不确定性分析或不确定性传播可以简单分为三类方 法：非概率方法、隶属函数和概率方法。 非概率方法包含区间分析、灵敏度分析和模糊逻辑。本文所应用的理论是区间分析，所要解决的不确 定问题是：已知有界不确定参数或未确知参数向量所在范围的上界向量和下界向量，确定结构响应的上界 和下界。在文中我们将求解具有有界不确定参数或未确知参数结构响应问题的集合理论凸方法的数学理论 模型应用在模型验证与确认活动中。利用具有有界不确定参数或未确知参数结构静力位移问题的集合理论 凸方法模型和具有有界不确定参数或未确知参数结构振动固有频率问题的集合理论凸方法模型的求解方 法给出模型确认中带有不确定区间的精确解或者解析解。从而为模型确认与验证中用到的带不确定性的结 构响应的基准解提供依据。有界不确定结构响应上下界的求解问题是全局最优化问题。只有结构响应是凸 函数时，在理论上才能保证有界不确定参数或未确知参数结构响应的上下界结构响应是全局最优解。 关键词：模型验证与确认 区间方法 不确定性 引言 模型验证和模型确认的概念自1967 年被明确提出以来受到了越来越多的关注。20 世纪90 年代末，模 型验证与确认的发展进入了一个新阶段，围绕解决预测难度评估问题，模型验证与确认的现代思想开始萌 芽。系统化的模型验证与确认方法最早在90 年代末出现于美国能源部的一项研究计划中。后来，在力学 界，特别是结构动力学领域得到了许多研究人员的关注。圣地亚实验室发表了许多关于现代模型验证与确 [1-3] 认方法学的报告 ，基本思想是，辨识和量化数值模拟实施过程每一环节 （建模、离散、算法选择、编程、 参数选取、计算和分析等，包括相关实验）产生的误差和不确定性，从而对全系统数值模拟的预测置信度 做出推断。现代模型验证与确认有三个显著特点：预测性、系统分层和不确定性量化。不确定性量化是模 1 国家自然科学基金委员会- 中国工程物理研究院联合基金资助项目 型验证与确认的关键也是难点所在。当具有足够的信息描述不确定参量的概率特征时，概率可靠性是一种 较为理想的模型。然而，由于其自身难以克服的种种局限性限制了概率可靠性理论在工程实际中的推广和 应用。基于此，区间分析、凸模型等非概率方法也在很好的发展。 集合理论凸方法是用一个集合来代表实际工程结构系统所有模型的。或者说，根据实际工程结构系统 的输入和输出所提供的信息，来确定一个总能包含真实工程结构系统 （这个工程结构系统可能不同于我们 抽象出来的数学模型，或者说，随着时间的推移，这个工程结构系统发生了变化。如结构参数的变化或误 差，没有考虑到的外界干扰以及载荷的变化，误差或不确定性）的数学模型。这个数学模型显然是一个结 构系统族或模型集。因此集合理论凸方法所建立的模型不是一个单一的结构系统模型，而是一个系统族或 模型集。在这个集合中，按问题的要求可以求出或确定结构系统响应所在的集合或结构系统响应的界限， 确切地说，可以求出结构系统的最大可能响应和最小可能响应。[4] 我们所要解决的不确定工程问题是：已知有界不确定参数或未确知参数向量 α ( ) ( , , , )T 所在范围的上界向量α ( ) ( , , , )T 和下界向量 i m 1 2 m i m 1 2 m α ( ) ( , , , )T ，如何确定结构响应的上界和下界。 i m 1 2 m 有界不确定结构响应上下界，即结构最大可能响应和最小可能响应的求解问题是全局优化问题。这是 一个目前在数学领域里还没有得到很好解决的问题。在优化数学理论中，对全局优化问题，只有目标函数 和约束函数均为可行集上的凸函数时，目标函数的任一局部极小值或极大值都是目标函数在非空可行集上 的全局极小值和极大值。在优化问题的求解过程中，在一般情况下，往往只能求出它的一个局部最优解。 所以，必须分析集合理论凸方法