本資料では Flat imlator Vr.5.. の下記改良成果についてご報告します iss iscolastic modl を利用した Film castig simlatio 機能の実装 iss iscolastic modl を利用した Matrial charactrizatio 機能の実

Size: px
Start display at page:

Download "本資料では Flat imlator Vr.5.. の下記改良成果についてご報告します iss iscolastic modl を利用した Film castig simlatio 機能の実装 iss iscolastic modl を利用した Matrial charactrizatio 機能の実"

Transcription

1 Flat imlatorvr.5.. 改良成果資料 発表用ダイジェスト版 5//5 株式会社 HAL Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

2 本資料では Flat imlator Vr.5.. の下記改良成果についてご報告します iss iscolastic modl を利用した Film castig simlatio 機能の実装 iss iscolastic modl を利用した Matrial charactrizatio 機能の実装 当機能は Matrial fit 改良成果資料をご参照下さい Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

3 iss iscolastic modl を利用した Film castig simlatio 機能の実装 解析対象モデル : 図 フィルムキャストプロセス Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

4 成形現象の支配方程式 質量保存方程式 : t 運動量保存方程式 : σ エネルギー保存方程式 : p. htot : T T T ρ C 3 p htot T Tair Tair : t Nomclatr : : : σ p : T : ρ : C p : : Vlocit compot Film thicss Etra strss compot trss compot Prssr Tmpratr Dsit Hat capacit Hat trasfr cofficit Ambit air tmpratr Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 3

5 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 4 構成方程式 : Nwtoia modl. D D 4 5 t t t t iss modl 未知量 未知量 モード数

6 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 5 支配方程式の数値解析法 6 Psdo bod forc 疑似体積力 の導入粘弾性応力粘性応力 粘弾性応力 ー粘性応力疑似体積力と呼び 上付き添え字 を付けて表現 7

7 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 6 有限要素離散化方程式 4 β β β β β β 4 β β β β β β. d d d d d d β β β β β β β β ξ 3 4 :-- :- 3: 4: ξ ξ ξ ξ 8 9 図 補間関数

8 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 7 粘弾性応力 体積力は要素重心に定義する 粘弾性応力未知量は コントロールボリューム法を利用して評価する 例えば に対するコントロールボリューム離散化方程式は以下に示す通り t d t d t コントロールボリューム法を利用して離散化対流項の差分近似として Upstram diffrc schm を採用当方程式の解法として OR 法に基づく陰解法を採用

9 テスト解析 mm/s 66 mm T ダイ流出口肉厚 3mm 均一 mm 4 mm/s Draw Ratio : 4 図 3 解析モデル Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 8

10 表 ケーススタディ一覧 Cas 緩和弾性率 Pa 緩和時間 s iss パラメータ Cas は粘性解析.Pa s を採用 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 9

11 解析結果 緩和時間 弾性 の増加に伴って フィルムネックイン量は減少します また iss モデルの パラメータを増加させると過渡一軸伸張粘度のひずみ硬化性が低下します Cas45 の結果を比較すると ひずみ硬化性の低下に伴って ネックイン量が増加する傾向が示されています Film width mm Film c i Cas Cas Cas Cas3 Cas4 Cas Dow stram distac mm 図 4 フィルム幅計算結果の比較 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

12 緩和時間 弾性 やひずみ硬化性が増すとフィルムの均一性が向上する傾向が示されます.3.5 Film thicss mm mm : Nc i の場合のフィルム肉厚 Cas Cas Cas Cas3 Cas4 Cas Film half width mm 図 5 Chill roll 位置フィルム肉厚計算結果の比較 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

13 粘弾性解析で評価されるセンターライン上の流速は 粘性解析結果と比較すると 全体的に低速になり 下流側で急激に加速される傾向を示します 4 Flow locit o th ctr li cm/s Dow stram distac mm Cas Cas Cas Cas3 Cas4 Cas5 図 6 センターライン上流速分布の比較 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd

14 -4.96~.64 KPa -.~.6 KPa -.488~3.93 KPa -.88~.39 KPa -.5~6.58 KPa.37~8.563 KPa Cas : Nwtoia flidpa s Cas3 : iss flid Pas 図 7 各種応力分布状態の比較 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 3

15 粘弾性解析の運用方法 4 自由表面計算の反復計算回数を指定します 粘弾性応力は各反復計算ステップ毎に正確に評価されます Film castig 計算パネルの実行ラジオボタンを ON とします 5 必要に応じて粘弾性応力解析の OR 反復計算回数と緩和係数を変更します 3 iss 粘弾性パラメータファイル Matrial fit で作成 を選択します iss 粘弾性解析のチェックボックスを ON とします 図 8 粘弾性フィルムキャスト解析用に解析実行タブメニューに設けられた入力項目 赤色表示 : 追加項目 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 4

