JSME-JT
|
|
- ちかこ こやぎ
- 4 years ago
- Views:
Transcription
1 Advance Publication by J-SAGE 日本機械学会論文集 ransactions of the JSME (in Japanese) DOI: /transjsme Received date : 13 November, 2017 Accepted date : 24 April, 2018 J-SAGE Advance Publication date : 7 May, 2018
2 柔軟マルチボディモデルを用いた移動式クレーンのロープの慣性の影響の検証 CHALERMPONG Kolawach *1, 山浦弘 *2, 原謙介 *2 Investigation of influence from inertia of hoisting rope on mobile crane using flexible multibody model Kolawach CHALERMPONG *1, Hiroshi YAMAURA *2 and Kensuke HARA *2 *1,*2 Department of Mechanical Engineering, School of Engineering, okyo Institute of echnology Ookayama, Meguro-ku, okyo , Japan Abstract his research studied an influence of the inertia of the hoisting rope on the dynamics of a mobile crane system. wo dynamics models of the mobile crane were formulated with the inertia included rope model and the massless rope model, and their dynamics response were compared. he floating frame of reference formulation was used in formulating the model of the deformable components. In the simulations, we variated the length of the rope and the mass of the hoisted load then applied the slewing motion to the boom for 10 seconds and let the system underwent free vibration for another 10 seconds while the boom angle was fixed. he deformation of the rope and the comparison of the hoisted load motion between two models were considered. he results illustrated that the rope model with included inertia was able to illustrated its bending during the crane operation and this deformation increased with the length of rope. It appeared that the deformation of the boom tip had influences on the deformation of the rope. In the motion of the hoisted load, the time history data illustrated that the inertia included rope model had slightly higher frequency than the massless rope model during the free motion, and the trajectory also showed the slightly larger sway. In addition, the total calculation time of the inertia included rope model was around 13 times longer. he difference in the motion of the hoisted load from both models increased with the length of rope but decreased when the mass of load increased. In conclusion, the model with included inertia of rope showed different response to the massless rope model at long rope and light load condition, and this difference was decreased at short rope and heavy load condition. Keywords : Flexible multibody dynamics, Floating frame of reference formulation, Mobile crane, Hoisting rope, Dynamics model 1. 背景クレーンの吊荷の振動制御手法の検討やリアルタイムシミュレーションの実施のため, 実機の動特性を詳細に表現でき, かつ高速に計算機シミュレーションが可能なクレーンの動力学モデルの開発が求められている. 従来の研究では, アウトリガー, ブームおよびロープの柔軟性を考慮した移動式クレーンの動力学モデルが提案されてきた (Kiliçaslan et al., 1999)(Sun and Kleeberger, 2003) (Maczyski and Wojciech, 2003) (Posiadala and Cekus, 2008) (rąbka, 2014). これらの多くはロープの慣性を省略し, 無質量柔軟体とみなしているため, ロープの変形は軸方向ひずみしか考慮されていない. しかし, 実際のロープには曲げやねじり変形が生じる場合もある. また, ロープの慣性を考慮した動力学モデルも存在している (rąbka, 2014, 2016) が, ロープ慣性の影響は十分に示されていない. そのため, ロープの慣性が移動式クレーンの動力学モデルに与える影響を調べるとともに, ロープの慣性を考慮したモデルと考慮しないモデルの動特性の差を検証することを, 本研究の目的とする. No [DOI: /transjsme ] 本論文は,Dynamics and Design Conference 2017 講演論文集, No.718 の掲載内容に基づいた論文である. *1 学生員, 東京工業大学大学院工学院機械系 ( 東京都目黒区大岡山 ) *2 正員, 東京工業大学大学院工学院機械系 of corresponding author: chalermpong.k.aa@m.titech.ac.jp
3 本研究ではロープの慣性を考慮する条件と考慮しない条件の移動式クレーンの動力学モデルを作成し, 二つのモデルの振動応答を比較する. 慣性を考慮したロープモデル, およびクレーンの柔軟体モデルは Floating Frame of Reference formulation( 浮動基準枠法 ) によって作成する (De Veubeke, 1976)(Shabana, 1997). この方式は Park の海上の浮体式クレーンの柔軟ブームの解析にも利用された実績がある (Park et al., 2011). 二つのモデルのブームに回転運動を入力し, シミュレーションを行う. シミュレーション結果として慣性を考慮したロープの曲げ変形, 二つのロープの振り子運動モードの比較, および二つのロープのモデルで計算した吊荷運動を示す. 最後に, 移動式クレーンのロープの慣性を考慮した動力学モデルと考慮しない動力学モデルの得失について示す. 2. ロープのモデル 2 1 慣性を考慮するロープのモデル本研究では浮動基準枠法 (Floating Frame of Reference formulation (FFRF) )(De Veubeke, 1976)(Shabana, 1997) と有限要素法を用いてロープのモデルを作成する. ロープやケーブルの変形に関する研究では, これらをはりとして考え, 柔軟モデルを作成する手法が用いられる. 例えば,erumichi のテザー (akehara et al., 2011),Ganguli のバイオリンのストリングス (Vinod Kumar and Ganguli, 2011) が挙げられる. 本研究では Euler Bernoulli のはり理論と線形ねじり変形により図 1 に示すような 5 つの要素からなる 3 次元はりモデルを作成した. ロープをシステムの i 番目の物体とし,r P を全体座標系の原点から見たロープ上の任意点 P までの位置ベクトルとすると以下の式で表現できる. Fig. 1 Schematic of the inertia included hoisting rope model r P = R i + A i (u o,ij + u f,ij ) = R i + A i u ij (1) R i は全体座標系の原点から見た i 番目の物体の座標系の原点の位置ベクトル,R i = [x i y i z i] である.A i は タイト ブライヤン角の回転行列,A i = A i (θ i ).θ i は体座標系から見られた i 番目の物体の角度であり,θ i = [φ i θ i ψ i ],φ i はロール角,θ i はピッチ角,ψ i はヨー角である.u o,ij は前変形要素 j 番上の任意点 P まで の i 番目の物体の座標系の原点から見た位置ベクトル,u f,ij は i 番目の物体の座標系から見た j 番目の要素上の 任意点 P の変形ベクトルであり.u f,ij は式 (2) で計算できる. u f,ij = S ij B ij q f,i (2) q f,i は i 番目の物体の座標系の一般化柔軟座標, B ij は j 番目の要素の境界条件で決めた柔軟座標変換行列であ る. 慣性を考慮するロープの一般化柔軟柔軟座標, 又は各要素の柔軟座標変換行列は付録 A に示す.S ij は形状 関数行列であり,Euler Bernoulli のはり理論と線形ねじり変形の形状関数行列は式 (3) のように示す.
4 1 ξ ij 0 0 6[ξ ij ξ ij 2 ]η ij 1 3ξ ij 2 + 2ξ ij 3 0 6[ξ ij ξ ij 2 ]ζ ij 0 1 3ξ ij 2 + 2ξ ij 3 0 (1 ξ ij )l ij ζ ij (1 ξ ij )l ij η ij S ij = [1 4ξ ij + 3ξ 2 ij ]l ij ζ ij 0 [ ξ ij + 2ξ 2 ij ξ 3 ij ]l ij [ 1 + 4ξ ij 3ξ 2 ij ]l ij η ij [ξ ij 2ξ 2 ij + ξ 3 ij ]l ij 0 ξ ij 0 0 6[ ξ ij + ξ ij 2 ]η ij 3ξ ij 2 2ξ ij 3 0 6[ ξ ij + ξ 2 ij ]ζ ij 0 3ξ 2 3 ij 2ξ ij 0 l ij ξ ij ζ ij l ij ξ ij η ij (3) [ [ 2ξ ij + 3ξ ij 2 ]l ij ζ ij 0 [ξ ij 2 ξ ij 3 ]l ij [2ξ ij 3ξ ij 2 ]l ij η ij [ ξ ij 2 + ξ ij 3 ]l ij 0 ] l ij は要素の長さであり, ξ ij = x ij /l ij, η ij = y ij /l ij, ζ ij = z ij /l ij である.x ij, y ij,z ij は j 番目の要素の座標系 の原点からの位置. 位置ベクトル式 (1) および運動エネルギーの式と弾性エネルギーの式より, ロープのモデ ルの j 番目の要素の慣性行列 M ij と剛性行列 K ff,ij は以下のように計算できる. I 3 3 M ij = ρ ij [ (A i u ij ) ] [I 3 3 A i u ij A i S ijb ij] dv ij (4) V ij (A i S ij B ij ) K ff,ij = (D ij S ij B ij ) E ij D ij S ij B ij dv ij (5) V ij ただし,ρ ij はロープの密度,I 3 3 は 3 3 の単位行列,V ij は要素 j 番の体積,u ij は i 番目の物体の座標系の原点 から見たの位置ベクトル,u ij, の対称行列である.D ij は微分演算子行列,E ij は材料弾性行列であり, 式 (6) と (7) に示す. 0 0 x ij D ij = 0 y ij x ij 0 [ z ij x ij] (6) E rope 0 0 E ij = [ 0 G rope 0 ] 0 0 G rope (7) 式 (7) 中の E rope はロープのヤング率,G rope はロープのせん断弾性率である. 従って, ロープの全体慣性行列 と剛性行列は式 (8) 及び (9) のように示す. 5 M i = M ij j=1 (8)
5 5 K i = K ij j=1 (9) 2 2 無質量ロープのモデル慣性の無いロープのモデルは軸方向にだけ変形する. そのため, このモデルでは座標系を省略し, ばねとして考える. 図 2 に無質量ロープによりブーム先端に掛けられた吊荷を示す. ばねの剛性率 k rope は式 (10) で計算する. Fig. 2 Schematic of the massless hoisting rope model k rope = E ropea rope L rope (10) ただし,A rope はロープの断面積,L rope はロープの長さである. 3. 移動式クレーンの動力学モデル 3 1 移動式クレーンの構造と仮定移動式クレーンのモデルをマルチボディダイナミクスの手法を用いて作成する. システムはシャシー (1 番目の物体 ), 旋回体 (2 番目の物体 ), 伸縮ブーム (3 番目の物体 ), ロープ (4 番目の物体 ) および吊荷 (5 番目の物体 ) からなる. 本研究の仮定では以下に明記する. (1) シャシーと旋回体は剛体である. (2) シャシーの回転は十分小さい, シャシーの回転行列を線変換して用いる. (3) ブームは柔軟体とし 3 つ要素の Euler Bernoulli のはりと線形ねじり変形として考え,FFRF によって作成し, 付録 B によりモデル化する. (4) 吊荷は質点である. (5) アウトリガーはばね ダンパーとして考える. ロープの慣性を考慮する多体系モデルを モデル 1 と呼び, 全部で 5 つの物体から構成する. 一方, ロープの慣性を考慮しない多体系モデルを モデル 2 と呼び, ロープを省略し, 全部で 4 つの物体から構成する. 二つのモデルを図 4 に示す.
6 (a) Chassis (b) Slewing Structure (c) Boom (d) Hoisting rope (e) Hoisted load Fig. 3 he components in a mobile crane system Fig. 4 Mobile crane system model 3 2 システムの拘束式 (1) 旋回体はシャシーの上に乗せて回転ジョイントで拘束する, 式 (11) と (12). (2) ブームと旋回体は回転ジョイントで拘束する, 式 (13) と (14). (3) ロープはブームの先端に球ジョイントで拘束する, 式 (15). (4) 吊荷はロープのもう一端に固定する, 式 (16) と (17). C 1 = R 1 + A 1 u 1,2 R 2 A 2 u 2,1 (11) C 2 = [φ 1 θ 1 ] [φ 2 θ 2 ] (12) C 3 = R 2 + A 2 u 2,3 R 3 (13) C 4 = [φ 2 ψ 2 ] [φ 3 ψ 3 ] (14)
7 C 5 = R 3 + A 3 u 3,4 R 4 (15) C 6 = R 4 + A 4 u 4,5 R 5 (16) C 7 = [e 4,26 e 4,27] (17) 式中の u i,k は i 番目の物体の座標系の原点から見た i 物体と k 番物体の拘束位置の位置ベクトルである. 従って, 全体システムの拘束式を式 (18) に示す. C(q, t) = 0 (18) 式 (18) 中の C の要素はモデルによって異なる. モデル 1 では C = [C 1 C 2 C 3 C 4 C 5 C 6 C 7 ] である. モデル 2 はロープの部分の拘束式がない代わりに, ばねの一般化力でブーム先端と吊荷を拘束するため, モデル 2 の拘束式は C = [C 1 C 2 C 3 C 4 ] である. モデル 2 のロープの一般化力は付録 C に導出する. 3 3 システムの運動方程式 ラグランジュの未定乗数法によって作成した運動方程式を式 (19) に示す. また, システムの拘束式の式 (18) を 2 階微分し, 式 (20) のように変形する. Mq + Kq + C q λ = Q (19) C q q = Q c (20) 式中の M はシステムの慣性行列,K はブームとロープ部分の柔軟性の剛性行列,C q は拘束式のヤコビ行列,λ はラグランジュ乗数,Q はシステムの一般化力を示す.Q には外力, 重力, コリオリ力ベクトルおよびアウトリ ガーの力を含む. また, モデル 2 の場合はロープのばね力も含む.Q c は拘束式を微分したヤコビ行列以外の残 り分. 最後に q はシステムの一般化座標である. モデル 1 の一般化座標は q = [R 1 θ 1 R 2 θ 2 R 3 θ 3 q f,3 R 4 θ 4 q f,4 R 5 ] となる. モデル 2 の一般化座標は q = [R 1 θ 1 R 2 θ 2 R 3 θ 3 q f,3 R 5 ] である. 本研究では拡大法及びバウムガルテの安定化法 (Baumgarte, 1972) を採用し式 (19) 及び (20) で, ラグランジュ 乗数を式 (21) のように示す. λ = (C q M 1 C q ) 1 (C q M 1 (Q Kq) γ) (21) 式 (21) 中の γ は式 (22) のように示す. γ = Q c αc β 2 C (22) 式 (22) 中の α 及び β は安定化に関するパラメータであり, 本研究では α = β = 15 を用いた. 4. 計算機シミュレーションの条件計算機シミュレーションの際にはブームの角度を 60 度で一定とし, 旋回体に回転運動を入力する. 入力回転速度を図 6 に示す. 移動式クレーンのモデルの物理パラメータは表 1 に示す通りである. 本研究の計算機シミュレーションでは式 (21) を式 (19) に入力し,MALab 2016b の ode23t を使い, 式 (19) を解いた.
8 Fig. 5 Input angular velocity of the slewing drive able 1 Specification of the mobile crane model Parameter Value Chassis mass [kg] Chassis dimension (width, length, thickness) 5, 2.5, 0.4 [m] Position of slewing drive joint from the center of chassis (x, y, z direction) 1.8, 0, 0.2 [m] Outrigger stiffness coefficient (x, y, z direction) (Maczyski and Wojciech, 2003) 1000, 1000, 7500 [kn/m] Outrigger damping coefficient (x, y, z direction) (Maczyski and Wojciech, 2003) 20, 20, 16 [kn s/m] Slewing drive mass 800 [kg] Slewing drive dimension (diameter, heights) 4, 1 [m] Position of boom joint from the center of slewing drive (x, y, z direction) 1, 0, 0.5 [m] Boom materials density 8000 [kg/m 3 ] Boom heights and width (H 31 W 31-H 32 W 32-H 33 W 33) [m 2 ] Boom thickness 0.01 [m] elescopic length (L 31-L 32-L 33) 20 [m] (8 [m] 6 [m] 6 [m]) Boom Young s modulus 210 [GPa] Boom Poisson s ratio 0.3 Rope length 5, 10, 15 [m] Rope diameter 30 [mm] Rope materials density (Feyrer, 2007) 7800 [kg/m 3 ] Rope Young s modulus (Feyrer, 2007) 196 [GPa] Rope Poisson s ratio (Feyrer, 2007) 0.3 Mass of load 100, 300, 500 [kg] 5. 計算機シミュレーションの結果および比較 5 1 振り子運動の固有振動数まず, 二つのロープのモデルに吊荷を掛けた場合の振り子運動の固有振動数を比較する. モデル 1 の場合, ロープと吊荷は荷物が付いた細い棒振り子であり, モデル 2 の場合は単振り子である.
9 (a) Influence of the length of rope (b) Enlarged section of 9-10[m] Fig. 6 Natural frequency of the pendulum motion ロープの長さ L に対する固有振動数 ω n はよく知られているようにロープが長くなると振動数が下がるが, 図 7 に示すようにモデル 1 の振動数はモデル 2 の振動数より高く, また, 吊り荷が重くなるとモデル 1 の固有振動数はモデル 2 の固有振動数に近づいている. これはロープの質量を考慮した場合, 吊り荷とロープを合わせた物体の質量中心のブーム先端からの位置が固有振動数を決定しているためである. 次に二つのロープのモデルの振り子固有振動数とモデル 1 のロープの曲げ変形の固有振動数, またブームの曲げ固有振動数を比較した結果を図 8 に示す. 図 8 に示すように振り子の振動数はロープの 1 次モードの曲げ変形より低いため, 振り子の運動はロープの変形に影響しない. しかし, ブームの縦方向 z と横方向 y の曲げ振動数は 7.89 m,9.66 m のロープ長さおいて 1 次モードのロープの曲げ変形の振動数と交わる. つまり, ロープ長さが 7 m を超える場合にはブームの柔軟性がロープの曲げ変形に影響すると考えられる. さらにロープが長くなるとロープの 1 次及び 2 次曲げ振動数が低くなり, ロープが 2 次モードの曲げ変形する可能性もあり, 振り子運動の固有振動数に近づくため, 吊荷の運動にロープの曲げ変形の影響が大きくなる可能性がある.
10 Fig. 7 Natural frequency of the pendulum motion and the deformation of structure 5 2 ロープの変形まず, 図 8 は各ロープの 3 番目の要素にある 3 番目のノードの変形ベクトルのクレーンの運動中の時間変化の大きさを示す. またはクレーンの回転入力終了後の 10 秒間の自由振動に着目し, 三つの吊荷条件と三つのロープ長さのモデル 1 のロープの形状を図 9 から 11 に示す. 図 8 により, ロープ長が長くなるとノードの変形量及び振動数が変化し, 変形値が大きくなり, 振動数が低くなる. 一方, 吊荷の質量の変化によっても変形の量及び振動数が変化するが, ロープ長の影響と比較するとその影響は小さい. また, 図 9 から 11 のロープの形状からも同様な傾向を見て取れる. これらの結果から, ロープ長が長くなると変形, 特に曲げ変形が増大するが, 吊荷の質量はロープの変形にほとんど影響しないということができる. Fig. 8 Magnitude of deformation vector of hoisting rope at node 3
11 (a) Hoisting rope in xyz plane (b) Hoisting rope in xz plane Fig. 9 Shape of the hoisting rope during residual vibration when L=5m (a) Hoisting rope in xyz plane (b) Hoisting rope in xz plane Fig. 10 Shape of the hoisting rope during residual vibration when L=10m
12 (a) Hoisting rope in xyz plane (b) Hoisting rope in xz plane Fig. 11 Shape of the hoisting rope during residual vibration when L=15m 5 3 吊荷の運動の時間変化二つのクレーンモデルを用いてシミュレーションした x と z 方向の吊荷位置の時間変化を図 12 から 17 に示す. ロープが長くなると二つのモデルの結果は異なり, モデル 1 の運動の方が少し高い振動数で低い振幅であることが分かる. しかし, 吊荷の重さが増えるとその差異が減少し, 同じような運動になる. これは固有振動数が近づくことと符合している. Fig. 12 ime history of the position in x direction of the hoisted load when L=5 m
13 Fig. 13 ime history of the position in x direction of the hoisted load when L=10 m Fig. 14 ime history of the position in x direction of the hoisted load when L=15 m Fig. 15 ime history of the position in z direction of the hoisted load when L=5 m
14 Fig. 16 ime history of the position in z direction of the hoisted load when L=10 m Fig. 17 ime history of the position in z direction of the hoisted load when L=15 m 5 4 吊荷の軌道 二つのモデルで計算した吊荷の軌道を比較し,xy 平面の軌道と xz 平面の軌道を図 18 から 23 に示す. 二つの 吊荷の軌道の差はロープが長い場合は大きく, 吊荷の質量が大きい場合は小さい. Fig. 18 rajectory of the hoisted load in xy plane when L=5 m
15 Fig. 19 rajectory of the hoisted load in xy plane when L=10 m Fig. 20 rajectory of the hoisted load in xy plane when L=15 m Fig. 21 rajectory of the hoisted load in xz plane when L=5 m Fig. 22 rajectory of the hoisted load in xz plane when L=10 m
16 Fig. 23 rajectory of the hoisted load in xz plane when L=15 m 最後に, モデル 1 およびモデル 2 で使ったシミュレーションの平均計算時間は 889 s および 70 s であった. この計算時間の違いは, モデル 1 はモデル 2 に比べ条件数が大きく, また, 一般化座標の数が約 2 倍となっているためと考えられる. これらの計算結果の比較により, モデル 2 はモデル 1 に比べ 13 分の 1 の時間で計算でき, その計算結果の差異が大きく出るロープが長く吊荷の質量が軽い場合でも差異は僅かであることがわかった. そのため, 特にロープの変形に注目する場合にはモデル 1 を用いるべきであるが, ロープの変形を詳細に必要としない場合にはモデル 2 を用いることにより, 実用的に十分な精度で高速に計算機シミュレーションが可能であることが分かった. 特に, リアルタイムシミュレーションを行うためには, モデル 2 の使用が必須である. 6. 結言本研究では移動式クレーンの動力学モデルのロープの慣性の影響を検証することを目的とした. ロープの慣性を考慮する条件と考慮しない条件の移動式クレーンの動力学モデルを柔軟多体系によって作成し, シミュレーションを行った. 得られた結果を以下に示す. (1) 慣性を考慮するロープのモデルを用いればクレーンの運転中のロープの曲げが表示でき. この変形はブームの柔軟性にも影響を受ける. また, ロープが長くなると吊荷の運動がロープの曲げ変形に影響する可能性もある. (2) 二つのモデルで計算した吊荷の軌道の差はロープが長くなると存在するが, 差異の大きさは僅かであり, 吊荷の質量が増加するにつれて軌道の差が減少する. (3) 慣性を考慮しないロープのモデルの計算時間は慣性を考慮するロープのモデルの約 13 分の 1 である. 謝 辞 本研究は株式会社タダノの助成を受けて行われたものである. 文 献 Baumgarte, J., Stabilization of constraints and integrals of motion in dynamical systems, Elsevier, Computer Methods in Applied Mechanics and Engineering, Vol.1, No.1 (1972), pp De Veubeke, B. F., he dynamics of flexible bodies, Elsevier, International Journal of Engineering Science, Vol.14, No.10 (1976), pp Feyrer, K., Wire ropes: tension, endurance, reliability (2007), pp.1-69, Springer-Verlag. Kiliçaslan, S., Balkan,. and Ider, S. K., ipping load of mobile cranes with flexible boom, Academic press, Journal of Sound and Vibration, Vol.223, No.4 (1999), pp Maczyski, A. and Wojciech, S., Dynamics of a mobile crane and optimisation of the slewing motion of its upper structure, Kluwer Academic, International Journal of Nonlinear Dynamics, Vol.32, No.3 (2003), pp
17 Park, K.-P., Cha, J.-H. and Lee, K.-Y., Dynamics factor analysis considering elastic boom effects in heavy lifting operations, Elsevier, Ocean Engineering, Vol.38, No.10 (2011), pp Posiadala, B. and Cekus, D., Discrete model of vibration of truck crane telescopic boom with consideration of the hydraulic cylinder of crane radius change in the rotary plane, Elsevier, Automation in Construction, Vol.17, No.3 (2008), pp Shabana, A. A., Flexible multibody dynamics: review of past and recent development, Kluwer Academic, Multibody System Dynamics, Vol.1, No.2 (1997), pp Sun, G. and Kleeberger, M., Dynamics responses of hydraulic mobile crane with consideration of the drive system, Elsevier, Mechanism and Machine heory, Vol.38, No.12 (2003), pp akehara, S., erumichi, Y. and Sogabe, K., Motion of a submerged tether subject to large deformations and displacements, he Japan Society of Mechanical Engineers, Journal of System Design and Dynamics, Vol.5, No.2 (2011), pp 日本機械学会編, マルチボディダイナミクス 1 基礎理論, コンピュータダイナミクス 3, コロナ社 (2006). 日本機械学会編, マルチボディダイナミクス 2 数値解析と実際, コンピュータダイナミクス 4, コロナ社 (2007). rąbka, A., Dynamics of telescopic cranes with flexible structural components, Elsevier, International Journal of Mechanical Sciences, Vol.88 (2014), pp rąbka, A., Influence of flexibilities of cranes structural components on load trajectory, Springer, Journal of Mechanical Science and echnology, Vol.30, No.1 (2016), pp Vinod Kumar, A. S. and Ganguli, R., Violin string shape functions for finite element analysis of rotating imoshenko beams, Elsevier, Journal of Finite Elements in Analysis and Design, Vol.47, No.9 (2011), pp References Baumgarte, J., Stabilization of constraints and integrals of motion in dynamical systems, Elsevier, Computer Methods in Applied Mechanics and Engineering, Vol.1, No.1 (1972), pp De Veubeke, B. F., he dynamics of flexible bodies, Elsevier, International Journal of Engineering Science, Vol.14, No.10 (1976), pp Feyrer, K., Wire ropes: tension, endurance, reliability (2007), pp.1-69, Springer-Verlag. Kiliçaslan, S., Balkan,. and Ider, S. K., ipping load of mobile cranes with flexible boom, Academic press, Journal of Sound and Vibration, Vol.223, No.4 (1999), pp Maczyski, A. and Wojciech, S., Dynamics of a mobile crane and optimisation of the slewing motion of its upper structure, Kluwer Academic, International Journal of Nonlinear Dynamics, Vol.32, No.3 (2003), pp Park, K.-P., Cha, J.-H. and Lee, K.-Y., Dynamics factor analysis considering elastic boom effects in heavy lifting operations, Elsevier, Ocean Engineering, Vol.38, No.10 (2011), pp Posiadala, B. and Cekus, D., Discrete model of vibration of truck crane telescopic boom with consideration of the hydraulic cylinder of crane radius change in the rotary plane, Elsevier, Automation in Construction, Vol.17, No.3 (2008), pp Shabana, A. A., Flexible multibody dynamics: review of past and recent development, Kluwer Academic, Multibody System Dynamics, Vol.1, No.2 (1997), pp Sun, G. and Kleeberger, M., Dynamics responses of hydraulic mobile crane with consideration of the drive system, Elsevier, Mechanism and Machine heory, Vol.38, No.12 (2003), pp akehara, S., erumichi, Y. and Sogabe, K., Motion of a submerged tether subject to large deformations and displacements, he Japan Society of Mechanical Engineers, Journal of System Design and Dynamics, Vol.5, No.2 (2011), pp he Japan Society of Mechanical Engineers ed., Multibody dynamics 1 fundamental theory, computer dynamics series 3 (2006), Corona Publishing (in Japanese). he Japan Society of Mechanical Engineers ed., Multibody dynamics 2 numerical analysis and problems, computer dynamics series 4 (2007), Corona Publishing (in Japanese).
18 rąbka, A., Dynamics of telescopic cranes with flexible structural components, Elsevier, International Journal of Mechanical Sciences, Vol.88 (2014), pp rąbka, A., Influence of flexibilities of cranes structural components on load trajectory, Springer, Journal of Mechanical Science and echnology, Vol.30, No.1 (2016), pp Vinod Kumar, A. S. and Ganguli, R., Violin string shape functions for finite element analysis of rotating imoshenko beams, Elsevier, Journal of Finite Elements in Analysis and Design, Vol.47, No.9 (2011), pp 付録 A 慣性を考慮するロープのモデルの有限要素式 本研究には慣性を考慮するロープのモデルは 4 番目の物体になり,5 つの要素で作成した. 慣性を考慮するロ ープの一般化柔軟座標は式 (A-1) に示す. q f,4 = [e 4,1 e 4,2 e 4,3 e 4,4 e 4,5 e 4,6 e 4,29 e 4,30] (A-1) 式 (A-1) 中の e i,6j 5 は要素 j 番の軸方向変形,e i,6j 4 は要素 j 番の横方向 (y) の曲げ変形,e i,6j 3 は要素 j 番の縦方向 (z) の曲げ変形,e i,6j 2 は要素 j 番のねじり変形,e i,6j 1 は要素 j 番の縦方向 (z) の曲げ変形の傾 き,e i,6j は要素 j 番の横方向 (y) の曲げ変形の傾きである ( ロープの場合,i = 4,j = 1,2,3,4,5). 式 (2) の各要素の変換行列,B ij, は式 (A-2) から (A-6) に示す. B 41 = [ I ] (A-2) B 42 = [I ] (A-3) B 43 = [ I ] (A-4) B 44 = [ I ] (A-5) B 45 = [ I ] (A-6) ただし,0 n m は n m の零行列である. 付録 B ブームのモデルの有限要素式 ブームのモデルは 3 番目の物体になり,3 つの要素で慣性を考慮するロープのモデルと同法によって作成した. 従ってブームの一般化柔軟座標は式 (B-1) に示す. q f,3 = [e 3,1 e 3,2 e 3,3 e 3,4 e 3,5 e 3,6 e 3,17 e 3,18] (B-1) ブームのモデルの各要素の変換行列は式 (B-2) から (B-4) に示す. B 31 = [ I ] (B-2) B 32 = [I ] (B-3) B 33 = [ I ] (B-4)
19 全体座標系の原点から見たブーム上の任意点 P までの位置ベクトル及び慣性行列と剛性行列は式 (1) から (9) まで, 又は式 (B-1) から (B-4) までで計算できる. 付録 C 慣性を考慮しないロープのモデルの一般化力 モデル 2 のロープのモデルはばねとして考え, 全体座標系から見たロープのベクトル,l rope, は図 C-1 に示す. Fig. C-1 Rope vector in global coordinate system 又は全体座標系の原点から見たブーム先端の位置ベクトルは式 (C-1) のように示す. r 3,tip = R 3 + A 3 (u o,33,tip + S 33,tip B 33 q f,3 ) = R 3 + A 3 u 33,tip (C-1) u o,33,tip はブームの物体座標系の原点から見た先端の重心の位置ベクトル. 従ってロープのベクトルは式 (C-2) のように示す. 又は S 33,tip = S 33 (x 33 = L 33, y 33 = 0, z 33 = 0) である. l rope = R 3 + A 3 u 33,tip R 5 (C-2) 式 (C-2) を用いてロープ全体座標系の原点見た弾力ベクトルは式 (C-3) に示す. F rope = k rope (l rope l rope,0 )l rope (C-3) ただし,l rope はロープベクトルの長さ,l rope,0 はロープの初期条件の長さ,l rope はロープの単位ベクトル. 全体座標系の原点見たロープのベクトル式 (C-2), ロープの弾力ベクトル式 (C-3) 及び仮想仕事式により, モ デル 2 のロープの一般化力は式 (C-4) のように示す. Q rope = F rope (F rope A 3 u 33,tip ) (F rope A 3 S 33,tip B 33 ) F rope ] [ (C-4)
JFE.dvi
,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho
More informationStudy on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yu
Study on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yuichiro KITAGAWA Department of Human and Mechanical
More informationThe Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh
The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Shigeru Kitamura, Member Masaaki Sakuma Genya Aoki, Member
More informationFig. 2 Signal plane divided into cell of DWT Fig. 1 Schematic diagram for the monitoring system
Study of Health Monitoring of Vehicle Structure by Using Feature Extraction based on Discrete Wavelet Transform Akihisa TABATA *4, Yoshio AOKI, Kazutaka ANDO and Masataka KATO Department of Precision Machinery
More informationThe Evaluation of LBB Behavior and Crack Opening Displacement on Statically Indeterminate Piping System Subjected to Monotonic Load The plastic collap
The Evaluation of LBB Behavior and Crack Opening Displacement on Statically Indeterminate Piping System Subjected to Monotonic Load The plastic collapse and LBB behavior of statically indeterminate piping
More informationInstability of Aerostatic Journal Bearings with Porous Floating Bush at High Speeds Masaaki MIYATAKE *4, Shigeka YOSHIMOTO, Tomoaki CHIBA and Akira CH
Instability of Aerostatic Journal Bearings with Porous Floating Bush at High Speeds Masaaki MIYATAKE *4, Shigeka YOSHIMOTO, Tomoaki CHIBA and Akira CHIBA Department of Mechanical Engineering, Tokyo University
More information第62巻 第1号 平成24年4月/石こうを用いた木材ペレット
Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting
More informationa) Extraction of Similarities and Differences in Human Behavior Using Singular Value Decomposition Kenichi MISHIMA, Sayaka KANATA, Hiroaki NAKANISHI a
a) Extraction of Similarities and Differences in Human Behavior Using Singular Value Decomposition Kenichi MISHIMA, Sayaka KANATA, Hiroaki NAKANISHI a), Tetsuo SAWARAGI, and Yukio HORIGUCHI 1. Johansson
More informationStudy on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Enginee
Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Engineering, Kanazawa University, Kakuma-machi, Kanazawa-shi,
More informationOCW-iダランベールの原理
講義名連続体力学配布資料 OCW- 第 2 回ダランベールの原理 無機材料工学科准教授安田公一 1 はじめに今回の講義では, まず, 前半でダランベールの原理について説明する これを用いると, 動力学の問題を静力学の問題として解くことができ, さらに, 前回の仮想仕事の原理を適用すると動力学問題も簡単に解くことができるようになる また, 後半では, ダランベールの原理の応用として ラグランジュ方程式の導出を示す
More information日本感性工学会論文誌
Vol.11 No.3 pp.397-406 2012 The Effect of Pleat Number and Skirt Length on the Visual Image of a Pleated Skirt by the Age Group Junghwa LEE*, Naomi MARUTA** and Taeko HIROKAWA*** * Bunka Gakuen University,
More information1..FEM FEM 3. 4.
008 stress behavior at the joint of stringer to cross beam of the steel railway bridge 1115117 1..FEM FEM 3. 4. ABSTRACT 1. BackgroundPurpose The occurrence of fatigue crack is reported in the joint of
More informationA Study on Throw Simulation for Baseball Pitching Machine with Rollers and Its Optimization Shinobu SAKAI*5, Yuichiro KITAGAWA, Ryo KANAI and Juhachi
A Study on Throw Simulation for Baseball Pitching Machine with Rollers and Its Optimization Shinobu SAKAI*5, Yuichiro KITAGAWA, Ryo KANAI and Juhachi ODA Department of Human and Mechanical Systems Engineering,
More informationVol. 62A ( ) Fundamental consideration on the redundancy of a steel Langer bridge based on nonlinear finite element analysis Isao Saiki, Kota Ka
Vol. 6A (6 3 ) Fundamental consideration on the redundanc of a steel Langer bridge based on nonlinear finite element analsis Isao Saiki, Kota Kawamura, Kaname Iwatsubo, Tetsuo Iwakuma -57 6-6-6 63- --
More informationMicrosoft Word - thesis.doc
剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル
More information第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える. 5.1 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f = l 2pL である. ただし, L [ 単位 m] は棒の長さ, [ 2 N / m ] 3 r[ 単位 Kg / m ] E r
第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える 5 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f l pl である ただし, L [ 単位 m] は棒の長さ, [ N / m ] [ 単位 Kg / m ] E は (5) E 単位は棒の材料の縦弾性係数 ( ヤング率 ) は棒の材料の単位体積当りの質量である l は境界条件と振動モードによって決まる無
More information2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように
3 章 Web に Link 解説 連続式 微分表示 の誘導.64 *4. 連続式連続式は ある領域の内部にある流体の質量の収支が その表面からの流入出の合計と等しくなることを定式化したものであり 流体における質量保存則を示したものである 2. 連続式 微分表示 の誘導図のような微小要素 コントロールボリューム の領域内の流体の増減と外部からの流体の流入出を考えることで定式化できる 微小要素 流入
More informationFig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).
Fig. 1 The scheme of glottal area as a function of time Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels). Fig, 4 Parametric representation
More information第6章 実験モード解析
第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法
More informationStudies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth
Studies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth and Foot Breadth Akiko Yamamoto Fukuoka Women's University,
More information<95DB8C9288E397C389C88A E696E6462>
2011 Vol.60 No.2 p.138 147 Performance of the Japanese long-term care benefit: An International comparison based on OECD health data Mie MORIKAWA[1] Takako TSUTSUI[2] [1]National Institute of Public Health,
More informationVisual Evaluation of Polka-dot Patterns Yoojin LEE and Nobuko NARUSE * Granduate School of Bunka Women's University, and * Faculty of Fashion Science,
Visual Evaluation of Polka-dot Patterns Yoojin LEE and Nobuko NARUSE * Granduate School of Bunka Women's University, and * Faculty of Fashion Science, Bunka Women's University, Shibuya-ku, Tokyo 151-8523
More informationOptical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd:
*1 *2 *3 PIV Measurement of Field of the Wind Turbine with a med Diffuser Kazuhiko TOSHIMITSU *4, Koutarou NISHIKAWA and Yuji OHYA *4 Department of Mechanical Engineering, Matsue National Collage of Technology,
More information日立金属技報 Vol.34
Influence of Misorientation Angle between Adjacent Grains on Magnetization Reversal in Nd-Fe-B Sintered Magnet Tomohito Maki Rintaro Ishii Mitsutoshi Natsumeda Takeshi Nishiuchi Ryo Uchikoshi Masaaki Takezawa
More information* Meso- -scale Features of the Tokai Heavy Rainfall in September 2000 Shin-ichi SUZUKI Disaster Prevention Research Group, National R
38 2002 7 2000 9 * Meso- -scale Features of the Tokai Heavy Rainfall in September 2000 Shin-ichi SUZUKI Disaster Prevention Research Group, National Research Institute for Earth Science and Disaster Prevention,
More informationr 0 r 45 r 90 F 0 n
Evaluation of Fatigue and Noise-and-vibration Properties of Automobile Partial Models Abstract Application of high strength steel sheets to automotive bodies requires evaluation technologies of fatigue
More information特-11.indd
CCH3000LJ Development of Luffing Jib Crawler Crane CCH3000LJ IHI IHI IHI IHI IHI IHI IHI In recent years, the need for the construction of social infrastructure, particularly in developing nations such
More informationSynthesis and Development of Electric Active Stabilizer Suspension System Shuuichi BUMA*6, Yasuhiro OOKUMA, Akiya TANEDA, Katsumi SUZUKI, Jae-Sung CHO
Synthesis and Development of Electric Active Stabilizer Suspension System Shuuichi BUMA*6, Yasuhiro OOKUMA, Akiya TANEDA, Katsumi SUZUKI, Jae-Sung CHO and Masaru KOBAYASHI Chassis Engineering Management
More information149 (Newell [5]) Newell [5], [1], [1], [11] Li,Ryu, and Song [2], [11] Li,Ryu, and Song [2], [1] 1) 2) ( ) ( ) 3) T : 2 a : 3 a 1 :
Transactions of the Operations Research Society of Japan Vol. 58, 215, pp. 148 165 c ( 215 1 2 ; 215 9 3 ) 1) 2) :,,,,, 1. [9] 3 12 Darroch,Newell, and Morris [1] Mcneil [3] Miller [4] Newell [5, 6], [1]
More information.I.v e pmd
Structural Design for Curved Panels by Laminated Composite Materials (Identification of Lamination Parameters Using Modal Testing Method ) Tetsuya NARISAWA, Shohei IWATA Abstract - Using a modal testing
More informationA Navigation Algorithm for Avoidance of Moving and Stationary Obstacles for Mobile Robot Masaaki TOMITA*3 and Motoji YAMAMOTO Department of Production
A Navigation Algorithm for Avoidance of Moving and Stationary Obstacles for Mobile Robot Masaaki TOMITA*3 and Motoji YAMAMOTO Department of Production System Engineering, Kyushu Polytecnic College, 1665-1
More informationA Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical
A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical calculation method of the gradient as a differential
More informationJournal of Textile Engineering, Vol.53, No.5, pp
ORIGINAL PAPER Journal of Textile Engineering (2007), Vol.53, No.5, 203-210 2007 The Textile Machinery Society of Japan Analyzing the Path and the Tension of a Yarn under a False-twist Process Using a
More informationsakigake1.dvi
(Zin ARAI) arai@cris.hokudai.ac.jp http://www.cris.hokudai.ac.jp/arai/ 1 dynamical systems ( mechanics ) dynamical systems 3 G X Ψ:G X X, (g, x) Ψ(g, x) =:Ψ g (x) Ψ id (x) =x, Ψ gh (x) =Ψ h (Ψ g (x)) (
More informationJAMSTEC Rep. Res. Dev., Volume 12, March 2011, 27 _ 35 1,2* Pb 210 Pb 214 Pb MCA 210 Pb MCA MCA 210 Pb 214 Pb * 2
JAMSTEC Rep. Res. Dev., Volume 12, March 2011, 27 _ 35 1,2* 1 1 1 1 210 Pb 210 Pb 214 Pb MCA 210 Pb MCA MCA 210 Pb 214 Pb 2010 10 4 2010 12 10 1 2 * 237-0061 2-15 046-867-9794 ogurik@jamstec.go.jp 27 210
More information原稿.indd
OTEC 18(2013),5968 小型の波浪発電を想定した浮体運動の最大化を目的とした浮体形状に関する研究 59 *1 *1 *2 Studies on the floating body shape to maximize the kinetic energy that are intended to be small wave power generator Shunya NISHIZAWA
More informationTitle 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue
Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Date 2017 URL http://hdl.handle.net/2433/229150 Right
More information国土技術政策総合研究所 研究資料
ISSN TECHNICAL NOTE of National Institute for Land and Infrastructure Management No255 September 2005 A Result of Investigation into the Actual Conditions of Japanese Ports Pavement and Cargo Handling
More informationStress Singularity Analysis at an Interfacial Corner Between Anisotropic Bimaterials Under Thermal Stress Yoshiaki NOMURA, Toru IKEDA*4 and Noriyuki M
Stress Singularity Analysis at an Interfacial Corner Between Anisotropic Bimaterials Under Thermal Stress Yoshiaki NOMURA, Toru IKEDA*4 and Noriyuki MIYAZAKI Department of Mechanical Engineering and Science,
More informationSafety Performance of Steel Deck Plate (Flat Decks) Used for Concrete Slab Moulding CONTENTS 1. Introduction ---------------------------------------------------------------- (2) 2. Flat Decks ------------------------------------------------------------------
More informationT05_Nd-Fe-B磁石.indd
Influence of Intergranular Grain Boundary Phases on Coercivity in Nd-Fe-B-based Magnets Takeshi Nishiuchi Teruo Kohashi Isao Kitagawa Akira Sugawara Hiroyuki Yamamoto To determine how to increase the coercivity
More information領域シンポ発表
1 次元の減衰運動の中の強制振動 ) ( f d d d d d e f e ce ) ( si ) ( 1 ) ( cos ω =ω -γ とおくと 一般解は 外力 f()=f siω の場合 f d d d d si f ce f ce si ) cos( cos si ) cos( この一般解は 1 φ は外力と変位との間の位相差で a 時間が経つと 第 1 項は無視できる この場合の振幅を
More information電子部品はんだ接合部の熱疲労寿命解析
43 Evaluation for Thermal Fatigue Life of Solder Joints in Electronic Components Haruhiko Yamada, Kazuyoshi Ogawa 2 63Sn- 37Pb 95Pb-5Sn Si Cu Si 63Sn-37Pb Since automotive electronic components are used
More information建築構造力学 I ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.
建築構造力学 I ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/050043 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 1 38 2 15 2 1 2 2 1 2 2 1977 2007 2015 10 ii F P = mα g =
More informationRubin Rubin
Yamagata Journal of Health Science, Vol. 7, 2004 2 2 Sumi MISAWA, Chizu KATAGIRI, Ryoko KOMATSU, Yohko FUJISAWA The first report focused on maternal role behaviors of women in childbirth, especially the
More informationIPSJ SIG Technical Report Vol.2012-CG-148 No /8/29 3DCG 1,a) On rigid body animation taking into account the 3D computer graphics came
3DCG 1,a) 2 2 2 2 3 On rigid body animation taking into account the 3D computer graphics camera viewpoint Abstract: In using computer graphics for making games or motion pictures, physics simulation is
More information1 Fig. 1 Extraction of motion,.,,, 4,,, 3., 1, 2. 2.,. CHLAC,. 2.1,. (256 ).,., CHLAC. CHLAC, HLAC. 2.3 (HLAC ) r,.,. HLAC. N. 2 HLAC Fig. 2
CHLAC 1 2 3 3,. (CHLAC), 1).,.,, CHLAC,.,. Suspicious Behavior Detection based on CHLAC Method Hideaki Imanishi, 1 Toyohiro Hayashi, 2 Shuichi Enokida 3 and Toshiaki Ejima 3 We have proposed a method for
More information食品工学.indb
, Vol. 11, No. 1, pp. 51-58, Mar. 2010 IH The Visualization and Quantification of the Flow in the Pan During Induction Heating and Gas Heating Haruna KAWAKAMI, Zensyu TOU, Mika FUKUOKA, and Noboru SAKAI
More informationFig. 1 KAMOME50-2 Table 1 Principal dimensions Fig.2 Configuration of the hydrofoils (Endurance and sprint foil) Fig. 3 Schematic view of the vortex l
A study on the simulation of the motion of personal full-submerged hydrofoil craft by Yutaka Terao, Member Summary A new energy utilization project developed by Tokai University was started in 1991. It
More informationAutumn II III Zon and Muysken 2005 Zon and Muysken 2005 IV II 障害者への所得移転の経済効果 分析に用いるデータ
212 Vol. 44 No. 2 I はじめに 2008 1 2 Autumn 08 213 II III Zon and Muysken 2005 Zon and Muysken 2005 IV II 障害者への所得移転の経済効果 17 18 1 分析に用いるデータ 1 2005 10 12 200 2 2006 9 12 1 1 2 129 35 113 3 1 2 6 1 2 3 4 4 1
More informationパソコンシミュレータの現状
第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に
More informationCore Ethics Vol.
Core Ethics Vol. - Core Ethics Vol. Core Ethics Vol. ,,,,,, Core Ethics Vol. - Core Ethics Vol. p. p. p. p. p. pp. - p. " pp. - 煞 ,,..,,....,.......,.,...,..,.,.......,..,,,,,,,,,,,,........,.,.....,...,.,,,,,,,,..,.....
More information06_学術_関節単純X線画像における_1c_梅木様.indd
Arts and Sciences X The formulation of femoral heard measurement corrected enlargement ratio using hip joints X-ray Imaging 1 2 1 1 1 2 Key words: Bipolar Hip Arthroplasty (BHA) Preoperative planning Enlargement
More information202
201 Presenteeism 202 203 204 Table 1. Name Elements of Work Productivity Targeted Populations Measurement items of Presenteeism (Number of Items) Reliability Validity α α 205 α ä 206 Table 2. Factors of
More information最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.
最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/052093 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 10 3 2000 2007 26 8 2 SI SI 20 1996 2000 SI 15 3 ii 1 56 6
More information& Vol.5 No (Oct. 2015) TV 1,2,a) , Augmented TV TV AR Augmented Reality 3DCG TV Estimation of TV Screen Position and Ro
TV 1,2,a) 1 2 2015 1 26, 2015 5 21 Augmented TV TV AR Augmented Reality 3DCG TV Estimation of TV Screen Position and Rotation Using Mobile Device Hiroyuki Kawakita 1,2,a) Toshio Nakagawa 1 Makoto Sato
More information44 2012 2013 3 35 48 法人化後の国立大学の収入変動 37 法人化後の国立大学の収入変動 2009 2005 2010 2012 2012 2008 2009a 2010 16 18 17 20 2 4 2012 38 44 2012 17 22 (1) (2) 2012 5 GP COE 30 WPI 1 2012 17 22 16 17 22 17 17 19 2012 2012
More information5D1 SY0004/14/ SICE 1, 2 Dynamically Consistent Motion Design of Humanoid Robots even at the Limit of Kinematics Kenya TANAKA 1 and Tomo
5D1 SY4/14/-485 214 SICE 1, 2 Dynamically Consistent Motion Design of Humanoid Robots even at the Limit of Kinematics Kenya TANAKA 1 and Tomomichi SUGIHARA 2 1 School of Engineering, Osaka University 2-1
More information<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1>
人工環境設計解析工学構造力学と有限要素法 ( 第 回 ) 東京大学新領域創成科学研究科 鈴木克幸 固体力学の基礎方程式 変位 - ひずみの関係 適合条件式 ひずみ - 応力の関係 構成方程式 応力 - 外力の関係 平衡方程式 境界条件 変位規定境界 反力規定境界 境界条件 荷重応力ひずみ変形 場の方程式 Γ t Γ t 平衡方程式構成方程式適合条件式 構造力学の基礎式 ひずみ 一軸 荷重応力ひずみ変形
More informationMotivation and Purpose There is no definition about whether seatbelt anchorage should be fixed or not. We tested the same test conditions except for t
Review of Seatbelt Anchorage and Dimensions of Test Bench Seat Cushion JASIC Motivation and Purpose There is no definition about whether seatbelt anchorage should be fixed or not. We tested the same test
More informationIPSJ SIG Technical Report NetMAS NetMAS NetMAS One-dimensional Pedestrian Model for Fast Evacuation Simulator Shunsuke Soeda, 1 Tomohisa Yam
1 1 1 1 1 NetMAS NetMAS NetMAS One-dimensional Model for Fast Evacuation Simulator Shunsuke Soeda, 1 Tomohisa Yamashita, 1 Masaki Onishi, 1 Ikushi Yoda 1 and Itsuki Noda 1 We propose the one-dimentional
More information6) , 3) L60m h=4m 4m φ19 SS400 σ y = kn/mm 2 E = 205.8kN/mm 2 Table1 4) 7 Fig.1 5 7) S S 2 5 (Fig.2 ) ( No.1, No.2, No.3, No.4)
Damages and Earthquake Resistant Performance of Steel Frame Structures with Self Strain Stress (Faculty of Architecture and Structural Engineering) Yutaka NIHO, Masaru TERAOKA (Professor Emeritus of KNCT)
More informationThe Japanese Journal of Health Psychology, 29(S): (2017)
Journal of Health Psychology Research 2017, Vol. 29, Special issue, 139 149Journal of Health Psychology Research 2016, J-STAGE Vol. Advance 29, Special publication issue, 139 149 date : 5 December, 2016
More informationON STRENGTH AND DEFORMATION OF REINFORCED CONCRETE SHEAR WALLS By Shigeru Mochizuki Concrete Journal, Vol. 18, No. 4, April 1980, pp. 1 `13 Synopsis A
ON STRENGTH AND DEFORMATION OF REINFORCED CONCRETE SHEAR WALLS By Shigeru Mochizuki Concrete Journal, Vol. 18, No. 4, April 1980, pp. 1 `13 Synopsis After Tokachioki Earthquake of 1968, the importance
More information: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs
15 C11-4 Numerical analysis of flame propagation in a combustor of an aircraft gas turbine, 4-6-1 E-mail: tominaga@icebeer.iis.u-tokyo.ac.jp, 2-11-16 E-mail: ntani@iis.u-tokyo.ac.jp, 4-6-1 E-mail: itoh@icebeer.iis.u-tokyo.ac.jp,
More informationTable 1. Reluctance equalization design. Fig. 2. Voltage vector of LSynRM. Fig. 4. Analytical model. Table 2. Specifications of analytical models. Fig
Mover Design and Performance Analysis of Linear Synchronous Reluctance Motor with Multi-flux Barrier Masayuki Sanada, Member, Mitsutoshi Asano, Student Member, Shigeo Morimoto, Member, Yoji Takeda, Member
More informationVol.7 No.2 ( ) in mm m/s 40 m/s 20 m/s m/s 20m/s 1999 US ,
Vol.7 No.2 (136-151) KAWAZOE Yoshihiko 13 1991; 1993; 1994; 2001; 2002e; 2003-136 - Vol.7 No.2 (136-151) 27 27 137 in 2 32 812mm 1997 2000 29 2030 m/s 40 m/s 20 m/s 0.48 14 25 m/s 20m/s 1999 US 2000 2001
More informationteionkogaku43_527
特集 : 振動流によるエネルギー変換 熱輸送現象と応用技術 * Oscillatory Flow in a Thermoacoustic Sound-wave Generator - Flow around the Resonance Tube Outlet - Masayasu HATAZAWA * Synopsis: This research describes the oscillatory
More informationJ. Jpn. Inst. Light Met. 65(6): 224-228 (2015)
65 62015 224 228 ** Journal of The Japan Institute of Light Metals, Vol. 65, No. 6 (2015), 224 228 2015 The Japan Institute of Light Metals Investigation of heat flow behavior on die-casting core pin with
More informationSICE東北支部研究集会資料(2012年)
77 (..3) 77- A study on disturbance compensation control of a wheeled inverted pendulum robot during arm manipulation using Extended State Observer Luis Canete Takuma Sato, Kenta Nagano,Luis Canete,Takayuki
More informationIPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of
GPU 1,a) 1,b) 1,c) 1,d) GPU 1 GPU Structure Of Array Array Of Structure 1. MPS(Moving Particle Semi-Implicit) [1] SPH(Smoothed Particle Hydrodynamics) [] DEM(Distinct Element Method)[] [] 1 Tokyo Institute
More informationStudy on Deformation of Fixation Devices of RED System in Maxillary Distraction Osteogenesis KoUJI OHTA, HIDETOSHI TOHMORI, MASARU SUGIYAMA, TAKESHI N
Study on Deformation of Fixation Devices of RED System in Maxillary Distraction Osteogenesis KoUJI OHTA, HIDETOSHI TOHMORI, MASARU SUGIYAMA, TAKESHI NOBUMORI, SHINGO INOUE, KOTARO TANIMOTO*, KAZUO TANNE
More informationuntitled
2010 58 1 39 59 c 2010 20 2009 11 30 2010 6 24 6 25 1 1953 12 2008 III 1. 5, 1961, 1970, 1975, 1982, 1992 12 2008 2008 226 0015 32 40 58 1 2010 III 2., 2009 3 #3.xx #3.1 #3.2 1 1953 2 1958 12 2008 1 2
More informationMicrosoft PowerPoint - シミュレーション工学-2010-第1回.ppt
シミュレーション工学 ( 後半 ) 東京大学人工物工学研究センター 鈴木克幸 CA( Compter Aded geerg ) r. Jaso Lemo (SC, 98) 設計者が解析ツールを使いこなすことにより 設計の評価 設計の質の向上を図る geerg の本質の 計算機による支援 (CA CAM などより広い名前 ) 様々な汎用ソフトの登場 工業製品の設計に不可欠のツール 構造解析 流体解析
More informationIPSJ SIG Technical Report Vol.2016-CE-137 No /12/ e β /α α β β / α A judgment method of difficulty of task for a learner using simple
1 2 3 4 5 e β /α α β β / α A judgment method of difficulty of task for a learner using simple electroencephalograph Katsuyuki Umezawa 1 Takashi Ishida 2 Tomohiko Saito 3 Makoto Nakazawa 4 Shigeichi Hirasawa
More informationKeywords: corotational method, Rigid-Bodies-Spring model, accuracy, geometrical nonlinearity
Keywords: corotational method, Rigid-Bodies-Spring model, accuracy, geometrical nonlinearity vo2=voi+(sina1+sina2)l/2+{f1(sin2a1+sin2a2) +Fz1(sina1cosai+sina2cosa2)}l/(2EA) woe=w01-l+(cosai+cosa2)l/2+{fy1(sinaicosai
More informationFig. 3 Coordinate system and notation Fig. 1 The hydrodynamic force and wave measured system Fig. 2 Apparatus of model testing
The Hydrodynamic Force Acting on the Ship in a Following Sea (1 St Report) Summary by Yutaka Terao, Member Broaching phenomena are most likely to occur in a following sea to relative small and fast craft
More informationテレビ番組による相互交流
Abstract: Import and Export of TV programs in Japan The TV programs for which contracts were concluded for export during fiscal 2001 included: Number of programs (number of titles) = 1,675 programs (series
More informationTHE JAPANESE JOURNAL OF PERSONALITY 2007, Vol. 15 No. 2, 217–227
2007 15 2 217 227 2007 1 2 445 238 207 1 7 7 (2005) (2005) (2006a) 2006a (2004) (assumed-competence) 218 15 2 1 2005 2004 (2004) (2004) (2004) 16PF (2006) 4 Hayamizu, 2002 2006a (2006a) 1 Hayamizu, Kino,
More informationDevelopment of Induction and Exhaust Systems for Third-Era Honda Formula One Engines Induction and exhaust systems determine the amount of air intake
Development of Induction and Exhaust Systems for Third-Era Honda Formula One Engines Induction and exhaust systems determine the amount of air intake supplied to the engine, and as such are critical elements
More information2 点吊り物理振子の振動解析: 弾性紐の効果 木ノ内 智貴*1 舟田 敏雄*1 桜井 賢人*1 大庭 勝久*1 青木 悠祐*1 宮内 太積*2 望月 孔二*3 An Analysis of Mode Coupling in Three Modes of Bifilar Suspension Phys
点吊り物理振子の振動解析: 弾性紐の効果 木ノ内 智貴* 舟田 敏雄* 桜井 賢人* 大庭 勝久* 青木 悠祐* 宮内 太積* 望月 孔二*3 An Analysis of Mode Coupling in Three Modes of Bifilar Suspension Physical-Pendulum: Effects of Elastic Strings Toshiki KINOUCHI*
More information国民年金保険料における未納 免除 猶予 追納の分析 Analysis of People's Decision-Making for the Absence of Contribution Payments, the Exemption, the Contribution Postponement
Powered by TCPDF (www.tcpdf.org) Title Sub Title Author Publisher 国民年金保険料における未納 免除 猶予 追納の分析 Analysis of people's decision-making for the absence of contribution payments, the exemption, the contribution
More information日本感性工学会論文誌
pp.389-402 2017 doi: 10.5057/jjske.TJSKE-D-17-00019 SKEL Fundamental Analysis on Designer s Inference Process Framework and Its Visualization Proposal of Inference Mapping Method to Assist Meta-cognition
More informationMicrosoft PowerPoint - ロボットの運動学forUpload'C5Q [互換モード]
ロボットの運動学 順運動学とは 座標系の回転と並進 同次座標変換行列 Denavit-Hartenberg の表記法 多関節ロボットの順運動学 レポート課題 & 中間試験について 逆運動学とは ヤコビアン行列 運動方程式 ( 微分方程式 ) ロボットの運動学 動力学 Equation of motion f ( ( t), ( t), ( t)) τ( t) 姿勢 ( 関節角の組合せ ) Posture
More informationDevelopement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Sum
Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Summary Previously, the authors developed the plastic
More information1 [1, 2, 3, 4, 5, 8, 9, 10, 12, 15] The Boston Public Schools system, BPS (Deferred Acceptance system, DA) (Top Trading Cycles system, TTC) cf. [13] [
Vol.2, No.x, April 2015, pp.xx-xx ISSN xxxx-xxxx 2015 4 30 2015 5 25 253-8550 1100 Tel 0467-53-2111( ) Fax 0467-54-3734 http://www.bunkyo.ac.jp/faculty/business/ 1 [1, 2, 3, 4, 5, 8, 9, 10, 12, 15] The
More informationVol.8 No (July 2015) 2/ [3] stratification / *1 2 J-REIT *2 *1 *2 J-REIT % J-REIT J-REIT 6 J-REIT J-REIT 10 J-REIT *3 J-
Vol.8 No.2 1 9 (July 2015) 1,a) 2 3 2012 1 5 2012 3 24, 2013 12 12 2 1 2 A Factor Model for Measuring Market Risk in Real Estate Investment Hiroshi Ishijima 1,a) Akira Maeda 2 Tomohiko Taniyama 3 Received:
More information(Microsoft PowerPoint - \221\34613\211\361)
計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.
More informationIPSJ SIG Technical Report Vol.2014-IOT-27 No.14 Vol.2014-SPT-11 No /10/10 1,a) 2 zabbix Consideration of a system to support understanding of f
1,a) 2 zabbix Consideration of a system to support understanding of fault occurrences based on the similarity of the time series Miyaza Nao 1,a) Masuda Hideo 2 Abstract: With the development of network
More informationTable 1 Type of polymeric coating materials Fig. 2 Results of suppressive effects of polymeric coating materials on the progress of neutralization of concrete. Table 2 Evaluation of the suppressive effects
More information空力騒音シミュレータの開発
41 COSMOS-V, an Aerodynamic Noise Simulator Nariaki Horinouchi COSMOS-V COSMOS-V COSMOS-V 3 The present and future computational problems of the aerodynamic noise analysis using COSMOS-V, our in-house
More informationI-2 (100 ) (1) y(x) y dy dx y d2 y dx 2 (a) y + 2y 3y = 9e 2x (b) x 2 y 6y = 5x 4 (2) Bernoulli B n (n = 0, 1, 2,...) x e x 1 = n=0 B 0 B 1 B 2 (3) co
16 I ( ) (1) I-1 I-2 I-3 (2) I-1 ( ) (100 ) 2l x x = 0 y t y(x, t) y(±l, t) = 0 m T g y(x, t) l y(x, t) c = 2 y(x, t) c 2 2 y(x, t) = g (A) t 2 x 2 T/m (1) y 0 (x) y 0 (x) = g c 2 (l2 x 2 ) (B) (2) (1)
More information1
半剛節が部材上の任意点にある部材剛性方程式 米子高専 川端康洋 稲田祐二. ピン半剛節を有する部材の解析の歴史 ()940 二見秀雄材の途中にピン接合点を有するラーメン材の算式とその応用建築学会論文集 つのピン節を含む部材の撓角法基本式と荷重項ピン節を含む部材の撓角法基本式と荷重項が求められている 以降 固定モーメント法や異形ラーメンの解法への応用が研究された 戦後には 関連する論文は見当たらない
More information1 Web [2] Web [3] [4] [5], [6] [7] [8] S.W. [9] 3. MeetingShelf Web MeetingShelf MeetingShelf (1) (2) (3) (4) (5) Web MeetingShelf
1,a) 2,b) 4,c) 3,d) 4,e) Web A Review Supporting System for Whiteboard Logging Movies Based on Notes Timeline Taniguchi Yoshihide 1,a) Horiguchi Satoshi 2,b) Inoue Akifumi 4,c) Igaki Hiroshi 3,d) Hoshi
More informationIR0036_62-3.indb
62 3 2016 253 272 1921 25 : 27 8 19 : 28 6 3 1921 25 1921 25 1952 27 1954 291960 35 1921 25 Ⅰ 0 5 1 5 10 14 21 25 34 36 59 61 6 8 9 11 12 16 1921 25 4 8 1 5 254 62 3 2016 1 1938.8 1926 30 1938.6.23 1939.9
More information013858,繊維学会誌ファイバー1月/報文-02-古金谷
Development of Non-Contact Measuring Method for Final Twist Number of Double Ply Staple Yarn Keizo Koganeya 1, Youichi Yukishita 1, Hirotaka Fujisaki 1, Yasunori Jintoku 2, Hironori Okuno 2, and Motoharu
More information75 unit: mm Fig. Structure of model three-phase stacked transformer cores (a) Alternate-lap joint (b) Step-lap joint 3 4)
3 * 35 (3), 7 Analysis of Local Magnetic Properties and Acoustic Noise in Three-Phase Stacked Transformer Core Model Masayoshi Ishida Kenichi Sadahiro Seiji Okabe 3.7 T 5 Hz..4 3 Synopsis: Methods of local
More informationNumerical Simulation for Abrupt Contraction Flow of Fiber Suspensions in Polymeric Fluid Kazunori Yasuda, Taro Nishimura* and Kiyoji Nakamura Departme
Numerical Simulation for Abrupt Contraction Flow of Fiber Suspensions in Polymeric Fluid Kazunori Yasuda, Taro Nishimura* and Kiyoji Nakamura Department of Mechanical Engineering, Osaka University, Suita,
More information新製品開発プロジェクトの評価手法
CIRJE-J-60 2001 8 A note on new product project selection model: Empirical analysis in chemical industry Kenichi KuwashimaUniversity of Tokyo Junichi TomitaUniversity of Tokyo August, 2001 Abstract By
More information今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未
力学 III GA 工業力学演習 X5 解析力学 5X 5 週目 立命館大学機械システム系 8 年度後期 今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未定乗数法
More information