JP2003512677A

JP2003512677A - 構造解析のための装置及び方法

Info

Publication number: JP2003512677A
Application number: JP2001532434A
Authority: JP
Inventors: ファン，ツィー，リアン; ロングツェング; ケネディー，ピーター，ケネス; ユ，フア，ガン
Original assignee: モールドフロー・ピーティーワイ・リミテッド
Priority date: 1999-10-15
Filing date: 2000-10-13
Publication date: 2003-04-02
Also published as: US6704693B1; KR100805970B1; CN1399750A; JP2011103138A; BRPI0015006B1; EP1247217A1; EP1247217A4; HK1050408A1; KR20030004299A; CA2387055A1; DE60045640D1; BR0015006A; ATE498868T1; JP5272025B2; WO2001029712A1; AU7893200A; AU773165B2; CN1211748C; NZ518409A; EP1247217B1

Abstract

(57)【要約】物体の応力は、表面メッシュに確定される対向するシェル要素（２０、２２）を介して、有限要素分析によってモデリングされる。表面メッシュは、物体の外部表面の表示上に生成される。複数の偏心シェル要素（２０、２２）は、シェル要素（２０、２２）における物体の厚さの半分に等しい厚さで、外部表面上のノード（ａ，ｂ，ｃ，ｄ，ｅ，ｆ，ｇ，ｈ）によって画定される。多点制約は、対向するシェル要素（２０、２２）に対応するノード間に確立される。動的、静的、あるいは熱的負荷のような境界条件及び荷重条件が特定される。そして、対応する移動、ひずみ、及び応力が計算される。その方法は、例えば、薄い壁を有する注入成形部分の構造的及びそり解析に適用できる。

Description

【発明の詳細な説明】

【０００１】発明の技術分野本発明は、負荷を受ける物体内の変形及び応力分布を決定における、特に、外
的又は内的負荷下のそれらの変形及び応力レベルを決定するための注入成形部分
の解析における専用の適用が決してないが、構成要素の構造解析のための方法及
び装置に関する。

【０００２】この発明は、薄壁構造、すなわち、１つの次元である厚さが少なくとも他の２
つの次元よりも数倍小さい最も重要な幾何学的特徴の構造解析に最も特別に適用
される。そのような構造の例は、金属、セラミック又は高分子材料の注入成形部
分、及びシートメタル（金属薄板）から形成される構造である。

【０００３】発明の背景技術特に注入成形による、構成要素及び部分の製造者は、設計コスト及び製品化ま
での時間を減らすために、設計プロセスにおける上流の部品の構造解析を動かす
ことを好む。設計技師の要求を満足するために、既存の商品は、有限要素解析（
ＦＥＡ）及びコンピュータ支援製図（ＣＡＤ）に統合される。Pro-Engineer（登
録商標）、CATIA（登録商標）、IDEAS（登録商標）、Solid Works（登録商標）
及びSolid Edge（登録商標）ブランドソリッドモデリングパッケージは、機械設
計及びラフティングに一般に用いられる。これらのパッケージは、構成要素幾何
学の３次元の写実的描写（「ソリッドモデル」として知られる）を生成するため
に用いられ得る。現在、ソリッドモデルに直接的に基づく構造解析パッケージは
、四面体及び六面体のような立体要素を用いる。

【０００４】ソリッドモデルの構造解析のために、ソリッドモデルによって画定された領域
は、ソリッド要素（ソリッドエレメント）と呼ばれる複数の小さい要素に分割さ
れる。このプロセスはメッシング（meshing）と呼ばれ、ソリッド要素の結果と
して生じる収集物はソリッドメッシュと呼ばれる。ソリッド要素は、通常、四面
体や六面体のような単純な幾何学的ソリッドである。ソリッドメッシュの生成は
、複雑な部分のためにめったに自動化されないけれども、近年改良されていた。
しばしば、ユーザは、メッシュがうまく生成されることを可能にするために、ソ
リッドモデルから機能を取り除く必要があろう。これは、非常に時間を浪費する
ことであり得、極端な場合には、構成要素の改造又は構成要素のある領域を必要
とするかもしれない。

【０００５】ソリッド要素の使用は、少なくとも負荷時の構成要素の構造応答を決定するこ
とにおいて、薄壁構造のためにシェル要素の使用に対して理論的利点を有しない
。しかしながら、大多数の構成要素モデリングは、ソリッドモデリングシステム
で行われ、ソリッド要素の使用は、より自然であり、幾何学的ソリッドモデルと
解析に用いられるメッシュとの間のより良いインターフェースを可能にする。特
定の問題が薄壁で囲まれた構成要素で発生する。この場合、正確な結果を達成す
るために、厚さ方向にいくつかの壁形状のソリッド要素があることを保証する必
要があった。これは、モデル内に多くの要素を導き、それ故、計算時間が長く、
大きなメモリが要求される。厚さを通して要素数を減らすために、より高いオー
ダの要素を使用できるが、そのようなメッシュの自動生成は未だに困難である。
大きいソリッド要素モデルのサイズを減らすために、ユーザは、特有の要素次元
を増やし、幾何学パターンを再メッシュしてもよい。自動メッシュ生成器は、わ
ずかな要素しか生成しないが、結果として生じる有限要素メッシュは、あまりに
も少ない要素のために、実際の応力分布をモデリングすることができないかもし
れない。さらに、厚さを通して不十分な要素を持つソリッド要素メッシュは、質
の悪い剛性マトリクス、単純曲げ（pure bending）のせん断ロッキング（shear
locking）及び不十分な擬態、曲げ支配の構造応答のような他の問題も有する。
これらは、有限要素解析の信頼性に深刻に影響を与え得る。