16 粘弾性解析時には コンター図作画項目として粘弾性応力 tttt の 4 項目が追加されます 図 9 粘弾性フィルムキャスト解析時に追加されたポスト処理項目 Copright Hpr Adacd imlatio Laborator Co. Ltd. All Rights Rsrd 5

FEM原理講座 (サンプルテキスト)

FEM原理講座 (サンプルテキスト) サンプルテキスト FEM 原理講座 サイバネットシステム株式会社 8 年 月 9 日作成 サンプルテキストについて 各講師が 講義の内容が伝わりやすいページ を選びました テキストのページは必ずしも連続していません 一部を抜粋しています 幾何光学講座については 実物のテキストではなくガイダンスを掲載いたします 対象とする構造系 物理モデル 連続体 固体 弾性体 / 弾塑性体 / 粘弾性体 / 固体

More information

Microsoft PowerPoint - 夏の学校(CFD).pptx

Microsoft PowerPoint - 夏の学校(CFD).pptx /9/5 FD( 計算流体力学 ) の基礎理論 性能 運動分野 夏の学校 神戸大学大学院海事科学研究科勝井辰博 流体の質量保存 流体要素内の質量の増加率 [ 単位時間当たりの増加量 ] 単位時間に流体要素に流入する質量 流体要素 Fl lm (orol olm) v ( ) ガウスの定理 v( ) /9/5 = =( ) b=b =(b b b ) b= b = b + b + b アインシュタイン表記

More information

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際 Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際に 収束判定に関するデフォルトの設定をそのまま使うか 修正をします 応力解析ソルバーでは計算の終了を判断するときにこの設定を使います

More information

<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1>

<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1> 人工環境設計解析工学構造力学と有限要素法 ( 第 回 ) 東京大学新領域創成科学研究科 鈴木克幸 固体力学の基礎方程式 変位 - ひずみの関係 適合条件式 ひずみ - 応力の関係 構成方程式 応力 - 外力の関係 平衡方程式 境界条件 変位規定境界 反力規定境界 境界条件 荷重応力ひずみ変形 場の方程式 Γ t Γ t 平衡方程式構成方程式適合条件式 構造力学の基礎式 ひずみ 一軸 荷重応力ひずみ変形

More information

(Microsoft PowerPoint - \221\34613\211\361)

(Microsoft PowerPoint - \221\34613\211\361) 計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.

More information

弥生会計/やよいの青色申告

弥生会計/やよいの青色申告 c c c c c 1 c c c c c c c c c c c 2 3 c c 4 a a a a a a a a a

More information

パソコンシミュレータの現状

パソコンシミュレータの現状 第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に

More information

Microsoft PowerPoint - 2_FrontISTRと利用可能なソフトウェア.pptx

Microsoft PowerPoint - 2_FrontISTRと利用可能なソフトウェア.pptx 東京大学本郷キャンパス 工学部8号館2階222中会議室 13:30-14:00 FrontISTRと利用可能なソフトウェア 2017年4月28日 第35回FrontISTR研究会 FrontISTRの並列計算ハンズオン 精度検証から並列性能評価まで 観測された物理現象 物理モデル ( 支配方程式 ) 連続体の運動を支配する偏微分方程式 離散化手法 ( 有限要素法, 差分法など ) 代数的な数理モデル

More information

Fuji Xerox Co., Ltd. All rights reserved.

Fuji Xerox Co., Ltd. All rights reserved. 2011 Fuji Xerox Co., Ltd. All rights reserved. 2 2011 Fuji Xerox Co., Ltd. All rights reserved. 2011 Fuji Xerox Co., Ltd. All rights reserved. 2011 Fuji Xerox Co., Ltd. All rights reserved. 2011 Fuji Xerox

More information

untitled

untitled http://www.riskdatabank.co.jp The of Japan, Ltd. All rights reserved. 2 The of Japan, Ltd. All rights reserved. 3 The of Japan, Ltd. All rights reserved. 4 The of Japan, Ltd. All rights reserved. 5 The

More information

(c) PIXTA Co. Ltd. All Rights Reserved.

(c) PIXTA Co. Ltd. All Rights Reserved. 2006-2011 (c) PIXTA Co. Ltd. All Rights Reserved. 2006-2011 (c) PIXTA Co. Ltd. All Rights Reserved. 2006-2011 (c) PIXTA Co. Ltd. All Rights Reserved. 2006-2011 (c) PIXTA Co. Ltd. All Rights Reserved. 2006-2011

More information

PPTテンプレート集 ver.1.0

PPTテンプレート集 ver.1.0 Copyright 2012 Radishbo-ya co., Ltd All Rights Reserved. 1 !!!!!!!!!! Copyright 2012 Radishbo-ya co., Ltd All Rights Reserved. 2 Copyright 2012 Radishbo-ya co., Ltd All Rights Reserved. Copyright 2012

More information

