本日の説明内容本日の説明内容本日の説明内容本日の説明内容 超音波ソフトウェア 超音波ソフトウェア 超音波ソフトウェア 超音波ソフトウェア ComWAVE ComWAVE ComWAVE ComWAVE 概要紹介概要紹介概要紹介概要紹介 車両用超音波センサ 車両用超音波センサ 車両用超音波センサ 車両

Transcription

1 自動車業界における超音波 シミュレーションの現状 (2) 科学システムサポートチーム 2015 年 2 月 19 日 Copyright (c)2006 ITOCHU Techno-Solutions Corporation

2 本日の説明内容本日の説明内容本日の説明内容本日の説明内容 超音波ソフトウェア 超音波ソフトウェア 超音波ソフトウェア 超音波ソフトウェア ComWAVE ComWAVE ComWAVE ComWAVE 概要紹介概要紹介概要紹介概要紹介 車両用超音波センサ 車両用超音波センサ 車両用超音波センサ 車両用超音波センサ 各種パーツの超音波による検査計測 各種パーツの超音波による検査計測 各種パーツの超音波による検査計測 各種パーツの超音波による検査計測 自動車業界に 自動車業界に 自動車業界に 自動車業界におけるおけるおけるおける超音波技術の利用超音波技術の利用超音波技術の利用超音波技術の利用 1 モデリング機能 モデリング機能 モデリング機能 モデリング機能 超音波とは 超音波とは 超音波とは 超音波とは 超音波シミュレータの活用方法超音波シミュレータの活用方法超音波シミュレータの活用方法超音波シミュレータの活用方法 自動車業界に 自動車業界に 自動車業界に 自動車業界におけるおけるおけるおける ComWAVE ComWAVE ComWAVE ComWAVE の活用活用活用活用

3 超音波とは 2

4 超音波とは ( 固体 気体中を伝搬する可聴域より高周波の波 ) 超音波デバイス 超音波検査 計測 超音波診断 治療 大腿骨 ヒト大腿部シミュレーション 周波数 :1GHz 100kHz~10MHz 100kHz~10MHz 波長 : 1µm 1mm 10mm 超音波接合 洗浄 可聴域 地震波 1.4E+00 周波数 : 10~ 100kHz 20Hz~20kHz 10Hz~0.1Hz 波長 : 50mm 100mm 10m 超音波領域 3

5 超音波の工学利用における 3 つのフェイズ 1. 超音波発信 - 主に圧電材を含む素子の特性を反映 : 電気特性 + 機械特性 - 材料値 構造 サイズの最適化による広帯域化 不要モード除去 指向性改善 出力増大 2. 超音波伝搬 - 試験体の特性を反映 : 機械特性 - 対象の材質 形状 サイズ 傷による伝搬経路 励起モードへの影響 超音波利用 3. 超音波受信 - 素子の特性を反映 : 電気特性 + 機械特性 発信 伝搬 受信 これらの課題解決のために超音波シミュレーションを援用することが 超音波現象解明の為の標準的な手法になりつつある 4

6 自動車業界における超音波技術の利用 5

7 自動車分野における超音波技術の利用 アルミ鋳造品検査 エンジンブロック クランクシャフト その他鋳造物の鋳巣の検査やスポット溶接部の検査 ハイブリット車電池電極接合 塗装剥離検査 各種部品の超音波洗浄 超音波検査 ケーブル接合 ( ワイヤハーネス ) 電子機器接合 車室内センサ 開閉部安全センサー 後方監視センサ 各種溶接部探傷 前方監視センサ 車軸のき裂探傷 ブレーキパッドのクラック 剥離検査 6

8 超音波ソフトウェア ComWAVE 概要紹介 7

9 超音波シミュレーションソフト ComWAVE TM ComWAVEは 超音波の見える化により 製品 設備 各種構造物の超音波検査 計測の最適化 各種超音波装置の設計 超音波計測ノイズの評価など 様々な超音波現象を評価できます 溶接部探傷解析 レールの探傷解析 (B スコーフ 表示 ) 3 次元音場分布図 超音波洗浄装置ワイヤーボンダ装置超音波流量計 8

10 有限要素法を採用し 超音波伝搬現象を忠実に再現 超音波探傷で使用する全ての超音波モードを取り扱うことができます クリーピング波 横波 縦波 モード変換波 レイリー波 ガイド波 回折波 9

