• 検索結果がありません。

連続 体構造物に対する汎 用低次元非線形振動解析

N/A
N/A
Protected

Academic year: 2021

シェア "連続 体構造物に対する汎 用低次元非線形振動解析"

Copied!
4
0
0

読み込み中.... (全文を見る)

全文

(1)

     博 士 ( 工 学 ) 原 田 学 位 論 文 題 名

連続体構造物に対する汎用低次元非線形振動解析 学位論文内容の要旨

近年の工業分野において,高速化,軽量化,エネルギーコスト低減など様々な理由,目 的から,機械の構造部材は厚肉から薄肉へと変わってきており,大振幅振動が生じやすい など,振動が機械に与える影響はますます大きくなってきている.そのため,動力学の応 用による薄肉構造物の定性的,定量的な振動特性を理解することの重要性が更に増してき ている.

  薄肉構造物の曲げ振動において,曲げ振幅が厚さに比べて十分に小さい微小変形を対象 とする場合,面内変形の影響を無視した線形理論を用いて,工学的に十分な精度の結果を 得ることができる.しかし,曲げ振幅が厚さと同程度かあるいはそれ以上になると,面内 変形と面外変形の連成を無視できなくなり,固有振動数が振幅とともに変化することとな るので線形理論では不十分になり,大たわみによる幾何学的非線形性を考慮した解析が必 要となる.分数調波共振,内部共振,結合調波共振,さらにはカオス振動など,線形理論 では予測不可能な非線形挙動が現実の振動問題ではたびたび生じる.故に,構造物の振動 特 性 を 正 確 に 把 握 す る た め に は , 非 線 形 理 論 の 適 用 が 必 須 と な る .   従来,幾何学的非線形性を考慮した連続体の非線形振動に関する研究は数多くなされて いる.解析手法も多数あるが,弱非線形性を有する構造物に対しては,古くから摂動法が 用いられ,強非線形性を有する構造物の動的挙動解析には,調和バランス法が広く用いら れている.しかし,これら解析的を方法は体系的で明確な解を得ることができ,対象の振 動特性を理解しやすく有意義である反面,境界条件や形状によっては,計算が複雑であっ たり解けない場合もあるなどの問題がある.一方,幾何学的非線形性を考慮に入れた有限 要素法(FEM)を用いた 解析も盛 んに行われている.FEMは,どのような境界条件でも解 けるという利点を有する反面,数値計算における時間やコストの面で最良の方法であると は言えない.また,数値解法であるために,解の構造を系統的に把握することは多くの場 合 困難 で ある.最 近では, 形状関数の 改良によ り得られ た低次元FEMモデルやFEMモ デルをべースにして得られた低次元モデルに,摂動法や調和バランス法などの解析的手法 を適用 する手法も提案されてきている.しかし,解析対象によっては,FEMよりは低自 由度だが,それなりに大きい自由度のモデルになる,適用可能な対象が限られているなど の問題も存在する.

  そこで本研究では,有限要素法と理論解析の両者の問題点を互いに補い,かつ解析対象 に制限をもたない方法の提案をめざす,

  本 論 文 は全4章 で 構成 さ れ てい る .以 下にそれ ぞれの章 の概要を 具体的に 示す,

935

(2)

  第1章 は 緒 論 で あ り ,本 研 究 の目 的 と意 義 お よび 各 章の 概 要 につ い て 述べ た .   第2章では,両端固定はりおよび中間支持を有する両端固定はりをモデルに,はりの1 次元非線形モード方程式を導出し,その妥当性を検証した.解析は面外1次モードにっい てのみ行われた.解析の手順は,まず幾何学的非線形性を考慮した変位―ひずみ関係式を 用い,仮想仕事の原理より非線形有限要素法方程式を導いた.得られた質量マトリクス,

線形剛性マトリクスからなる標準固有値問題を解くことにより,線形固有振動数,線形面 内モード,線形面外モードを得た.線形面外1次モードとその他の線形モードの中のある 1つのモードを用いて,非線形有限要素法方程式中の節点量ベクトルをモード座標に展開 することにより2自由度モデルを得た.両端固定はりにっいては十分な精度を持ったモデ ルが得られたが,中間支持を有する両端固定はりについては十分な精度を持ったモデルを 得られなかった.そこで,より汎用性のある導出方法をめざし,非線形有限要素法方程式 にNewton‑Raphson法 を適用し て得られ る変位(NLFEM解)の利用を試みた.その結果,