【０００６】薄壁構造は、典型的に、プレート及びシェル構成要素からなる。プレート及び
シェルのための古典理論がいくつか存在する。キルヒホッフ（Kirchhoff）理論
及びMindlin−Reissnerの理論が特に周知である。キルヒホッフ理論では、変形
前の中央表面（midsurface）への標準が垂線のままであるが、必ずしも変形後の
プレートに垂直ではない。中央表面への応力垂線は、どちらの理論でも無視され
る。多種のプレート及びシェル要素は、過去３５年にわたって異なるプレート及
びシェル理論に基づいて確立されていた。これらは、薄壁構造の正確な有限要素
解析を可能にするが、ＣＡＤシステムにおけるソリッド幾何学から得られなけれ
ばならないモデルを必要とする。解析のためのシェル要素モデルは、二次元の格
子又は曲線シェル要素からなる。一般に、要素の形状は、三角形又は四辺形のよ
うな単純な幾何学的形状である。それが要素の特性であるけれども、要素厚さは
、要素上に明確に示されない。シェル要素モデルは、ソリッドモデルの外壁間に
横たわる仮想断面（imaginary surface）上にシェル要素のメッシュを形成する
ことによって、ソリッド幾何学から生成され得る。この表面は、しばしばソリッ
ドモデルの中央平面（mid-plane）と呼ばれる。すべての場合に自動的に中央平
面表面を画定することはできず、シェル要素モデルの生成は、しばしば解析のた
めの別のモデルの構造を含む困難な仕事である。

【０００７】従って、薄壁構成要素の構造解析へのソリッド要素アプローチは、ソリッド幾
何学に容易に調和するという利点を有するが、シェル要素アプローチは、良い構
造的パフォーマンス、少ない計算時間、少ないメモリ要求、及びメッシュ生成の
容易さという利点を有する。しかしながら、低次数の要素の厚さを通して要素が
不十分ならば、ソリッド要素アプローチは、困難なメッシュ生成、高い要素数、
長い計算時間、高いメモリ要求、及び不十分な結果という不利を有するが、シェ
ル要素アプローチでは、シェルメッシュを生成する中央平面を得ることが難しい
。

【０００８】既存の境界要素方法は、ソリッド幾何学の表面で生成されるメッシュを用いる
ことによって、構成要素の構造解析を可能にするが、伝統的な境界要素方法は、
物質が等方性及び線形であることを要求する。また、境界要素方法は、方程式の
大きい縞のないシステムに至り、その解答は、多くのメモリ量を必要とする。

【０００９】上述のように、中央平面モデルが利用可能であるならば、シェル要素は、（一
般に薄壁）構造の構造解析に適している。確立したプレートシェル理論は、次数
が３から２に、すなわち、立体から平面に目立って減らされるように、シェル要
素に用いられる。一方、有限要素解析のためのＣＡＤパッケージからソリッドモ
デルを直接的に用いることが望ましい。

【００１０】そのようなシェルモデルは、一般に、三角形又は四辺形であり、平らであって
も曲がっていてもよい。各ノードには５又は６の自由度（ｄｏｆ）がある。最も
一般的な場合、自由度は、３つの並進移動と３つの回転からなる。図１は、それ
に付与される局所座標系を有する三角形のシェル要素を示し、その自由度は、そ
の座標系に参照される。局所ｘ、ｙ及びｚ方向におけるノードｎ（ｎ＝１，２又
は３）の並進移動自由度は、それぞれＵ_ｘｎ、Ｕ_ｙｎ及びＵ_ｚｎによって示され
る。同様に、局所ｘ、ｙ及びｚ軸についての回転は、それぞれθ_ｘｎ、θ_ｙｎ及
びθ_ｚｎによって示される。ノードが位置される要素を通す表面は、基準表面と
呼ばれる。通常、シェル要素は、基準表面として中央表面で定式化される。要素
基準表面が中央表面上にないならば、その要素は、基準表面が偏心率ε（図２を
見よ、そこでは、基準表面１０、中央表面１２及びノード１４が示される）と呼
ばれる中央表面から移動される距離を持つ偏心シェル要素であるといわれる。偏
心シェル要素の定式化は、標準のシェル要素の定式化を拡張することによって確
立され得る。中央表面上のひずみと基準表面の間、あるいは中央表面と基準表面
との間のノードの自由度のいずれかの関係は、標準のシェル要素の定式化を偏心
シェル要素の定式化への拡張のために用いられ得る。

【００１１】図３Ａ、３Ｂ及び３Ｃは、３つの可能な位置、すなわち、それぞれ要素の底面
、中央表面及び頂面に置かれる基準表面を持つ厚さｔの３つの平らな三角形のシ
ェル要素を示す。これらの図において、黒点はノード位置を示す。

【００１２】本発明の目的は、薄壁構造のためのシェル要素の利点の少なくともいくつかを
保持しながら、薄壁構造の構造解析のためにソリッド及びシェル要素を用いる制
限を対処することである。

【００１３】発明の概要それ故、本発明によれば、複数の表面部分を含む外部表面を有する物体の構造
応答を解析するための方法を提供する。その方法は、前記物体の３次元モデルを形成するステップであって、前記モデルは、前記外
部表面の表面メッシュ表現（見本、標本：representative）と複数の偏心シェル
要素を含み、前記各要素は、前記表面上の複数のノードによって画定され、前記
各ノードは、１以上の自由度を有する、形成ステップと、前記各要素における前記物体の厚さの半分を示す厚さを前記各要素に割り当て
るステップと、前記各要素のために、前記各要素のノードを含み、前記各要素における前記表
面メッシュと同時に起こる基準表面を画定するステップと、前記表面の対向する部分の一対のために、前記対向する部分の第１の部分にお
ける前記各ノードの前記自由度と、前記対向する部分の第２の部分における１以
上のノードの前記自由度との間の制約関係を確立するステップと、前記物体の有限要素構造解析を実行するステップと、前記構造解析の１以上の結果を出力するステップとを含む。

【００１４】従って、本発明によれば、シェルの基準表面は並置され、（ノード間の）多数
の制約が用いられ、厚さは、対応する表面部分の要素の収集が、ソリッドモデル
の中央平面に位置され、その位置でソリッドモデルの局所厚さを有するシェル要
素のメッシュとして、同一の構造的特徴を有することを保証するために割り当て
られる。有限要素構造解析を実行するステップは、境界条件と、外部及び／又は
内部負荷とを画定することを含む。従って、本発明は、ソリッド構成要素の構造
解析を実行するために表面メッシュ上に画定された修正シェル要素を用いる。