11 ComWAVE TM 主な機能 1. 有限要素法を採用し 超音波伝搬現象を忠実に再現 2. 異方性 不均質性材料を忠実にモデル化 3. 数十億要素の大規模問題への対応によりフル3D 超音波解析を実現 4.Aスコープ Bスコープ表示機能により 実験結果と容易に比較可能 5. 拡張リスタート機能 6. プローブ遅延時間自動計算機能 7. レイリー減衰機能 8. FEM- 外装法ハイブリッド解析機能 9. 超音波速度入力による異方性材料定義機能 10. き裂等接触解析機能 11. 流れ場との連成解析機能 12. JIS 試験片プローブテンプレート機能 13. ComWAVEの並列化効果とGPGPU 対応 14. 入出力機能 10

12 自動車業界における ComWAVE の活用 車両用超音波センサ 11

13 自動車業界における ComWAVE の活用 エンジンブロック クランクシャフト 車軸のき裂探傷 その他鋳造物の鋳巣の検査やスポット溶接部の検査 ケーブル接合 ( ワイヤハーネス ) 電子機器の電極電極接合 車室内センサ メッシュレス領域 後方監視 開閉部安全センサー 塗装剥離検査 角度 30 外挿元 前方監視センサ 各種部品の超音波洗浄 流量計測 12

14 超音波センサ解析事例 3 次元超音波センサ音場解析 解析モデル 40kHz パルス 要素サイズ :0.2[mm] 解析対象領域 V:60[mm] 60[mm] 63[mm] 要素数 :2,835 万要素 計算時間 :2e-4 [sec] 時間ステップ刻み幅 : e-007 [sec] ステップ数 :T/ t ~ 424 必要メモリ :~ 2.0GB 計算時間 :4 コア版 約 2,194 秒 センサ形状による音場の広がり方の違いを 音場分布の可視化により確認 センサ 1: 円形振動子 :Φ10mm センサ 2: 方形振動子 :5mm 10mm 13

15 センサ解析 : 空中超音波センサの音場解析 円形振動子 :Φ10mm および方形振動子 :5mmx10mm 円形振動子は 上下 左右どちらにも同じ拡がりを持って伝搬する 方形振動子は 上下方法より 左右に拡がりを持って伝搬する 左図 : 振動子変形図 右図 : 空中への超音波伝搬図 車両設置時の超音波伝搬解析バンパー等設置環境による影響を考慮要素数 :3.5 億要素 2000 ステップ 14

16 素子振動解析と ComWAVE の結合解析例 ( 超音波センサー ) 適切なセンサ形状 周波数 センサ適用エリア検討等に貢献素子振動振動解析または実験により取得圧電体の変位振動波形を得る 40kHz の電圧波形入力 ComWAVE による超音波伝搬解析 圧電体上に 得られた変位振動波形を入力 1.0 m 超音波伝搬解析 15

17 ハイブリッド解析 (FEM- キルヒホッフ法 ) 機能 フル FEM モデル 200 解析領域 ( 要素サイズ 0.25mm 約 8 億要素 ) 40 khz ウェーブレット 3 波 λ = 8.6 mm 300mm( 約 35λ) 200 φ10mm 縦波音速 :343e3 [mm/s] 円形振動子 密度 :1.24e-9 [kg/mm 3 ] 310 観測点 Step Step1 (FEM) Step2 ( 外挿 ) 使用メモリ 計算時間 47 GB 29 時間 (8 並列 ) - - FEM+ 外挿計算モデル X Y 単位 :mm Z Step 1:FEM 計算 ( 送信 ) FEM 領域 ( mm) 要素サイズ :0.25mm 要素数:1638 万 伝搬距離 50mm 外挿元領域積分面 :50 50mm 要素サイズ:0.25mm Source points:40000 点 (=50/ /0.25) Step 2: 外挿計算計算量 :40000 点 1 点 伝搬距離 300mm 観測点 Step Step1 (FEM) Step2 ( 外挿 ) 使用メモリ 計算時間 880 MB 8 分 (4 並列 ) 160 MB 1 分未満 (1 並列 ) 計算環境 (CPU:Intel Xeon X5675@3.07GHz) 計算規模 計算速度 計算精度すべての面で改善を見込める 伝搬距離が長くなるにつれて フルFEMの波形は崩れていく

