JP2011159213A - 溶接変形解析方法 - Google Patents
溶接変形解析方法 Download PDFInfo
- Publication number
- JP2011159213A JP2011159213A JP2010022232A JP2010022232A JP2011159213A JP 2011159213 A JP2011159213 A JP 2011159213A JP 2010022232 A JP2010022232 A JP 2010022232A JP 2010022232 A JP2010022232 A JP 2010022232A JP 2011159213 A JP2011159213 A JP 2011159213A
- Authority
- JP
- Japan
- Prior art keywords
- model
- analysis
- local
- local model
- plastic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
【課題】大型溶接構造物の溶接変形予測における計算精度の向上および計算時間の短縮を両立する。
【解決手段】メッシュを生成して解析対象となる溶接構造物をモデル化し熱弾塑性解析を行う方法において、溶接構造物の全体モデルを構築し(S1)、その全体モデルから、溶接部を含む局部モデルを抽出する(S2)。次に、抽出した局部モデルの、全体モデルの残り部分との境界部を拘束して(S3)熱弾塑性解析を行い(S4)、熱弾塑性解析の解析結果を含む局部モデルを全体モデルの残り部分に貼り付けて全体モデルを再構築する(S5)。その後、境界部の拘束を解放し(S6)、全体モデルの弾性解析を行うことで(S7)、溶接構造物の変形を計算する。
【選択図】図2
【解決手段】メッシュを生成して解析対象となる溶接構造物をモデル化し熱弾塑性解析を行う方法において、溶接構造物の全体モデルを構築し(S1)、その全体モデルから、溶接部を含む局部モデルを抽出する(S2)。次に、抽出した局部モデルの、全体モデルの残り部分との境界部を拘束して(S3)熱弾塑性解析を行い(S4)、熱弾塑性解析の解析結果を含む局部モデルを全体モデルの残り部分に貼り付けて全体モデルを再構築する(S5)。その後、境界部の拘束を解放し(S6)、全体モデルの弾性解析を行うことで(S7)、溶接構造物の変形を計算する。
【選択図】図2
Description
本発明は、数値解析による溶接構造物(単に、「構造物」と称する場合がある)の変形、残留応力などを計算する技術に関する。
溶接施工により大型構造物を製造する際に、溶接部近傍への熱累積とその後の冷却により溶接変形が発生する。その溶接変形を低減するために、通常では拘束冶具の取り付けや、溶接後の矯正作業などを行っている。このような状況において、有限要素法などの数値解析により変形を予測し、変形対策の適正化を図ることは、生産効率の向上およびコストの低減のために極めて重要である。
有限要素法による溶接変形の解析には、大きく分けると熱弾塑性解析と固有ひずみ法の二つの方法がある。熱弾塑性解析による溶接変形の解析は、解析対象とする溶接構造物において、溶接中の熱履歴を非定常熱伝導解析から求め、次に非線形解析である熱弾塑性解析により溶接中の変位、ひずみおよび応力の履歴を解析する方法である。
一方、固有ひずみ法による溶接変形の解析は、溶接部およびその近傍に生じた固有ひずみを溶接構造物に与え、線形解析である(熱)弾性解析により溶接変形を計算する方法である。溶接により生じる固有ひずみは、見かけのひずみから弾性ひずみを引いた値であり、固体力学の範囲では残留塑性ひずみと等しい。
この固有ひずみ法を用いて溶接構造物の溶接変形や残留応力を推定するには以下の手順をとる。まず、実際の溶接により発生した固有の溶接変形や残留応力を計測し、或いは、溶接構造物に用いられる溶接法、継手形状、溶接条件などをパターン化した単元モデルを用い、熱弾塑性解析により溶接変形や残留応力を求め、これらをひずみの適正値としてデータベース化しておく。次に、データベース化したひずみの適正値を溶接構造物のモデルに与え、弾性解析法を用いて溶接変形や残留応力を計算する。
特許文献1には、3次元モデルを用いた熱弾性解析を行うことにより得られた熱変形の結果を変形拘束条件として、2次元モデルを用いた熱弾塑性解析を行うことにより、解析精度が高く、且つ計算時間の削減可能な溶接残留応力解析方法が開示されている。
特許文献2には、溶接構造物全体のモデルを線形解析の固有ひずみ法で解析した個々の固有ひずみ成分を局所座標系に変換して溶接変形、残留応力を算出する解析方法が開示されている。
特許文献3には、溶接構造物の残留応力有限要素法で解析する際、溶接線から解析モデル境界までの軸方向距離を制約することにより解析時間を短縮する解析方法が開示されている。
特許文献4には、拘束条件で得られた溶接固有ひずみを用いて、溶接後拘束条件をなくしたときの溶接変形を求める解析方法が開示されている。
特許文献5には、溶接近傍のみ非線形解析を行い、被溶接物につき線形解析を行うことにより両領域の限界面における反力を計算し、その反力の差が所定の範囲内なるように収束演算することで、溶接変形を短時間に算出する解析方法が開示されている。
特許文献2には、溶接構造物全体のモデルを線形解析の固有ひずみ法で解析した個々の固有ひずみ成分を局所座標系に変換して溶接変形、残留応力を算出する解析方法が開示されている。
特許文献3には、溶接構造物の残留応力有限要素法で解析する際、溶接線から解析モデル境界までの軸方向距離を制約することにより解析時間を短縮する解析方法が開示されている。
特許文献4には、拘束条件で得られた溶接固有ひずみを用いて、溶接後拘束条件をなくしたときの溶接変形を求める解析方法が開示されている。
特許文献5には、溶接近傍のみ非線形解析を行い、被溶接物につき線形解析を行うことにより両領域の限界面における反力を計算し、その反力の差が所定の範囲内なるように収束演算することで、溶接変形を短時間に算出する解析方法が開示されている。
熱弾塑性解析法による変形解析は、溶接開始から完了までの過程を細かい時間間隔で解析するため、実溶接に近い変形履歴を模擬することができる。しかしながら、大型溶接構造物の場合、解析モデルには膨大な要素数が必要であり、熱弾塑性解析に長時間の計算時間を必要とするため、現実的には不可能なケースも多い。
特許文献1に示した解析方法であれば、熱弾塑性解析を行うものが2次元モデルであるため、計算時間の短縮はできる。しかし、複雑な大型構造物の溶接変形に対し、2次元モデルを適用できないケースが多く、適用が制限されているという問題がある。
特許文献1に示した解析方法であれば、熱弾塑性解析を行うものが2次元モデルであるため、計算時間の短縮はできる。しかし、複雑な大型構造物の溶接変形に対し、2次元モデルを適用できないケースが多く、適用が制限されているという問題がある。
固有ひずみ法変形解析は、構造物の弾性解析のみで解析することによる変形算出ができるため、計算時間を大幅に短縮できる。しかし、計算により得られた固有ひずみ分布を溶接構造物全体に与えることは容易ではない。
即ち、各単元モデルの熱弾塑性解析の際に用いたモデルのメッシュは、実際に計算される溶接構造物のモデルのメッシュと異なるため、単元モデルからデータベース化した固有ひずみを溶接構造物のモデルに与える際に、座標間の換算により誤差が生じ易い。その結果、場合によっては実溶接と異なる固有ひずみ分布を構造物モデルに付与してしまうという問題が発生する。
さらに、固有ひずみを求める際に用いられた単元モデルの拘束条件は、一般的には、溶接構造物の拘束条件を反映できていない。このため、構造物モデルへ付与することになる、単元モデルから得られた固有ひずみの分布が、実溶接により生じた実構造物の固有ひずみと異なるという問題もある。このような実構造物と異なった拘束条件で得られた固有ひずみを構造物モデルへ付与することにより計算された変形が、実溶接構造物の変形傾向に合わないケースが多い。
特許文献2、3、4には様々な固有ひずみ法を開示しており、座標系の換算や、解析モデルの簡略化又は計算範囲の縮小により、計算精度の向上や計算時間の短縮を図ることができる。しかし、構造物全体の変形に生じるひずみまたは応力分布などの計算ができないという問題がある。
特許文献5に開示された解析方法は、大型構造物に対しては収束計算に時間がかかるという問題がある。
そこで前記事情を鑑みて、本発明の目的は、大型溶接構造物の溶接変形予測における計算精度の向上および計算時間の短縮を両立することである。
前記目的を達成する本願発明の特徴は、固有ひずみの計算や構造物への付与(マッピング)を行わず、且つ溶接部およびその近傍のみに熱弾塑性解析を行う、いわば部分熱弾塑性解析(「溶接変形解析」と称する場合もある。)である。
具体的には、本発明の溶接変形解析方法は、メッシュ作成工程による解析対象となる構造物を全体モデルにし、部分熱弾塑性解析が必要な溶接部を含む局部モデルを、変形解析対象となる弾性解析が必要な溶接構造物の全体モデルから抽出する工程と、前記局部モデルと全体モデルの残り部分の境界部に拘束を与えて、前記局部モデルに対して熱弾塑性解析を行う工程と、前記局部モデルの熱弾塑性解析の工程が完了した後、前記熱弾塑性解析により得られた溶接終了後の最終のひずみ分布と応力分布とを含有する前記局部モデルを、前記全体モデルの残り部分に貼付けることにより全体モデルを再構築する工程と、全体モデルの弾性解析を行うことにより構造物の溶接変形と応力を求める工程を含んでいる。
前記局部モデルは前記全体モデルから直接抽出したものであるため、前記局部モデルのメッシュに構成された要素および節点の構成と相対座標関係が、前記全体モデルの一部と一致している。
前記局部モデルの熱解析を行う際に、実際の構造物に存在していない、前記局部モデルを抽出する工程により形成された、局部モデルと全体モデルの残り部分との解析上の境界部に、適当な拘束条件を与える。その際、境界部の拘束条件は、実構造物の拘束状況を考慮して設定する。
前記全体モデルの再構築は、前記局部モデルの熱弾塑性解析を行ったときに与えられた拘束条件を満たし、熱弾塑性解析により得られた溶接終了後の最終のひずみ分布と応力分布とを含有する前記局部モデルを、前記全体モデルに貼り付けてもよい。或いは、前記局部モデルを抽出した後の全体モデルの残り部分を、熱弾塑性解析により得られた溶接終了後の最終のひずみ分布と応力分布とを含有する前記局部モデルに貼り付けてもよい。
全体モデルの弾性解析による構造物の溶接変形を計算するときには、前記局部モデルの熱弾塑性解析を行ったときに与えられた拘束条件を解放し、実構造物に近い拘束条件を設定して弾性解析を行い、構造物全体の溶接変形と応力を求める。
詳細は、後記する。
詳細は、後記する。
本発明によれば、大型溶接構造物の溶接変形予測における計算精度の向上および計算時間の短縮を両立することができる。
次に、本発明を実施するための形態(以下、「実施形態」という。)について、適宜図面を参照しながら説明する。
≪構成≫
図1は、本実施形態の解析装置のハードウェア構成およびソフトウェア構成の構成図である。図1に示す解析装置100は、入力部110、出力部120、制御部130および記憶部140をハードウェア構成として備えるコンピュータである。この解析装置100は、記憶部140にCAD(Computer-Aided Design)データ141を記憶するとともに、ソフトウェア構成としてのモデル生成・加工部142、解析部143およびパラメータ取得部144を実現する溶接変形解析方法用のプログラムを記憶する。
図1は、本実施形態の解析装置のハードウェア構成およびソフトウェア構成の構成図である。図1に示す解析装置100は、入力部110、出力部120、制御部130および記憶部140をハードウェア構成として備えるコンピュータである。この解析装置100は、記憶部140にCAD(Computer-Aided Design)データ141を記憶するとともに、ソフトウェア構成としてのモデル生成・加工部142、解析部143およびパラメータ取得部144を実現する溶接変形解析方法用のプログラムを記憶する。
入力部110は、例えば、ユーザの操作を受け付けるマウスやキーボードであり、ユーザから入力された信号を制御部130に伝送する入力インタフェースも含む。
出力部120は、例えば、解析対象となる溶接構造物のCADデータやそのモデル、数値解析の計算結果などを表示するディスプレイであり、制御部130からの命令(描画命令を含む)に従い、所定の画像を表示する出力インタフェースも含む。
制御部130は、例えば、CPU(Central Processing Unit)であり、記憶部140が記憶するプログラムに記述されたコードを読み出し、相応の情報処理(数値解析の演算処理を含む)を実行する。
記憶部140は、例えば、ROM(Read Only Memory:「記憶媒体」ともいう)、RAM(Random Access Memory)、HDD(Hard Disc Drive)であり、前記したように、CADデータ141、モデル生成・加工部142、解析部143およびパラメータ取得部144として機能させるプログラムを記憶する。
CADデータ141は、解析対象となる溶接構造物を、CADによる設計で生成した、例えば、3次元のデータである。CADデータ141には溶接構造物の物性を示す物性値も含まれる。CADデータ141は、不図示の通信インタフェースを備え、ネットワークに接続した解析装置100が、そのネットワークに接続したデータベースサーバなどの外部から取得してもよい。
モデル生成・加工部142は、CADデータ141から、例えば溶接構造物の3次元の解析モデル(単に、「モデル」と称する場合がある)を生成したり、入力部110からの入力に従って加工したりする。モデルの生成には、CADデータ141から解析対象の形状を取得するとともに、入力部110からの入力により指定されたメッシュ構成に従って、メッシュの生成も行われる。形状を取得すれば、例えば当該溶接構造物の占める位置は、モデル内に予め設定した三次元座標(本実施形態では座標軸は直交しているものとして説明する。)の座標値によって定まるとともに、さらにその座標値に位置する部品を構成する物質の物性値が定まり、記憶部140に記憶される。前記物性値は、CADデータ141から取得される。
また、モデルの加工には、例えば、後記する局部モデルの抽出や、メッシュの節点の拘束、その解放、全体モデルの再構築などの処理が含まれる。
生成または加工されたモデル(全体モデルや局部モデルなど)は記憶部140に記憶される。
生成または加工されたモデル(全体モデルや局部モデルなど)は記憶部140に記憶される。
解析部143は、生成または加工したモデルに対し、数値解析を行い、所定の解析結果を求める。本実施形態で行われる数値解析には、例えば、非線形解析である熱弾塑性解析や、線形解析である弾性解析が含まれる。もし、本実施形態で固有ひずみ法を用いるのであれば、それは解析部143が行う。
パラメータ取得部144は、例えば、入力部110の入力によって指定されたパラメータを取得する。このパラメータは、モデル生成・加工部142によるモデルの生成、加工や解析部143による数値解析を実行するのに必要な値、条件などであり、具体的には、メッシュ構成に関する初期条件(例:使用するメッシュの数、サイズ、形状)、全体モデルから局部モデルを抽出するときに定める局部モデルと全体モデルの残り部分との境界部、または、メッシュの節点の拘束条件(例:拘束対象のメッシュの節点、その節点の拘束する方向(X方向、Y方向、Z方向の少なくとも一つ))である。
なお、メッシュの節点の「拘束」とは、数値解析を行ってもその節点の変位量をゼロにすることである。この拘束により、解析モデルの数値解析における位置の基準が定まる。3次元解析モデルにおいては、数値解析を行ううえで少なくとも6つの自由度の拘束が必要となる。ただし、拘束する節点の位置や方向は問わない。
以上で、本実施形態の解析装置の構成に関する説明を終了する。
以上で、本実施形態の解析装置の構成に関する説明を終了する。
≪処理≫
次に、本実施形態の解析装置100で実行される処理について説明する。この処理、つまり、主に溶接変形解析に関する処理は、制御部130による制御のもと、記憶部140に記憶されたプログラムを、RAMなどの記憶領域に読み出すことによって実現される。
図2は、本実施形態に係る有限要素法による溶接変形解析の処理を示すフローチャートである。この制御の主体は、制御部130である。
次に、本実施形態の解析装置100で実行される処理について説明する。この処理、つまり、主に溶接変形解析に関する処理は、制御部130による制御のもと、記憶部140に記憶されたプログラムを、RAMなどの記憶領域に読み出すことによって実現される。
図2は、本実施形態に係る有限要素法による溶接変形解析の処理を示すフローチャートである。この制御の主体は、制御部130である。
この処理では、まず、解析対象とした溶接構造物を3次元解析モデル化し、全体モデルを構築する(S1)。この全体モデルを構築する際に、熱弾塑性解析による十分な解析精度が得られるようにするために、溶接構造物に含まれる溶接部およびその近傍を含む熱弾塑性解析が必要な範囲に対しては、細かい(密な)メッシュを設定する。これ以外の範囲に対しては、弾性解析のみを行うため、粗いメッシュを設定する。
上記弾性解析のみを行う範囲のメッシュのサイズは、上記熱弾塑性解析が必要な範囲のメッシュサイズの数倍から数十倍である。その結果、全体モデルの要素数が減少し、計算時間の短縮ができる。これらのメッシュのサイズや要素数は、パラメータ取得部144のパラメータとして取得したものである。
次に、構築された全体モデルから、溶接部およびその近傍を含む、熱弾塑性解析が必要となる範囲の要素を抽出し、熱弾塑性解析用の局部モデルとする(S2)。この局部モデルは、全体モデルからそのまま抽出したものであるため、モデルおよびメッシュの作成が1回に限られ、メッシュの作成工程を短縮できる。局部モデルの抽出においては、パラメータ取得部144のパラメータとして、局部モデルの境界部と全体モデルの残り部分の境界部とが取得され、設定される。
次に、局部モデルを抽出した後、局部モデルの境界部に適当な拘束条件を与えて境界部を拘束し(S3)、局部モデルの熱弾塑性解析を行う(S4)。上記局部モデルの境界部は、上記全体モデルの残り部分との境界面である。また、境界部の拘束条件は、実構造物の拘束状況を考慮して設定する。前記拘束は、全体モデルの残り部分に形成された境界部についても行う。この拘束は、パラメータ取得部144のパラメータとして設定される。
なお、局部モデルの境界部の節点が拘束される状態で熱弾塑性解析を行うことが好ましいが、境界部の節点を拘束していない状態で局部モデルの熱弾塑性解析を実施してもよい。また、局部モデルの境界部(および全体モデルの残り部分の境界部)の節点の拘束は、すべての節点について行ってもよいし、一部の節点について行ってもよい。
なお、局部モデルの境界部の節点が拘束される状態で熱弾塑性解析を行うことが好ましいが、境界部の節点を拘束していない状態で局部モデルの熱弾塑性解析を実施してもよい。また、局部モデルの境界部(および全体モデルの残り部分の境界部)の節点の拘束は、すべての節点について行ってもよいし、一部の節点について行ってもよい。
次に、局部モデルの熱弾塑性解析を終了した後、解析の結果となる熱分布、応力分布およびひずみ分布を含む局部モデルを、全体モデルの残り部分に貼り付け(合わせ)、或いは、局部モデルを除く全体モデルの残り部分を、熱弾塑性解析の結果を含んでいる局部モデルに貼り付ける(合わせる)ことにより、弾性解析用の全体モデルを再構築する(S5)。
なお、局部モデルの熱弾塑性解析時に局部モデルの境界部に拘束を設定しない場合、弾性解析の全体モデルを再構築する際、全体モデルの残り部分の境界部の節点を、局部モデルの境界部の節点と合わせるように移動させ、全体モデルのメッシュを局部モデルに貼り付けることにする。
なお、局部モデルの熱弾塑性解析時に局部モデルの境界部に拘束を設定しない場合、弾性解析の全体モデルを再構築する際、全体モデルの残り部分の境界部の節点を、局部モデルの境界部の節点と合わせるように移動させ、全体モデルのメッシュを局部モデルに貼り付けることにする。
全体モデルを再構築した後、上記局部モデルの境界部の拘束条件を外して境界部の拘束を解放した後(S6)、再構築された全体モデルを用い、実際の構造物と同じ拘束条件(例:構造物の端部の拘束。詳細は後記。)を設定し、弾性解析を実施することにより、構造物の変形を計算する(S7)。
以上で、本実施形態の解析装置の処理に関する説明を終了する。
以上で、本実施形態の解析装置の処理に関する説明を終了する。
次に、本実施形態の解析方法を具体的な形状が定められた溶接構造物に適用した場合を実施例1として説明する。
図3は、解析対象としての溶接構造物の構成を示す図である。この溶接構造物は、パイプ状構造物であり、被溶接材Aと被溶接材Bとの周溶接Cにより構成される。溶接後の構造物は、長さが約5m、直径が約200mm、厚さが7mm〜13mmである。
図3は、解析対象としての溶接構造物の構成を示す図である。この溶接構造物は、パイプ状構造物であり、被溶接材Aと被溶接材Bとの周溶接Cにより構成される。溶接後の構造物は、長さが約5m、直径が約200mm、厚さが7mm〜13mmである。
この図3に示した溶接構造物を有限要素モデル化することにより3次元モデルを生成し、メッシュ作成工程(モデル生成・加工部142による処理の一つ)を実施することにより全体モデルを構築する。
図4は、図3の溶接構造物の全体モデルを示す図である。この全体モデルは、被溶接材Aの一部に対応するモデルである被溶接部1、被溶接材Bに対応するモデルである被溶接部2、および周溶接(溶接箇所)Cに対応するモデルである溶接部3により構成された局部モデル4(網掛け表示)と、被溶接材Aの残り部分に対応するモデルである被溶接部を含む全体モデルの残り部分5から構成される。
上記全体モデルのメッシュを作成する際に、溶接部3およびその近傍を含む熱弾塑性解析が必要な範囲を局部モデル4とし、細かいメッシュを設定する。これ以外の範囲、つまり、全体モデルの残り部分5に対しては弾性解析のみを行うため、粗いメッシュを作成する。上記熱弾塑性解析が必要な局部モデル4のメッシュのサイズは、上記弾性解析が必要な全体モデルの残り部分5のメッシュサイズの数分の一〜数十分の一である。ただし、説明の便宜上、図4においては、設定したメッシュの図示は省略してある。
上記熱弾塑性解析が必要な局部モデル4のメッシュは、弾性解析のみが行われる全体モデルの残り部分5のメッシュよりかなり小さいが、全体モデルの残り部分5と局部モデル4と共に同程度サイズのメッシュを設定しても構わない。これは、同程度の要素数である場合でも、弾性解析に係る時間が熱弾塑性解析に係る時間の数千分の一以下であるためである。即ち、全体モデルの残り部分5のメッシュを局部モデル4のそれと同程度に細かく分割すれば要素数が増加するものの、増加分のメッシュは弾性解析しか行われず、かつ、全体の解析時間は主に局部モデル4の熱弾塑性解析に支配されるため、全体の解析時間の増加分は、熱弾塑性解析の時間に比べ無視できる程度に小さくなる。
本実施例のモデルのメッシュ分割は実溶接条件と溶接ビード断面形状を考慮して決定する。要素タイプは、解析精度の高い六面体要素を用いるのが好ましい。また、四面体要素、三角柱要素を用いてもよいし、或いはこれらの混合要素を用いてもよい。
六面体メッシュを作成した全体モデルの規模は約100,000要素、150,000節点程度である。そのうち、溶接近傍を含む熱弾塑性解析が必要な局部モデル部分の規模は約25,000要素、40,000節点程度である。溶接は、レーザ・アークハイブリッド溶接により1パス溶接で行なったため、熱源モデルは、線状Gaussian熱源(レーザ熱源に相当)と点状Gaussian熱源(アーク熱源に相当)の複合移動熱源を用いた。解析の入熱条件は、実溶接条件および溶接ビード断面形状を考慮して決定した。溶接時間は実溶接と同じ90秒とした。
細かいメッシュが設定された局部モデルを含む全体モデルを作成した後、作成された全体モデルにおいて、溶接部およびその近傍を含む、熱弾塑性解析が必要な範囲を局部モデル4とし、局部モデルの抽出工程(モデル生成・加工部142による処理の一つ)を行う。
図5は、抽出された局部モデルを示す図である。この局部モデル4のメッシュは、図4に示す全体モデルの一部である局部モデル4のメッシュ構成と完全に一致している。例えば、抽出したメッシュの節点の座標値は、抽出前のそれと同一の値を示す。なお、抽出する際、局部モデル4には、全体モデルの残り部分5との境界面である境界部6が形成され、境界部6の位置は、パラメータ取得部144により設定される。
次に、抽出された局部モデル4を用い、熱弾塑性解析を行う。この局部モデル4の熱弾塑性解析には三つのステップ1〜3が含まれる。
局部モデル4の熱弾塑性解析のステップ1では、熱解析を実施する。具体的な流れとして、まず、実溶接条件を模擬する入熱条件を設定し、次に、溶接開始から完了までの過程を細かい時間間隔で熱解析を行い、実溶接に近い熱履歴を計算する。その結果、溶接部3およびその近傍の温度分布の履歴が得られる。
なお、上記局部モデル4の熱履歴を計算する際には、局部モデル4を用いたが、全体モデルを用いて熱解析を実施しても構わない。しかし、この場合、計算により得られた溶接部3およびその近傍の温度分布履歴(熱履歴)は、実溶接の熱履歴を模擬できる必要がある。
なお、上記局部モデル4の熱履歴を計算する際には、局部モデル4を用いたが、全体モデルを用いて熱解析を実施しても構わない。しかし、この場合、計算により得られた溶接部3およびその近傍の温度分布履歴(熱履歴)は、実溶接の熱履歴を模擬できる必要がある。
局部モデル4の熱弾塑性解析のステップ2では、次のステップ3に実施予定の熱弾塑性解析際に必要な条件となる境界部6の拘束条件を設定する。この境界部6は、上記局部モデル4の抽出工程により形成された、局部モデル4と全体モデルの残り部分5との解析上の境界面であり、実際の構造物にはこのような境界面が存在していない。なお、全体モデルの残り部分5の境界面の節点も拘束する。
図6は、局部モデルの境界部の拡大図である。図6を参照して、局部モデルの境界部の拘束条件の設定について説明する。本実施例の解析モデルは、符号10、11で示すような六面体要素により構成される。一つ一つの六面体要素毎に、符号20〜23で示すような節点は8個ある。ここで、節点20は局部モデル4の境界部6の境界面に所属せず、節点21〜23は境界部6の境界面に所属する。従って、局部モデル4の境界部6の拘束条件の設定は、節点21〜23のような境界部6の境界面にある節点に拘束条件を与える。
境界部6の拘束条件は、実構造物の拘束状況を考慮して設定する。本実施例の場合、図6に示す節点21〜23のような境界部6のすべての節点において、同図に示す座標方向であるX方向、Y方向、Z方向の変位量を0とする。
局部モデルの熱弾塑性解析のステップ3では、上記の境界部6の拘束条件を設定した後、熱弾塑性解析のステップ1にて得られた局部モデル4の熱履歴の解析結果をインプット条件として、非線形の熱弾塑性解析を実施する。ステップ3の解析により、局部モデル4のひずみ履歴および応力履歴が得られる。
次に、上記の局部モデル4の熱弾塑性解析の結果を用い、構造物の全体モデルを再構築し、この全体モデルに対し弾性解析を行い、構造物の変形を計算する。この全体モデルの弾性解析にも三つのステップ1〜3が含まれる。
なお、図7は、再構築された全体モデルを示す図である。
なお、図7は、再構築された全体モデルを示す図である。
全体モデルの弾性解析のステップ1では、熱弾塑性解析結果が含まれる局部モデル4を全体モデルに貼り付け、図7に示すような全体モデルを再構築する。局部モデル4の熱弾塑性解析を実施した際に、局部モデル4と全体モデルの残り部分5との境界面の節点は拘束されているため、局部モデル4の全体モデルの残り部分5への貼り付けは簡単に実現できる。また、局部モデル4の熱弾塑性解析により得られたひずみ分布および応力分布も全体モデルの残り部分5に貼り付けられる。その結果、再構築した全体モデルには、局部モデル4の熱弾塑性解析により得られる溶接近傍のひずみ分布、応力分布などの情報が含まれている。
上記の全体モデルを再構築する際に、局部モデル4を全体モデルの残り部分5に貼り付けるが、逆に局部モデル4を抽出することにより残った全体モデルの残り部分5のメッシュを、前記熱弾塑性解析結果を含有する局部モデル4に貼り付けても構わない。
全体モデルの弾性解析のステップ2では、局部モデル4の熱弾塑性解析を行う際に設定された局部モデルの境界部6の拘束条件を解放し、実構造物の拘束を模擬できる拘束条件を設定する。本実施例においては、図7に示すように、全体モデルの端部7に所属する節点41にX方向、Y方向、Z方向三つの拘束を設定し、節点42にY方向とZ方向の拘束を与え、節点43にZ方向の拘束を設定する。
全体モデルの弾性解析のステップ3では、局部モデル4のひずみ分布と応力分布が含まれ、新たな拘束条件(端部7での拘束)が設定された全体モデルを用い、弾性解析を行い、構造物の変形を計算する。
〔従来との比較〕
比較として、図3に示すような本実施例の解析対象を用い、従来の熱弾塑性解析法、従来の固有ひずみ法、本発明の部分熱弾塑性解析法により解析を行い、解析結果を得るまでの所有時間および解析精度について三つの解析方法を比較した。
比較として、図3に示すような本実施例の解析対象を用い、従来の熱弾塑性解析法、従来の固有ひずみ法、本発明の部分熱弾塑性解析法により解析を行い、解析結果を得るまでの所有時間および解析精度について三つの解析方法を比較した。
従来の熱弾塑性解析を行った際に、本実施例の全体モデルを用いて実施した。そのモデルのメッシュ構成は図7に示す本実施例に用いられた全体モデルと同じである。ただ、従来の熱弾塑性解析は全体モデルに対し熱弾塑性解析を実施するため、この全体モデルには温度、応力、ひずみ、変形量などの熱弾塑性解析の結果となる情報が一斉に含まれる。溶接は、レーザ・アークハイブリッド溶接により1パス溶接で行なったため、熱源モデルは、線状Gaussian熱源(レーザ熱源に相当)と点状Gaussian熱源(アーク熱源に相当)の複合移動熱源を用いた。解析の入熱条件は、実溶接条件および溶接ビード断面形状を考慮して決定した。溶接時間は実溶接と同じ90秒とした。
上記の解析モデルおよび溶接条件(入熱条件)を設定した後、本実施例に記載された前記局部モデル4の熱弾塑性解析法と同様な解析法を用い、全体モデルの熱弾塑性解析を行った(下記ステップ1〜3)。
熱弾塑性解析のステップ1では、熱解析を実施する。実溶接条件を模擬する入熱条件を設定し、溶接開始から完了までの過程を細かい時間間隔を設定して熱解析を行い、実溶接に近い熱履歴を計算する。その結果、溶接部3およびその近傍を含む全体モデルの温度分布の履歴が解析結果として得られる。
ステップ2では、実構造物の拘束を模擬できる、溶接変形解析に必要な境界条件となる拘束条件を設定する。本実施例に実施した全体モデルの弾性解析の際に設定された拘束条件と同様に、図7に示すような全体モデルの端部7に所属する節点41にX方向、Y方向、Z方向三つの拘束を設定し、節点42にY方向とZ方向の拘束を与え、節点43にZ方向の拘束を設定した。
ステップ3では、上記端部7の拘束条件を設定した後、熱弾塑性解析のステップ1に得られた全体モデルの熱履歴の解析結果をインプット条件として、非線形の熱弾塑性解析を実施する。実施することで、全体モデルのひずみ履歴および応力履歴、または各節点の変位履歴が得られる。その変位履歴から、溶接構造物の全体の変形を計算することができる。
また、従来の固有ひずみ法を実施する際に、まず熱弾塑性解析法により溶接部3およびその近傍の固有ひずみ分布を計算する必要がある。そのため、本実施例に用いられた熱弾塑性解析の局部モデル4のメッシュ構成を用いた。
図8は、局部モデルを示す図である。この図は、図5に示した局部モデル4の視点を変え、かつメッシュ構成が異なるものに等しく、局部モデル4の端部8が図示されている。
解析に必要な入熱条件も本実施例と同じに設定したが、拘束条件については図8に示すように、局部モデル4と全体モデル5の境界部6の反対側の端部8に適当な拘束を設定する。即ち、端部8の節点81にX方向、Y方向、Z方向三つの拘束を設定し、端部8の節点82にY方向とZ方向の拘束を与え、端部8の節点83にZ方向の拘束を設定する。
解析に必要な入熱条件も本実施例と同じに設定したが、拘束条件については図8に示すように、局部モデル4と全体モデル5の境界部6の反対側の端部8に適当な拘束を設定する。即ち、端部8の節点81にX方向、Y方向、Z方向三つの拘束を設定し、端部8の節点82にY方向とZ方向の拘束を与え、端部8の節点83にZ方向の拘束を設定する。
拘束条件を設定した後、前記局部モデルの熱弾塑性解析法と同様な解析法を用い、局部モデル4の熱弾塑性解析を行った。その結果、溶接部3およびその近傍の温度履歴、ひずみ履歴、応力履歴などが固有ひずみとして得られた。
次に、上記の結果から溶接近傍の固有ひずみを抽出し、図7に示す本実施例に用いられた全体モデルと同様なメッシュ構成を含み、現時点ではすべての要素および節点にひずみ、応力、温度などが一切にゼロの状態である全体モデルに、その固有ひずみを与える。
上記の固有ひずみを全体モデルに与えた後、全体モデルの拘束条件を設定し、弾性解析を行った。拘束条件の設定は、本実施例と同様に、図7に示すように、全体モデルの端部7に所属する節点41にX方向、Y方向、Z方向三つの拘束を設定し、節点42にY方向とZ方向の拘束を与え、節点43にZ方向の拘束を設定した。その後、全体モデルの弾性解析を行い、全体の変形を計算した。
図9は、本発明の溶接変形解析方法と、従来の熱弾塑性法、従来の固有ひずみ法の解析時間および解析精度の比較表である。変形量の解析精度が最も高いのは従来熱弾塑性解析法であるが(「◎」で表記)、本発明の解析方法でも、従来の固有ひずみ法(「△」で表記)よりも優れた、十分な解析精度が得られる(「○」で表記)。
一方、解析時間については、従来の熱弾塑性解析法が最も長く(「×」で表記)、解析対象によっては解析不可能のケースも多い。これに対し、本発明の解析方法は、固有ひずみ法(「△」で表記)よりも優れた、解析時間の大幅な短縮が実現できる(「○」で表記)。また、本発明の解析方法は、固有ひずみの計算と全体モデルへの付与工程とが不要となるため、解析時間が最も短い。
次に、本実施形態の解析方法を具体的な形状が定められた溶接構造物に適用した別の場合を実施例2として説明する。実施例2は、実施例1の境界部の拘束条件の設定を変更した場合に関する例である。
本実施例の解析対象は、実施例1に用いられた対象と同じ、図3に示すパイプ状溶接構造物である。また、解析方法は、境界部の拘束条件の設定に関係する処理以外は実施例1と同様である。つまり、実施例2の解析方法は、まずメッシュを生成する工程による解析対象となる構造物を全体モデルにし、熱弾塑性解析が必要な溶接部を含む局部モデルを、解析対象となる弾性解析が必要な溶接構造物の全体モデルから抽出する工程と、前記局部モデルの境界部と全体モデルの残り部分の境界部に拘束を与えることなく、前記局部モデルに対して熱弾塑性解析を行う工程と、前記局部モデルの熱弾塑性解析工程が完了した後、熱弾塑性解析により得られた溶接終了後の最終のひずみ分布と応力分布を含有する前記局部モデルを、熱弾塑性解析に起因する変位量を考慮しつつ、前記全体モデル残り部分へ貼り付けることにより全体モデルを再構築する工程と、全体モデルの弾性解析による構造物の溶接変形と応力を求める工程を含んでいる。
全体モデルのメッシュを生成する工程、熱弾塑性解析が必要な局部モデルを抽出する工程、局部モデルの熱弾塑性解析を行う工程、および全体モデルの弾性解析を行う工程、といった4つの解析ステップは、実施例1に示すステップと一致している。しかし、局部モデルの熱弾塑性解析の際に必要な境界部の拘束条件の設定や、全体モデルの再構築については、実施例1と異なる。
本実施例では、局部モデルの境界部を拘束することはないが、局部モデルの熱弾塑性解析を行うときには、局部モデルの端部は拘束する。局部モデルの端部の拘束の詳細について、図8を参照して説明する。本実施例においては、図8に示す局部モデル4は、図3に示す全体モデルから抽出したモデルであり、メッシュ構成はその全体モデルの一部と同一となる。
本実施例の局部モデル4の端部の拘束条件を設定する際には、図8に示すように、局部モデル4の境界部6の反対側となる端面である端部8に適当な拘束を設定する。例えば、節点81にX方向、Y方向、Z方向三つの拘束を設定し、節点82にY方向とZ方向の拘束を与え、節点83にZ方向の拘束を設定する。
また、弾性解析が必要な全体モデルを再構築する際に、全体モデルの残り部分5を、熱弾塑性解析結果を含んでいる局部モデル4に貼り付けることにする。これは、局部モデル4の熱弾塑性解析の際に、局部モデル4と全体モデルの残り部分5それぞれの境界部に拘束条件を設定してなかったため、局部モデル4の境界部6の節点には変位が発生し、全体モデルの残り部分の境界部の節点の座標と合わない、という現象に対応するためである。具体的には、熱弾塑性解析により発生した変位量の分だけ(変位量を相殺するように)全体モデルの残り部分5を移動させて貼り付ける。
図10は、実施例2における全体モデルの再構築の様子を示す図である。図10に示すように、局部モデル4の熱弾塑性解析を実施する前の全体モデルの残り部分5の境界部61の節点を移動させ、熱弾塑性解析を実施した局部モデル4の境界部6の節点と合わせる。
上記の全体モデルの残り部分5の境界部61の節点を、局部モデル4の境界部6の節点に合わせることには相応の時間がかかるものの、この節点合わせは境界部61の節点に限られるため、従来の固有ひずみ法のような溶接部近傍の膨大な節点の固有ひずみを全体モデルの残り部分5へ付与する工程と比べ、大幅に短縮できる。
なお、前記した貼り付けとは逆に、局部モデル4を全体モデルの残り部分5に貼り付けるようにして全体モデルを再構築してもよい。
なお、前記した貼り付けとは逆に、局部モデル4を全体モデルの残り部分5に貼り付けるようにして全体モデルを再構築してもよい。
全体モデルの残り部分5の境界部61の節点を局部モデル4の境界部6の節点と合わせることにより、全体モデルを再構築した後、再構築された全体モデルを用い、弾性解析を行い、溶接変形を計算する。
実施例2による解析方法の解析時間および解析精度は、図9に示した比較表と同様となる。
実施例2による解析方法の解析時間および解析精度は、図9に示した比較表と同様となる。
≪まとめ≫
本実施形態により、以下の効果が生じる。つまり、本実施形態の溶接変形解析によれば、溶接近傍の熱弾塑性解析により残留応力と残留ひずみの計算ができ、且つ、固有ひずみの計算や測定、更に固有ひずみの全体モデルへの付与などの工程を行わずに、構造物全体の弾性解析もできる。その結果、大型溶接構造物の変形予測が短時間かつ高精度でできる。
また、実施例2に示したように、局部モデルを拘束しないときの溶接変形解析であっても、大型溶接構造物の変形予測が短時間かつ高精度でできる。
本実施形態により、以下の効果が生じる。つまり、本実施形態の溶接変形解析によれば、溶接近傍の熱弾塑性解析により残留応力と残留ひずみの計算ができ、且つ、固有ひずみの計算や測定、更に固有ひずみの全体モデルへの付与などの工程を行わずに、構造物全体の弾性解析もできる。その結果、大型溶接構造物の変形予測が短時間かつ高精度でできる。
また、実施例2に示したように、局部モデルを拘束しないときの溶接変形解析であっても、大型溶接構造物の変形予測が短時間かつ高精度でできる。
≪その他≫
なお、前記した実施形態は、本発明を実施するための好適なものであるが、その実施形式はこれらに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形することが可能である。
なお、前記した実施形態は、本発明を実施するための好適なものであるが、その実施形式はこれらに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形することが可能である。
例えば、本実施形態では、3次元の解析モデルを用いた数値解析について説明した。しかし、1次元や2次元の解析モデルを用いてもよい。
また、本実施形態では、局部モデルを抽出するときに形成される局部モデルの境界部は、2以上形成してもよいし、メッシュの面に沿うことなく、メッシュの面と平行または非平行にメッシュを切断するように形成してもよい。メッシュを切断するように形成する場合は、その切断により形成された切断面を新たなメッシュの面としてもよい。
また、本実施形態では、局部モデルの境界部は、入力部から入力されたパラメータとして、ユーザが指定するようにした。しかし、局部モデルを抽出するときに、所望の解析時間および所望の解析精度を達成するように前記境界部を設定することができるプログラムを用いるようにしてもよい。
その他ハードウェア、ソフトウェア、等の具体的な構成について、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。
本発明の溶接変形解析方法は、形状を問わず、大型の溶接構造物の溶接変形の解析に有効である。
1 被溶接部
2 被溶接部
3 溶接部
4 局部モデル
5 全体モデルの残り部分
6 境界部
100 解析装置
110 入力部
120 出力部
130 制御部
140 記憶部
141 CADデータ
142 モデル生成・加工部
143 解析部
144 パラメータ取得部
2 被溶接部
3 溶接部
4 局部モデル
5 全体モデルの残り部分
6 境界部
100 解析装置
110 入力部
120 出力部
130 制御部
140 記憶部
141 CADデータ
142 モデル生成・加工部
143 解析部
144 パラメータ取得部
Claims (3)
- 溶接構造物をモデル化して数値解析を行う解析装置の溶接変形解析方法において、
前記解析装置の記憶部は、前記数値解析を行うためのメッシュを用いて生成した前記溶接構造物の全体モデルを記憶しており、
前記解析装置の制御部は、
前記全体モデルから、前記溶接構造物の溶接箇所に対応するモデルを含む局部モデルを抽出するステップと、
前記抽出した局部モデルの、前記全体モデルの残り部分との境界部を構成するメッシュの節点を拘束するステップと、
前記局部モデルにおいて、前記数値解析としての熱弾塑性解析を行うステップと、
前記局部モデルと前記全体モデルの残り部分とを合わせ、全体モデルを再構築するステップと、
前記再構築した全体モデルにおいて、前記拘束を解放し、前記数値解析としての弾性解析を行うことで、前記溶接構造物の変形を計算するステップと、を実行する
ことを特徴とする溶接変形解析方法。 - 前記全体モデルの残り部分に用いたメッシュは、前記局部モデルに用いたメッシュよりも粗である
ことを特徴とする請求項1に記載の溶接変形解析方法。 - 溶接構造物をモデル化して数値解析を行う解析装置の溶接変形解析方法において、
前記解析装置の記憶部は、前記数値解析を行うためのメッシュを用いて生成した前記溶接構造物の全体モデルを記憶しており、
前記解析装置の制御部は、
前記全体モデルから、前記溶接構造物の溶接箇所に対応するモデルを含む局部モデルを抽出するステップと、
前記局部モデルにおいて、前記数値解析としての熱弾塑性解析を行うステップと、
前記熱弾塑性解析により前記局部モデルに生じた変位量を相殺するように前記局部モデルまたは前記全体モデルの残り部分を移動して、前記局部モデルと前記全体モデルの残り部分とを合わせ、全体モデルを再構築するステップと、
前記再構築した全体モデルにおいて、前記数値解析としての弾性解析を行うことで、前記溶接構造物の変形を計算するステップと、を実行する
ことを特徴とする溶接変形解析方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010022232A JP2011159213A (ja) | 2010-02-03 | 2010-02-03 | 溶接変形解析方法 |
CN2011100366048A CN102152016A (zh) | 2010-02-03 | 2011-01-31 | 焊接变形分析方法 |
EP11152998A EP2363819A1 (en) | 2010-02-03 | 2011-02-02 | Method for simulation of welding distortion |
US13/020,075 US20110213594A1 (en) | 2010-02-03 | 2011-02-03 | Method for simulation of welding distortion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010022232A JP2011159213A (ja) | 2010-02-03 | 2010-02-03 | 溶接変形解析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2011159213A true JP2011159213A (ja) | 2011-08-18 |
Family
ID=43982376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010022232A Pending JP2011159213A (ja) | 2010-02-03 | 2010-02-03 | 溶接変形解析方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20110213594A1 (ja) |
EP (1) | EP2363819A1 (ja) |
JP (1) | JP2011159213A (ja) |
CN (1) | CN102152016A (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101291058B1 (ko) | 2012-01-18 | 2013-08-07 | 삼성중공업 주식회사 | 변형예측시스템 및 변형예측방법 |
JP2017161981A (ja) * | 2016-03-07 | 2017-09-14 | 株式会社東芝 | 解析装置、解析方法および解析プログラム |
KR20170110126A (ko) * | 2015-03-05 | 2017-10-10 | 가부시키가이샤 고베 세이코쇼 | 잔류 응력 추정 방법 및 잔류 응력 추정 장치 |
KR102060901B1 (ko) | 2018-08-14 | 2019-12-30 | 인포겟시스템 주식회사 | 고유 변형도 테이블을 이용한 메시 데이터 생성 방법 |
WO2023195454A1 (ja) * | 2022-04-05 | 2023-10-12 | 公立大学法人大阪 | 変形又は残留応力の予測方法及びプログラム |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5947565B2 (ja) * | 2012-03-02 | 2016-07-06 | 株式会社Jsol | 固有変形データの計算システム、計算プログラム及び計算方法 |
US9098648B2 (en) * | 2012-06-04 | 2015-08-04 | Dassault Systemes Simulia Corp. | Numerical modeling of lineweld connections |
CN103177149A (zh) * | 2012-10-29 | 2013-06-26 | 中国二十二冶集团有限公司 | 钢结构焊接变形分析方法 |
GB201222883D0 (en) | 2012-12-19 | 2013-01-30 | Airbus Operations Ltd | A method, apparatus and computer program product for modelling the non-linear structural response of a component |
CN108039093A (zh) * | 2013-03-11 | 2018-05-15 | 林肯环球股份有限公司 | 使用虚拟现实焊接***导入和分析外部数据 |
CN104057210B (zh) * | 2013-03-20 | 2017-02-08 | 上海通用汽车有限公司 | 基于强度理论的焊接设计方法 |
CN106662439B (zh) * | 2014-06-05 | 2019-04-09 | 联邦科学和工业研究组织 | 增材制造中变形的预测及最小化 |
CN104036150B (zh) * | 2014-06-27 | 2017-06-23 | 北京航空航天大学 | 一种在欧拉坐标系下计算一维理想弹塑性固体的技术 |
JP6345618B2 (ja) * | 2015-03-05 | 2018-06-20 | 株式会社神戸製鋼所 | 残留応力推定方法及び残留応力推定装置 |
JP6345617B2 (ja) * | 2015-03-05 | 2018-06-20 | 株式会社神戸製鋼所 | 残留応力推定方法及び残留応力推定装置 |
CN107330181A (zh) * | 2017-06-28 | 2017-11-07 | 上海交通大学 | 预测激光焊接变形量的实现方法 |
US20200331102A1 (en) * | 2017-09-08 | 2020-10-22 | University Public Corporation Osaka | Method and Apparatus for Analyzing Additively Manufactured Object, and Method and Apparatus for Additively Manufacturing an Object |
CN108460190B (zh) * | 2018-02-06 | 2022-03-04 | 北京卫星制造厂有限公司 | 一种预测蒙皮法兰结构焊缝预置应力焊接的方法 |
CN109299554A (zh) * | 2018-09-30 | 2019-02-01 | 西南石油大学 | 一种激光-电弧复合焊接热源模型的建立方法 |
CN109543333B (zh) * | 2018-12-02 | 2022-02-22 | 湖南大学 | 核电工程用六角管的焊接残余应力与变形的预测方法 |
CN111375930A (zh) * | 2018-12-27 | 2020-07-07 | 中车唐山机车车辆有限公司 | 一种焊接变形的控制方法 |
CN110598358A (zh) * | 2019-09-26 | 2019-12-20 | 华中科技大学 | 一种增材制造应力变形仿真方法、装置、设备及存储介质 |
CN110717291B (zh) * | 2019-09-26 | 2021-05-04 | 华中科技大学 | 一种焊接结构件变形仿真方法、装置、设备及存储介质 |
CN111307055B (zh) * | 2020-03-03 | 2024-03-05 | 成都管力安科技有限公司 | 一种管道数字孪生体系的设计方法 |
CN112784338A (zh) * | 2021-01-19 | 2021-05-11 | 上海跃影科技有限公司 | 一种模型尺寸构建方法及*** |
CN114619161B (zh) * | 2022-02-16 | 2023-04-18 | 江苏科技大学 | 一种薄板焊接变形的模型建构及其矫平方法 |
CN118106642A (zh) * | 2024-04-26 | 2024-05-31 | 华中科技大学 | 基于焊接纵向收缩力评估的焊接变形预测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1242211A1 (en) * | 1999-10-27 | 2002-09-25 | Caterpillar Inc. | Method and apparatus for providing a simulation of a welding process using integrated models |
JP4079638B2 (ja) | 2001-12-26 | 2008-04-23 | 株式会社東芝 | 有限要素法による残留応力解析方法 |
JP2004330212A (ja) | 2003-04-30 | 2004-11-25 | Toshiba Corp | 溶接構造物の解析方法および溶接構造物の解析装置 |
US7734450B2 (en) | 2004-03-29 | 2010-06-08 | Osaka Industrial Promotion Organization | Welding deformation computing method, welding deformation computing device, and computer program product |
JP5241991B2 (ja) | 2004-06-16 | 2013-07-17 | 川崎重工業株式会社 | 変形推定方法、プログラムおよび記録媒体 |
JP2009036669A (ja) | 2007-08-02 | 2009-02-19 | Toshiba Corp | 溶接残留応力解析方法および溶接残留応力解析システム |
CN101598654B (zh) * | 2008-11-21 | 2010-12-01 | 西北工业大学 | 一种摩擦焊接头塑性变形的物理模拟方法 |
-
2010
- 2010-02-03 JP JP2010022232A patent/JP2011159213A/ja active Pending
-
2011
- 2011-01-31 CN CN2011100366048A patent/CN102152016A/zh active Pending
- 2011-02-02 EP EP11152998A patent/EP2363819A1/en not_active Withdrawn
- 2011-02-03 US US13/020,075 patent/US20110213594A1/en not_active Abandoned
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101291058B1 (ko) | 2012-01-18 | 2013-08-07 | 삼성중공업 주식회사 | 변형예측시스템 및 변형예측방법 |
KR20170110126A (ko) * | 2015-03-05 | 2017-10-10 | 가부시키가이샤 고베 세이코쇼 | 잔류 응력 추정 방법 및 잔류 응력 추정 장치 |
KR101941811B1 (ko) | 2015-03-05 | 2019-01-23 | 가부시키가이샤 고베 세이코쇼 | 잔류 응력 추정 방법 및 잔류 응력 추정 장치 |
JP2017161981A (ja) * | 2016-03-07 | 2017-09-14 | 株式会社東芝 | 解析装置、解析方法および解析プログラム |
KR102060901B1 (ko) | 2018-08-14 | 2019-12-30 | 인포겟시스템 주식회사 | 고유 변형도 테이블을 이용한 메시 데이터 생성 방법 |
WO2023195454A1 (ja) * | 2022-04-05 | 2023-10-12 | 公立大学法人大阪 | 変形又は残留応力の予測方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
US20110213594A1 (en) | 2011-09-01 |
EP2363819A1 (en) | 2011-09-07 |
CN102152016A (zh) | 2011-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2011159213A (ja) | 溶接変形解析方法 | |
Samper et al. | Form defects tolerancing by natural modes analysis | |
Rai et al. | Finite element method based machining simulation environment for analyzing part errors induced during milling of thin-walled components | |
Navangul et al. | Error minimization in layered manufacturing parts by stereolithography file modification using a vertex translation algorithm | |
Ratchev et al. | FEA-based methodology for the prediction of part–fixture behaviour and its applications | |
EP3259691B1 (en) | Additive manufacturing method involving computer-aided simulation of additive manufacturing processes | |
EP3118817B1 (en) | Post-processing system for finite element analysis | |
US11947333B2 (en) | Dual lattice representation for crash simulation and manufacturing | |
CN103455651A (zh) | 使用有限元分析模型数值模拟具有热影响区的结构 | |
WO2020166346A1 (ja) | 製造監視支援装置、製造監視支援方法及び製造監視支援プログラム | |
Zongo et al. | Geometric deviations of laser powder bed–fused AlSi10Mg components: numerical predictions versus experimental measurements | |
Ulu et al. | Manufacturability oriented model correction and build direction optimization for additive manufacturing | |
Abarca et al. | Multi-scale modeling for prediction of residual stress and distortion in Ti–6Al–4V semi-circular thin-walled parts additively manufactured by laser powder bed fusion (LPBF) | |
US20210110087A1 (en) | Analysis mesh generation method, recording medium, and analysis mesh generation device | |
CN113010960B (zh) | 零部件制造控制数据生成方法、装置和计算机设备 | |
JP4809683B2 (ja) | 解析モデル生成装置,解析モデル生成方法および解析モデル生成プログラム | |
JP2007102595A (ja) | 解析メッシュ生成装置 | |
JP5343052B2 (ja) | 構造物解析方法、プログラムおよび解析装置 | |
Gopalakrishnan et al. | A digital engineering framework to facilitate automated data exchange between geometric inspection and structural analysis | |
JP2010176573A (ja) | 金型設計装置およびその方法 | |
JP7166102B2 (ja) | 造形データ作成装置、造形データ作成装置用プログラム、造形データ作成方法および造形方法 | |
JP7474382B2 (ja) | 点ごとの重ね合わせ手順に基づくスケーリング方法及びそのシステム | |
JP2002207777A (ja) | 中立面モデルの生成方法 | |
Kukreja et al. | Voxel-based adaptive toolpath planning using graphics processing unit for freeform surface machining | |
KR20100049729A (ko) | 선박 건조시 모델 정보를 이용한 고정정반의 조립품 자동배치 방법 |