【００１５】解析の結果は、予想される変形と、いくらかの内部又は外部負荷条件下におけ
る物体の応力状態とを示すデータ又はイメージ（画像）を含んでもよい。

【００１６】表面部分は、平らでなくともよい。

【００１７】好ましくは、前記対向する部分の前記第１の部分の前記各ノードのために、前
記方法は、前記対向するノードが存在するとき、前記ノードの前記自由度と、前
記ノードの反対側に対向する部分対の第２の部分におけるノードである対向する
ノードの自由度との間の制約関係、あるいは、対向するノードが存在しないとき
、前記ノードの前記自由度と、前記ノードの反対側に対向する部分の要素である
対向する要素のノードの自由度との間の制約関係を確立するステップを含む。

【００１８】表面の部分は、シート又は端（エッジ）として、それに応じて、要素は、シー
ト要素又は端（エッジ）要素として分類されてもよい。

【００１９】好ましくは、前記物体の厚さの半分を示す厚さは、そのような厚さが画定でき
るならば、前記物体の実際の厚さから得られる。前記物体の厚さを画定できない
ところでは、前記要素の厚さは、隣接した要素の厚さ、あるいは隣接した要素の
厚さに比例しているとみなされてもよい。

【００２０】従って、物体の端における、あるいはその端部におけるような、重要な厚さを
画定することが不可能であるところでさえ、その厚さを示す厚さがなお画定され
得る。

【００２１】好ましくは、その方法は、その間で前記物体が最も薄い、前記物体の前記表面
の対向する部分である前記対向する部分対を選択するステップを含む。

【００２２】従って、その方法は、薄い物体のための最適な結果を提供することが期待され
る。その場合、制約関係は、好ましくは、一般にその間で前記物体が最も薄い表
面部分のノード間に画定される。

【００２３】メッシュは、三角形、四辺形又は（多角形を含む）他の単純形状の格子からな
り、いずれも平らかあるいは曲がっていてもよい。これらの形状の多くは、ＣＡ
Ｄシステムによって容易に生成され得る。

【００２４】構造解析は、その形成中物体に誘導される内部応力又は負荷の影響からなって
もよく、その場合、構造解析は、これらの内部負荷によってもたらされる物体の
結果として生じる変形又はそりからなるであろう。そのような変形は、以下「そ
り（warpage）」と称される。

【００２５】その方法は、物体のステレオリトグラフ表示を作り、あるいは組み込むことに
よって、表面メッシュを作り出してもよい。より好ましくは、その方法は、より
小さい要素を持つより細かいメッシュを生成することによるような）ステレオリ
トグラフ表示を改善するステップを含んでもよい。

【００２６】また、本発明は、複数の表面部分を含む外部表面を有する物体の構造を解析す
るための装置を提供する。その装置は、前記物体の３次元モデルを形成するためのモデリング手段であって、前記モデ
ルは、前記外部表面の表面メッシュ表現（見本、標本：representative）と複数
の偏心シェル要素を含み、前記各要素は、前記表面上の複数のノードによって画
定され、前記各ノードは、１以上の自由度を有する、モデリング手段と、前記各要素における前記物体の厚さの半分を示す厚さを前記各要素に割り当て
るための手段と、前記各要素のために、前記各要素のノードを含み、前記各要素における前記表
面メッシュと同時に起こる基準表面を画定するための手段と、前記表面の対向する部分の一対のために、前記対向する部分の第１の部分にお
ける前記各ノードの前記自由度と、前記対向する部分の第２の部分における１以
上のノードの前記自由度との間の制約関係を確立するための制約手段と、前記物体の有限要素構造解析を実行する解析手段と、前記構造解析の１以上の結果を出力する出力手段とを有する。

【００２７】これらの各手段は、コンピュータプログラムプロダクト又は部分を含んでもよ
く、その装置は、好ましくは、そのようなコンピュータプログラム部分を実行す
るためのコンピュータを含む。

【００２８】解析の結果は、予想される変形と、いくらかの内部又は外部負荷条件下におけ
る物体の応力状態とを示す（例えば、プリントアウト、ディスプレイ又はコンピ
ュータファイルの形式で）データ又はイメージを含んでもよい。

【００２９】好ましくは、前記対向する部分の前記第１の部分の前記各ノードのために、前
記制約手段は、前記対向するノードが存在するとき、前記ノードの前記自由度と
、前記ノードの反対側に対向する部分対の第２の部分におけるノードである対向
するノードの自由度との間の制約関係、あるいは、対向するノードが存在しない
とき、前記ノードの前記自由度と、前記ノードの反対側に対向する部分の要素で
ある対向する要素のノードの自由度との間の制約関係を確立するために、実施可
能である。

【００３０】好ましくは、ある厚さ（すなわち、各要素における物体の厚さの半分を示す厚
さ）を前記各要素に割り当てるための前記手段は、前記各要素における前記物体
の全厚さを示す厚さ、そのような厚さを画定できるならば、前記物体の実際の全
厚さを用い、前記物体の厚さを画定できない場合には、隣接する要素の厚さに等
しい、あるいは隣接する要素の厚さに比例する厚さを前記要素に割り当てるため
に、実施可能である。

【００３１】従って、物体の厚さの半分を示す厚さを計算することにおいて、その装置は、
好ましくは、可能であれば、実際の全厚さから半分の厚さを得る。

【００３２】好ましくは、その装置は、その間で前記物体が最も薄い、前記物体の前記表面
の対向する部分を前記対向する部分対として選択するように構成される。

【００３３】メッシュは、平らな又は曲がっている、三角形、四辺形あるいは他の単純形状
の格子からなってもよい。これらの形状の多くは、ＣＡＤシステムによって容易
に生成され得る。

【００３４】解析手段は、形成中物体に誘導される内部応力又は負荷の影響を解析すること
が実施可能であるか、あるいはそのように構成されてもよい。その場合、構造解
析は、これらの内部負荷によってもたらされる物体の結果として生じる変形又は
そりからなるであろう。