18 超音波センサ解析事例 超音波アレイセンサ解析例 3 次元大空間超音波センサの送信 伝搬 受信解析 センサ配置検討 見逃し防止等に貢献 メッシュレス領域 縦波音速 :340e3 [mm/s] 密度 :1.24e-9 [kg/mm 3 ] 320 mm 円形振動子 (φ10mm) をアレー状に配置 角度 30 アレイ開口長 100mm 外挿元 40 khz ウェーブレット 3 波 λ = 8.6 mm 角度に応じて 遅延時間をそれぞれ入力 圧力コンター表示 ( 水平面 ) 解析モデル (3 つの外挿元を指定 ) 参考文献 : 羽多野裕之, 山里敬也, 片山正昭, 自動車用近距離レーダネットワークのための超音波アレイエミッタ, 自動車技術, Vol.62, No.7, pp (2008) 17

19 超音波センサ解析 2 次元車室内超音波伝搬解析 解析モデル 40kHz 連続波 要素サイズ :1.0[mm] 解析対象領域 V:6000[mm] 2100[mm] 要素数 :1260 万要素 計算時間 :2e-2 [sec] 時間ステップ刻み幅 : e-006 [sec] ステップ数 :T/ t ~ 8500 必要メモリ :~ 1.2GB 計算時間 :4 コア版 約 1.6 時間 送信子傾き 45 度 周波数 40kHz Φ10mm 2300mm 受信子傾き - 45 度 Φ10mm 反射板 6000 mm 反射板と車体との間に隙間あり 車室内に取り付けられた送信センサから照射された超音波の伝搬 反射板による反射について車室内の音場分布の可視化により確認し 受信センサでの信号強度の評価を行う mm

20 超音波センサ解析 2 次元車室内超音波伝搬解析 解析結果 (Z1: 空気の音響インピーダンス Z2: 反射板材料の音響インピーダンス ) (a) Z1に比較して Z2がかなり大きい場合 (Z1<<Z2) 反射板で全反射 送信子 受信子 反射板 全反射 反射率 =(Z2-Z1)/(Z1+Z2)=1.0 (b) Z1に比較して Z2が5 倍程度の場合 (Z1:Z2=1:5) 反射板の一部で吸収 一部で吸収 送信子 受信子 反射板 反射率 =(Z2-Z1)/(Z1+Z2)=4/6=0.667 受信波形音響インピーダンス ( 材料の密度と音速の積で算出 ) を指定することで 材料の吸収率の考慮が可能 全反射の場合に比べて 振幅値が小さくなっているのを確認できる 19

21 超音波センサ分野 : 3 次元車室内超音波伝搬解析 車室内超音波センサ解析 解析規模 : 要素数 :2.4 億要素ステップ数 :5,000 必要メモリ :6GB 計算時間 : 約 25 時間 (1CPU コア, 圧力波計算 ) 計算時間 : 約 1.7 時間 (1GPU GPU(K40), 圧力波計算 ) GPU 計算は CPU と比較して約 15 倍高速! 車室内の様々な反射源からの反射 吸収を考慮した 超音波伝搬状況を把握可能 3 次元音場分布図 20

22 自動車業界における ComWAVE の活用 各種パーツの超音波による検査計測 21

23 各種パーツの超音波による検査計測 エンジンブロック クランクシャフト 車軸のき裂探傷 その他鋳造物の鋳巣の検査やスポット溶接部の検査 ケーブル接合 ( ワイヤハーネス ) 電子機器の電極電極接合 車室内センサ メッシュレス領域 塗装剥離検査 後方監視 開閉部安全センサー 角度 30 外挿元 前方監視センサ 各種部品の超音波洗浄 流量計測 22

24 ワイヤーボンディング解析例 ( 超音波接合 ) 3 次元モデル全要素数 : 約 6 億要素 メッシュ幅 :0.01mm 入力周波数 :100kHz 金線 φ30μm 15.5mm 3mm 11.2mm 先端部半透明表示 23

25 3 次元試計算結果 ( 金線とチップとの接触がある場合金線とチップの接触解析機能を使用 ) 超音波による配線接合 電池等の電極部接合の効率化で利用下図は 14μsでの変形図 ( 変位 :30 倍に拡大 青 ~ 赤は0~30μm) 全要素数 :208,435,744 ステップ数 :18,473 (25μs) 要素大きさ :10μm 計算時間 :10.5 時間 (4GPU) 73.5 時間 (4CPU コア ) 24

