(57)【要約】【課題】構造体の溶接変形を比較的簡素な計算で短時間に精度良く計算する。【解決手段】溶接変形の計算システムは、構造体の初期形状、初期応力、剛性、溶接線に対応する応力解放領域、応力解放領域における応力解放率、及び溶接線に対応する固有変形データをそれぞれ取得する初期形状取得手段、初期応力取得手段、剛性取得手段、応力解放領域取得手段、応力解放率取得手段、及び固有変形データ取得手段と、構造体の初期形状の応力解放領域に初期応力と剛性をそれぞれ応力解放率に基づいて補正した補正初期応力と補正剛性を適用し、構造体の初期形状の応力解放領域を除く領域に初期応力と剛性とを適用し、弾性ＦＥＭ解析によって構造体の応力解放変形を算出する応力解放変形算出手段と、算出した構造体の応力解放変形後の解析モデルに固有変形データを適用し、弾性ＦＥＭ解析によって構造体の溶接変形を算出する溶接変形算出手段とを有する。【選択図】図９
2014160008.tif 000002
【特許請求の範囲】
【請求項１】
複数の板状部材を溶接してなる構造体の溶接変形を弾性ＦＥＭ解析によって計算する溶接変形の計算システムであって、
前記構造体の初期形状を取得する初期形状取得手段と、
前記構造体の初期応力を取得する初期応力取得手段と、
前記構造体の剛性を取得する剛性取得手段と、
前記構造体の溶接線に対応する応力解放領域を取得する応力解放領域取得手段と、
前記構造体の応力解放領域における応力解放率を取得する応力解放率取得手段と、
前記構造体の溶接線に対応する固有変形データを取得する固有変形データ取得手段と、
前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出する応力解放変形算出手段と、
前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出する溶接変形算出手段と、
を有することを特徴とする溶接変形の計算システム。
【請求項２】
前記応力解放変形算出手段は、前記応力解放領域に、前記初期応力取得手段によって取得した初期応力をゼロに補正した補正初期応力と前記剛性取得手段によって取得した剛性をゼロに補正した補正剛性とを適用し、前記応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出する、
ことを特徴とする請求項１に記載の溶接変形の計算システム。
【請求項３】
前記初期形状取得手段は、弾性ＦＥＭ解析を用いて算出した構造体の初期形状を取得し、
前記初期応力取得手段は、弾性ＦＥＭ解析を用いて算出した構造体の初期応力を取得する、
ことを特徴とする請求項１又は請求項２に記載の溶接変形の計算システム。
【請求項４】
前記構造体が複数の溶接線を有するときに、
前記応力解放領域取得手段は、前記構造体の各溶接線に対応する応力解放領域を取得し、
前記応力解放率取得手段は、前記構造体の各応力解放領域における応力解放率を取得し、
前記固有変形データ取得手段は、前記構造体の各溶接線に対応する固有変形データを取得し、
前記応力解放変形算出手段は、
最初の溶接線における応力解放変形算出時には、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、
最初の溶接線以降の溶接線における応力解放変形算出時には、当該溶接線の前の溶接線について算出した前記構造体の溶接変形後の形状及び応力を用い、前記構造体の溶接変形後の形状の解析モデルの応力解放領域に、前記構造体の溶接変形後の応力を応力解放率に基づいて補正した補正応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、前記構造体の溶接変形後の応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、
前記溶接変形算出手段は、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出する、
ことを特徴とする請求項１から請求項３の何れか１項に記載の溶接変形の計算システム。
【請求項５】
複数の板状部材を溶接してなる構造体の溶接変形を弾性ＦＥＭ解析によって計算する溶接変形の計算プログラムであって、
コンピュータを、前記構造体の初期形状を取得する初期形状取得手段、
前記構造体の初期応力を取得する初期応力取得手段、
前記構造体の剛性を取得する剛性取得手段、
前記構造体の溶接線に対応する応力解放領域を取得する応力解放領域取得手段、
前記構造体の応力解放領域における応力解放率を取得する応力解放率取得手段、
前記構造体の溶接線に対応する固有変形データを取得する固有変形データ取得手段、
前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出する応力解放変形算出手段、及び、
前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出する溶接変形算出手段として機能させる、
ことを特徴とする溶接変形の計算プログラム。
【請求項６】
コンピュータを、前記応力解放変形算出手段として機能させるときは、前記応力解放領域に、前記初期応力取得手段によって取得した初期応力をゼロに補正した補正初期応力と前記剛性取得手段によって取得した剛性をゼロに補正した補正剛性とを適用し、前記応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出するように機能させる、
ことを特徴とする請求項５に記載の溶接変形の計算プログラム。
【請求項７】
コンピュータを、前記初期形状取得手段として機能させるときは、弾性ＦＥＭ解析を用いて算出した構造体の初期形状を取得するように機能させ、前記初期応力取得手段として機能させるときは、弾性ＦＥＭ解析を用いて算出した構造体の初期応力を取得するように機能させる、
ことを特徴とする請求項５又は請求項６に記載の溶接変形の計算プログラム。
【請求項８】
コンピュータを、前記構造体が複数の溶接線を有するときに、
前記応力解放領域取得手段として機能させるときは、前記構造体の各溶接線に対応する応力解放領域を取得するように機能させ、
前記応力解放率取得手段として機能させるときは、前記構造体の各応力解放領域における応力解放率を取得するように機能させ、
前記固有変形データ取得手段として機能させるときは、前記構造体の各溶接線に対応する固有変形データを取得するように機能させ、
前記応力解放変形算出手段として機能させるときは、
最初の溶接線における応力解放変形算出時には、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、
最初の溶接線以降の溶接線における応力解放変形算出時には、当該溶接線の前の溶接線について算出した前記構造体の溶接変形後の形状及び応力を用い、前記構造体の溶接変形後の形状の解析モデルの応力解放領域に、前記構造体の溶接変形後の応力を応力解放率に基づいて補正した補正応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、前記構造体の溶接変形後の応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出するように機能させ、
前記溶接変形算出手段として機能させるときは、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出するように機能させる、
ことを特徴とする請求項５から請求項７の何れか１項に記載の溶接変形の計算プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の板状部材を溶接してなる構造体の溶接変形を有限要素法を用いた解析（ＦＥＭ解析）によって計算する溶接変形の計算システム及び計算プログラムに関し、コンピュータを用いた計算の技術分野に属する。
【背景技術】
【０００２】
複数の板状部材を溶接して構造体を製造する場合、溶接時の熱によって構造体が全体としてどのように変形するかを予め予測しておくことが重要であり、これを、コンピュータを用いて行おうとする試みがなされている。例えば特許文献１には、溶接構造体を基本的な溶接継手の集合体とみなし、各継手における溶接線の固有変形を溶接順序に従って構造体に与えることにより、構造体全体としての溶接変形を予測する方法が開示されている。
【０００３】
このようにして溶接線の固有変形を用いて構造体の溶接変形を予測することが可能であるが、予測した構造体の溶接変形と実際に溶接した構造体の溶接変形とが相違する場合がある。この１つの原因として、構造体を構成する各部材について溶接前の成形プロセスに応じて生じる初期応力の影響が考えられる。
【０００４】
したがって、構造体の初期応力を考慮して構造体の溶接変形を予測することで、構造体の溶接変形を精度良く予測することができるものと考えられる。構造体の初期応力を考慮したものとして、例えば非特許文献１には、プレス加工によって生じた初期応力を考慮して熱弾塑性ＦＥＭ解析を用いて構造体の溶接変形を予測するものが開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平０６−１８０２７１号公報
【非特許文献】
【０００６】
【非特許文献１】Ninshu Ma, etc, ''Simulation of Welding Deformation of Automobile Part Considering Stamping Results using Explicit Solver in LS-DYNA'', Proceedings of the 8th International Welding Symposium organized by Japan Welding Society, (2008), pp161
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、非特許文献１に記載されるように、構造体の初期応力を考慮して熱弾塑性ＦＥＭ解析を用いて構造体の溶接変形を計算する場合、構造体の溶接変形を精度良く予測することが可能であるが、構造体を構成する各部材について熱依存性を有する材料物性データを用意する必要がある。また、熱弾塑性ＦＥＭ解析を用いて構造体の溶接変形を計算する場合、熱依存性を考慮して計算することから計算量が非常に多く計算時間が長くなる。
【０００８】
そこで、本発明は、構造体の溶接変形を比較的簡素な計算で短時間に精度良く計算することができる溶接変形の計算システム及び計算プログラムを提供することを課題とする。
【課題を解決するための手段】
【０００９】
前記課題を解決するため、本発明は、次のように構成したことを特徴とする。
【００１０】
まず、本願の請求項１に記載の発明は、複数の板状部材を溶接してなる構造体の溶接変形を弾性ＦＥＭ解析によって計算する溶接変形の計算システムであって、前記構造体の初期形状を取得する初期形状取得手段と、前記構造体の初期応力を取得する初期応力取得手段と、前記構造体の剛性を取得する剛性取得手段と、前記構造体の溶接線に対応する応力解放領域を取得する応力解放領域取得手段と、前記構造体の応力解放領域における応力解放率を取得する応力解放率取得手段と、前記構造体の溶接線に対応する固有変形データを取得する固有変形データ取得手段と、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出する応力解放変形算出手段と、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出する溶接変形算出手段と、を有することを特徴とする。
【００１１】
また、請求項２に記載の発明は、請求項１に記載の発明において、前記応力解放変形算出手段は、前記応力解放領域に、前記初期応力取得手段によって取得した初期応力をゼロに補正した補正初期応力と前記剛性取得手段によって取得した剛性をゼロに補正した補正剛性とを適用し、前記応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出することを特徴とする。
【００１２】
また、請求項３に記載の発明は、請求項１又は請求項２に記載の発明において、前記初期形状取得手段は、弾性ＦＥＭ解析を用いて算出した構造体の初期形状を取得し、前記初期応力取得手段は、弾性ＦＥＭ解析を用いて算出した構造体の初期応力を取得することを特徴とする。
【００１３】
また、請求項４に記載の発明は、請求項１から請求項３の何れか１項に記載の発明において、前記構造体が複数の溶接線を有するときに、前記応力解放領域取得手段は、前記構造体の各溶接線に対応する応力解放領域を取得し、前記応力解放率取得手段は、前記構造体の各応力解放領域における応力解放率を取得し、前記固有変形データ取得手段は、前記構造体の各溶接線に対応する固有変形データを取得し、前記応力解放変形算出手段は、最初の溶接線における応力解放変形算出時には、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、最初の溶接線以降の溶接線における応力解放変形算出時には、当該溶接線の前の溶接線について算出した前記構造体の溶接変形後の形状及び応力を用い、前記構造体の溶接変形後の形状の解析モデルの応力解放領域に、前記構造体の溶接変形後の応力を応力解放率に基づいて補正した補正応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、前記構造体の溶接変形後の応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、前記溶接変形算出手段は、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出することを特徴とする。
【００１４】
さらに、請求項５に記載の発明は、複数の板状部材を溶接してなる構造体の溶接変形を弾性ＦＥＭ解析によって計算する溶接変形の計算プログラムであって、コンピュータを、前記構造体の初期形状を取得する初期形状取得手段、前記構造体の初期応力を取得する初期応力取得手段、前記構造体の剛性を取得する剛性取得手段、前記構造体の溶接線に対応する応力解放領域を取得する応力解放領域取得手段、前記構造体の応力解放領域における応力解放率を取得する応力解放率取得手段、前記構造体の溶接線に対応する固有変形データを取得する固有変形データ取得手段、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出する応力解放変形算出手段、及び、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出する溶接変形算出手段として機能させることを特徴とする。
【００１５】
また、請求項６に記載の発明は、請求項５に記載の発明において、コンピュータを、前記応力解放変形算出手段として機能させるときは、前記応力解放領域に、前記初期応力取得手段によって取得した初期応力をゼロに補正した補正初期応力と前記剛性取得手段によって取得した剛性をゼロに補正した補正剛性とを適用し、前記応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出するように機能させることを特徴とする。
【００１６】
また、請求項７に記載の発明は、請求項５又は請求項６に記載の発明において、コンピュータを、前記初期形状取得手段として機能させるときは、弾性ＦＥＭ解析を用いて算出した構造体の初期形状を取得するように機能させ、前記初期応力取得手段として機能させるときは、弾性ＦＥＭ解析を用いて算出した構造体の初期応力を取得するように機能させることを特徴とする。
【００１７】
また、請求項８に記載の発明は、請求項５から請求項７の何れか１項に記載の発明において、コンピュータを、前記構造体が複数の溶接線を有するときに、前記応力解放領域取得手段として機能させるときは、前記構造体の各溶接線に対応する応力解放領域を取得するように機能させ、前記応力解放率取得手段として機能させるときは、前記構造体の各応力解放領域における応力解放率を取得するように機能させ、前記固有変形データ取得手段として機能させるときは、前記構造体の各溶接線に対応する固有変形データを取得するように機能させ、前記応力解放変形算出手段として機能させるときは、最初の溶接線における応力解放変形算出時には、前記構造体の初期形状の解析モデルの応力解放領域に、前記初期応力取得手段によって取得した初期応力を応力解放率に基づいて補正した補正初期応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の初期形状の解析モデルの応力解放領域を除く領域に、前記初期応力取得手段によって取得した初期応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出し、最初の溶接線以降の溶接線における応力解放変形算出時には、当該溶接線の前の溶接線について算出した前記構造体の溶接変形後の形状及び応力を用い、前記構造体の溶接変形後の形状の解析モデルの応力解放領域に、前記構造体の溶接変形後の応力を応力解放率に基づいて補正した補正応力と前記剛性取得手段によって取得した剛性を応力解放率に基づいて補正した補正剛性とを適用し、前記構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、前記構造体の溶接変形後の応力と前記剛性取得手段によって取得した剛性とを適用し、弾性ＦＥＭ解析によって前記構造体の応力解放変形を算出するように機能させ、前記溶接変形算出手段として機能させるときは、前記応力解放変形算出手段によって算出した前記構造体の応力解放変形後の解析モデルに、固有変形データを適用し、弾性ＦＥＭ解析によって前記構造体の溶接変形を算出するように機能させることを特徴とする。
【発明の効果】
【００１８】
以上の構成により、本願各請求項の発明によれば、次の効果が得られる。
【００１９】
まず、本願の請求項１に記載の発明によれば、構造体の初期形状、初期応力、剛性、溶接線に対応する応力解放領域、応力解放領域における応力解放率、及び溶接線に対応する固有変形データが取得され、構造体の初期形状の応力解放領域に初期応力と剛性を応力解放率に基づいて補正した補正初期応力と補正剛性を適用すると共に応力解放領域を除く領域に初期応力と剛性を適用して、弾性ＦＥＭ解析によって構造体の応力解放変形が算出され、算出した構造体の応力解放変形後の解析モデルに固有変形データを適用して、弾性ＦＥＭ解析によって構造体の溶接変形が算出される。
【００２０】
これにより、構造体の初期応力を考慮して構造体の溶接変形を計算することができるので、構造体の溶接変形を精度良く計算することができる。また、熱弾塑性ＦＥＭ解析を用いることなく、固有変形データを用いた弾性ＦＥＭ解析によって構造体の溶接変形を計算することができるので、熱弾塑性ＦＥＭ解析を用いた場合に比して比較的簡素な計算で短時間に計算することができる。加えて、構造体の溶接時における溶接部の状態に応じて応力解放率を適宜設定することができ、構造体の溶接変形をより精度良く計算することができる。
【００２１】
また、請求項２に記載の発明によれば、構造体の初期形状の応力解放領域に初期応力をゼロに補正した補正初期応力と剛性をゼロに補正した補正剛性とを適用し、応力解放領域を除く領域に初期応力と剛性とを適用し、弾性ＦＥＭ解析によって構造体の応力解放変形が算出されることにより、構造体の溶接時に、溶接部において初期応力及び剛性がゼロとなる場合について構造体の溶接変形を精度良く短時間に計算することができる。
【００２２】
また、請求項３に記載の発明によれば、弾性ＦＥＭ解析を用いて算出した構造体の初期形状と初期応力が取得されることにより、構造体を構成する各板状部材についてプレス成形シミュレーションなどの弾性ＦＥＭ解析を用いたシミュレーションを行う場合、そのシミュレーション結果を用いることができるので、構造体の初期形状及び初期応力を容易に取得することができ、前記効果をより有効に奏することができる。
【００２３】
また、請求項４に記載の発明によれば、構造体が複数の溶接線を有するとき、最初の溶接線における応力解放変形算出時には、構造体の初期形状の応力解放領域に初期応力と剛性を応力解放率に基づいて補正した補正初期応力と補正剛性とを適用し、応力解放領域を除く領域に初期応力と剛性とを適用し、最初の溶接線以降の溶接線における応力解放変形算出時には、前の溶接線について算出した構造体の溶接変形後の形状及び応力を用い、構造体の溶接変形後の形状の応力解放領域に、構造体の溶接変形後の応力と剛性を応力解放率に基づいて補正した補正応力と補正剛性とを適用し、応力解放領域を除く領域に、構造体の溶接変形後の応力と剛性とを適用し、弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。これにより、複数の溶接線を有する構造体の溶接変形を予測する際に前の溶接線についての溶接変形結果を用いて溶接線ごとに応力解放変形解析と溶接変形解析を順次行うことができ、構造体の溶接変形を精度良く短時間に計算することができる。
【００２４】
そして、請求項５〜８に記載の溶接変形の計算プログラムに関する発明によれば、これをコンピュータで実行することにより、溶接変形の計算システムに関する請求項１〜４に記載の発明と同様の効果を奏することができる。
【図面の簡単な説明】
【００２５】
【図１】本発明の第１実施形態に係る溶接変形の計算を説明するための説明図である。
【図２】応力解放変形の計算を説明するための説明図である。
【図３】固有変形を説明するための説明図である。
【図４】本発明の第１実施形態に係る溶接変形の計算に用いる構造物を説明するための説明図である。
【図５】前記構造物の板厚方向における外面及び内面の初期応力を示す図である。
【図６】前記構造物の中心点について板厚方向における初期応力を示すグラフである。
【図７】本発明の第１実施形態に係るシステムの全体構成を示す図である。
【図８】図７に示す記憶装置の構成を示す図である。
【図９】前記システムにおいて溶接変形を計算する動作を示すフローチャートである。
【図１０】応力解放変形解析の動作を示すフローチャートである。
【図１１】剛性データ入力画面を示す図である。
【図１２】応力解放領域・応力解放率データ入力画面を示す図である。
【図１３】固有変形データ入力画面を示す図である。
【図１４】初期応力を考慮した場合に弾性ＦＥＭ解析による変形図である。
【図１５】初期応力を考慮しない場合に弾性ＦＥＭ解析による変形図である。
【図１６】熱弾塑性ＦＥＭ解析に用いた材料物性データを示すグラフである。
【図１７】初期応力を考慮した場合に熱弾塑性ＦＥＭ解析による変形図である。
【図１８】初期応力を考慮しない場合に熱弾塑性ＦＥＭ解析による変形図である。
【図１９】応力解放領域・応力解放率データ入力画面の別の画面を示す図である。
【図２０】本発明の第２実施形態に係る溶接変形の計算を説明するための説明図である。
【図２１】本発明の第２実施形態に係るシステムにおいて溶接変形を計算する動作を示すフローチャートである。
【発明を実施するための形態】
【００２６】
以下、本発明の実施形態について添付図面を参照しながら説明する。
図１は、本発明の第１実施形態に係る溶接変形の計算を説明するための説明図である。図１に示すように、本発明の第１実施形態に係る溶接変形の計算では、複数の板状部材を溶接してなる構造体に生じる溶接時の変形を、溶接時に入力される熱による溶接部の軟化に伴って溶接部の初期応力が解放されて構造体が変形することを模擬した応力解放変形と、応力解放変形後の構造体に適用する固有変形を用いた溶接変形とを合わせたものとして計算する。
【００２７】
同一形状及び同一物性を有する板状の第１部材Ｍ１と第２部材Ｍ２とを突合せ溶接によって接合して構造体１を形成する場合、図１（ａ）に示すように、第１部材Ｍ１と第２部材Ｍ２とを突き合わせた構造体１ａに溶接線Ｌ１で突合せ溶接すると、構造体１ａから構造体１ｂに溶接時に変形する。
【００２８】
この溶接時の変形を、図１（ｂ）に示すように、溶接線Ｌ１に対応する溶接部を模擬した応力解放領域Ｒ１の初期応力が解放されて構造体１ａから構造体１ｃに変形する応力解放変形と、図１（ｃ）に示すように、応力解放変形後の構造体１ｃに対して溶接線Ｌ１で突合せ溶接することにより構造体１ｃから構造体１ｂに変形する固有変形を用いた溶接変形とからなるものとして計算する。
【００２９】
図２は、応力解放変形の計算を説明するための説明図である。図２（ａ）に示すように、第１部材Ｍ１と第２部材Ｍ２とを突き合わせた構造体１ａは、溶接前の成形プロセスなどに応じて初期応力σ１２を有している。また、構造体１ａは、剛性、具体的にはヤング率Ｅ１２を有している。なお、図２では、第１部材Ｍ１と第２部材Ｍ２とは、これに限定するものではないが、同一形状及び同一物性を有するものとして示している。
【００３０】
応力解放変形の計算では先ず、図２（ｂ）に示すように、溶接線に対応する溶接部を模擬した応力解放領域Ｒ１の初期応力σ１２及び剛性Ｅ１２を所定の応力解放率に基づいて補正して補正初期応力σ１２´及び補正剛性Ｅ１２´を算出し、応力解放領域Ｒ１の初期応力σ１２及び剛性Ｅ１２を補正初期応力σ１２´及び補正剛性Ｅ１２´に代える。応力解放率として、例えば５０％や１００％などの種々の値を設定することができる。応力解放率して１００％を設定した場合は、補正初期応力σ１２´及び補正剛性Ｅ１２´は共にゼロのときの変形が算出される。
【００３１】
そして、構造体１ａの応力解放領域Ｒ１に補正初期応力σ１２´及び補正剛性Ｅ１２´を適用し、構造体１ａの応力解放領域Ｒ１を除く領域に初期応力σ１２及び剛性Ｅ１２を適用し、既知の弾性ＦＥＭ解析を用いて、図２（ｃ）に示すように、構造体１ａから構造体１ｃに変形する構造体の応力解放変形を算出する。応力解放変形後の構造体１ｃは、算出された応力σ１２´´と剛性Ｅ１２とを有している。
【００３２】
本実施形態では、このようにして算出された応力解放変形後の構造体１ｃに対して、前述したように、溶接線に対応する固有変形を適用し、既知の弾性ＦＥＭ解析を用いて構造体の溶接変形を算出する。
【００３３】
図３は、固有変形を説明するための説明図であり、突合せ溶接時の固有変形を模式的に示している。複数の板状部材を溶接する際には溶接時の熱によって溶接線近傍の所定の領域に固有ひずみが生じることとなるが、固有変形とは、この固有ひずみを固有ひずみが生じる領域で積分したものであり、縦収縮、横収縮、縦曲げ、横曲げの４つの基本変形によって規定したものである。
【００３４】
図３では、第１部材Ｍ１と第２部材Ｍ２を溶接線Ｌ１で突合せ溶接によって接合する場合について固有変形ひずみ領域の溶接線方向に単位長さを有する部分を部分的に取り出して示している。縦収縮は、図３（ａ）に示すように、溶接線方向の収縮であり、横収縮は、図３（ｂ）に示すように、溶接線方向と直交する幅方向の収縮であり、縦曲げは、図３（ｃ）に示すように、溶接線方向の曲げであり、横曲げは、図３（ｄ）に示すように、溶接線方向と直交する方向の曲げである。
【００３５】
以下、本発明の実施形態に係る溶接変形の計算について説明する。
図４は、本発明の第１実施形態に係る溶接変形の計算に用いる構造物を説明するための説明図である。図４に示すように、本実施形態に係る溶接変形の計算では、構造体２として、半円筒状に成形された板状部材Ｍ１１を用いた。構造体２について、板状部材Ｍ１１の外面に溶接線Ｌ１１で板状部材Ｍ１２が溶接されることを模擬して、板状部材Ｍ１１の外面に溶接線Ｌ１１の溶接を行うときの溶接変形を計算した。
【００３６】
構造体２、具体的には板状部材Ｍ１１は初期応力を有するものを用いた。構造体２の初期形状及び初期応力としては、プレス成形金型を用いて鋼板（ブランク）を半円筒状に成形するプレス成形を既知の弾性ＦＥＭ解析を用いてプレス成形シミュレーションして算出したプレス成形後の形状及び応力を用いた。
【００３７】
図５は、前記構造物の板厚方向における外面及び内面の初期応力を示す図である。図５（ａ）は、構造物２としての板状部材Ｍ１１の外面の初期応力を示し、図５（ｂ）は、板状部材Ｍ１１の内面の初期応力を示している。図５では、斜線ハッチングを用いて初期応力を表示している。
【００３８】
また、図６は、前記構造物の中心点について板厚方向における初期応力を示すグラフである。図６では、図４に示す構造物２の中心点Ｐ１について板厚方向における初期応力を実線で示し、構造体２をプレス成形金型を用いて成形するときにプレス成形金型の下死点における構造物２の中心点Ｐ１についての板厚方向における応力を二点鎖線で示している。なお、図５及び図６に示す構造体２の初期応力は、プレス成形シミュレーションを用いて算出したものであるため、図６では中心点Ｐ１の板厚方向位置を積分点を用いて示している。
【００３９】
溶接変形の計算に用いる構造体２は、図５及び図６に示すような初期応力を有し、図５（ａ）に示すように、板状部材Ｍ１１の外面では初期応力として圧縮応力を有し、板状部材Ｍ１１の内面では初期応力として引張応力を有している。また、構造体２の中心点Ｐ１では、外面において圧縮応力を有し、内面において引張方向を有し、図６の実線に示す初期応力を有している。本実施形態では、このように初期応力を有する構造体の溶接変形を計算する。
【００４０】
図７は、本発明の第１実施形態に係るシステムの全体構成を示す図である。図７に示すように、本発明の第１実施形態に係るシステムは、コンピュータ１０を中心として構成され、コンピュータ１０は、中央処理装置１１と、溶接変形の計算に必要なデータなどを入力するためのキーボードなどの入力装置１２と、溶接変形の計算結果などを表示するためのディスプレイなどの表示装置１３と、溶接変形を計算するためのプログラムなどを記憶するメモリなどの記憶装置１４と、溶接変形の計算結果などを出力するプリンタなどの出力装置１５とを有し、固有変形データを保存する固有変形データベース２０に接続されている。
【００４１】
中央処理装置１１は、入力装置１２、表示装置１３及び出力装置１５を制御すると共に、記憶装置１４及び固有変形データベース２０にアクセス可能に構成されている。中央処理装置１１は、入力装置１２を介して入力された情報と、記憶装置１４に記録されているプログラムやデータと、固有データベース２０に保存されている固有変形データとを用いて、溶接変形の計算をすると共に、溶接変形の計算結果を記憶装置１４に保存するように構成されている。
【００４２】
図８は、図７に示す記憶装置の構成を示す図である。図８に示すように、記憶装置１４は、プログラム記憶部とデータ記憶部とを有しており、プログラム記憶部には、構造体の形状データを有限要素分割して解析モデルを作成するための解析モデル作成プログラムと、構造体の解析モデルの応力解放領域に構造体の初期応力及び剛性を応力解放率に基づいて補正した補正初期応力及び補正剛性を適用すると共に構造体の解析モデルの応力解放領域以外の領域に構造体の初期応力及び剛性を適用し、既知の弾性ＦＥＭ解析によって構造体の応力解放変形を計算するための応力解放変形計算プログラムと、構造体の応力解放変形後の解析モデルに固有変形データを適用し、既知の弾性ＦＥＭ解析によって構造体の溶接変形を計算するための溶接変形計算プログラムと、入力画面や解析結果などを表示するための表示プログラムとが記憶されている。
【００４３】
一方、データ記憶部には、構造体の初期形状が記録される形状データファイルと、構造体の初期応力が記録される応力データファイルと、構造体の剛性が記録される剛性データファイルと、構造体の溶接線に応じて設定された応力解放領域及び該応力解放領域における応力解放率が記録される応力解放領域・応力解放率データファイルと、構造体に対する応力解放変形の計算結果が記録される応力解放変形結果ファイルと、応力解放変形後の構造体に対する固有変形を用いた溶接変形の計算結果が記録される溶接変形結果ファイルとが備えられている。
【００４４】
図９は、前記システムにおいて溶接変形を計算する動作を示すフローチャートである。構造体の溶接変形を計算する前に、コンピュータ１０には先ず、ユーザによって入力装置１２を介して構造体の初期形状の形状データ、及び構造体の初期応力データが登録され、データ記憶部の形状データファイル、及び応力データファイルに登録される。
【００４５】
構造体の初期形状及び初期応力として、図５及び図６に示すプレス成形シミュレーションの計算結果のように、弾性ＦＥＭ解析を用いて算出した構造体の初期形状及び初期応力を用いることができる。また、構造体の初期形状及び初期応力として、実際に構造体の形状及び応力を測定した実験結果を用いることも可能である。
【００４６】
また、コンピュータ１０には、ユーザによって入力装置１２を介して構造体の剛性データ、構造体の溶接線に対応する応力解放領域及び応力解放領域における応力解放率が登録され、データ記憶部の剛性データファイル、応力解放領域・応力解放率データファイルに登録される。
【００４７】
図１１は、剛性データ入力画面を示す図である。図１１に示すように、構造体の剛性データを入力する際には、表示装置１３に剛性データ入力画面３０が表示され、該画面３０には、溶接線番号入力部３１と剛性データ入力部３２とが設けられると共に、登録ボタン３３が設けられている。
【００４８】
溶接線番号入力部３１には、構造体の溶接線番号を入力することができるようになっている。剛性データ入力部３２には、剛性データとしてのヤング率を入力することができるようになっている。そして、剛性データ入力画面３０において、剛性データが入力された後に、ユーザによって登録ボタン３３が押されると、剛性データが、データ記憶部の剛性データファイルに記録される。
【００４９】
図１２は、応力解放領域・応力解放率データ入力画面を示す図である。図１２に示すように、構造体の溶接線に対応する応力解放領域と応力解放領域における応力解放率を入力する際には、表示装置１３に応力解放領域・応力解放率データ入力画面４０が表示され、該画面４０には、溶接線番号入力部４１と応力解放領域・応力解放率データ入力部４２とが設けられると共に、登録ボタン４３が設けられている。
【００５０】
溶接線番号入力部４１には、構造体の溶接線番号を入力することができるようになっている。応力解放領域・応力解放率データ入力部４２には、構造体の溶接線に対応する初期応力を解放させる応力解放領域の幅と深さを入力することができると共に、構造体の応力解放領域における初期応力を解放させる応力解放率を入力することができるようになっている。
【００５１】
図１２に示す画面４０では、応力解放領域として幅２０ｍｍ、深さ１．６ｍｍである断面矩形状の応力解放領域を設定し、応力解放率として１００％、すなわち初期応力をすべて解放させて補正初期応力をゼロにするように設定したものが示されている。応力解放領域としては、断面矩形状に限らず、種々の断面形状を設定することが可能である。応力解放率についても、種々の応力解放率を設定することが可能である。
【００５２】
応力解放領域・応力解放率データ入力画面４０において、応力解放領域データと応力解放率データが入力された後に、ユーザによって登録ボタン４３が押されると、応力解放領域データ及び応力解放率データが、データ記憶部の応力解放領域・応力解放率データファイルに記録される。
【００５３】
また、コンピュータ１０に接続される固有変形データベース２０には、予め種々の溶接線に対応する固有変形データが保存されているが、固有変形データベース２０に保存されていない溶接線に対応する固有変形データを記録する場合は、ユーザによって入力装置１２を介して構造体の溶接線に対応する固有変形データが登録され、固有変形データベース１２に登録される。
【００５４】
図１３は、固有変形データ入力画面を示す図である。構造体の溶接線に対応する固有変形データを入力する際には、表示装置１３に固有変形データ入力画面５０が表示され、該画面５０には、溶接線入力部５１と材料物性入力部５２と固有変形入力部５３とが設けられると共に、登録ボタン５４が設けられている。
【００５５】
溶接線入力部５１には、アーク溶接などの溶接方法を入力することができると共に突合せ溶接などの溶接タイプを入力することができるようになっている。溶接線入力部５１にはまた、溶接条件としての電流、電圧、速度及び熱効率を入力することができるようになっている。
【００５６】
材料物性入力部５２には、鉄などの材料を入力することができると共に、板厚、密度、比熱及び熱膨張係数を入力することができるようになっている。また、固有変形入力部５３には、構造体の溶接線に対応する固有変形データ、具体的には縦収縮量、横収縮量、縦曲げ量、横曲げ量を入力することができるようになっている。そして、固有変形データ入力画面５０において、溶接線、材料物性及び固有変形のデータが入力された後に、ユーザによって登録ボタン５４が押されると、溶接線に対応する固有変形が固有変形データベース２０に記録される。
【００５７】
このようにして、構造体の初期形状、初期応力、剛性、溶接線に対応する応力解放領域、応力解放領域における応力解放率、及び溶接線に対応する固有変形データが登録された状態で、図９に示すように、構造体の溶接変形の計算が行われる。
【００５８】
構造体の溶接変形を計算する際には先ず、構造体の溶接変形を計算するための各種データが取得される。具体的には、構造体の初期形状データが取得され（ステップＳ１）、構造体の初期応力データが取得され（ステップＳ２）、構造体の剛性データが取得され（ステップＳ３）、構造体の溶接線に対応する応力解放領域データ及び応力解放領域における応力解放率データが取得され（ステップＳ４）、構造体の溶接線に対応する固有変形データが取得される（ステップＳ５）。そして、構造体の応力解放変形解析が行われる（ステップＳ６）。
【００５９】
図１０は、応力解放変形解析の動作を示すフローチャートである。ステップＳ６における構造体の応力解放変形解析では、図１０に示すように、ステップＳ４において取得した構造体の応力解放領域についてステップＳ２において取得した構造体の初期応力をステップＳ４において取得した応力解放率に基づいて補正して補正初期応力が算出されると共に、ステップＳ４において取得した構造体の応力解放領域についてステップＳ３において取得した構造体の剛性をステップＳ４において取得した応力解放率に基づいて補正して補正剛性が算出される（ステップＳ１１）。
【００６０】
図１２に示すように、応力解放領域における応力解放率として１００％が設定されている場合、ステップＳ１１では、構造体の応力解放領域における初期応力をゼロに補正した補正初期応力が算出され、構造体の応力解放領域における剛性をゼロに補正した補正剛性が算出される。
【００６１】
そして、ステップＳ１において取得した構造体の初期形状の解析モデルの応力解放領域に、ステップＳ１１において算出した補正初期応力及び補正剛性が適用され、構造体の初期形状の解析モデルの応力解放領域を除く領域に、ステップＳ２及びステップＳ３において取得した構造体の初期応力及び剛性が適用され、既知の弾性ＦＥＭ解析によって構造体の応力解放変形が算出される（ステップＳ１２）。
【００６２】
図１２に示すように、応力解放領域における応力解放率として１００％が設定されている場合、ステップＳ１２では、構造体の初期形状の応力解放領域に初期応力をゼロに補正した補正初期応力と剛性をゼロに補正した補正剛性とが適用され、構造体の初期形状の応力解放領域を除く領域に構造体の初期応力及び剛性が適用され、既知の弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。
【００６３】
ステップＳ１２において構造体の応力解放変形が算出されると、算出した構造体の応力解放変形が応力解放変形結果ファイルに記録される。構造体の応力解放変形の計算結果として、構造体の応力解放変形後の形状及び応力が算出され、これらが応力解放変形結果ファイルに記録される。
【００６４】
ステップＳ６において応力解放変形解析が行われると次に、固有変形を用いた溶接変形解析が行われる（ステップＳ７）。ステップ７では、ステップＳ６において算出した構造体の応力解放変形後の解析モデルに、具体的には応力解放変形後の形状及び応力を有する解析モデルに、ステップＳ５において取得した溶接線に対応する固有変形データが適用され、既知の弾性ＦＥＭ解析によって構造体の溶接変形が算出される。
【００６５】
ステップＳ７において固有変形を用いた構造体の溶接変形が算出されると、算出した構造体の溶接変形の解析結果が、溶接変形結果ファイルに記録されると共に表示装置１３に出力される（ステップＳ８）。溶接変形結果ファイルには、構造体の溶接変形の解析結果として算出した構造体の溶接変形後の形状及び応力が記録される。一方、表示装置１３には、構造体の溶接変形後の形状が構造体の初期形状からのたわみ量と共に表示される。
【００６６】
図１４は、初期応力を考慮した場合に弾性ＦＥＭ解析による変形図である。図１４に示すように、解析結果を表示する際には、表示装置１３に解析結果画面６０が表示され、該画面６０には、解析結果表示部６１が設けられると共に、終了ボタン６２が設けられている。
【００６７】
解析結果表示部６１には、構造体の溶接変形の解析結果が表示され、構造体の溶接変形後の形状が構造体の初期形状からの変形量であるたわみ量と共に表示されるようになっている。構造体のたわみ量は、たわみ量の大小が色の濃淡で表示されると共にたわみ量の最大値が表示される。なお、構造体のたわみ量は、構造体の板厚方向の各積分点におけるたわみ量を平均化して示される。
【００６８】
図１４は、図５及び図６に示す初期応力を有する構造体２について、図４に示す溶接線Ｌ１１で溶接する場合について、図１１に示す剛性データ、図１２に示す応力解放領域及び応力解放率などを適用し、応力解放変形解析及び固有変形を用いた溶接変形解析を行って算出した構造体２の溶接変形の解析結果を示している。
【００６９】
また、構造体２について、応力解放変形解析を行うことなく、固有変形を用いた溶接変形解析のみによる構造体２の溶接変形についても算出した。具体的には、構造体２の初期形状の解析モデルに、溶接線Ｌ１１に対応する固有変形データを適用して既知の弾性ＦＥＭ解析によって構造体の溶接変形を算出した。
【００７０】
図１５は、初期応力を考慮しない場合に弾性ＦＥＭ解析による変形図であり、図１５は、構造体２の初期形状の解析モデルに、溶接線Ｌ１１に対応する固有変形データを適用して既知の弾性ＦＥＭ解析によって構造体の溶接変形を算出した解析結果を示している。図１４及び図１５に示す解析結果から、構造体の初期応力を考慮しない構造体の溶接変形の計算では、構造体の初期応力を考慮した構造体の溶接変形の計算に比してたわみ量が大きく算出されることが分かる。
【００７１】
本実施形態ではまた、本実施形態に係る構造体の溶接変形の計算結果を、熱弾塑性ＦＥＭ解析を用いた構造体の溶接変形の計算結果と比較するために、構造体２について、初期応力を考慮する場合と初期応力を考慮しない場合について熱弾塑性ＦＥＭ解析を用いて構造体の溶接変形を算出した。
【００７２】
熱弾塑性ＦＥＭ解析を用いた構造体２の溶接変形の計算では、溶接時に加えられる熱によって、構造体２における溶接線Ｌ１１について、幅２０ｍｍで板厚方向において外面から三分の二の深さまで１８００度の温度（初期温度）になるものとして、熱弾塑性ＦＥＭ解析を用いて構造体の溶接変形を算出した。
【００７３】
図１６は、熱弾塑性ＦＥＭ解析に用いた材料物性データを示すグラフである。熱弾塑性ＦＥＭ解析を用いた溶接変形の計算では、熱依存性を有する材料物性データとして、図１６に示す材料物性データを用いた。構造体２について、図１６に示すポアソン比、降伏応力、ヤング率、熱膨張係数及び塑性硬化係数を用いた。また、構造体２の初期形状及び初期応力としては、本実施形態に係る構造体２の溶接変形の計算と同様の初期形状及び初期応力を用いた。
【００７４】
図１７は、初期応力を考慮した場合に熱弾塑性ＦＥＭ解析による変形図である。図１８は、初期応力を考慮しない場合に熱弾塑性ＦＥＭ解析による変形図である。
【００７５】
図１７及び図１８についても、構造体の溶接変形後の形状を構造体の初期形状からの変形量であるたわみ量と共に表示している。構造体のたわみ量は、たわみ量の大小が色の濃淡で表示すると共にたわみ量の最大値を表示し、構造体のたわみ量は、構造体の板厚方向の各積分点におけるたわみ量を平均化して示している。
【００７６】
図１７に示すように、初期応力を考慮した場合について熱弾塑性ＦＥＭ解析を用いた溶接変形の解析結果は、図１４に示す溶接変形の解析結果とたわみ量の最大値が略等しく略同様の解析結果が得られた。図１８に示すように、初期応力を考慮しない場合についても熱弾塑性ＦＥＭ解析を用いた溶接変形の解析結果は、図１５に示す溶接変形の解析結果とたわみ量の最大値が略等しく略同様の解析結果が得られた。これら解析結果から分かるように、本実施形態に係る溶接変形の計算では、熱弾塑性ＦＥＭ解析を用いた溶接変形の計算と同様に、初期応力を考慮して構造体の溶接変形を精度良く計算することができる。
【００７７】
このように、本実施形態に係る溶接変形の計算では、構造体の初期形状、初期応力、剛性、溶接線に対応する応力解放領域、応力解放領域における応力解放率、及び溶接線に対応する固有変形データが取得され、構造体の初期形状の応力解放領域に初期応力と剛性を応力解放率に基づいて補正した補正初期応力と補正剛性を適用すると共に応力解放領域を除く領域に初期応力と剛性を適用して、弾性ＦＥＭ解析によって構造体の応力解放変形が算出され、算出した構造体の応力解放変形後の解析モデルに固有変形データを適用して、弾性ＦＥＭ解析によって構造体の溶接変形が算出される。
【００７８】
これにより、構造体の初期応力を考慮して構造体の溶接変形を計算することができるので、構造体の溶接変形を精度良く計算することができる。また、熱弾塑性ＦＥＭ解析を用いることなく、固有変形データを用いた弾性ＦＥＭ解析によって構造体の溶接変形を計算することができるので、熱弾塑性ＦＥＭ解析を用いた場合に比して比較的簡素な計算で短時間に計算することができる。加えて、構造体の溶接時における溶接部の状態に応じて応力解放率を適宜設定することができ、構造体の溶接変形をより精度良く計算することができる。
【００７９】
また、構造体の初期形状と初期応力として、弾性ＦＥＭ解析を用いて算出した構造体の初期形状と初期応力が取得されることにより、構造体を構成する各板状部材についてプレス成形シミュレーションなどの弾性ＦＥＭ解析を用いたシミュレーションを行う場合、そのシミュレーション結果を用いることができるので、構造体の初期形状及び初期応力を容易に取得することができる。
【００８０】
前述した実施形態では、図１２に示すように、応力解放領域として幅と深さを入力する際に、応力解放領域の深さは、寸法を用いて入力するようになっているが、応力解放領域の深さを、深さの割合として入力するようにすることも可能である。
【００８１】
図１９は、応力解放領域・応力解放率データ入力画面の別の画面を示す図である。図１９に示すように、応力解放領域・応力解放率データ入力画面４０における応力解放領域・応力解放率データ入力部４２において、応力解放領域の深さを、深さの割合（％）で入力するようにすることも可能である。
【００８２】
また、前述した実施形態では、構造体の初期形状及び初期応力として、図５及び図６に示すプレス成形シミュレーションの計算結果を用いているが、構造体が複数の溶接線を有する場合、１つの溶接線について溶接変形を計算するために応力解放変形解析及び固有変形を用いた溶接変形解析を行った後に、その解析結果を、次の溶接線についての溶接変形の計算に用いることも可能である。
【００８３】
図２０は、本発明の第２実施形態に係る溶接変形の計算を説明するための説明図である。図２０に示すように、本発明の第２実施形態に係る溶接変形の計算では、板状部材Ｍ２１と板状部材Ｍ２２とを溶接線Ｌ２１で溶接すると共に板状部材Ｍ２１と板状部材Ｍ２３とを溶接線Ｌ２２で溶接して製造される構造体３のように複数の溶接線を有する構造体の溶接変形を計算する。
【００８４】
本発明の第２実施形態に係る溶接変形の計算においても、前述した第１実施形態に係る溶接変形と同様のシステムが用いられるが、第２実施形態に係る溶接変形の計算では、複数の溶接線に対応して、各溶接線に対応する応力解放領域及び各応力解放領域における応力解放率が登録されるようになっている。
【００８５】
図２１は、本発明の第２実施形態に係るシステムにおいて溶接変形を計算する動作を示すフローチャートである。本発明の第２実施形態に係るシステムにおいても、本発明の第１実施形態に係るシステムにおける溶接変形の計算動作と同様に、構造体の溶接変形を計算する際には先ず、構造体の溶接変形を計算するための各種データが取得される。
【００８６】
構造体の初期形状データが取得され（ステップＳ２１）、構造体の初期応力データが取得され（ステップＳ２２）、構造体の剛性データが取得され（ステップＳ２３）、構造体の各溶接線に対応する応力解放領域データ及び応力解放領域における応力解放率データが取得され（ステップＳ２４）、構造体の各溶接線に対応する固有変形データが取得される（ステップＳ２５）。そして、構造体の応力解放変形解析が行われる（ステップＳ２６）。
【００８７】
ステップＳ２６における構造体の応力解放変形解析では、最初の溶接線における応力解放変形算出時には、構造体の初期形状の解析モデルの応力解放領域に、ステップＳ２２において取得した構造体の初期応力をステップ２４において取得した応力解放率に基づいて補正した補正初期応力とステップ２３において取得した構造体の剛性をステップ２４において取得した応力解放率に基づいて補正した補正剛性とが適用され、構造体の初期形状の解析モデルの応力解放領域を除く領域に、ステップＳ２２及びステップＳ２３において取得した構造体の初期応力及び剛性が適用され、既知の弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。
【００８８】
ステップＳ２６において応力解放変形解析が行われると次に、固有変形を用いた溶接変形解析が行われる（ステップＳ２７）。ステップ２７では、ステップＳ２６において算出した構造体の応力解放変形後の解析モデルに、具体的には応力解放変形後の形状及び応力を有する解析モデルに、ステップＳ２５において取得した溶接線に対応する固有変形データが適用され、既知の弾性ＦＥＭ解析によって構造体の溶接変形が算出される。
【００８９】
そして、複数の溶接線のすべてについて溶接変形の解析が終了したか否かが判定される（ステップＳ２８）。ステップＳ２８での判定結果がノー（ＮＯ）の場合、すなわち、全溶接線の解析が終了していない場合は、次の溶接線についてステップＳ２６及びステップＳ２８が繰り返される。
【００９０】
次の溶接線についてのステップＳ２６における構造体の応力解放変形解析では、前の溶接線についてステップＳ２７において算出した溶接変形後の形状及び応力を用い、構造体の溶接変形後の形状の解析モデルの応力解放領域に、構造体の溶接変形後の応力をステップ２４において取得した応力解放率に基づいて補正した補正応力とステップＳ２３において取得した構造体の剛性をステップ２４において取得した応力解放率に基づいて補正した補正剛性とが適用され、構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、構造体の溶接変形後の応力とステップＳ２３において取得した剛性が適用され、既知の弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。
【００９１】
このように、最初の溶接線以降の溶接線における応力解放変形算出時には、当該溶接線の前の溶接線について算出した構造体の溶接変形後の形状及び応力を用い、構造体の溶接変形後の形状の解析モデルの応力解放領域に、構造体の溶接変形後の応力を応力解放率に基づいて補正した補正応力と構造体の剛性を応力解放率に基づいて補正した補正剛性とが適用され、構造体の溶接変形後の形状の解析モデルの応力解放領域を除く領域に、構造体の溶接変形後の応力と構造体の剛性とが適用され、弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。
【００９２】
次の溶接線についても、ステップＳ２６において応力解放変形解析が行われると次に、固有変形を用いた溶接変形解析が行われ（ステップＳ２７）、ステップＳ２６において算出した構造体の応力解放変形後の解析モデルに、具体的には応力解放変形後の形状及び応力を有する解析モデルに、ステップＳ２５において取得した溶接線に対応する固有変形データが適用され、既知の弾性ＦＥＭ解析によって構造体の溶接変形が算出される。
【００９３】
そして、複数の溶接線のすべてについて溶接変形の解析が終了するまで、ステップＳ２６及びステップＳ２７が繰り返され、ステップＳ２８での判定結果がイエス（ＹＥＳ）になると、すなわち、全溶接線の解析が終了すると、算出した構造体の溶接変形の解析結果が、溶接変形結果ファイルに記録されると共に表示装置１３に出力される（ステップＳ２９）。
【００９４】
このように、本発明の第２実施形態に係る溶接変形の計算では、構造体が複数の溶接線を有する場合、最初の溶接線における応力解放変形算出時には、構造体の初期形状の応力解放領域に初期応力と剛性を応力解放率に基づいて補正した補正初期応力と補正剛性とを適用し、応力解放領域を除く領域に初期応力と剛性とが適用され、最初の溶接線以降の溶接線における応力解放変形算出時には、前の溶接線について算出した構造体の溶接変形後の形状及び応力を用い、構造体の溶接変形後の形状の応力解放領域に、構造体の溶接変形後の応力と剛性を応力解放率に基づいて補正した補正応力と補正剛性とが適用され、応力解放領域を除く領域に、構造体の溶接変形後の応力と剛性とを適用し、弾性ＦＥＭ解析によって構造体の応力解放変形が算出される。
【００９５】
これにより、複数の溶接線を有する構造体の溶接変形を予測する際に前の溶接線についての溶接変形結果を用いて溶接線ごとに応力解放変形解析と溶接変形解析を順次行うことができ、構造体の溶接変形を精度良く短時間に計算することができる。
【００９６】
本発明は、例示された実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において、種々の改良及び設計上の変更が可能である。
【産業上の利用可能性】
【００９７】
以上のように、本発明によれば、構造体の溶接変形を比較的簡素な計算で短時間に精度良く計算することができるので、複数の板状部材を溶接してなる構造体の溶接変形を有限要素法を用いて予測する場合に、好適に利用される可能性がある。
【符号の説明】
【００９８】
１０  コンピュータ
１１  中央処理装置
１２  入力装置
１３  表示装置
１４  記憶装置
１５  出力装置
２０  固有変形データベース  
【図１】
2014160008.tif 000003
【図２】
2014160008.tif 000004
【図３】
2014160008.tif 000005
【図４】
2014160008.tif 000006
【図５】
2014160008.tif 000007
【図６】
2014160008.tif 000008
【図７】
2014160008.tif 000009
【図８】
2014160008.tif 000010
【図９】
2014160008.tif 000011
【図１０】
2014160008.tif 000012
【図１１】
2014160008.tif 000013
【図１２】
2014160008.tif 000014
【図１３】
2014160008.tif 000015
【図１４】
2014160008.tif 000016
【図１５】
2014160008.tif 000017
【図１６】
2014160008.tif 000018
【図１７】
2014160008.tif 000019
【図１８】
2014160008.tif 000020
【図１９】
2014160008.tif 000021
【図２０】
2014160008.tif 000022
【図２１】
2014160008.tif 000023