【００３５】好ましい実施の形態の詳細な記述本発明によれば、３次元の本体における構造解析を実行するために３次元の物
体の境界において画定されるシェル要素のメッシュを用いる方法は、以下に記述
される。

【００３６】多くの概念を導入するために、一例として、図４厚さｔの薄い方形棒（バー）
が描写される。バーの構造応答は、本体の外部上にメッシュとして画定されるシ
ェルモデルによって表される。図４において、方形バーを画定するメッシュは、
次のように示される頂点にノードを置くことによって画定される四辺形のシェル
要素のセットになるように画定される：

【表１】

【００３７】明瞭のために、端部（要素３、４、５及び６）は、今のところ無視されるが、
後ほど論じられる。（バーの頂面と底面を表す）番号１及び２のシェル要素がバ
ーの構造応答を画定するためにどのように用いられ得るかがまず論じられる。

【００３８】余談として、それぞれノードａｄｇｈ及びｂｃｆｅによって画定される基準表
面で、厚さｔの従来のシェル要素が要素１及び２のために用いられ、その要素が
層状のシェルモデルを形成するために互いに結合されたとみなされるならば、結
果として生じるモデルは、シェルモデルの領域がバーの外方に延びるような（図
５を見よ、そこでは、符号１６が要素１の基準表面を示し、符号１８が要素２の
基準表面を示す）、より不自然なモデルの応答を提供するだろう。実際に、厚さ
ｔのバーのために、このシェル要素モデルは、厚さ２ｔのバーを表す。正確な構
造応答を得るために、本発明によれば、０．５ｔの厚さは、要素１及び２に帰さ
れ、図６に示されるように、要素１のための基準表面２０は、（ａｄｇｈのノー
ドによって画定される）バーの頂面に移動され、要素２のための基準表面２２は
、（ノードｂｃｆｅによって画定される）バーの底面に移動される。

【００３９】上述では、端部（エッジ）の効果が無視された、これは以下で扱われる。要素
は互いに結合されたとみなされる。実際には、これは、本発明によれば、他点制
約を用いることによって達成され、詳細は以下で論じられる。

【００４０】さらに、平らなプレートのソリッドモデルの単純な場合が図７に示される。そ
れは、頂部表面と底部表面を有する二重に覆われて現れる。幾何学的見地から、
平らなプレートは、それぞれ半分の厚さである２つのプレートの結合とみなされ
る。頂部プレートは、基準表面として頂部表面を持つ偏心シェル要素を用いてモ
デリングされ、そこでは、プレート仮定が用いられる。同様に、底部プレートは
、基準表面として底部表面を持つ偏心シェル要素を用いてモデリングされる。頂
部表面と底部表面の結合は、頂部表面の変形と底部表面の変形の間の位置関係か
ら生じるとみなされる。この関係は、古典プレート理論の仮定を用いて確立され
る。

【００４１】本発明における構造解析の好ましい方法は、次のステップを含み、それぞれは
詳細に論じられる：１．外部境界上のメッシュの創造２．厚さと偏心の割り当て：適合及び不適合要素３．ノードの表面へのマッチング４．ノードｎの自由度とそのマッチング要素ｅの自由度との間の制約方程式の確
率５．内部及び外部負荷、並びに境界条件の指定６．解析実行及び設計改善

【００４２】・外部境界上のメッシュの創造ソリッドモデルから始め、中央平面表面よりむしろ計算領域を作り出すために
、３次元の物体を画定する外部表面を利用せよ。三角形又は四辺形の要素を用い
て外部表面をメッシュせよ。計算領域は、一部の一般に対向する第１及び第２の
表面のメッシュ表示からなる。例えば、Ｔ形状部のメッシュ又は計算領域が図８
に示される。

【００４３】・厚さと偏心の割り当て：適合及び不適合要素第１の表面を第１の半厚部の基準表面と、第１の表面を第２の半厚部の基準表
面とみなしなさい。２つの表面の要素は、そのような要素間の厚さを識別する能
力に基づいて適合される。これらのシェル要素の厚さを適合要素間の距離の半分
に割り当てなさい。また、第１及び第２の表面上の要素の偏心が決定され、それ
は要素厚の半分である。要素の偏心の（正又は負の）記号は、選択されたシェル
要素定式化の予め定められたルールに依存する。

【００４４】好ましくは、適合できない要素である、前記第１及び第２の表面のあらゆる不
適合要素は、そのような隣接する適合要素が存在するところに近接する適合要素
、あるいは、そのような隣接する適合要素が存在せず、前記隣接する不適合要素
が厚さに割り当てられたところに隣接する不適合要素の平均厚さである厚さに割
り当てられる。これらの要素の偏心を適合要素として割り当てなさい。従って、
結局、第１及び第２の表面のすべての適合及び不適合要素が厚さ及び偏心を割り
当てられ得る。

【００４５】前記第１及び第２の表面間の表面であり、前記第１又は第２の表面のいずれか
に隣接するエッジ表面の各要素は、エッジ表面の要素が隣接する第１及び第２の
表面の要素の厚さの０．０１〜１００．０（好ましくは、１／３）倍の厚さに割
り当てられる。

【００４６】好ましくは、前記第１又は第２の表面に隣接しないエッジ表面の各要素は、厚
さを割り当てられた前記エッジ表面の要素に隣接する平均厚さである厚さに割り
当てられる。

【００４７】薄壁部分のエッジ変形特性をシミュレートするために、エッジ表面の各要素は
、人工的異方性物質特性を割り当てられる。主な物質軸方向は、部分厚方向に沿
う。部分厚方向に平行なヤング率（Ｅ１）は、前記エッジ表面の前記要素が隣接
する前記第１及び第２の表面の物質の平均ヤング率の０．００１〜１０００（好
ましくは、１５）倍に割り当てられ得る。部分厚方向に垂直なヤング率（Ｅ２）
は、前記エッジ表面の前記要素が隣接する前記第１及び第２の表面の物質の平均
ヤング率の０．００１〜１００（好ましくは、０．１）倍に割り当てられ得る。

