感 度 解 析 を用 い た振動 系 のパ ラ メー タ同 定* (収束性の改薯に関する検討) 鞍 谷 文 保航 岩 壼 卓 三艦 沖 田 耕 三寧1 PammeterIde魏ti且cationofVibrationSy8tem8U8ingSen8itMtyAnaly8b (Con8i"e凱ionoHmprovemεntofConvergenee) F㎜iyasuKURATANI,Taku301WATSUBOandKozoOKITA P陀vi。u叩aper5discussedam舳odforidentifyingmodeiparamεtεr5。f匙v量bration5ystemby mεansofanon-linearleast5quarεminimi2ationpro6εdurεofanevaluationfunctionwhi¢hconsistg ofthedi鉦 ετen`eofnaturalfτequendεs,mode5hapecoεf五cient50fnatur置1mode5betw㏄nth¢ meagureddataandthεdata仁alculated`romth巴e8tima把dmoddparamete霜.ThemεthDdusin醒the p艶viou8eΨaluationfunctiontendston㏄es5itatgalargenumberofiteτation&Thi5paperpropose菖 amethodu51nganew¢ Ψaluatio鯛functionforthepu叩oseofdeαeasin騒thenumb¢rofiteratlonsand steadyin8cDnver区ence5.1tco皿tain35calarΨalueswhiぐheΨaluatethe5imilarityofmode5hape betw㎝themεa5uredn呂tufaimodεandthe£alcul飢edoneof田ohorder,in3tεadofεachmode由ape oo¢ 伍cientofnatu囮1modε5.Thee旅 詑tsofusingthe訊bovε ・de飴{きd5亡alarvaluesa陀di5cusscd. Numericalmodelsa祀examin已dtoverifytheefnciεncyofthemethodu5inganewevaluation function. IdEntifiしation,Sensi‡ivity,N邑tu岨1ModqConverg¢nc¢ κ εy研or幽:Vibration,P融rameter lmprovement 1.ま え が き た 同 定 法 は,実 験 で 得 られ た固 有 振 動 数,固 有 モ ー ド と物 理 モ デ ルの 固有 振 動 数,固 有 モ ー ドを比 較 し,両 機械 構造 物 の振 動 挙動 の 予測 や 特性 改 善 を精度 よ く 者 の残差 二 乗和 が 最小 にな る よ うに,感 度 解 析 を用 い 行 うた め には,解 析 に 用 い る数学 モ デ ルの 信頼 性 を向 て物 理 モ デル の物理 的 特性 を繰返 し修 正 す る反復 改 善 上 させ る こ とが 重 要 とな る.そ こで,振 動 実験 ヂー タ 計算 法 で あ る.と ころで,検 討 した 同 定法 の よ うな 反 を も とに.対 象物 の 動特 性 を的確 に表 現 す る数 学 モ デ 復改 善 計 算 法 を実 用 化 す る上 で は,収 束速 度 収 束 安 ル を決 定 し よ うとす る,い わ ゆ る 系の 実験 的 同 定 に 関 定性 が重 大 な課 題`5絢 とな っ て い る,著 者 の 方 法 に お す る研 究mが 行 わ れ て い る,系 の実 験 的 同定 に は動 特 い て も,同 定 計 算過 程 に お い て多 くの反 復 計 算 を要 し 性 を固有 振 動 数,モ ー ド減 衰比,固 有 モー ドで 表 現 す た り,収 束 安 定性 が悪 くな る場 合 が あ り,改 善 の余 地 る実 験 的 モー ド解 析 法 と,質 量 、剛性,減 衰 係 数 な ど が認 め られ た. の物 理 的特 性 をパ ラ メー タ として物 理 モ デ ルで表 現 す そ こ で,本 報 で は,同 定 の評 価 に 固 有 振 動 数,固 有 る方法 が あ る.理 論解 析 との 結合,構 造 変 更 に よる最 モー ドを用 い る場 合 に問題 を限 定 した 上 で,同 定計 算 適化 な どへの拡 張 性 を考 えれ ば,系 を物 理 的 特 性 で 表 の収 束性 を改 善 させ る方法 に っ いて 検討 した.そ の結 現 す る方法 が 望 ま しい. 果,従 来 の 同 定法 に おい て,固 有 モ ー ドの各 成 分 を評 著 者 ら も,こ れ まで に,系 を物 理 モ デ ル で 表 現 し, 価関 数 とす る こ とが,収 束性 を悪 くす る要 因 の一 つで その物 理 的特 性 を実験 デー タか ら同 定 す る方法 につ い あ る こ とが推 測 され た.そ こで,収 束 性 を改 善 す る た て検討 して きた.そ して,ば ね ・質 量 系の 同定9}、弾性 め に、従 来 とは異 な る固 有 モ ー ドの評 価 方 法 を 用 い, 体 で支 持 され た剛 体 系 の 同 定1鮒,を行 い,実 験 デ ー タ その評 価 方法 で得 られ た評価 量 を評価 関 数 とす る同定 か ら物 理 的特 性 が同 定 で きる こ とを確認 した.検 討 し 法 を提 案 す る.そ して,従 来 の同 定法 に比 べ て.収 束 性 噌 平 成2年 暮 月30日 第65期 通 常 総 会 講 演 会 に お い て 講 演 , 原稿受付 平成 言年4月25日, 朝 正員 、兵 庫 県 立 工 業技 術 セ ン タ ー1◎654神 戸 市須磨区行 平 町31-12〕, 岡 正 員 ,神 声 大 学 工学 部1蜘57神 戸 市 灘 区 六 甲 台町1). が大 幅 に向 上 す る こ とを示 す. 2.従 来 の 同定 法 の問 題 点 と改善 策 2・1従 来 の 同定 法 従 来 の 同 定 法 掛の 計 算 手 順 を 説 明 す る.実 験 で 得 られ た 固 有 振 動 数 力¢お よ び 固 通 常 は,1回 の 修正 作 業 で 最良 な モ デ ルパ ラ メー タ 有 モー ド仏 碍 の 各成 分 を並 べ た 同定 の 目標 とす る次 の 値 が 得 ら れ る 場 合 は 少 な く,最 良 な 値 が 得 ら れ る ま の ベ ク トル{尺 り を 設 定 す る. で 修 正 作 業 を 繰 返 す こ と に な る. 伍 ド}=伍 ε,…,ノノ,{〃 暉 り1,…,〔 〃。 曜}7}ア … …(1〕 し か し,式(2)の 評 価 関 数 を 用 い た 場 合,初 期 値 と 次 に,初 期 モ デ ル の 固 有 値 解 析,感 度 解 析 を 行 い,初 し て 設 定 し た モ デ ル パ ラ メ ー タの 値 が 真 値 に 近 い 場 合 期 モ デ ル の 固 有 振 動 数,固 有 モ ー ドお よ び モ デ ル パ ラ に も,各 モ デ ル パ ラ メ ー タ の 値 の 組 合 せ に よ っ て は, メ ー タ に 対 す る固 有 振 動 数,固 有 モ ー ドの 一 次 感 度 を 収 束 性 が 悪 く な っ た り,発 散 す る場 合 を経 験 した. 算 出 す る.算 出 した感度 を用 い て モ デルパ ラメ ー タ変 2・1・1従 来の評価関数の 問題点 こ こ で は,収 更 後 の 固 有 振 動 数,固 有 モ ー ドの 変 化 予 測 値{1ぞ}を 予 束 性 が 悪 くな る 原 因 を 考 察 す る.式(2)の 測 し,そ の 予 測 値{R}と は,各 目 標 値{刃`}と す る 次 式 の 評 価 関 数8が の 相 違 を評 価 最 小 とな る よ う に モ デ ル パ ラ メ ー タ を 修 正 す る. モ ー ド次 数 の 固 有 モ ー ドの 各 成 分 の 大 き さ を評 価 量 と し て い る.し た が っ て,式{4)の よ う に感 度 を 用 い て 固 有 モ ー ドの 各 成 分 の 変 化 予 測 を 行 う 際 に は, モ デ ル パ ラ メ ー タ の 変 更 に伴 う固 有 モ ー ドの 各 成 分 の E二告({1ぞ8}一{1ぞ})ア[師 ノ1({1i3つ 一{12})・ ・ … 一・ 響(2) こ こ で、lW}は 評価 関数 で 目標 値 の 相 対 的 重 要 度 を 決 め る 対 角 マ 変 化 は 連 続 的 で な け れ ば な ら な い.そ こ で,物 理 モ デ ト リ ッ ク ス で あ り,そ の 要 素 蹴,は 評 価 関 数 を 無 次 元 ル の 固 有 モ ー ド{以 後,計 算 値 と呼 ぶ)を 規 準 化 す る場 化 す る ため 合 に は,各 モ ー ド次 数 ご と に 目 標 値 と対 応 す る測 定 位 翫 冨 の,!擢2{3) と な る.こ 置 の ど こ か の 成 分 の 値 が,常 こ で,",は ゴ番 め の 目標 値 の 相 対 釣 重 み 係 に 目標 値 と同 じ にな る よ う に 計 算 値 を 規 準 化 す る こ と で 連 続 関 数 と し て い る. 数 で あ る.ま た,初 期 モ デ ル の 固 有 振 動 数,固 有 モ ー ド し か し,固 有 モ ー ドを 規 準 化 す る と き の 基 準 点 の 選 び を 要 素 と す る ベ ク トル を{R氏1},固 有 振 動 数,固 有 モ ー 方 に よ っ て,固 有 モ ー ドの 各 成 分 の 変 化 の 様 子 が 異 な ドの 感 度 を 要 素 と す る 感 度 マ ト リ ッ ク ス を 【Sl,各 モ る.特 に,固 有 振 動 数 の 小 さ い 順 に 目標 値 と計 算 値 の デ ル パ ラ メ ー タ の 変 更 率 ハ を 要 素 と す るベ ク トル を 固 有 モ ー ドを 対 応 させ よ う と して も同 じ 変 形 形 態 の 固 {γ}と す る と,予 測 値{π}は 有 モ ー ドが 対 応 し な い,あ る い は 固 有 モ ー ドの 形 状 は 次 式 で 求 め ら れ る. {κ}={1{㌔}十[s】{γ}(4) 式(2)の 評 価 関 数Eを 近 い が 成 分 が 最 大 と な る 測 定 位 置 が 異 な る 場 合 に は, 最 小 に す る 変 更 率{γ}は,式 代 入 し,Eを 基 準 点 の 選 び 方 に よ っ て 固 有 モ ー ドの 各 成 分 の 変 化 の (4)を 式(2)に 最 小 に す る 条 件 δEノ δ{γ} =Oか ら 導 出 され た 次 式 に よ り求 め る こ とが で き る . {γ}=(〔5]ア 【哨[s〕)-1 た. そ の 場 合 の 問 題 点 を,図1の5自 X[S】7{晒 〆 】{{1そ つ 一{ノ守o})(5} を 用 い て 説 明 す る,例 ま た,修 正 後 の モ デ ル パ ラ メ ー タ の 値 は,式{5}で 算 出 され た γ趣を 次 式 に 代 入 す る こ と で 算 出 で き る. ゴ轟=ゴ幽o(1+γ^)(6) こ こ で,4岨,4凸 様 子 が 大 き く異 な り,収 束 性 が 悪 くな る場 合 が 見 られ は同 定 し よ う とす る 々番 め の モ デ ル パ ラ メ ー タ の 初 期 値 お よ び 修 正 後 の 値 を 表 す. 更 し た と き の4次 由 度 ば ね ・質 量 系 と し て,図1の 質 量5だ け を変 固 有 モ ー ドの 各 成 分 の 変 化 の 様 子 を 図2に 示 す.た 1,3の 振 幅 が 目 標 値 と 同 じに な る よ う に 計 算 値 を 規 準 だ し,図2(a},(b}は,そ 化 し た 場 合 で あ る.図1の れ ぞれ 質 量 モ デ ル で は,質 量5を1.5 倍 に す る と 振 幅 が 最 大 と な る位 置 が 質 量5か ら質 量 弓 に 変 わ る た め,基 準 点 の 選 び 方 に よ り各 成 分 の 変 化 の ノ〃 司0 2000N滞 o 10-1D 01. 多 状 況 が 大 き く 異 な る.こ れ に 伴 い,一 次 感 度 だ け で 予 測 で き る 各 成 分 の 範 囲 も異 な っ て く る.そ 2(a)あ 3。。。N設h°kg る い は(b)の こ で 、図 固 有 モ ー ドの 変 化 の 状 況 を も と に 固 有 モ ー ドの 一 次 感 度 を 用 い て モ デ ル パ ラ メ ー タ を 4㎝N,凝 5(脚 ゜°kq 浩ミ ㈹ ㎏ 》⇒ ζ 〉 1413・kg< 4000N'mく 繍 磯 .1騨 !/! 〉 く 質世5をL5倍 図15自 に変更 由 度 ば ね ・質 量系 と4次 固 有 モ ー ド 修 正 す る 場 合 に,質 量1,3の 用 す れ ば,よ ど ち ら を 基 準 点 と して 採 り適 正 な 質 量5の す る こ と は 困 難 で あ る.な て も,質 量1,3を 修 正 が 行 え る か を判 断 お,4次 固 有 モ ー ドに 関 し 基 準 点 に 選 ぶ 場 合 だ け で な く,他 の 質 量 を 基 準 点 に選 ぶ こ と も 可 能 で あ り,他 の モ ー ド次 数 に お い て も基 準 点 を 選 ぶ 場 合 に は 同様 な 問 題 が 生 じ る. した が って,固 有 モ ー ドの各 成 分 の 大 きさ を評 価 関 に よ っ て 定 ま る.し た が づ て,両 固 有 モ ー ドの 成 分 比 数 とす る場 合 に は,一 次 感度 を用 い て モ デ ルパ ラ メー の 違 い は ベ ク トル の 方 向 の 違 い と な り,両 固 有 モ ー ド タの適正 な 修正 を行 うため の基 準点 の 選 び方,す なわ の 相 違 は 両 ベ ク トル の な す 角 の で 評 価 で き る こ と に ち規 準化 の方 法 に優 昧 な点 が あ り,基 準 点 の選 び方 が な る.角 悪 い と収 束 性 が悪 くな る可能 性 が 大 きい.し か し,収 度 が 零 とな れ ば 両 固 有 モ ー ドが 一 致 す る. こ の よ う な 固 有 モ ー ドの 評 価 量 と し て,任 意 に 規 準 束性 を改 善 す るた め の適 正 な規 準化 の 方 法 を定 め る基 化 さ れ た 二 つ の ベ ク トル{A腸,{莇}間 準 は見 つ け に くい. よ り一 義 的 に 定 ま る次 式 のCOSユ6レ を 用 い る. 2・2新 た な固 有 モ ー ドの評 価 量 を評 価 閲 数 と した 同 定法 の なす角 傷 に cos2畠=わ =({ハ 弓`}7{1嚇)21({1W}711鴇 りx{1嚇 従来 の 同定 法 で 用 い た評 価 関 数 で は,固 有 「鵬}) モー ドの 各成 分 の大 きさ を評 価 す るた め に規 拳化 の問 (7) 題 が 生 ず る.そ こで,各 成 分 ご とに 目標 値 と計 算値 の 評 価 関 数 £'は,式(7}で 定 義 し た 固 有 モ ー ドの 評 価 値 を比 較 せ ず に 目標 値 と計算 値 の固 有 モ ー ドの相 違 が 量 を 用 い て,目 擦 値{Pう と 計 算 値{P}の 評価 で きれ ば,こ れ までの 問題 が解 消 され,収 束 性 が か ら な る次 式 とす る, 、E'=({P"レ{PD7[明({、P`}一{PD… 改 善 で きる と考 え られ る.以 下 で は,規 準 化 の 方法 に とらわれ な い 固有 モー ドの評価 量 を評 価 関 数 とす る同 こ こで,{P"}と{P}は 残 差二乗和 … …(8) 次 の よ う な ベ ク トル で あ る. 定 法 を提 案 す る. 捗に訟1:::鴛::読ll:1器}一{9) 固有 モー ドは,成 分 の 大 き さが定 ま らず と も、成 分 の比 が定 まれ ば一義 的 に決 定 で きる,そ こで,成 分 比 た だ し,カ は 固 有 振 動 数,η だ けを用 い た固有 モー ドの 評価 方 法 を考 える.例 えば、 評 価 量 で あ る。 同 定 計 算 手 順 は 従 来 と同 じ方 法 を 用 い 図3左 に 示 す3自 由度 系 の ノ次 固有 モ ー ドに お い て、 る こ と が で き る.す 実験 値 の固 有 モ ー ド{1鴨り,計 算 値 の固 有 モ ー ド{劫} るco52θ 」の 一 次 感 度 を 用 い て,{1鴨,{ハ と も,質 量1の 振 幅 をX方 向 に,質 量2の 振 幅 を γ れ るco52畠 方 向 に.質 量3の 振幅 をZ方 向 に取 った結 果,両 固 有 は 式{7)の 固 有 モ ー ドの な わ ち,モ デ ル パ ラ メ ー タ に 対 す の 値 が1に ㌧ ・}から 得 ら 近 づ くよ うに モ デ ル パ ラメ ー タの 値 を 修 正 す る こ と に な る. モ ー ドが図3右 に示 す ベ ク トル とな った とす る.こ こ 2・2・1新 で,ペ タ トル の 方向 は それ ぞ れの 固 有 モ ー ドの 成 分 比 に及 ぼ す影響 た な 固 有 モ ー ドの 評 価 量 が 収 束 性 の 改 善 式(7)の 固 商 モ ー ドの 評 価 量 を 用 い る こ とが,収 束 性 の 改 善 に ど の よ う に 役 立 つ か を 検 討 す る.そ 芒 .曳to こ で,cos2θ 、の 値 が モ デ ル パ ラ メ ー タ の 変 更 に よ り どの よ うに 変 化 す る か を調 ぺ る.co52畠 をモ デ .2 ル パ ラ メ ー タ の 変 更 率 γ幽で 偏 微 分 す る と 次 式 と な = きoo る, 塞 δco83θ レ!δ γ凸;2ル 屯 爵 萄 xiBXD/C-831(〆 一1.08 1.0152-025 箋 こ こ で, 野しatio0 5/m瓢 儲)質 量1を 基準点 、4={劫 図2質 三 止3を 基準点 量5の 変 更 に 伴 う4次 固 有 モー ドの 変化 臼 ハる},8={ハ C・=〔 劫 り ア{ハ凋,D={ハ E={ハ {b}質 隻xC))(E】(10) 弓}7{δ{現}!δ γ霞} 弓り 『{私} 与り γ{δ{助 〕1δγ轟} な お,δ{1鴇・}1δ み は 既 報 吻で 用 い た モ デ ル パ ラ メ ー タ 正 す る 場 合 に も,必 ず 安 定 に モ デ ル パ ラ メ ー タ の 真 値 の 変 更 率 γ凸に 対 す る 固 有 モ ー ドの 一 次 感 度 で あ る, を 求 め る こ と が 可 能 と な る。 今,計 算 値 の 固 有 モ ー ドが モ デ ル パ ラ メ ー タ 変 更 後 も ま た,図4に お い て,質 量5の 変 更 がco♂ 島の変化 長 さ が 変 化 し な い よ う に 規 準 化 す る場 合 を 考 え る.そ に 及 ぼ す 影 響 が モ ー ド次 数 ご と に 異 な る,す な わ ち 質 の 場 合 に は,固 有 モ ー ドの 長 さ の 感 度 〔 式(1ω のE〕 が 量5の 零 と な る た め, の モ ー ド次 数 へ の 影 響 が 少 な い こ と が わ か る.し か し, B21(〆1xC)xE=O(11) と な る.ま xCの た,Schwa㍑ お い て も、Eを モ ー ド へ の 影 響 が 大 き く、他 他 の 質 量 で は 質 量5と の 不 等 式 か ら導 か れ る β ≦ 君 両 辺 を 変 更 率 為 で 偏 微 分 し た2(BxD)≦2(C Xβ)に 変 更1ま3次,4次 の 傾 向 は, モー 一 ド次 数 ご と に 異 な る 傾 向 と似 て い る.こ 零 と お く と, り,cos2島 (12) B>く 」0≦0 は 傾 向 が 異 な る,二 モ デ ルパ ラメ ー タが 固有 振動 数 の 変化 に与 え る影響 が の こ とよ を 目 標 と す れ ば,モ デ ル パ ラ メ ー タ 間 の 目 標 に 対 す る相 関 を,従 来 の 固 有 モ ー ドの 各 成 分 を 目 標 の 関 係 が 導 か れ る.し た が っ て,変 更 率 ハ が 零 か ら 増 と す る 場 合 に 比 べ て 小 さ くす る こ と が で き る と考 え ら 加,す れ る. な わ ち モ デ ルパ ラ メ ー タ の 値 が 真 値 か ら 大 き く な る 場 合 に は,式(10)に お い て 式(11),(12)を 考慮す る こ と に よ り、 量 を 目 標 と す る と, δco呂ユθ,1δ γ馬≦o(13) の 関 係 が 導 か れ る.式(13〕 るCOげ 以 上,本 報 で 同 定 の 評 価 に 用 い る 固 有 モ ー ドの 評 価 ① 目標 と する 固有 モ ー ドと計算 値 の 固有 モ ー ドの相違 は モ デ ルパ ラメー タ に対 す θ」の 感 度 が 常 に 負 と な り,co♂ のが単調 に減 を 評 価 す る と き の 曖 昧 さ が な くな る た め,固 有 モ ー ド の 評 価 が 一 義 的 に で き る, 少 す る こ と を意 味 す る.ま た,変 更 率 γ幽が 零 か ら負 の ② 評 価 量 がOか ほ う に 増 加,す 次 数 間 の 数 値 の ば ら つ きが 小 さ くな る. な わ ちモ デ ル パ ラ メ ー タ の 値 が 真 値 か ら 小 さ くな る 場 合 に も,式(1ω,(12)に 負 で あ る こ と を 考 慮 す れ ば,co♂ す る こ とが わ か る.し ③ モ デ ルパ ラメー タの変 化 に対 して、モデ ル パ ラ メー 劇 の 値 が 単調 に減 少 タ の 真 値 を頂 点 と す る 上 に凸 な 二 次 関 数 的 に 振 る 舞 う た が っ て,モ デ ル パ ラ メ ー タ の れ 島 は モ デ ル パ ラ メー タ の 真 値 を 頂 目標 値 を 設 定 で き る た め.個 々 の 目標 値 に 対 す る収 束 安 定 性 が 向 上 す る. ④ モ デ ル パ ラ メ ー タ間 の 目標 値 に 対 す る 相 関 が 小 さ く 点 と す る 上 に 凸 な 関 数 と な り,感 度 は 必 ず モ デ ル パ ラ な る た め,式(5)の メ ー タ の 正 し い 修 正 方 向 を 示 す こ とが わ か る.こ こ で, つ な が る. 図1の5自 な どの 利 点 が あ る.こ のcos呈aの 由 度 系 を例 に と り,質 量5を 変 化 の 状 況 を 図4に 範 囲 の 数 値 と な る た め 、モ ー ド お い て 、γ禽が 値 が 真 値 よ り小 さ い 場 合 の 感 度 は 常 に 正 とな る,こ ら の こ と よ り,co♂ ら1の 変 更 した と き 示 す.図4に 係 数 行 列[S}ア[朗 【SIの 対 角 化 に れ ら の 事 項 は,同 定 計 算 に お け お い て, る 収 束 の 不 安 定 性 を抑 制 す る 効 果 σ,があ り,従 来 の 同 各 モ ー ド と も上述 の よ う に 上 に凸 な 関 数 とな っ て お 定 法 に 比 べ て,収 束 性 を 向 上 さ せ る こ と が 見 込 ま れ り,変 化 の 状 況 は な め ら か で 、質 量 の 変 更 が2倍 る. 程度 ま で は 二 次 関 数 的 な 振 舞 い を し て い る こ とが わ か る. 翫 数値実験およσ考寮 し た が っ て,一 次 感 度 を 用 い て モ デ ル パ ラ メ ー タ を修 梶 案 し た 同 定 法 が,従 来 の 同 定 法 に 比 べ て 収 束 性 を 改 善 で き る こ と を 確 認 す る た め に,数 値 実 験 を 行 っ た. 3・1ば ね ・質 量 系 の 場 合 図1に 示 す5自 由度 の ば ね ・質 量 系 を 取 り上 げ,検 討 を行 う. 表1の 左 か ら2列 め に数 値 実 験 の 目 標 値 を 得 る た め の 固 有 値 解 析 に 用 い た 諸 元 を 示 す.図1の モ デ ルにお い て は 、5個 の 固 膏 振 動 数 と固 有 モ ー ドの 固 有 ペ ア が 得 ら れ る.そ こ で,そ れ らの 値 を 目標 値 と し て 同 定 計 算 を 行 っ た.同 定 計 算 の 初 期 モ デ ル は 、表1の 3列 め か ら5列 め に 示 す3通 す る モ デ ル パ ラ メ ー タ は,モ め に 既 知 と し た 質 量 規1(固 り と した.同 左か ら 定 しようと デル の大 き さを定 め るた 有値 解 析 に用 い た質 量 勘 の値 を用 い る〉 を除 い た残 りの 質 量 職 か ら 肋 お よ 違 い は あ る が 安 定 に 収 束 し て お り,従 来 法 に比 べ て も び ばね定 数 島 か ら 鳥 の計10個 少 な い 反 復 計 算 回 数 で 残 差 が 小 さ く な っ て い る こ とが 図5に,同 とした. 定 計 算 の反 復 計 算 回 数 と残 差 の関 係 を示 す.た だ し,縦 軸 は 目標 値 と計 算値 の 各成 分 の残 差 の 絶 対値 和 で あ る.図5に な評価 関 数E'を お い て,実 線 が 提 案 した新 た 用 い た同 定 法(以 後 、提案 法 と呼 ぶ) の 場 合 で,破 線 が 従 来 の評 価 関 数Eを わ か る. 3・2ば ね 要 素 で 支 持 され た 剛 体 系 の 場 合 モ デ ル と し て,図6に 次 の 示 す剛体 が ば ね で支 持 され た モ デ ル を取 り上 げ る.3・1節 で 扱 っ た ば ね ・質 量 系 モ デ ル 用 いた同定法 に お い て は,モ デ ル パ ラ メ ー タ の 初 期 値 が 真 値 と か な 似 後,従 来 法 と呼 ぶ)の 場 合 で あ る.な お,両 同定 法 り異 な っ て い て も,固 有 モ ー ドの 順 序 が 入 れ 替 わ っ た の違 い を評価 す るた め,提 案 法 の 残 差 の数 値 は式(7) りす る こ とが な く、同 定 計 算 に お い て は 比 較 的 取 扱 い の 固有 モー ドの評 価 量 で は な く,従 来 法 と同 じ固 有 モ が 容 易 な モ デ ル で あ る.し か し,図6の ー ドの各 成 分 の残 差 の絶 対 値和 を用 い た. 性 モ ー メ ン ト,重 心 位 置 の 座 標 な ど が モ デ ル パ ラ メ ー 図5に お い て,従 来 法 で は収 束 の状 況 が 初 期 モ デル モ デ ル で は,慣 タ に 含 まれ る た め,固 有 モ ー ドの 順 序 が 入 れ 替 わ り や に よ って大 き く異 な る こ とが わ か る.特 に,モ デ ルB す く.そ れ だ け 同 定 が 困 難 と な る.し た が っ て,初 期 モ の場 合 に は反復 回 数2か ら4に おい て残 差 が増 大 して デ ル の 値 の 設 定 に 制 約 が 加 わ る た め,初 期 モ デ ル の 値 お り,不 安 定 な 収束 状 況 とな って い る.し か し,提 案 法 に 依 存 し に く い 同 定 法 が 必 要 と な る, の場 合 に は,三 つ の モ デ ル に お い て,多 少 の収 束 性 の 図6の モ デ ル に お い て も,表2の 左 か ら2列 めの値 を 用 い て 固 有 値 解 析 を 行 い,同 定 計 算 の 目 標 値 を 得 表1ば 醜od¢1 ㎜ εtεr mll㎎ た.初 ね ・質 量 系 の モ デ ル パ ラ メ ー タ E翼 ● 轄t 7邑luo 示 す3通 1nit軸1服hに 閥od璽1Al岡odo1塾1閥odεlO lo』 期 モ デ ル は、表2の り と し た.な 表2剛 30.o 皿3k感 40.o 30.o 46』 」0』 麗odo1 E冨 邑ot m4㎞ 30』 15』 35』 40』 園lr四 已t¢r 鞄luo m51昭 20.o 40』 20』 50』 m㎞ k1眈 2』Xlo3 4』 翼103 2』XIO, 3』XIO3 1.δ 箕103 2.5xlo3 5』xloヨ 4.0区103 8』Mo3 看.5翼1窃3 4』X103 ■ k5賄 h囁 5.OXIO3 4』XIO, 團ノ■ コ』X163 めに モ デ ル に お い て もモ 体 系の モ デ ル パ ラメ ー タ 血2k匹 猛4即 お,図6の め か ら5列 10』 15』 擁 即騨 k3眈 左 か ら3列 2.5刈03 6.Ox103 1.5×103 35.o 20』 5.5×10コ 4.Ox103 5』X10, 自.OX103 2.OXIO3 田 ヨ 一 一 _一 一 1"it二 亀1V亀luε 質odo1Al鮪od¢i8i隠odolC 50』 50』 o.肥 0,065 o.罰 o,価 o.肥5 o.50 1.壮5 0,105 0,025 o.㎝25 o.10 o.伍o o』15 o.oo15 0.10 o.鮨o o.ユo km{四1■ o,030 L55×105 o.oo80 K昌y囲/■ k皿 艮!瞳 k財 眈 J翼1昭 ・ 圃守 J,1蟷12 Jz㎞ 皿ε 3.0刈03 1.o翼 一 ・ ■を 咀 y8咀 2ε 咀 o.岡 0.α50 155買104 LOOX105 1,70冨10ら 1.7似10・ 1.㎝ 5』OX105 5.00xlo噂 5』o翼lo, 瓦105 L75翼105 L75翼104 1.oo翼loら kbz軒/■ 4.50MIO筋 4.50xlo4 5.00)dO5 k㏄ 5』9xlo5 匠.00N104 5.㎝Xloら k〔h胆1● 1』D翼105 1.90×104 1.oo}doう hdz四 5.50刈05 5,50x104 5.oo刈05 世〆塁 〆置 ● 0トXYZ:Fi其edcoor〔 圭ir】at巳 6一π帽Y弓ヱ信:C¢ 匝t¢rofgraviし 図6剛 アcoordinate 体 系モデル デ ル の 大 き さ を定 め る た め に 全 体 の 質 量 餌 は 既 知 と し た.そ し て 、残 りのXσ,γ メ ン トみ,ん,ゐ,重 σ,2F°軸 回 りの 慣 性 モ ー 心 位 置 の 座 標 翫,弥,砺 持 ば ね の ば ね 定 数 々ロ,々鯉,々皿,々的,妬,妬,々 計M個 お よ び支 血,々由 の の モ デ ル パ ラ メ ー タ を 同 定 し た, 図7に,反 復 計 算回 数 と残差 の絶 対 値 和 の 関 係 を示 す.図7に お い て,従 来 法 で は 前 節 と同 様 に収 束 の 状 況 が 初 期 モ デ ル に よ っ て 大 き く異 な る こ と が わ か る, モ デ ルAの 初 期 値 は 表2に 示 す よ う に 真 値 の1110に 設 定 し て お り,こ の よ う に 全 体 が 同 率 で 縮 小,拡 大 さ れ た 場 合 に は,モ ー ドの 順 序 が 入 れ 替 わ り に く く,従 来 法 で も収 束 が 速 い.し か し,モ デ ルB,Cの よ う に, 初 期 値 を任 意 に 設 定 した 場 合 に は,収 束 性 が 悪 くな る 二 と が あ り,モ デ ルBの 場 合 に は20回 要 して も収 束 に は 至 ら な か っ た.し の反復計 算 を か し,提 案 法 を 用 い た 場 合 に は,三 つ の モ デ ル に お い て ,多 少 の 収 束 性 の 違 い は あ る が 安 定 に収 束 し て お り,初 期 モ デ ル の 影 響 を 受 け に くい こ とが わ か る.ま た,従 来 法 に 比 べ て も少 な い 反 復 計 算 回 数 で 残 差 が 小 さ く な っ て お り,提 4.結 関 数 を用 い た同定 法 を提 案 した. (3)提 案 法 の 有 効 性 が 立 証 で き る. 書 案 した同定 法 の 有効 性 を確 認 す るた め に数 値 実 験 を行 い,従 来 の 同 定 法 に比 べ て,収 束 速 度,収 束 安 定性 が 大幅 に向 上 す る こ とを示 した. 感 度 解 析 を用 い た 振 動 系 の モ デ ル パ ラ メ ー タ 同 定 に 文 お い て,反 復 計 算 の 収 束 性 を 改 善 さ せ る 方 法 に つ い て 検 討 し た.そ (1)同 の 内 容 を ま とめ る と次 の よ う に な る. 定 の 評 価 に 固 有 モ ー ドの 個 々 の 成 分 を用 い る こ と が,収 束 性 を 悪 く す る 要 因 の 一 つ で あ る こ と を 示 し た. (2)収 束 性 を 改 善 さ せ る た め に は,各 モ ー ド次 数 ご と に 固 有 モ ー ドペ ク トル の 方 向 を 評 価 す る 量 を 評 価 関 数 に用 い る ことが有 効 で あ る こ と を述 べ、その 評 価 〔1〕 長 舩,モ 献 ー ド解 析.{19呂5転1{}5培 風 館, 2} 藤FH・ 3〕 鞍 谷 ・韓 川 ・ 沖 田.磯 論,5牛506.C(1988}.脳 1} 粧 谷 ・藤h「 ・ 沖 田.機 論,55-512.C{1989},84⑪. 卸 中 嗣 ・高 畑 、機 論,岳4-51)7,Cq曾 61 大 熊 ・長 祉.文 献15)け2557ヘ に1 中 「11・小 柳,晶 小 12」.吏 新 田 ・沖 田 ・鞍 谷,機 論,与2-476,C〔19邑6},1224. 四. 呂61,3530. ー ジ, 乗`ゴ、 に よ る 実 験 干'一 ク解 極,硅997, 京 メく学 出 販 会.
© Copyright 2025 ExpyDoc