26 突合せ継ぎ手溶接部探傷シミュレーション例 鋳造物の鋳巣の検査やスポット溶接部の検査 超音波プローブ適用範囲決定等に利用 inspection with longitudinal wave 5,000 万要素 2,739 steps (30us) 2GPU の並列計算で約 13 分 (CPU 計算では約 10 時間 ) 探傷周波数 :2MHz( 横波入射角度 60 度 ) inspection with share wave

27 TOFD 法による超音波探傷解析 鋳造物の鋳巣検査や溶接部巣検査や溶接部の検査および内部欠陥検出装置開発で利用 送信 (2MHz, 横波 45 度 ) 受信 き裂上端エコー き裂下端エコー 差分 ( き裂ありーき裂なし ) き裂上端エコー き裂下端エコー 26

28 フェーズドアレイによる欠陥の画像化解析 画像化機能 ( 開口合成表示機能 ) アレイ素子 :32ch 2MHz ピッチ2.4mm 水鋼材 鋼材表面平面 合成開口表示 : 欠陥 電子部品等の欠陥検出で利用 複雑形状部分や微小欠陥の検出に有利 1 素子ずつ送信 全素子受信を 32ch 繰り返す 鋼材全領域に焦点を当てた計算を実施し合成開口画像を作成 鋼材表面曲面 合成開口表示 : 従来の B スキャン画像 : 5 素子の電子操作 焦点位置は欠陥中央部 27

29 板波伝搬解析モデル (3 次元解析 ) ボディー周りの傷 塗装剥離検出検査効率化で利用 航空機関係では 複合材検査で利用 時刻歴変位波形 : φ60 入力中心周波数 :1MHz 3mm 変位ベクトル図 : 15µSec 変位ベクトル図アニメーション : 18µSec 非対称モード 27µSec 対称モード 28

30 車軸の超音波探傷シミュレーション例 鉄道車両の車軸を模擬した小型輪軸での探傷のシミュレーション例 ( 鉄道総研様ご発表論文 1) 車輪と車軸の接合部に生じたき裂を探傷 車輪を外した検査と外さない検査での探傷波形の違いを確認 実験を行うことが困難な条件をシミュレーションで代用 適切な検査方法を素早く検討することが可能 10.7μs 31.0μs 48.2μs 縦波 横波 車輪のコーナで反射した横波 き裂で反射した横波 車軸 - 車輪間で反射した横波車輪に透過する横波 約 10.7 億要素 解析モデル (3 次元モデル ) 解析結果 ( 上図 : 超音波伝搬図 下図 : 探傷波形 ) 解析結果では 車輪を外さない場合は 外した場合に比べてエコー高さが10.6dB 低下した 実験では13.1dBであり よい一致を示している なお実験との違いは 解析では車輪装着時にき裂周辺の面圧や摩擦力の分布を考慮していないためと考えられる 1:Makino etc. Modeling and Simulation of Ultrasonic Testing on Miniature Wheelset. International Congress on Ultrasonics. AIP Conf Proc. 1433, (2012). 29

31 超音波洗浄機中の超音波伝搬解析 超音波により各種部品の洗浄効率化での利用 水中の圧力波伝搬 モデル 1 : ガラスを振動子に平行に設置モデル 2 : ガラスを傾けて設置 要素サイズ :0.5mm 全要素数 :6,060,000 ステップ数 :4,275 (300μs) 計算時間 :3.5 時間 Pentium D 2.8GHz 1 使用メモリ : 約 420MB OS: Linux 30

32 水の流れを考慮した超音波流量計の解析例 (3 次元解析 )(1/2) 配管内の流れ計測効率化での利用解析モデル管 : 内径 = 50mm 外径 = 56mm ( 板厚 = 3mm) 水の流れ送信子 :100kHz 振動子 :φ = 8mm 流れ速度 (km/h) 流れ場速度分布図 (y-z 面 ) 水 流れ速度 (km/h) 720 受信子 : 振動子 :φ = 8mm 流れ場速度分布図 (x-z 面 ) 0 31

33 水の流れを考慮した超音波流量計の解析例 (3 次元解析 )(2/2) Flow3ddistribution の実行結果 32