【００４８】好ましくは、前記第１又は第２の表面に隣接しないエッジ表面の各要素は、そ
れぞれ材料特性を割り当てられた、部分厚方向に平行な平均ヤング率と前記エッ
ジ表面の隣接する要素の部分厚方向に垂直な平均ヤング率とである、部分厚方向
（Ｅ１）に平行及び垂直なヤング率を割り当てられる。

【００４９】エッジ表面の各要素のための偏心率は０である。

【００５０】・表面上のポイントへのノードのマッチング図９において、ノードｎとポイントＰとの関連線が、ノードｎにおける前記第
１の表面２４の垂線、ポイントＰにおける前記第２の表面２６の垂線、又はこれ
ら２つの垂線のある種の平均ベクトルに平行であるというルールに基づいて、前
記第１の表面（示されるように、底部表面２４）上のノードｎと一般に対向する
第２の表面（示されるように、頂部表面２６）上のポイントＰとの間の適合関係
を確立せよ。ポイントＰは、一般に、特定の要素内にあるが、ノードであり得、
前記第２の表面上の要素の辺上にあり得る。ポイントＰは、そのノードｎの適合
ポイントとみなされる。ポイントＰが位置される要素は、同様に容易に決定され
、ノードの適合要素とみなされる。ポイントＰがいくつかの要素で共有されるな
らば、これらの要素の一つは、適合要素として選ばれる。この実施の形態では、
最も小さい領域の要素が適合要素として選ばれる。

【００５１】・ノードｎの自由度とそのマッチング要素ｅの自由度との間の制約方程式の確
率キルヒホッフ及びMindlin−Reissnerプレート理論の両方で、変形前の中央表
面への垂線が直線のままであるという仮定が適用できる。上述のように、前記第
１の表面の要素は第１の半厚部を表し、前記第２の表面の要素は第２の半厚部を
表す。第１の半厚部と第２の半厚部が一つとして負荷に耐え得るように、それら
を互いに結合させるために、上記仮定は、前記第１の表面メッシュの変形と第２
の表面メッシュの変形との間で実行されなければならない。数値的に、前記第１
の表面上のノードの自由度は、前記第２の表面上のノードの自由度に関連すべき
である。

【００５２】図１０において、ノードｎが対向する表面（要素ｅ）上に適合ポイントＰを有
するので、それらの連結線は、中央表面の垂線に、正確にあるいは非常に厳密に
平行である。変形前の中央表面への垂線が直線のままであるという仮定に基づい
て、ノードｎのすべての又はいくつかの自由度（転移配置及び回転）は、適合要
素の局所座標系の適合ポイントＰにおける配置変数（転移配置及び回転）の線形
関数として表現され得る。

【００５３】適合要素の適合ポイントの位置及び選択されたシェル要素の配置補間関数に基
づいて、適合ポイントにおける配置変数（転移配置及び回転）は、適合要素の局
所座標系における適合要素のノードの自由度の線形関数として表現され得る（図
１０）。

【００５４】それ故、上記２つの関係によれば、ノードのすべて又は一部の自由度は、適合
要素のノードの自由度の線形関数として表現され得る。これらの関係は、全体的
な有限要素行列式に組み込まれ得る典型的多点制約（ＭＰＣ）方程式である。Ｍ
ＰＣ方程式を扱うためのアルゴリズムは、よく確立され、有限要素解析における
種々の出版されたテキストに利用可能である。変換方法、ラグランジュ乗数及び
多点制約を実行するペナルティ関数方法は、現在の問題にすべて適用できる。

【００５５】三角形の切子面シェル要素がこの実施の形態で用いられる。これは、１８の自
由度（各ノードにおいて６）を持つ３ノード要素である。その要素は、Berganと
Felippa（50 (1985) 25の応用力学及び工学における計算方法）の局所薄膜定式
化をBatozとLardeur（28 (1989) 533の工学の数滴方法の国際ジャーナル）の曲
げ定式化に重ね合わせ、結合された方程式を大域座標系に変換することによって
構築される。

【００５６】局所基準表面垂線についてのドリリング回転自由度は、薄膜定式化において用
いられる。それは、以下によって局所要素系に定義される：

【式１】

【００５７】次に、ノードｎの自由度とその適合要素ｅの自由度との関係が定義される（図
１０を見よ）。キルヒホッフとMindlinの両プレート理論に用いられる、変形前
の中央表面への垂線が変形後も直線のままであるという仮定に基づいて、ノード
ｎの自由度と適合ポイントＰの配置及び回転との間に次の関係がある。

【００５８】

【式２】ここで、「ｈ」は、ノードｎとその適合ポイントＰとの間の距離であり、Ｕ_ｘｎ、Ｕ_ｙｎ、Ｕ_ｚｎ、θ_ｘｎ、θ_ｙｎ及びθ_ｚｎは、その適合要素ｅの局所要素座
標系に参照されるノードｎにおける局所自由度であり、Ｕ_ｘｐ、Ｕ_ｙｐ、Ｕ_ｚｐ、θ_ｘｐ、θ_ｙｐ及びθ_ｚｐは、ポイントＰにおけるとともに、適合要素ｅの局
所要素座標系における配置及び回転である。

【００５９】 θ_ｚｎを定義する式（２）の最後の関係が式（１）を用いて得られることに留
意されたい。

【００６０】式（２）は、行列形式に書き直され得る：

【式３】

【００６１】要素定式化から、ポイントＰにおける配置及び回転は、形状関数補完法によっ
て得られる。すなわち、

【式４】

【００６２】

【式５】

【００６３】式（３）、（４）及び（５）を結合して、ノードｎの自由度とその適合要素ｅ
の自由度との間の関係は、以下のように得られる：

【式６】

【００６４】これらの制約方程式は、大域座標系に変換され得、それによって、大域座標系
におけるノードｎの自由度とその適合要素ｅの自由度との間の最終的な制約方程
式を提供する。

