材料物性参数识别的梯度正则化方法Ξ.pdfVIP

第 17 卷第 1 期 计 算 力 学 学 报 V o l. 17 N o. 1 2000 年 2 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 2000 文章编号:(2000) 材料物性参数识别的梯度正则化方法 1 1 2 1 3 3 刘迎曦 , 　王登刚 , 　张家良 , 　李守巨 , 　路振刚 , 　于亚军 ( 1 大连理工大学工程力学系, 大连 116023; 2 辽宁工学院, 锦州 121001; 3 白山发电厂, 桦甸 132400) 摘　要: 本文对梯度正则化方法( ) 作了进一步的研究, Gradient R egularization M ethod 给出一种建立了梯度正则化迭代算法和选择正则参数的简明实用方法。文中椭圆算子 方程参数识别算例不仅说明了 GR 法具有广泛的适用性和一定的抗噪音能力, 而且收 敛速度较快, 具有较大的收敛范围。 关键词: 梯度正则化方法, F rechet 导数, 正则参数, 参数识别, 不适定, 椭圆算子方程 中图分类号: P 63142　　　文献标识码:A 1　引　　言 反问题普遍存在于自然现象和工程实际应用中, 而其中以待定未知模型参数较为常见。待 定未知参数有不同的途径可循, 最明显最直接的方法是测量这些参数在一些离散点处的值, 但 是这种方法在实际工程应用中往往是非常困难甚至是行不通的。这时人们不得不转去测量与 ( ) 待定参数有一定关系 一般通过微分关系联系着 的其他可观测参量在边界上的值, 来推断模 型参数的值, 这种根据给定的可测参量的观测值, 推断模型参数的值即是我们所说的参数识别 问题。 近年来反问题已成为各个科学领域的热门课题。但是反问题的不适定性仍然是反问题理 论研究和工程应用的一个瓶颈问题。由于反问题具有不适定性, 强非线性, 计算量大等难点, 而 且从实际问题中提出的反问题类型各异, 涉及的微分方程各不相同, 因此基本上没有统一的求 解方法, 但是我们所用的很多方法都可以用有限的原则来描述, 即可以为反演建立一种理 [ 1 ] 论 。总的来说反问题的求解方法可分为直接法和间接法, 直接法一般计算速度快, 但误差积 累严重, 抗噪音能力弱, 因此没有间接法那样的广泛适用性。间接法一般为数值迭代, 其中有些 方法具有比较广泛的适用性, 梯度正则化方法就是其中之一。 ( ) 梯度正则化方法首先由文献[2 ]和[3 ]提出, 并验证了其在一维双曲方程 一维波动方程 中的适用性。目前人们认为椭圆方程算子参数比识别双曲方程算子参数和抛物方程算子参数 [ 4 ] 困难的多 , 本文验证了采用该法识别二维椭圆算子参数的适用性, 给出一种建立了梯度正则 化迭代算法的简明方式, 对正则参数的选取作了初步的研究, 以希望扩大算法的收敛范围, 减 小算法对初值的严重依赖性, 降低其对先验信息的要求, 使其更具有实用性。本文工作是文献 [3 ] 的继续。 2　梯度正则化方法原理及算法 2 1　反问题的数学模型 从工程、物理中提出的反问题都有明确的数学形式和定量要求, 其一般的数学模式是偏微 收稿日期: 　基金项目: 国家自然科学基金资助, 项目编号: ( )