34 空中 ( 水中 ) 超音波解析 空中 ( 水中 ) 超音波 特長 例 対象物に非接触で超音波を送受信できる 減衰により空中は 水中より遠くまで伝わらない ( 空気中では 100m 程度 水中では 10km 程度伝搬 ) 空気を介した非破壊検査のエコーは微弱なため 大出力 高感度なセンサ開発が必須 コウモリ イルカのエコーロケーション 自動車等の超音波式バックセンサー 非破壊検査 : 接触媒質を使えない検査に適している 超音波診断装置 ( 医療 ) 等 33

35 空中 ( 水中 ) 伝搬 - 反射解析手順 全領域 FEM 計算では不可能な大規模空間の解析を可能にする手法 FEM 計算領域 音場引継ぎ領域 : 領域内の全要素に ある指定時刻の音場 ( スナップショット ) または境界波形を引き渡す 受信センサ 送信センサ FEM 計算領域 一様媒質 往路伝搬経路 解析手順 step1 : 送信計算 ( 有限要素法 ) step2 : 往路伝搬計算 ( 外挿法 ) step3 : 反射計算 ( 有限要素法 ) step4 : 復路伝搬計算 ( 外挿法 ) step5 : 受信計算 ( 有限要素法 ) FEM 計算領域 反射体 計算規模 計算速度 計算精度すべての面で改善を見込める

36 非接触空中超音波による探傷解析例 Step1(FEM) 100 x 100 x 100 mm 40kHz ウェーブレット3 波 積分面 λ = 8.6 mm Step2( 外挿 ) 1m (117λ 117λ) Step3(FEM) 240 x 240 x 100 mm リスタート領域ボイド 空気の泡 水など観測面貼り合わせたアクリル板 100 x 100 x 10 mm 水 検査体モデル 縦波音速 1480 m/s 密度 1.0e-6 kg/mm3 空気の泡ボイド縦波音速 343 m/s 密度 1.24e-9 kg/mm3 アクリル板 縦波音速 2700 m/s 横波音速 1400 m/s 密度 1.2e-6 kg/mm3 Step 計算タイプ ステップ数 要素数 計算量 Source points Target points 使用メモリ 計算時間 Step1 FEM 計算 万 GB 2 分 (4 並列 ) Step2 外挿計算 ( リスタートファイル出力 ) 万 1.4 GB 30 分 (4 並列 ) Step3 FEM 計算 億 GB 3.2 時間 (4 並列 ) 計算環境 (CPU:Intel Xeon X5675@3.07GHz) 参考 URL: 35

37 非接触空中超音波による探傷解析例電池など水を使えない材料の検査に有利 様々な材料を超音波が透過透過する際の特性 ( 反射 透過率 水や空気による屈折 回折効果等 ) の評価に貢献圧力コンター図 ( 上段 : 平面図 下段 : 側面図 ) 空気の泡 ボイドの箇所の圧力が大きくなっている ( 赤くなっている ) のが分かる 空気の泡 水 観測面 アクリル板 ボイド 平面図 : アクリル板の下半分の中央断面側面図 : アクリル板の中央断面 36

38 モデリング機能 37

39 基本形状組み合わせによるモデル化基本形状組み合わせによるモデル化基本形状組み合わせによるモデル化基本形状組み合わせによるモデル化立方体立方体立方体立方体円柱円柱円柱円柱だ円体だ円体だ円体だ円体四角錐四角錐四角錐四角錐回転回転回転回転だ円体円体円体円体トーラストーラストーラストーラス平面平面平面平面球形球形球形球形組合せモデル組合せモデル組合せモデル組合せモデル FEM FEM FEM FEM メッシュ生成メッシュ生成メッシュ生成メッシュ生成 38

40 汎用プリプロセッサ GiD によるモデル化 任意形状モデリング CAD データからのモデル取込 [ 対応フォーマット ] IGES ファイル,DXF ファイル, Parasolid,ACIS ファイル,VDA ファイル, Rhino ファイル,Shape ファイル NASTRAN メッシュ,STL メッシュ, VRML メッシュ,3D Studio, GID メッシュ 39

41 JIS 試験片テンプレート機能本テンプレートの特長 JIS Z 3060 に準拠した検査体および溶接部モデルテンプレートを実装 各種探触子と検査体等の組み合わせも自由自在な自動モデル結合機能 40