【００６５】・内部及び外部負荷、並びに境界条件の指定メッシュのための負荷、制約条件及び材料特性が指定される。負荷条件は、あ
らゆる種類の通常の外部及び／又は内部の負荷を含む。作用する負荷は、ポイン
ト、エッジ、パッチ又は容量で定義され得る。負荷タイプは、静的又は動的な、
実際の力又は熱負荷、残留応力／ひずみ及び初期応力／ひずみであり得る。材料
は、等方性又は異方性であり得る。

【００６６】・解析実行及び設計改善物体の構造の有限要素解析は、配置、ひずみ及び応力を得るために行われる。

【００６７】解析の結果は、構成要素の応力、ひずみ、偏差及び変形形状を含む。それ故、
本発明は、例えば、以下のようなものを含む、あらゆる構造解析の通常の結果を
決定するために用いられる：１．材料は、あまりにも高い応力又はひずみレベルを受ける、２．構成要素は、その意図的な設計機能を実行するために、十分に厳格である、
３．偏向部分は、アセンブリの他の部分に強い影響を与える。

【００６８】解析結果を考慮することによって、ユーザは、負荷条件下の部分のパフォーマ
ンスを改善するために、部分幾何学あるいは材料タイプを修正することができる
。

【００６９】この好ましい方法は、注入成形部分のそり解析のために用いられてもよい。ま
ず、内部負荷の一例として、成形中の残留応力分布は、Moldflow Corporationに
よって提供されるような、フロー及び冷却解析から決定される。用語「成形中の
残留応力分布（in-mold residual stress distribution）」は、部分が成形中で
ちょうど放出前の処理のために、部分の応力分布に言及する。この分布は、構成
要素の偏向形状及び結果として生じる残留応力分布を作るためにここで記述され
る装置への入力として用いられてもよい。本発明の好ましい方法のユーザは、部
分幾何学、処理条件あるいは材料タイプを修正することができ、設計が改善され
たことをチェックするために、さらなる解析を実行することができる。

【００７０】その代わりに、残留応力、残留ひずみ又は熱負荷は、構造解析で用いられると
き、変形形状が計算され得、それによって、処理の効果、部分幾何学及び材料タ
イプを明らかにするような、他の市販のソフトウェアから得られてもよい。

【００７１】注入成形部分のそり解析のために、厚さ変更は、以下のように、第１の表面の
ノードの局所偏向（Ｕ_ｚｎ）と、対向する表面の適合ポイントの局所偏向（Ｕ_ｚ _ｐ）との間の特別な制約方程式（２）を修正することによって強制されるべきで
ある：

【式７】ここで、ψは、厚さ収縮量を表し、その値は、部分の成形中の残留応力レベルに
基づく。この残留応力レベルが平面内の残留応力であることに留意されたい。

【００７２】この好ましい実施の形態では、主な初期ひずみは、エッジ表面上の要素を除く
すべての要素のために平面内の成形中の残留応力から計算され、ψは、これらの
主な初期ひずみを単純に平均することによって得られる。それ故、厚さ収縮量ψ
は、ちょうど部分の放出前の平面内応力分布から結果として生じる平均平面内応
力である。

【００７３】フロー及び冷却解析からエッジ表面の要素についての成形中の残留応力結果は
、捨てられる。特別な人工の成形中の残留応力は、部分変形のエッジ要素の制約
のための非現実的なそりを防止するために、これらのエッジ要素に割り当てられ
る。

【００７４】第１又は第２の表面の一つ上の要素ｅ_２と辺ａｂを共有する、図１１のｅ_１の
ようなエッジ要素のために、それらの初期の応力は、以下のように計算される。
第１に、要素ｅ_２の成形中の残留応力に基づいて、辺ａｂに沿うひずみが計算さ
れ、第２に、このひずみが部分厚方向に垂直な主ひずみとして要素ｅ_１にとられ
、第３に、厚さ収縮量ψが部分厚方向に平行な主ひずみとして要素ｅ_１にとられ
る。最終的に、これら２つの主ひずみから、同等の成形中の残留応力をエッジ要
素ｅ_１のために得ることができる。

【００７５】第１又は第２の表面上のあらゆる要素とどの辺も共有しない、図１１のｅ_３の
ようなエッジ要素のために、それらの成形中の残留応力は、人工的成形中残留応
力を割り当てられた前記エッジ表面の隣接するエッジ要素の平均成形中残留応力
に割り当てられる。

【００７６】結論従って、本発明の本質的概念は、構造解析を実行するためのソリッドモデル上
の表面メッシュの使用である。理想的に、本発明は、薄壁構造に適する。また、
薄壁構造に小さい塊になった領域があり、応力／ひずみ状態がこれらの塊になっ
た領域には重要でないならば、それは特に有用である。その解析は、線形、非線
形、バックリング及び形態頻度解析を含む、静的及び動的構造解析のために用い
られ得る。多くの異なるタイプの負荷が用いられてもよい。解析の結果は、応力
、ひずみ、部分の変形及び大域座標系における偏向を含む。

【００７７】また、この発明は、注入成形部分のそり解析のために用いられ得る。そり解析
に用いられるとき、その結果は、構造解析の場合に関するが（応力、ひずみ、部
分の変形及び大域座標系における偏向）、変形形状は、注入成形プロセス中発揮
される残留応力のためである。

【００７８】境界要素方法にあるような、負荷及び材料タイプの特別な制限はない。ソリッ
ドモデルは、ＳＴＬ、ＩＧＥＳ、ＡＣＩＳ、ＰＡＲＡＳＯＬＩＤなどのようない
くつかのファイル形式で入力され得る。実際に、本方法の使用を可能にするため
には、表面メッシュが部分モデルにおいて作られることのみが必要である。

【００７９】上記を考慮して、本発明の方法及び装置は、材料特性が処理によって影響を及
ぼされる注入成形構成要素の解析に特に適している。この点で、この出願に開示
される方法は、材料特性が異方性である注入成形構成要素における構造解析（線
形又は非線形）を実行するために用いられてもよい。そのような解析は、注入成
形材料の機械的特性を予測できる、Moldflow Corporationによって提供されるよ
うな、他形式の解析と調和することができる。機械的特性の偏向は、分子配向、
結晶形態又は短い繊維強化材の追加から生じてもよい。