線 形 面外1次 モ ード と 線 形面 外1次 固有 振 動 数に お けるNLFEM解 の 面 外成分 を0とし たもの(擬似非線形面内1次モード,とする)を用いるのが,最もよいモデルを得る節点 量ベクトルの展開であることが明らかになった.この手法を用いることにより,両端固定 はりについてはより精度の高いモデルを,中間支持を有する両端固定はりについても十分 な精度を持ったモデルを得ることができた.得られた2自由度モデルの面内慣性を無視す ることにより1自由度モデルが得られ,モデルの非線形性を解析するために,摂動法や調 和バランス法などの解析的手法を適用することが可能となった.得られた1自由度モデル に第一種完全楕円積分を適用して得られた背骨曲線と他の解析法によって得られた結果と を比較することにより,得られた1自由度モデルは十分な精度を有していることが確かめ られた.

  第3章で は, 内部共振 のない4辺 単純支持 された板 および4辺 固定された 板,1:3内 部共 振を有する4辺固定された板などを対象に,第2章と同様の手法を用い,板の1次元

(内部共振を有する板にっいては2次元)非線形モード方程式を導出しその妥当性を検証 した.内部共振のない板については,線形面外1次モードと擬似非線形面内1次モードを 用い十分な精度を持った1自由度モデルを得ることが可能であることが示された.内部共 振を有する板については,3C01〜〜cvsの関係にあるので,線形面外1次モードと擬似非線 形面 内1次モードの他に線形面外5次モードと擬似非線形面内5次モードを用い,十分な 精度を持った2自由度モデルを得ることが可能であることが示された.また,極端な辺長 比を 有する長 方形板お よび2種類の 支持条件が混在するL字状板に対して,本手法を適 用した場合においてもよい結果を得られ,本手法は高い汎用性を有していることが示され た,

  第4章 は 結 諭 で あ り , 本 研 究 を 通 じ て 得 ら れ た 成 果 を 取 り ま と め た .

936 ‑

(3)

学位論文審査の要旨 主査

副査 副査 副査

教授 教授 教授 助教授

山田 佐々木 三上 小林

学 位 論 文 題 名

    兀 一彰     隆 幸徳

連続 体構造物に対する汎 用低次元非線形振動解析

  近 年 の 工 業 分 野 に お い て 、 高 速 化 、 軽 量 化 、 工ネ ル ギ ーコ ス ト 低減 な ど 様々 な 理 由 、 目 的 か ら 、 機 械 の 構 造 部 材 は 厚 肉 か ら 薄 肉 へ と 変 わ っ て き て お り 、 大 振 幅 振 動 が 生 じ や す い な ど 、 振 動 が 機 械 に 与 え る 影 響 は ま す ま す 大 き く な っ て き て い る 。 そ の た め 、 薄 肉 構 造 物 の 振 動 特 性 を 、 動 力 学 の 応 用 に よ っ て 定 性 的 、 定 量 的 に 把 握 す る こ と の 重 要 性 が こ れ ま で 以 上 に 増 し て き て い る 。

  幾 何 学 的 非 線 形 性 を 考 慮 し た 連 続 体 の 非 線 形 振 動 に 関 す る 研 究 は 数 多 く な さ れ て い る 。 解 析 手 法 も 多 数 あ る が 、 弱 非 線 形 性 を 有 す る 構 造 物 に 対 し て は 、 古 く か ら 摂 動 法 が 用 い ら れ 、 強 非 線 形 性 を 有 す る 構 造 物 の 動 的 挙 動 解 析 に は 、 調 和 バ ラ ン ス 法 が 広 く 用 い ら れ て い る 。 し か し 、 こ れ ら 解 析 的 な方 法 は 体系 的 で 明確 な 解 を得 る こ とが で き 、 対 象 の 振 動 特 性 を 理 解 し や す く 有 意 義 で あ る 反 面 、 境 界 条 件 や 形 状 に よ っ て は 、 計 算 が 複 雑 で あ っ た り 解 け な い 場 合 も あ る な ど の 問 題 が あ る 。 一 方 、 幾 何 学 的 非 線 形 性 を 考 慮 に 入 れ た 有 限 要 素 法(FEM)を 用 い た 解 析 も 盛 ん に 行 わ れ て い る 。FEMは 、 ど の よ う な 境 界 条 件 で も 解 け る と い う 利 点 を 有 す る 反 面 、 数 値 計 算 に お け る 時 間 や コ ス ト の 面 で 最 良 の 方 法 で あ る と は 言 え な い 。 ま た 、 数 値 解 法 で あ る た め に 、 多 く の 場 合 、 解 の 構 造 を 系 統 的 に 把 握 す る こ と は 困 難 で あ る 。 最 近 で は 、 形 状 関 数 の 改 良 に よ り 得 ら れ た 低 次 元FEMモ デ ル やFEMモ デ ル を べ ー ス に し て 得 ら れ た 低 次 元 モ デ ル に 、 摂 動 法 や 調 和 バ ラ ン ス 法 な ど の 解 析 的 手 法 を 適 用 す る 手 法 も 提 案 さ れ て き て い る 。 し か し 、 解 析 対 象 に よ っ て は 、FEMよ り は 低 自 由 度 だ が 、 そ れ な り に 大 き い 自 由 度 の モ デ ル に な る 、 適 用 可 能 な 対 象 が 限 ら れ て い る 、 な ど の 問 題 も 存 在 す る 。   本 論 文 で は 、 有 限 要 素 法 と 理 論 解 析 の 両 者 の 問 題 点 を 互 い に 補 い 、 か つ 解 析 対 象 に 制 限 を も た な い 解 析 手 法 に 関 す る 基 礎 的 な 研 究 を 行 な っ た も の で あ り 、 そ の 主 要 な 成 果 は 次 の2点 に 要 約 さ れ る 。

  (1)両 端 固 定 は り お よ び 中 間 支 持 を 有 す る 両 端 固 定 は り を モ デ ル に 、 は り の1次 元 非 線 形 モ ー ド 方 程 式 を 導 出 し 、 そ の 妥 当 性 を 検 証 し た 。 解 析 は 面 外1次 固 有 振 動 数 付