42 写真取り込みによるモデル化き裂等の複雑形状を詳細にモデル化が可能 突合せ溶接部テンプレートによるモデル化 減肉設定用の写真 41 き裂設定用の写真 テンプレートへの写真の適用メッシュ生成 ( 減肉 き裂をモデル化 ) 41

43 写真取り込みによるモデル化機能道路面を忠実に再現 路面での反射状況解析に利用 コンクリート材料写真 ( 1) 材料 ID:1 材料 ID:2 材料 ID:3 骨材あり 6 種類の材料 8 諧調でモデル化 材料 ID:4 材料 ID:5 材料 ID:6 骨材なし 1: フリー写真素材サイト somephoto から引用 42

44 モデル作成 表示機能 複数プロジェクト同時表示 ツリー構造表示ツリー構造で各プロジェクトを表示でき 操作性が向上 Python 機能 Python で開発された他のコードも使用可能 43

45 モデル作成 表示機能 B スコープ可視化画面 44

46 超音波シミュレータの活用方法 45

47 超音波シミュレータの活用方法超音波シミュレータの活用方法超音波シミュレータの活用方法超音波シミュレータの活用方法 ComWAVE ComWAVE ComWAVE ComWAVE データベースデータベースデータベースデータベース超音波センサ設計時超音波センサ設計時超音波センサ設計時超音波センサ設計時目的にあった センサ形状 周波数 材質等を試験体を作らず 様々なパラメータについてシミュレーションで検討可能 これによりこれによりこれによりこれにより 設計工数の削減設計工数の削減設計工数の削減設計工数の削減 高度な最適設計度な最適設計度な最適設計度な最適設計 センサセンサセンサセンサ適用範囲適用範囲適用範囲適用範囲拡大を実現大を実現大を実現大を実現 センサ運用時センサ運用時センサ運用時センサ運用時センサ利用時の様々な状況 ( センサ音場やエコー ) を事前にシミュレーションしデータベース化 センサ運用時にシミュレーションデータベースと比較 分析実施 これによりこれによりこれによりこれにより 誤検知の低減 高精 誤検知の低減 高精 誤検知の低減 高精 誤検知の低減 高精度な検知 高信頼性を実現 度な検知 高信頼性を実現 度な検知 高信頼性を実現 度な検知 高信頼性を実現 センサセンサセンサセンサ技術者の育成時技術者の育成時技術者の育成時技術者の育成時超音波を取り扱うにあたり 注意すべき点や回折 干渉等の理解が難しい事例をシミュレーションによりデータベース化し育成時利用 これによりれによりれによりれにより 効率的な技術者育成効率的な技術者育成効率的な技術者育成効率的な技術者育成 超音波特性の深い理解超音波特性の深い理解超音波特性の深い理解超音波特性の深い理解 シミュシミュシミュシミュレーションと実験を比較レーションと実験を比較レーションと実験を比較レーションと実験を比較し 実験技実験技実験技実験技術向上を実現術向上を実現術向上を実現術向上を実現 2 つのセンサによる超音波の 2 つのセンサによる超音波の 2 つのセンサによる超音波の 2 つのセンサによる超音波の干渉現象の可視化干渉現象の可視化干渉現象の可視化干渉現象の可視化ネットワークネットワークネットワークネットワーク超音波アレイソナー解析によるセンサ音場解析超音波アレイソナー解析によるセンサ音場解析超音波アレイソナー解析によるセンサ音場解析超音波アレイソナー解析によるセンサ音場解析円形円形円形円形 方形センサの違いによる方形センサの違いによる方形センサの違いによる方形センサの違いによる 振動および超音波音場の変化振動および超音波音場の変化振動および超音波音場の変化振動および超音波音場の変化 46

48 お知らせ 非破壊評価総合展 に出展 メンテナンス レジリエンス TOKYO 2015 ( 旧 : ものづくり NEXT ) 公式ホームページ : 日時 :2015 年 7 月 22 日 ( 水 )~2015 年 7 月 24 日 ( 金 ) 開催場所 : 東京ビッグサイト 来場の際には弊社ブースへお立寄りください 47

49 お問合せ ~ アプリケーション導入からシミュレーションの自動化 最適化まで ~ あなたの業務をシミュレーションにより実現します 科学システム事業部 CAE ソリューション営業部 本社住所東京都千代田区霞が関 霞が関ビル 東京 大阪 TEL: /FAX: TEL: /FAX: URL: ctc comwave 検索