【００８０】また、本発明に開示される方法は、重合、金属充填、セラミック充填又は金属
材料から作られる注入成形構成要素のそり形状を決定するために用いられてもよ
い。処理効果のために、注入成形材料は、異なる方向に異なる量だけ縮み得る。
本発明は、注入成形部分の残留応力分布を予測できる他の解析製品と調和され得
る。この残留応力分布は、処理のために材料の熱加工特性の変化を考慮すること
によって決定される。この応力分布を構造解析の負荷として用いることによって
、そり形状及び最終的な構成要素の収縮量が決定され得る。そのような解析が過
度のそりレベルを示すならば、構成要素幾何学、材料タイプ又は処理条件が変更
され得る。もう一つの解析は、改良が注目されるまで実行され得る。

【００８１】本発明の精神及び範囲内の修正は、当業者にとって容易に果たされ得る。それ
故、この発明が上記に一例として記述される特定の実施の形態に制限されないこ
とを理解されたい。

【図面の簡単な説明】

好ましい実施の形態は、添付図面を参照して、一例として記述される。

【図１】定義された局所座標系を持つ先行技術の典型的な三角形シェル要素の描写であ
る。

【図２】要素基準表面が中央表面にない、先行技術の偏心シェル要素を描写する。

【図３】３つの可能な位置の一つに置かれる基準表面を持つ厚さｔの先行技術の平らな
三角形シェル要素をそれぞれ示す。

【図４】モデリングされた物体の一例として薄い長方形の棒の図である。

【図５】先行技術、従来のシェル要素でモデリングされた図４の棒の図である。

【図６】本発明の好ましい実施の形態における方法において、シェル要素でモデリング
された図４の棒の図である。

【図７】本発明によって、それぞれ厚さ０．５ｈの２つのプレートの接合としてどのよ
うに見せられ得るかを示す、厚さｈの単純なプレートのような物体のソリッドモ
デルの図である。

【図８】部分の第１及び第２の一般に対向する表面のメッシュ表示を含むＴ形状部分の
モデルの図である。

【図９】本発明における物体の対向する表面の要素に物体の一表面の要素のノードのマ
ッチングの図である。

【図１０】マッチング要素の局所座標系の図である。

【図１１】本発明における物体のモデルの端要素の図である。