(4)

近についてのみ行った。解析の手順は、まず幾何学的非線形性を考慮した変位一ひず み関係式を用い、仮想仕事の原理より非線形有限要素法方程式を導いた。得られた質 量マト1J クス、線形剛性マトリクスからなる標準固有値問題を解くことにより、線形 固有振動数、線形面内モード、線形面外モードを得た。線形面外 1 次モードとその他 の線形モードの中の特定の1 つのモードを用いて、非線形有限要素法方程式中の節点 量ベクトルをモード座標に展開することにより 2 自由度モデルを得たが、より汎用性 のある導出方法をめざし、非線形有限要素法方程式に Newton 一 Raphson 法を適用し て得られる変位 (NLFEM 解)の利用を試みた。その結果、線形面外1 次モードと線形 面 外 1 次 固有 振動 数にお ける NLFEM 解 の面 外成 分を0 と した もの( 非線 形面 内1 次 モードとする)を用いるのが、最もよいモデルを得る節点量ベクトルの展開であるこ とが明らかになった。この手法を用いることにより、両端固定はりについてはより精 度の高いモデルを、中間支持を有する両端固定はりについても十分な精度を持ったモ デルを得ることができた。また、得られた2 自由度モデルの面内慣性を無視すること により 1 自由度モデルが得られ、モデルの非線形特性を解析するために、摂動法や調 和 バ ラ ン ス 法 な ど の 解 析 的 手 法 を 適 用 す る こ と が 可 能 と な っ た 。   (2) 内部 共振の ない 4 辺単純支持された板および4 辺固定された板,1:3 内部共振 を有する4 辺固定された板などを対象に、板の 1 次元(内部共振を有する板について は 2 次元)非線形モード方程式を導出しその妥当性を検証した。すなわち、内部共振 のない板については、線形面外1 次モードと非線形面内 1 次モードを用い十分な精度 を持った1 自由度モデルを得ることが可能であることを示した。 3 り1 与の 5 の関係に ある内部共振を有する板については、線形面外1 次モードと非線形面内1 次モードの 他に線形面外5 次モードと非線形面内5 次モードを用いることにより、十分な精度を 持 っ た 2 自 由 度 モ デ ル を 得 る こ と が 可 能 で あ る こ と を 示 し た 。    これを要するに、著者は、連続体構造物に対する汎用低次元非線形振動解析法の提 案とその有用性を論じたものであって、非線形振動に関して有益な知見を与えており、

振動工学の進歩に貢献するところ大なるものがある。よって著者は、・北海道大学博士

   (工学)の学位を授与される資格あるものと認める。

参照

関連したドキュメント

This paper considers a possibility of decision whether the robot hand is having a correct work or not by using the analysis of the mechanical vibration of robot that is doing

This paper examines the diffuse bifurcation mode of a hollow circular cylinder specimen consisting of a non-coaxial Cam-clay model.. The specimen undergoes a

spread takes small values for fast time varying pole. p osition, and large values for slow time

振動流中および一様 流中に没水 した小口径の直立 円柱周辺の3次 元流体場 に関する数値解析 を行った.円 柱高 さの違いに よる流況および底面せん断力

Series of numerical analysis to estimate structural frequency and modal damping were conducted for a two-dof model using the simulated external forces induced by impulse force and

は ) 変調が激し $A\backslash$ ときに, 小さ $A\backslash$ スケールの砕波 (spilling breaker) や表面張力波が確認されて $A\backslash$ る [1].. Dysthe 方程式によると ,

せん断帯の数値解析は、材料の非線形性だけでなく初期形状の非対称性や材料の非均質性

しかしながら,式 (8) の Courant 条件による時間増分