BE-EPK732 2007 Panasonic Cycle Technology Co., Ltd. All Rights Reserved. 2 4 9 12 14 20 21 22 24 26 28 29 30 31 32 33 34 36 3 5 4 6 7 9 8 10 11 13 12 1. 5. 6. 3. 7. 2. 4. 15 14 16 17 18 19 20 21 1. 2.

More information

構造力学Ⅰ第12回

構造力学Ⅰ第12回 第 回材の座屈 (0 章 ) p.5~ ( 復習 ) モールの定理 ( 手順 ) 座屈とは 荷重により梁に生じた曲げモーメントをで除して仮想荷重と考える 座屈荷重 偏心荷重 ( 曲げと軸力 ) 断面の核 この仮想荷重に対するある点でのせん断力 たわみ角に相当する曲げモーメント たわみに相当する ( 例 ) 単純梁の支点のたわみ角 : は 図 を仮想荷重と考えたときの 点の支点反力 B は 図 を仮想荷重と考えたときのB

More information

Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1

Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1 Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1 Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 2 60 50 40 30 20 10 0 20173 20183 Copyright 2017 JAPAN POST BANK CO., LTD.

More information

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P P.11 P.14 P.15 P.16 P.16 P.17 P.19 P.20 P.22 P P P P P P P P P

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P P.11 P.14 P.15 P.16 P.16 P.17 P.19 P.20 P.22 P P P P P P P P P 201628 3 2016 5 13 P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P.10 2016 P.11 P.14 P.15 P.16 P.16 P.17 P.19 P.20 P.22 P.23 10 P.24 11 P.26 12 P.27 13 P.28 14 P.28 15 P.29 16 P.30 17 P.31 P.33 P.34 Copyright 2016 JAPAN

More information

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P.10 P.12 P.13 P.14 P.14 P.15 P.17 P.18 P.20 P P P P P.25 P.27 P.28 Copyright 2016 JAPAN POST BA

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P.10 P.12 P.13 P.14 P.14 P.15 P.17 P.18 P.20 P P P P P.25 P.27 P.28 Copyright 2016 JAPAN POST BA 201729 3 1 2016 8 12 P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P.10 P.12 P.13 P.14 P.14 P.15 P.17 P.18 P.20 P.21 10 P.22 11 P.23 12 P.24 13 P.25 P.27 P.28 Copyright 2016 JAPAN POST BANK CO., LTD. All Rights Reserved.

More information

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt シミュレーション工学 ( 後半 ) 東京大学人工物工学研究センター 鈴木克幸 CA( Compter Aded geerg ) r. Jaso Lemo (SC, 98) 設計者が解析ツールを使いこなすことにより 設計の評価 設計の質の向上を図る geerg の本質の 計算機による支援 (CA CAM などより広い名前 ) 様々な汎用ソフトの登場 工業製品の設計に不可欠のツール 構造解析 流体解析

More information

<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D>

<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D> 離散化手法とスキームの基礎 と選択法 007//6 宇宙航空研究開発機構情報 計算工学センター嶋英志 本講習の目的 基礎的な計算法の性質を述べ 各手法の持つ長所短所を理解することによって 手法の背景を理解した正しい選択に近づくこと クーラン数 風上差分 等の広い範囲の CFD 技術に共通の概念について その意味とイメージを把握すること 本講習の方針 様々な流体方程式の基礎となる移流方程式を用いて色々な計算法の特徴を計算例を示しながら解説する

More information

...Y..FEM.pm5

...Y..FEM.pm5 . 剛塑性有限要素法 名古屋大学大学院工学研究科. はじめに. 剛塑性体の構成式.. 降伏条件.. 構成方程式 ([D] マトリックス ). 節点速度 ひずみ速度関係..[B] マトリックス.. 四角形一次要素の [B] マトリックス.4 4 仮想仕事の原理 ( 剛性マトリックス ([K] マトリックス )).5 非線形方程式の解法.5. 直接代入法.5.wto-Raphso 法.6 非圧縮性の拘束と数値積分.7

More information

計算機シミュレーション

計算機シミュレーション . 運動方程式の数値解法.. ニュートン方程式の近似速度は, 位置座標 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます. 本来は が の極限をとらなければいけませんが, 有限の小さな値とすると 秒後の位置座標は速度を用いて, と近似できます. 同様にして, 加速度は, 速度 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます.

More information

NumericalProg09

