JP4809645B2 - メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム - Google Patents

メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム Download PDF

Info

Publication number
JP4809645B2
JP4809645B2 JP2005248421A JP2005248421A JP4809645B2 JP 4809645 B2 JP4809645 B2 JP 4809645B2 JP 2005248421 A JP2005248421 A JP 2005248421A JP 2005248421 A JP2005248421 A JP 2005248421A JP 4809645 B2 JP4809645 B2 JP 4809645B2
Authority
JP
Japan
Prior art keywords
node
calculating
predetermined value
ratio
line segment
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.)
Active
Application number
JP2005248421A
Other languages
English (en)
Other versions
JP2007065802A (ja
Inventor
優耶 山下
寛世 三村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Japan Research Institute Ltd
Original Assignee
Japan Research Institute Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Japan Research Institute Ltd filed Critical Japan Research Institute Ltd
Priority to JP2005248421A priority Critical patent/JP4809645B2/ja
Publication of JP2007065802A publication Critical patent/JP2007065802A/ja
Application granted granted Critical
Publication of JP4809645B2 publication Critical patent/JP4809645B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、有限要素法を用いて二次元又は三次元の形状モデルを有限要素法を用いて解析する場合に、形状モデルが互いに接する接点又は接線近傍であっても、有限要素解析を行うことが可能なメッシュに分割することができるメッシュ分割方法、該方法を用いた有限要素解析装置、及びコンピュータプログラムに関する。
コンピュータ技術の進展により、より詳細な現象を扱う非線形問題、過渡問題等の応答解析をコンピュータシミュレーションで行うことが実用上可能となってきている。斯かるコンピュータシミュレーションにおいては、様々な物理現象を扱うことができ、適用範囲が広いことから、解析手法として有限要素法（Finite Element Method：以下、ＦＥＭ）が良く用いられている。
有限要素法では、解析対象となる形状モデルをコンピュータへ入力し、有限要素（メッシュ）を生成し、数値解析を行う。有限要素法による解の妥当性は、有限要素の大きさ、品質（辺長比）等により大きく影響されることから、一般に要素寸法を解析対象の材料ごとに、あるいは形状モデルの要素面、要素辺、節点等に対して指定することにより、最適な有限要素解析モデルを作成する。
メッシュの自動生成を行う場合、「要素」を構成する多角形や多面体が「正」形状から外れてアスペクト比が大きい「要素」、すなわち辺長比が非常に大きな形状の要素が生成される。例えば要素が三角形である場合、細長くて尖った形状の三角形、換言すれば３辺のうち１辺が他の２辺に比べて極端に短い三角形が生成される。なお、辺長比とは、要素を構成する辺のうち、最長辺の長さを最短辺の長さで除した比を意味しており、正多角形では１になるのに対し、細長くて尖った多角形では著しく大きくなる。
辺長比が大きな形状の要素が含まれたメッシュを対象として解析計算を行う場合、解が発散して収束しない、あるいは、誤差の大きな解析結果を出力するという問題点があった。
特に複数の形状モデルが接する場合、例えば二次元モデルでは２つの円が接する場合、三次元モデルでは２つの円柱が接する場合等には、接点又は接線の近傍では辺長比が過大なメッシュに分割されやすく、上述したような問題が生じやすい。そこで例えば特許文献１では、メッシュに自動分割する場合に、接点又は接線近傍でメッシュ分割される要素を所定の大きさ以上になるよう制約条件を付与することにより、メッシュ分割の自動化を図っている。
また、特許文献２では、接点又は接線の近傍に生じやすい辺長比が過大なメッシュを解析対象から除外することにより、解を収束させ、出力される解析結果の誤差を抑制している。特に「辺長比が大きくて形状の悪い要素」については、例えば三角形の短辺の両端にある２節点のうちの一方の節点を他方の節点に重ね合わせる等の処理により解析対象から除外している。
特許第３６３２１２７号公報 特開２００４−０５４８６３号公報
しかし、特許文献１のように、メッシュの粗密に関する制約条件を、ユーザが指定することは困難であり、解析結果の誤差を一定範囲内に収束させるためには試行錯誤を繰り返す必要があり、解析結果を出力するまでの時間を短縮することが困難であるという問題点があった。
また、特許文献２のように、接点又は接線の近傍に生じやすい辺長比が過大なメッシュを解析対象から除外するためには、解析用のプログラム自体を変更する必要があり、特に節点の重ね合わせ等により解析対象から除外する方法では、解析対象の形状が大きく変化することから、解析結果に大きな誤差が含まれることになり好ましくない。
本発明は斯かる事情に鑑みてなされたものであり、粗密に関する制約条件を付与することなく、複数の形状モデルが接する接点又は接線近傍においても、精度良く解析結果を出力することができるようメッシュ分割するメッシュ分割方法、該方法を用いた有限要素解析装置、及びコンピュータプログラムを提供することを目的とする。
上記目的を達成するために第１発明に係るメッシュ分割方法は、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、前記有限要素解析装置は、前記プロセッサにより、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出するステップと、曲線上まで移動した点における曲線の曲率半径を算出するステップと、前記移動距離を前記曲率半径で除した比率を算出するステップと、算出した比率が所定値より大きいか否かを判断するステップと、算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、算出した比率が所定値以下であると判断した場合、メッシュ分割を停止するステップとを実行することを特徴とする。
また、第２発明に係るメッシュ分割方法は、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、前記有限要素解析装置は、前記プロセッサにより、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出するステップと、曲面上まで移動した点における曲面の曲率半径を算出するステップと、前記移動距離を前記曲率半径で除した比率を算出するステップと、算出した比率が所定値より大きいか否かを判断するステップと、算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、算出した比率が所定値以下であると判断した場合、メッシュ分割を停止するステップとを実行することを特徴とする。
また、第３発明に係るメッシュ分割方法は、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、前記有限要素解析装置は、前記プロセッサにより、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断するステップと、所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出するステップと、曲線上まで移動した点における曲線の曲率半径を算出するステップと、前記移動距離を前記曲率半径で除した比率を算出するステップと、算出した比率が所定値より大きいか否かを判断するステップと、算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止するステップとを実行することを特徴とする。
また、第４発明に係るメッシュ分割方法は、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、前記有限要素解析装置は、前記プロセッサにより、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断するステップと、所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出するステップと、曲面上まで移動した点における曲面の曲率半径を算出するステップと、前記移動距離を前記曲率半径で除した比率を算出するステップと、算出した比率が所定値より大きいか否かを判断するステップと、算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返するステップと、算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止するステップとを実行することを特徴とする。
また、第５発明に係る有限要素解析装置は、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段と、曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段と、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、算出した比率が所定値より大きいか否かを判断する判断手段と、該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段と、前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段とを備えることを特徴とする。
また、第６発明に係る有限要素解析装置は、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段と、曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段と、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、算出した比率が所定値より大きいか否かを判断する判断手段と、該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段と、前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段とを備えることを特徴とする。
また、第７発明に係る有限要素解析装置は、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段と、算出した線分長が所定値より大きいか否かを判断する第１の判断手段と、該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段と、曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段と、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、算出した比率が所定値より大きいか否かを判断する第２の判断手段と、該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段と、前記第２の判断手段で算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段と、前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段とを備えることを特徴とする。
また、第８発明に係る有限要素解析装置は、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段と、算出した線分長が所定値より大きいか否かを判断する第１の判断手段と、該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段と、曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段と、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、算出した比率が所定値より大きいか否かを判断する第２の判断手段と、該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段と、算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段と、前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段とを備えることを特徴とする。
また、第９発明に係るコンピュータプログラムは、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、前記コンピュータを、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段、曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、算出した比率が所定値より大きいか否かを判断する判断手段、該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動し、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段、及び前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段として機能させることを特徴とする。
また、第１０発明に係るコンピュータプログラムは、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、前記コンピュータを、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段、曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、算出した比率が所定値より大きいか否かを判断する判断手段、該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動し、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段、及び前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段として機能させることを特徴とする。
また、第１１発明に係るコンピュータプログラムは、複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、前記コンピュータを、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段、算出した線分長が所定値より大きいか否かを判断する第１の判断手段、該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段、曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、算出した比率が所定値より大きいか否かを判断する第２の判断手段、該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段、前記第２の判断手段で算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段、及び前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段として機能させることを特徴とする。
また、第１２発明に係るコンピュータプログラムは、複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、前記コンピュータを、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段、算出した線分長が所定値より大きいか否かを判断する第１の判断手段、該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段、曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段、前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、算出した比率が所定値より大きいか否かを判断する第２の判断手段、該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段、算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段、及び前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段として機能させることを特徴とする。
第１発明、第５発明、及び第９発明では、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定し、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出し、曲線上まで移動した点における曲線の曲率半径を算出する。算出した移動距離を曲率半径で除した比率が所定値より大きい場合、第３の節点を曲線上へ移動して、第１、第２、第３の節点を結ぶ有限要素にメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を再算出し、上述した処理を繰り返す。算出した比率が所定値以下である場合、それ以上メッシュ分割をせず、有限要素解析を実行する。これにより、２つの二次元形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑えることが可能となる。
第２発明、第６発明、及び第１０発明では、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定し、第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出し、曲面上まで移動した点における曲面の曲率半径を算出する。算出した移動距離を曲率半径で除した比率が所定値より大きい場合、第３の節点を曲面上へ移動して、第１、第２、第３の節点を結ぶ有限要素にメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を再算出し、上述した処理を繰り返す。算出した比率が所定値以下である場合、それ以上メッシュ分割をせず、有限要素解析を実行する。これにより、２つの三次元形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑制することが可能となる。
第３発明、第７発明、及び第１１発明では、複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断する。第１の節点と第２の節点とを結ぶ線分長が所定値より大きいと判断した場合、該線分上にメッシュ分割に用いる第３の節点を設定し、該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出し、曲線上まで移動した点における曲線の曲率半径を算出する。算出して移動距離を曲率半径で除した比率が所定値より大きい場合、第３の節点を曲線上へ移動して、第１、第２、第３の節点を結ぶ有限要素にメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を算出し、上述した処理を繰り返す。算出して移動距離を曲率半径で除した比率が所定値以下である場合、第３の節点に基づいてメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を算出し、上述した処理を繰り返す。第１の節点と第２の節点とを結ぶ線分長が所定値以下である場合、それ以上メッシュ分割をせず、有限要素解析を実行する。これにより、２つの二次元形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑えることが可能となる。
第４発明、第８発明、及び第１２発明では、複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断する。第１の節点と第２の節点とを結ぶ線分長が所定値より大きいと判断した場合、第１の節点と第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定し、該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出し、曲面上まで移動した点における曲面の曲率半径を算出する。移動距離を前記曲率半径で除した比率が所定値より大きい場合、第３の節点を曲面上へ移動して、第１、第２、第３の節点を結ぶ有限要素にメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を算出し、上述した処理を繰り返す。算出した比率が所定値以下である場合、第３の節点に基づいてメッシュ分割し、第３の節点を第２の節点とし、線分上に新たに設定した第４の節点を第３の節点として、移動距離及び曲率半径を算出し、上述した処理を繰り返す。第１の節点と第２の節点とを結ぶ線分長が所定値以下である場合、それ以上メッシュ分割をせず、有限要素解析を実行する。これにより、２つの三次元形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑制することが可能となる。
本発明によれば、複数の二次元形状モデル又は三次元形状モデルが接する接点又は接線近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑制することが可能となる。
以下、本発明をその実施の形態を示す図面に基づいて具体的に説明する。
（実施の形態１）
図１は、本発明の実施の形態１に係るメッシュ分割方法を適用した有限要素解析装置１の構成を示すブロック図である。図１において、有限要素解析装置１は、少なくとも、ＣＰＵ（中央演算装置）１１、記憶手段１２、ＲＯＭ１３、ＲＡＭ１４、通信回線に接続する通信手段１５、マウス及びキーボード等の入力手段１６、ディスプレイ等の出力手段１７及び補助記憶手段１８で構成される。
ＣＰＵ１１は、内部バス１９を介して有限要素解析装置１の上述したようなハードウェア各部と接続されており、上述したハードウェア各部を制御するとともに、ＲＯＭ１３に格納された制御プログラム又は補助記憶手段１８であるＣＤ−ＲＯＭ、ＤＶＤ等の（可搬型）記録媒体２を用いて記憶手段１２へ導入された制御プログラムに従って、種々のソフトウェア的機能を実行する。なお、記憶手段１２は、ハードディスク等の固定型記憶媒体であり、上述した制御プログラムの他、処理に必要なデータを事前に記憶しておく。
ＲＡＭ１４は、ＳＲＡＭ、フラッシュメモリ等で構成されソフトウェアの実行時に発生する一時的なデータを記憶する。通信手段１５は、内部バス１９に接続されており、外部からのデータの取得、外部装置の動作制御データ等を送受信する。
入力手段１６は、有限要素解析装置１を操作するために必要な文字キー、テンキー、各種のファンクションキー等を備えたキーボード、マウス等の入力媒体である。出力手段１７は、液晶表示装置、ＣＲＴディスプレイ等の表示装置であり、有限要素解析装置１の動作状態を表示したり、ユーザへ操作入力を促す画面を表示したり、解析結果をグラフィカルに表示するための画像データの表示等を行う。なお、出力手段１７をタッチパネル方式とすることにより、入力手段１６の各種のファンクションキーの内の一部又は全部を出力手段１７が代用することも可能である。
なお、有限要素解析装置１は、通信手段１５に接続されている外部のコンピュータから本発明に係るコンピュータプログラムをダウンロードし、ＣＰＵ１１にて処理を実行する形態であってもよい。
以下、上述した構成の有限要素解析装置１を用いた構造解析処理の手順について説明する。本実施の形態１では、２つの円柱状の三次元形状モデルが接する場合に、上面の両円が接する接点近傍の電磁場解析を例に説明する。図２は、本実施の形態１で用いる三次元形状モデルの斜視図である。図２の例では、２つの円柱状の三次元形状モデルが側面で接しており、説明を簡単にするために上面で接する両円の接点近傍でのメッシュ分割方法について説明する。
図３は、本発明の実施の形態１に係る有限要素解析装置１のＣＰＵ１１の処理手順を示すフローチャートである。本実施の形態１では、２つの円柱状の形状モデルが接する接点（接線）近傍でのメッシュ分割を例に挙げて説明する。図４は、受け付けた２つの円柱状の形状モデルの上面の一部を示す模式図である。
図４では、２つの円が接する接点を節点Ｐ１として、相隣り合う節点Ｐ２とを結ぶ線分により、メッシュ分割を開始する。有限要素解析装置１のＣＰＵ１１は、オペレータによって操作された入力装置１５から、又は通信手段１５を介して外部のコンピュータから、２つの円柱を表す三次元形状モデルを受け付ける(ステップＳ３０１)。そして、ＣＰＵ１１は、受け付けた形状モデルが接する接点の座標値を第１の節点Ｐ１の座標値として取得する(ステップＳ３０２)。
ＣＰＵ１１は、第１の節点Ｐ１に相隣り合い、円の周上の第２の節点Ｐ２を設定し、第２の節点Ｐ２の座標値を取得する（ステップＳ３０３）。ＣＰＵ１１は、第１の節点Ｐ１と第２の節点Ｐ２とを結ぶ線分上に、第３の節点Ｐ３（例えば線分Ｐ１Ｐ２の中点）を設定し、第３の節点の座標値を取得する（ステップＳ３０４）。
ＣＰＵ１１は、第３の節点Ｐ３を円の周上へ円の中心から周方向へ移動した場合の、周上の新たな節点Ｐ３’の座標値を取得する（ステップＳ３０５）。ＣＰＵ１１は、第３の節点Ｐ３が新たな節点Ｐ３’へ移動した場合の移動距離Ｄ（Ｐ３Ｐ３’）及び新たな節点Ｐ３’での曲線の曲率半径Ｒを算出する（ステップＳ３０６）。なお、図４の例では新たな節点Ｐ３’が存在する曲線が円であることから、該円の半径の長さを算出する。
ＣＰＵ１１は、算出した移動距離Ｄを曲率半径Ｒで除算して第３の節点Ｐ３の移動比率Ｄ／Ｒを算出し（ステップＳ３０７）、算出した移動比率Ｄ／Ｒが所定値より大きいか否かを判断する（ステップＳ３０８）。ＣＰＵ１１が、算出した移動比率が所定値より大きいと判断した場合（ステップＳ３０８：ＹＥＳ）、ＣＰＵ１１は、第１の節点Ｐ１と新たな節点Ｐ３’とを、及び第２の節点Ｐ２と新たな節点Ｐ３’とを、それぞれ結ぶ線分により、接点近傍のメッシュ分割を行う(ステップＳ３０９)。
ＣＰＵ１１は、新たな節点Ｐ３’を第２の節点Ｐ２とし、新たな節点Ｐ３’の座標値を第２の節点Ｐ２の座標値とする（ステップＳ３１０）。またＣＰＵ１１は、新たな節点Ｐ４を第１の節点Ｐ１と第２の節点Ｐ２とを結ぶ線分上に設定し、新たな節点Ｐ４を第３の節点Ｐ３とし、新たな節点Ｐ４の座標値を取得して、第３の節点Ｐ３の座標値とする（ステップＳ３１１）。ＣＰＵ１１は、ステップＳ３０５へ処理を戻し、上述した処理を繰り返す。
ＣＰＵ１１が、算出した移動比率が所定値以下であると判断した場合（ステップＳ３０８：ＮＯ）、ＣＰＵ１１は、これ以上メッシュ分割した場合、分割された有限要素の辺長比が大きく、一定の精度を保持して有限要素解析することが困難であると判断し、節点Ｐ３を移動することなくメッシュ分割を終了し（ステップＳ３１２）、有限要素解析を実行する（ステップＳ３１３）。
上述したＣＰＵ１１の判断処理において、例えば所定値が０．００１である場合、半径１の円弧モデルに対して０．００１のトレランスが存在するという相対トレランスを意味している。図５は、相対トレランスが０．００１である場合の円弧とトレランスとの関係を模式的に示す図である。図５（ａ）のように２つの点が離れていないとみなす最大距離、すなわち点と線、点と面、又は線と面とが離れていないとみなす最大距離がトレランス０．００１であり、実際には図５（ｂ）に示す図５（ａ）の破線で囲まれた部分の拡大図のようにメッシュを構成する線分と円弧とは離れているのに対して、離れていないとみなすものである。
このように、トレランスを超えてメッシュ分割しないように制約条件を付加することにより、２つの円柱状の形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑えることが可能となる。
なお、上述した実施の形態１では、ＣＰＵ１１が、算出した移動比率が所定値以下であると判断した場合（ステップＳ３０８：ＮＯ）、ＣＰＵ１１は、第３の節点Ｐ３を移動することなくメッシュ分割を終了（ステップＳ３１２）しているが、特にこれに限定されるものではなく、第３の節点Ｐ３を生成しなくても良いし、第３の節点Ｐ３を新たな節点Ｐ３’へ移動し、第１の節点Ｐ１と新たな節点Ｐ３’とを、及び第２の節点Ｐ２と新たな節点Ｐ３’とを、それぞれ結ぶ線分により、接点近傍のメッシュ分割を行ってから、メッシュ分割を終了するものであっても良い。
（実施の形態２）
以下、本発明の実施の形態２に係るメッシュ分割方法を適用した有限要素解析装置１について説明する。本実施の形態２に係るメッシュ分割方法を適用した有限要素解析装置１の構成は実施の形態１と同様であることから、同一の符号を付することで詳細な説明を省略する。本実施の形態２は、メッシュ分割する有限要素の１辺の長さの大小判断を実施の形態１に係るメッシュ分割方法に組み合わせた点に特徴を有する。
本実施の形態２でも、実施の形態１と同様、２つの円柱状の三次元形状モデルが接する場合に、上面の両円が接する接点近傍の電磁場解析を例に説明する。説明を簡単にするために上面で接する両円の接点近傍でのメッシュ分割方法について説明する。
図６は、本発明の実施の形態２に係る有限要素解析装置１のＣＰＵ１１の処理手順を示すフローチャートである。本実施の形態２は、２つの円柱状の形状モデルが接する接点（接線）近傍でのメッシュ分割を例に挙げて説明する。図７は、受け付けた２つの円柱状の形状モデルの上面の一部を示す模式図である。
図７では、２つの円が接する接点を節点Ｐ１として、相隣り合う節点Ｐ２とを結ぶ線分により、メッシュ分割を開始する。有限要素解析装置１のＣＰＵ１１は、オペレータによって操作された入力装置１５から、又は通信手段１５を介して外部のコンピュータから、２つの円柱を表す三次元形状モデルを受け付ける(ステップＳ６０１)。そして、ＣＰＵ１１は、受け付けた形状モデルが接する接点の座標値を第１の節点Ｐ１の座標値として取得する(ステップＳ６０２)。
ＣＰＵ１１は、第１の節点Ｐ１に相隣り合い、円の周上の第２の節点Ｐ２を設定し、第２の節点Ｐ２の座標値を取得する（ステップＳ６０３）。ＣＰＵ１１は、第１の節点Ｐ１と第２の節点Ｐ２との距離Ｐ１Ｐ２を算出して、距離Ｐ１Ｐ２が所定値より大きいか否かを判断する（ステップＳ６０４）。
ＣＰＵ１１が、距離Ｐ１Ｐ２が所定値より大きいと判断した場合（ステップＳ６０４：ＹＥＳ）、ＣＰＵ１１は、第１の節点Ｐ１と第２の節点Ｐ２とを結ぶ線分上に、第３の節点Ｐ３（例えば線分Ｐ１Ｐ２の中点）を設定し、第３の節点の座標値を取得する（ステップＳ６０５）。
ＣＰＵ１１は、第３の節点Ｐ３を円の周上へ円の中心から周方向へ移動した場合の、周上の新たな節点Ｐ３’の座標値を取得する（ステップＳ６０６）。ＣＰＵ１１は、第３の節点Ｐ３が新たな節点Ｐ３’へ移動した場合の移動距離Ｄ（Ｐ３Ｐ３’）及び新たな節点Ｐ３’での曲線の曲率半径Ｒを算出する（ステップＳ６０７）。なお、図７の例では新たな節点Ｐ３’が存在する曲線が円であることから、該円の半径の長さを算出する。
ＣＰＵ１１は、算出した移動距離Ｄを曲率半径Ｒで除算して第３の節点Ｐ３の移動比率Ｄ／Ｒを算出し（ステップＳ６０８）、算出した移動比率Ｄ／Ｒが所定値より大きいか否かを判断する（ステップＳ６０９）。ＣＰＵ１１が、算出した移動比率が所定値より大きいと判断した場合（ステップＳ６０９：ＹＥＳ）、ＣＰＵ１１は、第１の節点Ｐ１と新たな節点Ｐ３’とを、及び第２の節点Ｐ２と新たな節点Ｐ３’とを、それぞれ結ぶ線分により、接点近傍のメッシュ分割を行う(ステップＳ６１０)。
ＣＰＵ１１は、新たな節点Ｐ３’を第２の節点Ｐ２とし、新たな節点Ｐ３’の座標値を第２の節点Ｐ２の座標値とする（ステップＳ６１１）。またＣＰＵ１１は、新たな節点Ｐ４を第１の節点Ｐ１と第２の節点Ｐ２とを結ぶ線分上に設定し、新たな節点Ｐ４を第３の節点Ｐ３とし、新たな節点Ｐ４の座標値を取得して、第３の節点Ｐ３の座標値とする（ステップＳ６１２）。ＣＰＵ１１は、ステップＳ６０６へ処理を戻し、上述した処理を繰り返す。
ＣＰＵ１１が、算出した移動比率が所定値以下であると判断した場合（ステップＳ６０９：ＮＯ）、ＣＰＵ１１は、これ以上曲線形状に合わせたメッシュ分割（以下、フィッティング処理という）を行った場合、フィッティング処理により分割された有限要素の辺長比が大きく、一定の精度を保持して有限要素解析することが困難であると判断し、フィッティング処理を行うことなく第３の節点Ｐ３に基づいて辺長比を保持したメッシュ分割を行い（ステップＳ６１３）、第３の節点Ｐ３を第２の節点Ｐ２とし、第３の節点Ｐ３の座標値を第２の節点Ｐ２の座標値とする（ステップＳ６１４）。ＣＰＵ１１は、ステップＳ６１２へ処理を戻し、上述した処理を繰り返す。
ＣＰＵ１１が、距離Ｐ１Ｐ２が所定値以下であると判断した場合（ステップＳ６０４：ＮＯ）、ＣＰＵ１１は、節点Ｐ３を移動することなくメッシュ分割を終了し（ステップＳ６１５）、有限要素解析を実行する（ステップＳ６１６）。
このように、第１の節点Ｐ１と第２の節点Ｐ２との距離Ｐ１Ｐ２が所定値如何にならない範囲でメッシュ分割するように制約条件を付加することにより、２つの円柱状の形状モデルが接する接点近傍においても、辺長比を所定の比率以内に抑制することができることから「辺長比が大きくて形状の悪い要素」の発生を未然に防止することができ、解析対象である形状モデルの形状を大きく変えることなく有限要素解析を行うことができ、解析による誤差を最小限に抑えることが可能となる。
なお、上述した実施の形態１及び２では、２つの円柱状の三次元形状モデルの接線（接点）近傍でのメッシュ分割を例に挙げて説明しているが、複数の形状モデルが点又は線で接する場合であれば、二次元形状モデル同士、三次元形状モデル同士、二次元形状モデルと三次元形状モデルとの組み合わせ、のいずれであっても同様の効果が期待できる。また、電磁場解析に限定されるものではなく、有限要素法を用いるあらゆる有限要素解析に適用可能であることは言うまでも無い。
本発明の実施の形態１に係るメッシュ分割方法を適用した有限要素解析装置の構成を示すブロック図である。 本実施の形態１で用いる三次元形状モデルの斜視図である。 本発明の実施の形態１に係る有限要素解析装置のＣＰＵの処理手順を示すフローチャートである。 受け付けた２つの円柱状の形状モデルの上面の一部を示す模式図である。 相対トレランスが０．００１である場合の円弧とトレランスとの関係を模式的に示す図である。 本発明の実施の形態２に係る有限要素解析装置のＣＰＵの処理手順を示すフローチャートである。 受け付けた２つの円柱状の形状モデルの上面の一部を示す模式図である。
符号の説明
１ 有限要素解析装置
１１ ＣＰＵ
１２ 記憶手段
１３ ＲＯＭ
１４ ＲＡＭ
１５ 通信手段
１６ 入力手段
１７ 出力手段
１８ 補助記憶手段

Claims (12)

  1. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、
    前記有限要素解析装置は、前記プロセッサにより、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出するステップと、
    曲線上まで移動した点における曲線の曲率半径を算出するステップと、
    前記移動距離を前記曲率半径で除した比率を算出するステップと、
    算出した比率が所定値より大きいか否かを判断するステップと、
    算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、
    算出した比率が所定値以下であると判断した場合、メッシュ分割を停止するステップと
    を実行することを特徴とするメッシュ分割方法。
  2. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、
    前記有限要素解析装置は、前記プロセッサにより、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、
    該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出するステップと、
    曲面上まで移動した点における曲面の曲率半径を算出するステップと、
    前記移動距離を前記曲率半径で除した比率を算出するステップと、
    算出した比率が所定値より大きいか否かを判断するステップと、
    算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、
    算出した比率が所定値以下であると判断した場合、メッシュ分割を停止するステップと
    を実行することを特徴とするメッシュ分割方法。
  3. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、
    前記有限要素解析装置は、前記プロセッサにより、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断するステップと、
    所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出するステップと、
    曲線上まで移動した点における曲線の曲率半径を算出するステップと、
    前記移動距離を前記曲率半径で除した比率を算出するステップと、
    算出した比率が所定値より大きいか否かを判断するステップと、
    算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、
    算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、
    前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止するステップと
    を実行することを特徴とするメッシュ分割方法。
  4. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する処理を実行するプロセッサを備えた有限要素解析装置によるメッシュ分割方法であって、
    前記有限要素解析装置は、前記プロセッサにより、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長が所定値より大きいか否かを判断するステップと、
    所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定するステップと、
    該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出するステップと、
    曲面上まで移動した点における曲面の曲率半径を算出するステップと、
    前記移動距離を前記曲率半径で除した比率を算出するステップと、
    算出した比率が所定値より大きいか否かを判断するステップと、
    算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割するステップと、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返するステップと、
    算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離及び前記曲率半径を算出し、前記移動距離を前記曲率半径で除した比率が所定値より大きいか否かの判断を繰り返すステップと、
    前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止するステップと
    を実行することを特徴とするメッシュ分割方法。
  5. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段と、
    曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段と、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、
    算出した比率が所定値より大きいか否かを判断する判断手段と、
    該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段と、
    前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段と
    を備えることを特徴とする有限要素解析装置。
  6. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、
    該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段と、
    曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段と、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、
    算出した比率が所定値より大きいか否かを判断する判断手段と、
    該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段と、
    前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段と
    を備えることを特徴とする有限要素解析装置。
  7. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段と、
    算出した線分長が所定値より大きいか否かを判断する第１の判断手段と、
    該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段と、
    曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段と、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、
    算出した比率が所定値より大きいか否かを判断する第２の判断手段と、
    該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段と、
    前記第２の判断手段で算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段と、
    前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段と
    を備えることを特徴とする有限要素解析装置。
  8. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割する有限要素解析装置であって、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段と、
    算出した線分長が所定値より大きいか否かを判断する第１の判断手段と、
    該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段と、
    該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段と、
    曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段と、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段と、
    算出した比率が所定値より大きいか否かを判断する第２の判断手段と、
    該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段と、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段と、
    算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段と、
    前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段と
    を備えることを特徴とする有限要素解析装置。
  9. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、
    前記コンピュータを、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段、
    曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、
    算出した比率が所定値より大きいか否かを判断する判断手段、
    該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動し、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段、及び
    前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段
    として機能させることを特徴とするコンピュータプログラム。
  10. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、
    前記コンピュータを、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、
    該第３の節点を線分に対向する曲面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段、
    曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、
    算出した比率が所定値より大きいか否かを判断する判断手段、
    該判断手段で、比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動し、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、
    前記第３の節点を前記第２の節点とし、前記線分上に新たに設定した第４の節点を前記第３の節点として、前記移動距離算出手段へ戻る繰り返し手段、及び
    前記判断手段で、比率が所定値以下であると判断した場合、メッシュ分割を停止する分割停止手段
    として機能させることを特徴とするコンピュータプログラム。
  11. 複数の二次元形状モデルが点で接する場合に、接点近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、
    前記コンピュータを、
    複数の二次元形状モデルの接点を第１の節点とし、該第１の節点と一の二次元形状モデルの曲線上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段、
    算出した線分長が所定値より大きいか否かを判断する第１の判断手段、
    該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、
    該第３の節点を線分に直交する方向に曲線上まで移動した場合の移動距離を算出する移動距離算出手段、
    曲線上まで移動した点における曲線の曲率半径を算出する曲率半径算出手段、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、
    算出した比率が所定値より大きいか否かを判断する第２の判断手段、
    該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲線上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段、
    前記第２の判断手段で算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段、及び
    前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段
    として機能させることを特徴とするコンピュータプログラム。
  12. 複数の三次元形状モデルが点又は線で接する場合に、接点又は接線近傍を、有限要素解析を実行することが可能な複数のメッシュに分割するコンピュータで実行可能なコンピュータプログラムであって、
    前記コンピュータを、
    複数の三次元形状モデルの接点又は複数の三次元形状モデルが接触する接線上の１点を第１の節点とし、該第１の節点と一の三次元形状モデルの曲面上に位置する相隣る第２の節点とを結ぶ線分長を算出する節点間距離算出手段、
    算出した線分長が所定値より大きいか否かを判断する第１の判断手段、
    該第１の判断手段で所定値より大きいと判断した場合、前記第１の節点と前記第２の節点とを結ぶ線分上にメッシュ分割に用いる第３の節点を設定する節点設定手段、
    該第３の節点を線分に対向する局面の法線方向に曲面上まで移動した場合の移動距離を算出する移動距離算出手段、
    曲面上まで移動した点における曲面の曲率半径を算出する曲率半径算出手段、
    前記移動距離を前記曲率半径で除した比率を算出する比率算出手段、
    算出した比率が所定値より大きいか否かを判断する第２の判断手段、
    該第２の判断手段で算出した比率が所定値より大きいと判断した場合、前記第３の節点を曲面上へ移動して、前記第１、第２、第３の節点を結ぶ有限要素にメッシュ分割する分割手段、
    前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第１の繰り返し手段、
    算出した比率が所定値以下であると判断した場合、前記第３の節点に基づいてメッシュ分割し、前記第３の節点を第２の節点とし、前記第１の節点と該第２の節点とを結ぶ線分上に新たに設定した第４の節点を第３の節点として、前記移動距離算出手段へ戻る第２の繰り返し手段、及び
    前記第１の判断手段で前記第１の節点と前記第２の節点とを結ぶ線分長が所定値以下である場合、メッシュ分割を停止する分割停止手段
    として機能させることを特徴とするコンピュータプログラム。
JP2005248421A 2005-08-29 2005-08-29 メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム Active JP4809645B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005248421A JP4809645B2 (ja) 2005-08-29 2005-08-29 メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005248421A JP4809645B2 (ja) 2005-08-29 2005-08-29 メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム

Publications (2)

Publication Number Publication Date
JP2007065802A JP2007065802A (ja) 2007-03-15
JP4809645B2 true JP4809645B2 (ja) 2011-11-09

Family

ID=37927975

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005248421A Active JP4809645B2 (ja) 2005-08-29 2005-08-29 メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム

Country Status (1)

Country Link
JP (1) JP4809645B2 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20090005638A (ko) * 2007-07-09 2009-01-14 한국과학기술원 불일치 요소망을 해결하기 위한 다절점 천이 유한요소모델링 방법 및 기록매체
JP6295093B2 (ja) * 2014-02-04 2018-03-14 日本ユニシス株式会社 有限要素解析装置および板要素分割方法
CN104036126A (zh) * 2014-06-04 2014-09-10 广东电网公司电力科学研究院 肺泡收缩和扩展过程中颗粒物运动的数值模拟方法和系统
KR102094979B1 (ko) * 2018-03-30 2020-04-17 주식회사 마이다스아이티 1d 기반 설계데이터로부터 유한요소해석을 위한 2d 해석모델 생성 시스템 및 방법

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09282491A (ja) * 1996-04-18 1997-10-31 Nec Corp 有限要素法解析モデルにおける不適合メッシュ間の接合 方法
JP3807911B2 (ja) * 2000-09-07 2006-08-09 シャープ株式会社 解析装置、解析方法および解析プログラムを記録した記録媒体

Also Published As

Publication number Publication date
JP2007065802A (ja) 2007-03-15

Similar Documents

Publication Publication Date Title
KR101741402B1 (ko) 평행 지오데식 곡선들에 의해 모델링되는 부품의 설계
US10055896B1 (en) Interactive vertex manipulation system and methods for geometry repair
JP2006134305A (ja) 最小許容差スタックおよび最大許容差スタックの自動的計算を実行可能な方法、コンピュータ可読データ記憶媒体、およびコンピュータ支援設定システム
CN107066646B (zh) 受约束非对称细分网格中拓扑结构变化的方法、系统、介质
JP4809645B2 (ja) メッシュ分割方法、有限要素解析装置、及びコンピュータプログラム
Zegard et al. Advancing building engineering through structural and topology optimization
US20130226530A1 (en) Mesh generation system
EP3520084B1 (en) System and method for element quality improvement in 3d quadrilateral-dominant surface meshes
JP4620565B2 (ja) 解析メッシュ生成装置
JP2002283816A (ja) タイヤ有限要素モデルの作成方法、タイヤ有限要素モデルの作成装置およびプログラム
Remacle et al. Transient Mesh Adaptation Using Conforming and Non Conforming Mesh Modifications.
JP2007065803A (ja) 有限要素解析方法、有限要素解析装置、及びコンピュータプログラム
US7451048B2 (en) Method and apparatus for analyzing electromagnetic wave, and computer product
JPH06259404A (ja) 計算機を利用した解析方法
JP2007156604A (ja) 曲面作成方法及び曲面作成プログラム並びに３次元形状処理装置
JP2002207777A (ja) 中立面モデルの生成方法
JP2003132103A (ja) 設計支援方法及び設計支援装置
JP2014032473A (ja) 情報処理装置およびその方法
JP2006323714A (ja) 解析用メッシュデータ生成方法
JP3905373B2 (ja) メッシュ生成方法、メッシュ生成装置、コンピュータプログラム、及び記録媒体
JP4073732B2 (ja) タイムチャート作成装置、タイムチャート作成プログラム及び記録媒体
EP4195088A1 (en) 3d axis machining design
CN115719403B (zh) 构建刀具轨迹的方法和装置、电子设备和存储介质
JP2009129339A (ja) タイヤ性能シミュレーション方法、タイヤ性能シミュレーション装置、及びタイヤ性能シミュレーションプログラム
US20230168660A1 (en) Based on multiple manufacturing process variations, producing multiple contours representing predicted shapes of an ic design component

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080813

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110329

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110512

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110809

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110819

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140826

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4809645

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140826

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250