───────────────────────────────────────────────────── フロントページの続き (81)指定国ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＣＹ，ＤＥ，ＤＫ，ＥＳ，ＦＩ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＯＡ(ＢＦ，ＢＪ，ＣＦ，ＣＧ，ＣＩ，ＣＭ，ＧＡ，ＧＮ，ＧＷ，ＭＬ，ＭＲ，ＮＥ，ＳＮ，ＴＤ，ＴＧ)，ＡＰ(ＧＨ，ＧＭ，ＫＥ，ＬＳ，ＭＷ，ＭＺ，ＳＤ，ＳＬ，ＳＺ，ＴＺ，ＵＧ，ＺＷ)，ＥＡ(ＡＭ，ＡＺ，ＢＹ，ＫＧ，ＫＺ，ＭＤ，ＲＵ，ＴＪ，ＴＭ)，ＡＥ，ＡＧ，ＡＬ，ＡＭ，ＡＴ，ＡＵ，ＡＺ，ＢＡ，ＢＢ，ＢＧ，ＢＲ，ＢＹ，ＢＺ，ＣＡ，ＣＨ，ＣＮ，ＣＲ，ＣＵ，ＣＺ，ＤＥ，ＤＫ，ＤＭ，ＤＺ，ＥＥ，ＥＳ，ＦＩ，ＧＢ，ＧＤ，ＧＥ，ＧＨ，ＧＭ，ＨＲ，ＨＵ，ＩＤ，ＩＬ，ＩＮ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＣ，ＬＫ，ＬＲ，ＬＳ，ＬＴ，ＬＵ，ＬＶ，ＭＡ，ＭＤ，ＭＧ，ＭＫ，ＭＮ，ＭＷ，ＭＸ，ＭＺ，ＮＯ，ＮＺ，ＰＬ，ＰＴ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＳＬ，ＴＪ，ＴＭ，ＴＲ，ＴＴ，ＴＺ，ＵＡ，ＵＧ，ＵＺ，ＶＮ，ＹＵ，ＺＡ，ＺＷ (72)発明者ツェングロングオーストラリア 3131 ビクトリア州フォレストヒル、スティーブンスロード 16 (72)発明者ケネディー，ピーター，ケネスアメリカ合衆国 14850 ニューヨーク州イタカザパークウェイ 315 (72)発明者ユ，フア，ガンオーストラリア 3105 ビクトリア州ブーリーン、ヤラバレーブバード 113 Ｆターム(参考） 5B046 JA08

Claims

【特許請求の範囲】

【請求項１】複数の表面部分を含む外部表面を有する物体の構造応答を解析
するための方法であって、前記物体の３次元モデルを形成するステップであって、前記モデルは、前記外
部表面の表面メッシュ表現と複数の偏心シェル要素を含み、前記各要素は、前記
表面上の複数のノードによって画定され、前記各ノードは、１以上の自由度を有
する、形成ステップと、前記各要素における前記物体の厚さの半分を示す厚さを前記各要素に割り当て
るステップと、前記各要素のために、前記各要素のノードを含み、前記各要素における前記表
面メッシュと同時に起こる基準表面を画定するステップと、前記表面の対向する部分の一対のために、前記対向する部分の第１の部分にお
ける前記各ノードの前記自由度と、前記対向する部分の第２の部分における１以
上のノードの前記自由度との間の制約関係を確立するステップと、前記物体の有限要素構造解析を実行するステップと、前記構造解析の１以上の結果を出力するステップと、を含むことを特徴とする方法。
【請求項２】前記表面部分のいくらかあるいはすべてが平らでないことを特
徴とする請求項１に記載の方法。
【請求項３】前記対向する部分の前記第１の部分の前記各ノードのために、
前記方法は、前記対向するノードが存在するとき、前記ノードの前記自由度と、
前記ノードの反対側に対向する部分対の第２の部分におけるノードである対向す
るノードの自由度との間の制約関係、あるいは、対向するノードが存在しないと
き、前記ノードの前記自由度と、前記ノードの反対側に対向する部分の要素であ
る対向する要素のノードの自由度との間の制約関係を確立するステップを含むこ
とを特徴とする請求項１又は２に記載の方法。
【請求項４】前記物体の厚さの半分を示す厚さは、そのような厚さが画定で
きるならば、前記物体の実際の厚さから得られ、前記物体の厚さを画定できない
ところでは、前記要素の厚さは、隣接した要素の厚さ、あるいは隣接した要素の
厚さに比例しているとみなされることを特徴とする請求項１乃至３のいずれかに
記載の方法。
【請求項５】その間で前記物体が最も薄い、前記物体の前記表面の対向する
部分である前記対向する部分対を選択するステップを含むことを特徴とする請求
項１乃至４のいずれかに記載の方法。
【請求項６】前記メッシュは、平らであるか曲がっている、三角形、四辺形
又は（多角形を含む）他の単純形状の格子からなることを特徴とする請求項１乃
至５のいずれかに記載の方法。
【請求項７】その形成中、前記物体に誘導される内部応力又は負荷の影響を
解析するための請求項１乃至６のいずれかに記載の方法。
【請求項８】前記内部応力又は負荷によってもたらされる前記物体の変形又
はそりを解析するための請求項７に記載の方法。
【請求項９】前記方法は、前記物体のステレオリトグラフ表示を作り、ある
いは組み込むことによって、前記表面メッシュを作り出すステップを含むことを
特徴とする請求項１乃至８のいずれかに記載の方法。
【請求項１０】前記物体のステレオリトグラフ表示を作り、あるいは組み込
むこと、及び前記ステレオリトグラフ表示を改善することによって、前記表面メ
ッシュを作り出すステップを含むことを特徴とする請求項１乃至９のいずれかに
記載の方法。
【請求項１１】前記方法は、物体の構造応答を解析するために、前記方法を
符号化するコンピュータプログラムを処理するコンピュータによって実行される
ことを特徴とする請求項１乃至１０のいずれかに記載の方法。
【請求項１２】請求項１乃至１１のいずれかに記載の方法に従って、物体の
構造応答を解析することを含むことを特徴とする物体を形成する方法。
【請求項１３】請求項１２に記載の方法に従って形成される物体又は物体の
部分。
【請求項１４】複数の表面部分を含む外部表面を有する物体の構造を解析す
るための装置であって、前記物体の３次元モデルを形成するためのモデリング手段であって、前記モデ
ルは、前記外部表面の表面メッシュ表現と複数の偏心シェル要素を含み、前記各
要素は、前記表面上の複数のノードによって画定され、前記各ノードは、１以上
の自由度を有する、モデリング手段と、前記各要素における前記物体の厚さの半分を示す厚さを前記各要素に割り当て
るための手段と、前記各要素のために、前記各要素のノードを含み、前記各要素における前記表
面メッシュと同時に起こる基準表面を画定するための手段と、前記表面の対向する部分の一対のために、前記対向する部分の第１の部分にお
ける前記各ノードの前記自由度と、前記対向する部分の第２の部分における１以
上のノードの前記自由度との間の制約関係を確立するための制約手段と、前記物体の有限要素構造解析を実行する解析手段と、前記構造解析の１以上の結果を出力する出力手段と、を有することを特徴とする装置。
【請求項１５】前記各手段は、コンピュータプログラムプロダクト又は部分
からなることを特徴とする請求項１４記載の装置。
【請求項１６】そのようなコンピュータプログラム部分を実行するためのコ
ンピュータを含むことを特徴とする請求項１５記載の装置。
【請求項１７】前記解析の前記結果は、予想される変形と、いくらかの内部
又は外部負荷条件下における物体の応力状態とを示すデータ又はイメージを含む
ことを特徴とする請求項１４乃至１６のいずれかに記載の装置。
【請求項１８】前記対向する部分の前記第１の部分の前記各ノードのために
、前記制約手段は、前記対向するノードが存在するとき、前記ノードの前記自由
度と、前記ノードの反対側に対向する部分対の第２の部分におけるノードである
対向するノードの自由度との間の制約関係、あるいは、対向するノードが存在し
ないとき、前記ノードの前記自由度と、前記ノードの反対側に対向する部分の要
素である対向する要素のノードの自由度との間の制約関係を確立するために、実
施可能であることを特徴とする請求項１４乃至１７のいずれかに記載の装置。
【請求項１９】ある厚さを前記各要素に割り当てるための前記手段は、前記
各要素における前記物体の全厚さを示す厚さ、そのような厚さを画定できるなら
ば、前記物体の実際の全厚さを用い、前記物体の厚さを画定できない場合には、
隣接する要素の厚さに等しい、あるいは隣接する要素の厚さに比例する厚さを前
記要素に割り当てるために、実施可能であることを特徴とする請求項１４乃至１
８のいずれかに記載の装置。
【請求項２０】前記装置は、その間で前記物体が最も薄い、前記物体の前記
表面の対向する部分を前記対向する部分対として選択するように構成されること
を特徴とする請求項１４乃至１９のいずれかに記載の装置。
【請求項２１】前記メッシュは、平らであるか曲がっている、三角形、四辺
形あるいは他の単純形状の格子からなることを特徴とする請求項１４乃至２０の
いずれかに記載の装置。
【請求項２２】前記解析手段は、形成中物体に誘導される内部応力又は負荷
の影響を解析することが実施可能であるか、あるいはそのように構成されること
を特徴とする請求項１４乃至２０のいずれかに記載の装置。
【請求項２３】請求項１乃至１１のいずれかにおける、物体の構造応答を解
析する方法のために、該方法を符号化するコンピュータプログラムで提供され、
あるいはそれを処理するコンピュータ。
【請求項２４】請求項１乃至１１のいずれかにおける、物体の構造応答を解
析する方法のために、該方法を具現化するコンピュータプログラムで提供される
コンピュータに読み取り可能な記憶媒体。