NumericalProg09 数値解析および プログラミング演習 [08 第 9 回目 ] の解法 - 4. Ruge-Kua( ルンゲ クッタ 法 Ruge-Kua-Gill( ルンゲ クッタ ジル / ギル 法 5. 多段解法 解法の対象 常微分方程式 d( d 初期値条件 (, の変化に応じて変化する の値を求める. ( 0 ( 0 と 0 は,give 0 常微分方程式の初期値問題 と言う. 3 Ruge-Kua 法の導出

More information

OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二

OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二 OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 勉強会 @ 富山富山県立大学中川慎二 * OpenFOAM のソースコードでは, 基礎式を偏微分方程式の形で記述する.OpenFOAM 内部では, 有限体積法を使ってこの微分方程式を解いている. どのようにして, 有限体積法に基づく離散化が実現されているのか,

More information

Microsoft Word - NumericalComputation.docx

Microsoft Word - NumericalComputation.docx 数値計算入門 武尾英哉. 離散数学と数値計算 数学的解法の中には理論計算では求められないものもある. 例えば, 定積分は, まずは積分 ( 被積分関数の原始関数をみつけること できなければ値を得ることはできない. また, ある関数の所定の値における微分値を得るには, まずその関数の微分ができなければならない. さらに代数方程式の解を得るためには, 解析的に代数方程式を解く必要がある. ところが, これらは必ずしも解析的に導けるとは限らない.

More information

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P P P P P P P P.31

P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P P P P P P P P.31 201729 3 2017 5 15 P. 2 P. 4 P. 5 P. 6 P. 7 P. 9 P.10 2017 P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P.24 10 P.25 11 P.27 12 P.28 13 P.29 14 P.29 15 P.30 16 P.31 17 P.32 P.34 P.35 Copyright 2017

More information

スライド 1

スライド 1 H25 創造設計演習 ~ 振動設計演習 1~ 1 ゆれない片持ち梁の設計 振動設計演習全体 HP(2011 年度まで使用 今は閲覧のみ ): http://hockey.t.u-tokyo.ac.jp/shindousekkei/index.html M4 取付ネジ 2 Xin 加振器 50mm 幅 30mm 材料 :A2017または ABS 樹脂 計測点 :Xout 2mm? Hz CAD 所望の特性になるまで繰り返す?

More information

Microsoft PowerPoint - 知財報告会H20kobayakawa.ppt [互換モード]

Microsoft PowerPoint - 知財報告会H20kobayakawa.ppt [互換モード] 亀裂の変形特性を考慮した数値解析による岩盤物性評価法 地球工学研究所地圏科学領域小早川博亮 1 岩盤構造物の安定性評価 ( 斜面の例 ) 代表要素 代表要素の応力ひずみ関係 変形: 弾性体の場合 :E,ν 強度: モールクーロン破壊規準 :c,φ Rock Mech. Rock Engng. (2007) 40 (4), 363 382 原位置試験 せん断試験, 平板載荷試験 原位置三軸試験 室内試験

More information

Microsoft PowerPoint - suta.ppt [互換モード]

Microsoft PowerPoint - suta.ppt [互換モード] 弾塑性不飽和土構成モデルの一般化と土 / 水連成解析への適用 研究の背景 不飽和状態にある土構造物の弾塑性挙動 ロックフィルダム 道路盛土 長期的に正確な予測 不飽和土弾塑性構成モデル 水頭変動 雨水の浸潤 乾湿の繰り返し 土構造物の品質変化 不飽和土の特徴的な力学特性 不飽和土の特性 サクション サクション s w C 飽和度が低い状態 飽和度が高い状態 サクションの効果 空気侵入値 B. サクション増加

More information

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード]

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード] 圧密問題への逆問題の適用 一次元圧密と神戸空港の沈下予測 1. 一次元圧密の解析 2. 二次元圧密問題への適用 3. 神戸空港の沈下予測 1. 一次元圧密の解析 一次元圧密の実験 試験システムの概要 分割型圧密試験 逆解析の条件 未知量 ( 同定パラメータ ) 圧縮指数 :, 透水係数 :k 初期体積ひずみ速度 : 二次圧密係数 : 観測量沈下量 ( 計 4 点 ) 逆解析手法 粒子フィルタ (SIS)

More information

PowerPoint Presentation

PowerPoint Presentation Non-linea factue mechanics き裂先端付近の塑性変形 塑性域 R 破壊進行領域応カ特異場 Ω R R Hutchinson, Rice and Rosengen 全ひずみ塑性理論に基づいた解析 現段階のひずみは 除荷がないとすると現段階の応力で一義的に決まる 単純引張り時の応カーひずみ関係 ( 構成方程式 ): ( ) ( ) n () y y y ここで α,n 定数, /

More information

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように 3 章 Web に Link 解説 連続式 微分表示 の誘導.64 *4. 連続式連続式は ある領域の内部にある流体の質量の収支が その表面からの流入出の合計と等しくなることを定式化したものであり 流体における質量保存則を示したものである 2. 連続式 微分表示 の誘導図のような微小要素 コントロールボリューム の領域内の流体の増減と外部からの流体の流入出を考えることで定式化できる 微小要素 流入

More information

OpenCAE勉強会 公開用_pptx

OpenCAE勉強会 公開用_pptx OpenCAE 勉強会岐阜 2013/06/15 ABAQUS Student Edition を用い た XFEM き裂進展解析事例報告 OpenCAE 学会員 SH 発表内容 ABAQUS Student Edition とは? ABAQUS Student Edition 入手方法など - 入手方法 / インストール - 解析 Sample ファイルの入手方法 etc. XFEM について -XFEM

More information

85 4

85 4 85 4 86 Copright c 005 Kumanekosha 4.1 ( ) ( t ) t, t 4.1.1 t Step! (Step 1) (, 0) (Step ) ±V t (, t) I Check! P P V t π 54 t = 0 + V (, t) π θ : = θ : π ) θ = π ± sin ± cos t = 0 (, 0) = sin π V + t +V

More information

Microsoft Word 卒業論文2.doc

Microsoft Word 卒業論文2.doc 平成 6 年度 卒業論文 狭窄部を有する血管内の血流の有限要素解析 高知工科大学工学部知能機械システム工学科知能流体力学研究室 清水昌彦 目次 第 章緒言 - 本研究を行う背景と目的 - 血液の性質 -3 数値計算 - 有限要素法の概要 第 章基礎方程式 - 支配方程式 -- 連続の式 5 -- コーシーの運動方程式 6 --3 血液の構成方程式 6 - 無次元化 7 第 3 章解析手法 3- 有限要素解析

More information

位相最適化?

位相最適化? 均質化設計法 藤井大地 ( 東京大学 ) 位相最適化? 従来の考え方 境界形状を変化させて最適な形状 位相を求める Γ t Ω b Γ D 境界形状を変化させる問題点 解析が進むにつれて, 有限要素メッシュが異形になり, 再メッシュが必要になる 位相が変化する問題への適応が難しい Γ Γ t t Ω b Ω b Γ D Γ D 領域の拡張と特性関数の導入 χ Ω ( x) = f 0 f x Ω x

More information

FrontISTR による熱応力解析 東京大学新領域創成科学研究科人間環境学専攻橋本学 2014 年 10 月 31 日第 15 回 FrontISTR 研究会 < 機能 例題 定式化 プログラム解説編 熱応力解析 / 弾塑性解析 >

FrontISTR による熱応力解析 東京大学新領域創成科学研究科人間環境学専攻橋本学 2014 年 10 月 31 日第 15 回 FrontISTR 研究会 < 機能 例題 定式化 プログラム解説編 熱応力解析 / 弾塑性解析 > FronISR による熱応力解析 東京大学新領域創成科学研究科人間環境学専攻橋本学 214 年 1 月 31 日第 15 回 FronISR 研究会 < 機能 例題 定式化 プログラム解説編 熱応力解析 / 弾塑性解析 > FronISR に実装されている定式化を十分に理解し, 解きたい問題に対してソースコードを自由にカスタマイズ ( 要素タイプを追加, 材料の種類を追加, ユーザサブルーチンを追加

More information

国土技術政策総合研究所 研究資料

国土技術政策総合研究所 研究資料 3. 解析モデルの作成汎用ソフトFEMAP(Ver.9.0) を用いて, ダムおよび基礎岩盤の有限要素メッシュを8 節点要素により作成した また, 貯水池の基本寸法および分割数を規定し,UNIVERSE 2) により差分メッシュを作成した 3.1 メッシュサイズと時間刻みの設定基準解析結果の精度を確保するために, 堤体 基礎岩盤 貯水池を有限要素でモデル化する際に, 要素メッシュの最大サイズならびに解析時間刻みは,

More information

Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1

Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1 Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 1 Copyright 2017 JAPAN POST BANK CO., LTD. All Rights Reserved. 2 60 50 40 30 20 10 0 20173 20183 Copyright 2017 JAPAN POST BANK CO., LTD.

More information

第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える. 5.1 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f = l 2pL である. ただし, L [ 単位 m] は棒の長さ, [ 2 N / m ] 3 r[ 単位 Kg / m ] E r

第 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 information

A Luvens ICCG 未収束時にワーニングを出力するようにした A Luvens 非線形計算未収束時に計算をストップするようにした A Luvens 外部回路に電流源素子を追加 A Curie 浮き電極の境界条件を追加 A Hertz ポートの境界条件で差動ペアの設定が可能になった A Her

A Luvens ICCG 未収束時にワーニングを出力するようにした A Luvens 非線形計算未収束時に計算をストップするようにした A Luvens 外部回路に電流源素子を追加 A Curie 浮き電極の境界条件を追加 A Hertz ポートの境界条件で差動ペアの設定が可能になった A Her ********************************************************* Femtet2016.1.1 更新履歴 ********************************************************* [A] : 機能追加 [M] : 機能変更 [B] : バグ修正 A モデラ 応力解析 圧電解析の時 異方性材料の材料の座標軸を表示する

More information

第6章 実験モード解析

第6章 実験モード解析 第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法

More information

P. 2 P. 4 P. 5 P. 6 P. 7 P. 8 P. 9 P P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P P P P P P P.30 16

P. 2 P. 4 P. 5 P. 6 P. 7 P. 8 P. 9 P P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P P P P P P P.30 16 2018303 2018 5 15 P. 2 P. 4 P. 5 P. 6 P. 7 P. 8 P. 9 P.10 2018 P.11 P.13 P.15 P.16 P.17 P.17 P.18 P.20 P.21 P.23 P.24 10 P.25 11 P.27 12 P.28 13 P.29 14 P.29 15 P.30 16 P.32 1 2 20172016 Copyright 2018

More information

で通常 0.1mm 程度であるのに対し, 軸受内部の表面の大きさは通常 10mm 程度であり, 大きさのスケールが100 倍程度異なる. 例えば, 本研究で解析対象とした玉軸受について, すべての格子をEHLに用いる等間隔構造格子で作成したとすると, 総格子点数は10,000,000のオーダーとなる

で通常 0.1mm 程度であるのに対し, 軸受内部の表面の大きさは通常 10mm 程度であり, 大きさのスケールが100 倍程度異なる. 例えば, 本研究で解析対象とした玉軸受について, すべての格子をEHLに用いる等間隔構造格子で作成したとすると, 総格子点数は10,000,000のオーダーとなる 論文の内容の要旨 論文題目 転がり軸受における枯渇弾性流体潤滑とマクロ流れのマルチスケール連成解析手法の開発 氏名柴﨑健一 転がり軸受は, 転動体が, 外輪および内輪上の溝を転がることにより, 軸を回転自在に支持する機械要素であり, 長寿命化, 低摩擦化が強く求められている. 軸受の摩耗や焼付を防ぎ, 寿命を延ばすため, 通常は潤滑油またはグリースなどの潤滑剤が用いられる. 潤滑油は, 転がり接触する二表面間に表面粗さよりも厚い膜を形成し,

More information

大気環境シミュレーション

大気環境シミュレーション 第 3 回 (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.0 () 0 =.5 (3) 0 =.0 締切 04 年 月 6 日 ( 月 ) 夕方まで 提出先 347 室 オーバーフロー失敗ゴメンなさい (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.330 () 0 =.33 (3) 0

More information

LAN CATE TEL LAN CATE TEL 2 3

LAN CATE TEL LAN CATE TEL 2 3 C Matsushita Electric Industrial Co., Ltd. All Rights Reserved. TH-32LX80 TH-26LX80 TH-20LX80 LAN CATE TEL LAN CATE TEL 2 3 4 1 2 1 2 1 2 3 1 2 3 4 # $ % 6 7 1 1 2 3 2 3 4 8 9 10 1 3 2 11 4 4 6 4 6 1 2

More information

ブック 1.indb

ブック 1.indb 20 29 29 18 21 29 10 30 31 10 11 12 30 13 10 30 14 11 30 15 12 16 13 17 14 18 15 19 16 20 17 21 18 10 20 29 82 83 84 85 86 87 88 20 10 89 20 12 11 90 20 13 12 91 20 14 13 92 20 14 14 93 15 15 94 15 16

More information

問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた

問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた 問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた 測定データを図 1-2 に示す データから, オーステナイト系ステンレス鋼どうしの摩擦係数を推定せよ

More information

ハピタス のコピー.pages

ハピタス のコピー.pages Copyright (C) All Rights Reserved. 10 12,500 () ( ) ()() 1 : 2 : 3 : 2 4 : 5 : Copyright (C) All Rights Reserved. Copyright (C) All Rights Reserved. Copyright (C) All Rights Reserved. Copyright (C) All

More information

Copyright 2008 All Rights Reserved 2

Copyright 2008 All Rights Reserved 2 Copyright 2008 All Rights Reserved 1 Copyright 2008 All Rights Reserved 2 Copyright 2008 All Rights Reserved 3 Copyright 2008 All Rights Reserved 4 Copyright 2008 All Rights Reserved 5 Copyright 2008 All

More information

<94F E4F8EB25F >

<94F E4F8EB25F > JGS 5 土の三軸試験の供試体作製 設置 初期状態% 設)炉容器 No. 後供試体を用いる試験の基準番号と名称 JGS 51-9 土の繰返し非排水三軸試験 試 料 の 状 態 1) 乱さない 土粒子の密度 ρ s g/cm 供 試 体 の 作 製 ) トリミング 液 性 限 界 w L ) % 土 質 名 称 礫まじり粘土質砂 塑 性 限 界 w P ) % 1 5.1.96.98 質量 m i

More information

09.pptx

09.pptx 講義内容 数値解析 第 9 回 5 年 6 月 7 日 水 理学部物理学科情報理学コース. 非線形方程式の数値解法. はじめに. 分法. 補間法.4 ニュートン法.4. 多変数問題への応用.4. ニュートン法の収束性. 連立 次方程式の解法. 序論と行列計算の基礎. ガウスの消去法. 重対角行列の場合の解法項目を変更しました.4 LU 分解法.5 特異値分解法.6 共役勾配法.7 反復法.7. ヤコビ法.7.

More information

板バネの元は固定にします x[0] は常に0です : > x[0]:=t->0; (1.2) 初期値の設定をします 以降 for 文処理のため 空集合を生成しておきます : > init:={}: 30 番目 ( 端 ) 以外については 初期高さおよび初速は全て 0 にします 初期高さを x[j]

板バネの元は固定にします x[0] は常に0です : > x[0]:=t->0; (1.2) 初期値の設定をします 以降 for 文処理のため 空集合を生成しておきます : > init:={}: 30 番目 ( 端 ) 以外については 初期高さおよび初速は全て 0 にします 初期高さを x[j] 機械振動論固有振動と振動モード 本事例では 板バネを解析対象として 数値計算 ( シミュレーション ) と固有値問題を解くことにより振動解析を行っています 実際の振動は振動モードと呼ばれる特定パターンが複数組み合わされますが 各振動モードによる振動に分けて解析を行うことでその現象を捉え易くすることが出来ます そこで 本事例では アニメーションを活用した解析結果の可視化も取り入れています 板バネの振動

More information

施設・構造1-5b 京都大学原子炉実験所研究用原子炉(KUR)新耐震指針に照らした耐震安全性評価(中間報告)(原子炉建屋の耐震安全性評価) (その2)

施設・構造1-5b 京都大学原子炉実験所研究用原子炉(KUR)新耐震指針に照らした耐震安全性評価(中間報告)(原子炉建屋の耐震安全性評価) (その2) 原子炉建屋屋根版の水平地震応答解析モデル 境界条件 : 周辺固定 原子炉建屋屋根版の水平方向地震応答解析モデル 屋根版は有限要素 ( 板要素 ) を用い 建屋地震応答解析による最上階の応答波形を屋根版応答解析の入力とする 応答解析は弾性応答解析とする 原子炉建屋屋根版の上下地震応答解析モデル 7.E+7 6.E+7 実部虚部固有振動数 上下地盤ばね [kn/m] 5.E+7 4.E+7 3.E+7

More information

Microsoft Word - 4_構造特性係数の設定方法に関する検討.doc

Microsoft Word - 4_構造特性係数の設定方法に関する検討.doc 第 4 章 構造特性係数の設定方法に関する検討 4. はじめに 平成 年度 年度の時刻歴応答解析を実施した結果 課題として以下の点が指摘 された * ) 脆性壁の評価法の問題 時刻歴応答解析により 初期剛性が高く脆性的な壁については現在の構造特性係数 Ds 評価が危険であることが判明した 脆性壁では.5 倍程度必要保有耐力が大きくなる * ) 併用構造の Ds の設定の問題 異なる荷重変形関係を持つ壁の

More information

Microsoft PowerPoint - 第8章

Microsoft PowerPoint - 第8章 講義予定 案. 9/ 数値シミュレーションの手続き テキスト第 章. 9/ 9 偏微分方程式と解析解 テキスト第 章 3. 9/6 休講 4. 9/30 差分方程式とそのスキーム テキスト第 3 章 変換 テキスト第 4 章 5. 0/ 7 計算 テキスト第 5 章 連立一次方程式の解法 テキスト第 6 章 6. 0/ 流れ関数 ポテンシャルによる解法 テキスト第 7 章 7. 0/8 流速 圧力を用いた解法

More information

平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損

平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損 平成 3 年度 JAXA 航空プログラム公募型研究報告会資料集 (3 年度採用分 1 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損傷発生位置の推定発表内容 (1 荷重同定 1:11 点衝撃荷重同定 ( 荷重同定 : 分布荷重同定 (3 今後の予定

More information

高砂2013表紙

高砂2013表紙 Takasago Thermal Engineering Co.,Ltd. Copyright 2015 Takasago Thermal Engineering Co., Ltd. All Rights Reserved. 2015 CONTENTS 1 Green Air Technology Profile 2015 Green Air Technology Profile 2015 2 !

More information

静的弾性問題の有限要素法解析アルゴリズム

静的弾性問題の有限要素法解析アルゴリズム 概要 基礎理論. 応力とひずみおよび平衡方程式. 降伏条件式. 構成式 ( 応力 - ひずみ関係式 ) 有限要素法. 有限要素法の概要. 仮想仕事の原理式と変分原理. 平面ひずみ弾性有限要素法定式化 FEM の基礎方程式平衡方程式. G G G ひずみ - 変位関係式 w w w. kl jkl j D 構成式応力 - ひずみ関係式 ) (. 変位の境界条件力の境界条件境界条件式 t S on V

More information

NS NS Scalar turbulence 5 6 FEM NS Mesh (A )

NS NS Scalar turbulence 5 6 FEM NS Mesh (A ) 22 3 2 1 2 2 2 3 3 4 NS 4 4.1 NS............ 5 5 Scalar turbulence 5 6 FEM 5 6.1 NS.................................... 6 6.2 Mes A )................................... 6 6.3.....................................

More information

1 911 34/ 22 1012 2/ 20 69 3/ 22 69 1/ 22 69 3/ 22 69 1/ 22 68 3/ 22 68 1/ 3 8 D 0.0900.129mm 0.1300.179mm 0.1800.199mm 0.1000.139mm 0.1400.409mm 0.4101.199mm 0.0900.139mm 0.1400.269mm 0.2700.289mm

More information

液晶ディスプレイ取説TD-E432/TD-E502/TD-E552/TD-E652/TD-E432D/TD-E502D

液晶ディスプレイ取説TD-E432/TD-E502/TD-E552/TD-E652/TD-E432D/TD-E502D 1 2 3 4 5 6 7 1 2 3 4 5 6 7 2 2 2 1 1 2 9 10 11 12 13 14 15 16 17 1 8 2 3 4 5 6 7 1 2 3 4 5 6 7 8 9 10 9 11 12 13 13 14 15 16 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 1 2 3 4 5 6 7 8 9 11 12

More information

000-.\..

000-.\.. 1 1 1 2 3 4 5 6 7 8 9 e e 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 10mm 150mm 60mm 25mm 40mm 30mm 25 26 27 1 28 29 30 31 32 e e e e e e 33 e 34 35 35 e e e e 36 37 38 38 e e 39 e 1 40 e 41 e 42 43

More information

1 2 http://www.japan-shop.jp/ 3 4 http://www.japan-shop.jp/ 5 6 http://www.japan-shop.jp/ 7 2,930mm 2,700 mm 2,950mm 2,930mm 2,950mm 2,700mm 2,930mm 2,950mm 2,700mm 8 http://www.japan-shop.jp/ 9 10 http://www.japan-shop.jp/

More information

第18回海岸シンポジウム報告書

第18回海岸シンポジウム報告書 2011.6.25 2011.6.26 L1 2011.6.27 L2 2011.7.6 2011.12.7 2011.10-12 2011.9-10 2012.3.9 23 2012.4, 2013.8.30 2012.6.13 2013.9 2011.7-2011.12-2012.4 2011.12.27 2013.9 1m30 1 2 3 4 5 6 m 5.0m 2.0m -5.0m 1.0m

More information

1 1 36 223 42 14 92 4 3 2 1 4 3 4 3429 13536 5 6 7 8 9 2.4m/ (M) (M) (M) (M) (M) 6.67.3 6.57.2 6.97.6 7.27.8 8.4 5 6 5 6 5 5 74 1,239 0 30 21 ( ) 1,639 3,898 0 1,084 887 2 5 0 2 2 4 22 1 3 1 ( :) 426 1500

More information

1 C 2 C 3 C 4 C 1 C 2 C 3 C

1 C 2 C 3 C 4 C 1 C 2 C 3 C 1 e N >. C 40 41 2 >. C 3 >.. C 26 >.. C .mm 4 C 106 e A 107 1 C 2 C 3 C 4 C 1 C 2 C 3 C 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124

More information

(1519) () 1 ( ) () 1 ( ) - 1 - - 2 - (1531) (25) 5 25,000 (25) 5 30,000 25,000 174 3 323 174 3 323 (1532) () 2 () 2-3 - - 4 - (1533) () 1 (2267)204 () (1)(2) () 1 (2267)204 () (1)(2) (3) (3) 840,000 680,000

More information

平成24年財政投融資計画PDF出後8/016‐030

平成24年財政投融資計画PDF出後8/016‐030 24 23 28,707,866 2,317,737 26,390,129 29,289,794 2,899,665 24 23 19,084,525 21,036,598 1952,073 24 23 8,603,613 8,393,427 967,631 925,404 202,440 179,834 217,469 219,963 66,716 64,877 3,160,423 2,951,165

More information

[mm] [mm] [mm] 70 60 50 40 30 20 10 1H 0 18 19 20 21 22 23 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 1 2 3 4 5 6 7 8 9 10 11 12 60 50 40 30 20 10 0 18 19 20 21 22 23 24 1 2 3 4

More information

STC-IP10

STC-IP10 1 2 3 1 1 4 5 6 7 8 1 9 2 3 4 5 6 10 7 8 11 12 9 10 11 12 13 13 1 2 14 3 4 5 6 7 15 8 16 9 17 10 11 12 18 13 14 15 19 1 2 20 1 2 21 3 4 5 22 1 2 3 4 23 1 2 1 24 2 3 4 5 25 1 2 3 26 1 2 3 4 27 1 2 28 1

More information