JP7284402B2 - Evaluation device, evaluation method and program - Google Patents
Evaluation device, evaluation method and program Download PDFInfo
- Publication number
- JP7284402B2 JP7284402B2 JP2019190732A JP2019190732A JP7284402B2 JP 7284402 B2 JP7284402 B2 JP 7284402B2 JP 2019190732 A JP2019190732 A JP 2019190732A JP 2019190732 A JP2019190732 A JP 2019190732A JP 7284402 B2 JP7284402 B2 JP 7284402B2
- Authority
- JP
- Japan
- Prior art keywords
- node
- stress
- nodes
- load
- virtual
- 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
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T30/00—Transportation of goods or passengers via railways, e.g. energy recovery or reducing air resistance
Description
本発明は、構造部材の評価装置、評価方法およびプログラムに関する。 The present invention relates to a structural member evaluation device, evaluation method, and program.
鉄道車両用台車枠(以下、台車枠と略記する。)は、車体を支持すると共に、車輪、車軸、主電動機、駆動装置、およびサスペンション部品等の種々の台車部品を支持する構造部材である。鉄道車両の走行時には、上記の種々の台車部品から台車枠に荷重が伝達される。このため、台車枠には、これらの荷重に対して十分な強度および耐久性が要求される。 A railroad vehicle bogie frame (hereinafter abbreviated as a bogie frame) is a structural member that supports a vehicle body and various bogie parts such as wheels, axles, traction motors, driving devices, and suspension parts. When the railway vehicle is running, loads are transmitted from the above-described various bogie components to the bogie frame. Therefore, the bogie frame is required to have sufficient strength and durability to withstand these loads.
そこで、従来、台車枠の強度を評価するために種々の方法が提案されている。例えば、特許文献1に開示された台車枠の強度評価方法では、台車枠のFEMモデルを用いて台車枠の強度評価が行われる。 Therefore, conventionally, various methods have been proposed to evaluate the strength of the bogie frame. For example, in the bogie frame strength evaluation method disclosed in Patent Document 1, bogie frame strength evaluation is performed using an FEM model of the bogie frame.
具体的には、特許文献1の強度評価方法では、FEMモデルの表面上に、仮想的に歪みゲージが定義される。歪みゲージは、FEMモデルにおいて、強度評価の対象となる表面上の節点Aと、当該節点Aの近傍の表面上の2つの節点B,Cから内挿して求められる点Dとを結ぶように定義される。そして、A-D間における荷重負荷前後の距離の変化に基づいて、A-D方向の応力が算出される。特許文献1の強度評価方法では、上記のようにして算出された応力に基づいて、台車枠の強度が評価される。 Specifically, in the strength evaluation method of Patent Document 1, strain gauges are virtually defined on the surface of the FEM model. In the FEM model, the strain gauge is defined to connect a node A on the surface to be evaluated for strength and a point D obtained by interpolating from two nodes B and C on the surface near the node A. be done. Then, the stress in the AD direction is calculated based on the change in the distance between AD before and after the load is applied. In the strength evaluation method of Patent Document 1, the strength of the bogie frame is evaluated based on the stress calculated as described above.
しかしながら、本発明者らによる研究の結果、特許文献1に開示された方法では、内挿して求められる点Dが、FEMモデルの表面上ではなく、要素が存在しない空間上に位置する場合があることが分かった。この場合、実際には算出する必要がない方向の応力が算出されるため、台車枠の強度評価を適切に行うことができない。 However, as a result of research by the present inventors, in the method disclosed in Patent Document 1, the point D obtained by interpolation may be located not on the surface of the FEM model but in a space where no elements exist. I found out. In this case, since the stress is calculated in a direction that does not actually need to be calculated, the strength evaluation of the bogie frame cannot be performed appropriately.
そこで、本発明は、構造部材の強度評価を適切に行うことができる、評価装置、評価方法およびプログラムを提供することを目的としている。 SUMMARY OF THE INVENTION Accordingly, an object of the present invention is to provide an evaluation apparatus, an evaluation method, and a program capable of appropriately evaluating the strength of a structural member.
本発明は、下記の評価装置、評価方法およびプログラムを要旨とする。 The gist of the present invention is the following evaluation device, evaluation method, and program.
(1)評価対象となる構造部材のFEM解析用の解析モデルを取得する、モデル取得部と、
前記解析モデルの表面上の任意の節点を対象節点として、前記表面上において前記対象節点と共通の要素に属する複数の節点から、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせを基準節点ペアとして抽出する、抽出部と、
前記抽出部に抽出された前記基準節点ペアを構成する2つの節点の間に仮想節点を規定する、規定部と、
前記解析モデルに荷重を与える、荷重負荷部と、
前記荷重負荷部が荷重を与える前後における前記解析モデルの表面上の各節点の変位を取得する、変位取得部と、
前記対象節点および前記仮想節点の変位に基づいて、前記構造部材の前記対象節点に対応する部位の強度に関する値を算出する、算出部と、
を備え、
前記抽出部は、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない場合には、基準節点ペアとして抽出しない、構造部材の評価装置。
(1) a model acquisition unit that acquires an analysis model for FEM analysis of a structural member to be evaluated;
With an arbitrary node on the surface of the analysis model as a target node, a combination of two nodes that can form a triangle together with the target node is selected from a plurality of nodes on the surface that belong to a common element with the target node. an extraction unit for extracting as a reference node pair;
a defining unit that defines a virtual node between two nodes forming the reference node pair extracted by the extracting unit;
a load applying unit that applies a load to the analysis model;
a displacement acquisition unit that acquires the displacement of each node on the surface of the analysis model before and after the load applying unit applies the load;
a calculation unit that calculates a strength value of a portion of the structural member corresponding to the target node based on the displacement of the target node and the virtual node;
with
The extraction unit does not extract a reference node pair even if the combination of two nodes that can form a triangle with the target node does not have a common element to which both nodes belong. evaluation equipment.
(2)前記算出部は、前記強度に関する値として、ひずみ、応力および安全率のうちの少なくとも一つを出力する、上記(1)に記載の構造部材の評価装置。 (2) The structural member evaluation device according to (1), wherein the calculation unit outputs at least one of strain, stress, and safety factor as the value related to strength.
(3)前記荷重負荷部は、異なる複数の荷重条件ごとに、前記解析モデルに荷重を与え、
前記変位取得部は、前記複数の荷重条件ごとに各節点の変位を取得し、
前記算出部は、
前記複数の荷重条件ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出する、第1算出部と、
前記第1算出部によって算出された前記ひずみに基づいて、前記対象節点と前記仮想節点との間に生じる応力を前記荷重条件ごとに算出する、第2算出部と、
前記第2算出部によって荷重条件ごとに算出された応力に基づいて、それらの平均応力成分の和を取って前記対象節点と前記仮想節点との間の複合平均応力を算出するとともに、それらの変動応力成分の二乗和平方根を取って前記対象節点と前記仮想節点との間の複合変動応力を算出する、第3算出部と、
予め設定された複合平均応力、複合変動応力および応力限界の関係と、前記第3算出部によって算出された前記複合平均応力および前記複合変動応力とに基づいて、前記構造部材の前記対象節点に対応する部位の安全率を算出する、第4算出部と、
を含む、上記(1)または(2)に記載の構造部材の評価装置。
(3) the load application unit applies a load to the analysis model for each of a plurality of different load conditions;
The displacement acquisition unit acquires the displacement of each node for each of the plurality of load conditions,
The calculation unit
a first calculator that calculates the strain generated between the target node and the virtual node for each of the plurality of load conditions;
a second calculator that calculates the stress generated between the target node and the virtual node for each load condition based on the strain calculated by the first calculator;
Based on the stress calculated for each load condition by the second calculation unit, the sum of the average stress components is calculated to calculate the composite average stress between the target node and the virtual node, and their fluctuations a third calculation unit that calculates a composite fluctuating stress between the target node and the virtual node by taking the square root of the sum of squares of the stress components;
Corresponding to the target node of the structural member based on the preset relationship between the composite mean stress, the composite fluctuating stress, and the stress limit, and the composite mean stress and the composite fluctuating stress calculated by the third calculation unit A fourth calculation unit that calculates the safety factor of the part to
The structural member evaluation device according to (1) or (2) above, comprising:
(4)前記規定部は、前記基準節点ペアを構成する2つの節点のうち、一方を第1基準節点、他方を第2基準節点とし、前記第1基準節点と前記第2基準節点との距離に対する第1基準節点と仮想節点との距離の比を0~1の範囲で変化させることによって仮想節点の位置を変更し、
前記第1算出部は、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出し、
前記第2算出部は、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じる応力を算出し、
前記第3算出部は、前記仮想節点の位置ごとに、前記複合平均応力および前記複合変動応力を算出し、
前記第4算出部は、前記仮想節点の位置ごとに前記安全率を算出し、算出された安全率のうち最も低い安全率を出力する、上記(3)に記載の構造部材の評価装置。
(4) The defining part defines one of the two nodes forming the reference node pair as a first reference node and the other as a second reference node, and sets the distance between the first reference node and the second reference node. changing the position of the virtual node by changing the ratio of the distance between the first reference node and the virtual node in the range of 0 to 1,
The first calculation unit calculates a strain generated between the target node and the virtual node for each position of the virtual node,
The second calculation unit calculates a stress generated between the target node and the virtual node for each position of the virtual node,
The third calculation unit calculates the composite average stress and the composite fluctuating stress for each position of the virtual node,
The structural member evaluation device according to (3) above, wherein the fourth calculation unit calculates the safety factor for each position of the virtual node, and outputs the lowest safety factor among the calculated safety factors.
(5)前記構造部材は、鉄道車両用台車枠である、上記(1)から(4)のいずれかに記載の構造部材の評価装置。 (5) The structural member evaluation apparatus according to any one of (1) to (4) above, wherein the structural member is a bogie frame for a railway vehicle.
(6)コンピュータによって実行される評価方法であって、
評価対象となる構造部材のFEM解析用の解析モデルを取得する、モデル取得ステップと、
前記解析モデルの表面上の任意の節点を対象節点として、前記表面上において前記対象節点と共通の要素に属する複数の節点から、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせを基準節点ペアとして抽出する、抽出ステップと、
前記抽出ステップで抽出された前記基準節点ペアを構成する2つの節点の間に仮想節点を規定する、規定ステップと、
前記解析モデルに荷重を与える、荷重負荷ステップと、
前記荷重負荷ステップで荷重を与える前後における前記解析モデルの表面上の各節点の変位を取得する、変位取得ステップと、
前記対象節点および前記仮想節点の変位に基づいて、前記構造部材の前記対象節点に対応する部位の強度に関する値を算出する、算出ステップと、
を備え、
前記抽出ステップでは、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない場合には、基準節点ペアとして抽出しない、構造部材の評価方法。
(6) A computer implemented evaluation method comprising:
a model acquisition step of acquiring an analysis model for FEM analysis of a structural member to be evaluated;
With an arbitrary node on the surface of the analysis model as a target node, a combination of two nodes that can form a triangle together with the target node is selected from a plurality of nodes on the surface that belong to a common element with the target node. an extraction step of extracting as reference node pairs;
a defining step of defining a virtual node between two nodes constituting the reference node pair extracted in the extracting step;
a loading step that applies a load to the analytical model;
a displacement acquisition step of acquiring the displacement of each node on the surface of the analysis model before and after the load is applied in the load application step;
a calculation step of calculating a strength value of a portion of the structural member corresponding to the target node based on the displacements of the target node and the virtual node;
with
In the extraction step, even if a combination of two nodes that can form a triangle with the target node, if there is no common element to which both nodes belong, the structural member is not extracted as a reference node pair. evaluation method.
(7)前記算出ステップでは、前記強度に関する値として、ひずみ、応力および安全率のうちの少なくとも一つを出力する、上記(6)に記載の構造部材の評価方法。 (7) The structural member evaluation method according to (6) above, wherein in the calculating step, at least one of strain, stress, and safety factor is output as the value related to strength.
(8)前記荷重負荷ステップでは、異なる複数の荷重条件ごとに、前記解析モデルに荷重を与え、
前記変位取得ステップでは、前記複数の荷重条件ごとに各節点の変位を取得し、
前記算出ステップは、
前記複数の荷重条件ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出する、第1算出ステップと、
前記第1算出ステップで算出された前記ひずみに基づいて、前記対象節点と前記仮想節点との間に生じる応力を前記荷重条件ごとに算出する、第2算出ステップと、
前記第2算出ステップで荷重条件ごとに算出された応力に基づいて、それらの平均応力成分の和を取って前記対象節点と前記仮想節点との間の複合平均応力を算出するとともに、それらの変動応力成分の二乗和平方根を取って前記対象節点と前記仮想節点との間の複合変動応力を算出する、第3算出ステップと、
予め設定された複合平均応力、複合変動応力および応力限界の関係と、前記第3算出ステップで算出された前記複合平均応力および前記複合変動応力とに基づいて、前記構造部材の前記対象節点に対応する部位の安全率を算出する、第4算出ステップと、
を含む、上記(6)または(7)に記載の構造部材の評価方法。
(8) in the load applying step, a load is applied to the analysis model for each of a plurality of different load conditions;
In the displacement acquisition step, the displacement of each node is acquired for each of the plurality of load conditions;
The calculating step includes:
a first calculation step of calculating the strain generated between the target node and the virtual node for each of the plurality of load conditions;
a second calculation step of calculating the stress generated between the target node and the virtual node for each load condition based on the strain calculated in the first calculation step;
Based on the stresses calculated for each load condition in the second calculation step, the sum of the average stress components is calculated to calculate the composite average stress between the target node and the virtual node, and their fluctuations a third calculating step of taking the square root of the sum of squares of the stress components to calculate a composite fluctuating stress between the target node and the virtual node;
Corresponding to the target node of the structural member based on the preset relationship between the composite mean stress, the composite fluctuating stress, and the stress limit, and the composite mean stress and the composite fluctuating stress calculated in the third calculation step A fourth calculation step for calculating the safety factor of the part to be
The method for evaluating the structural member according to (6) or (7) above, comprising:
(9)前記規定ステップでは、前記基準節点ペアを構成する2つの節点のうち、一方を第1基準節点、他方を第2基準節点とし、前記第1基準節点と前記第2基準節点との距離に対する第1基準節点と仮想節点との距離の比を0~1の範囲で変化させることによって仮想節点の位置を変更し、
前記第1算出ステップでは、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出し、
前記第2算出ステップでは、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じる応力を算出し、
前記第3算出ステップでは、前記仮想節点の位置ごとに、前記複合平均応力および前記複合変動応力を算出し、
前記第4算出ステップでは、前記仮想節点の位置ごとに前記安全率を算出し、算出された安全率のうち最も低い安全率を出力する、上記(8)に記載の構造部材の評価方法。
(9) In the defining step, one of two nodes constituting the reference node pair is defined as a first reference node and the other is defined as a second reference node, and a distance between the first reference node and the second reference node is determined. changing the position of the virtual node by changing the ratio of the distance between the first reference node and the virtual node in the range of 0 to 1,
In the first calculation step, a strain generated between the target node and the virtual node is calculated for each position of the virtual node;
In the second calculation step, a stress generated between the target node and the virtual node is calculated for each position of the virtual node;
In the third calculation step, the composite mean stress and the composite fluctuating stress are calculated for each position of the virtual node;
The structural member evaluation method according to (8) above, wherein in the fourth calculation step, the safety factor is calculated for each position of the virtual node, and the lowest safety factor among the calculated safety factors is output.
(10)前記構造部材は、鉄道車両用台車枠である、上記(6)から(9)のいずれかに記載の構造部材の評価方法。 (10) The structural member evaluation method according to any one of (6) to (9) above, wherein the structural member is a bogie frame for a railway vehicle.
(11)上記(6)から(10)のいずれかに記載の評価方法をコンピュータに実行させるプログラム。 (11) A program that causes a computer to execute the evaluation method according to any one of (6) to (10) above.
本発明によれば、構造部材の強度を適切に評価することができる。 ADVANTAGE OF THE INVENTION According to this invention, the intensity|strength of a structural member can be evaluated appropriately.
以下、本発明の実施の形態に係る構造部材の評価装置、評価方法およびプログラムについて説明する。なお、以下においては、構造部材の一例として鉄道車両用台車枠を挙げ、鉄道車両用台車枠の各部の強度に関する値を評価装置によって算出する場合について説明する。 Hereinafter, a structural member evaluation apparatus, evaluation method, and program according to an embodiment of the present invention will be described. In the following description, a bogie frame for a railway vehicle is taken as an example of a structural member, and a case in which values relating to the strength of each portion of the bogie frame for a railway vehicle are calculated by an evaluation device will be described.
[装置構成]
図1は、本発明の一実施形態に係る評価装置の概略構成を示すブロック図である。図1に示すように、本実施形態に係る評価装置10は、モデル取得部12、抽出部14、規定部16、荷重負荷部18、変位取得部20、および算出部22を備える。
[Device configuration]
FIG. 1 is a block diagram showing a schematic configuration of an evaluation device according to one embodiment of the present invention. As shown in FIG. 1 , the
図2は、本実施形態に係る評価装置10の評価対象となる鉄道車両用台車枠の一例を示す斜視図である。図2に示すように、台車枠100は、一対の側梁102と、一対の側梁102を接続する横梁104とを備えている。横梁104は、一対のパイプ部材104aを有している。以下においては、評価装置10が、台車枠100の解析モデルを用いて、台車枠100の各部の強度に関する値を算出する場合について説明する。ただし、評価装置10によって強度が評価される台車枠の構成は図2に示す台車枠100に限定されず、評価装置10は、種々の構成の台車枠の強度評価を行うために用いられる。
FIG. 2 is a perspective view showing an example of a railway vehicle bogie frame to be evaluated by the
図1に示すように、評価装置10のモデル取得部12は、台車枠100の解析モデルを取得する。解析モデルは、複数の要素および複数の節点を有する。本実施形態では、複数の節点には、コーナー節点(一次節点:要素の角に配置された節点)に加えて、中間節点(二次節点:コーナー節点とコーナー節点との間に位置するように要素辺上に配置された節点)が含まれてもよい。なお、解析モデルとしては、有限要素法(FEM)による応力解析において用いられる公知の解析モデルを利用することができるので、解析モデルについての詳細な説明は省略する。
As shown in FIG. 1 , the
本実施形態では、モデル取得部12は、他の装置において作成された解析モデルを与えられることによって解析モデルを取得してもよいし、公知のFEM解析ソフトと同様に、ユーザの操作に基づいて解析モデルを作成することによって解析モデルを取得してもよい。
In this embodiment, the
抽出部14は、解析モデルの表面上の任意の節点を対象節点として、解析モデルの表面上において対象節点と共通の要素に属する複数の節点から、複数の基準節点ペアを抽出する。具体的には、抽出部14は、対象節点とともに三角形を形成することができる2つの節点の複数の組み合わせをそれぞれ基準節点ペアとして抽出する。なお、対象節点は、強度に関する値(例えば、安全率)の算出対象となる節点である。本実施形態では、抽出部14は、例えば、ユーザによって指定された領域に属する複数の節点を対象節点として選択してもよく、解析モデルの全ての節点を対象節点として選択してもよい。また、抽出部14は、コーナー節点のみを対象節点として選択してもよい。以下、図面を用いて基準節点ペアの抽出方法について説明する。
The extracting
図3は、抽出部14による基準節点ペアの抽出方法を説明するための図である。なお、図3には、解析モデルの表面上の複数の要素のうちの一部の要素e1,e2,e3,e4が示されている。また、図3には、複数の節点として、コーナー節点A,B,C,D,E,Fおよび中間節点G,H,I,J,K,L,M,N,Oが示されている。なお、図3において、コーナー節点Bとコーナー節点Fとを結ぶ要素辺は存在しないものとする。すなわち、要素e1と要素e4との間は、要素が存在しない空間である。
FIG. 3 is a diagram for explaining how the
例えば、図3において、要素e1~e4に属するコーナー節点Aを対象節点とする。この場合、抽出部14は、対象節点Aと共通の要素e1~e4に属する複数の節点B~O(対象節点Aを除く節点)から、基準節点ペアを抽出する。具体的には、抽出部14は、対象節点Aとともに三角形を形成することができる2つの節点の組み合わせを基準節点ペアとして抽出する。
For example, in FIG. 3, the corner node A belonging to elements e1 to e4 is set as the target node. In this case, the
例えば、図3に示す要素B~Oのうち、節点Bおよび節点Cは、対象節点Aと共通の要素e1に属し、かつ対象節点Aとともに三角形を形成することができる。したがって、抽出部14は、節点Bおよび節点Cの組み合わせを基準節点ペアとして抽出する。また、例えば、節点Bおよび節点Hも、対象節点Aと共通の要素e1に属し、かつ対象節点Aとともに三角形を形成することができる。したがって、抽出部14は、節点Bおよび節点Hの組み合わせも基準節点ペアとして抽出する。
For example, among the elements B to O shown in FIG. 3, the node B and the node C belong to the common element e1 with the object node A and can form a triangle with the object node A. Therefore, the
ただし、抽出部14は、対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない節点の組み合わせを、基準節点ペアから除外する。例えば、図3に示す要素B~Oのうち、節点Hおよび節点Kは、対象節点Aと共通の要素e1,e2,e3に属し、かつ対象節点Aとともに三角形を形成することができる。しかしながら、節点Hおよび節点Kが属する共通の要素は存在しない。具体的には、節点Hは、要素e1に属し、節点Kは要素e2および要素e3に属している。このため、節点H,Kの組み合わせは、双方の節点が属する共通の要素が存在しない組み合わせとなり、抽出部14は、節点H,Kの組み合わせを基準節点ペアとして抽出しない。同様に、例えば、節点Hおよび節点Oは、対象節点Aと共通の要素e1,e4に属し、かつ対象節点Aとともに三角形を形成することができる。しかしながら、節点Hおよび節点Oが属する共通の要素は存在しない。具体的には、節点Hは、要素e1に属し、節点Oは要素e4に属している。このため、節点H,Oの組み合わせも、双方の節点が属する共通の要素が存在しない組み合わせであり、抽出部14は、節点H,Oの組み合わせを基準節点ペアとして抽出しない。
However, the
本実施形態では、抽出部14は、上記のルールに従って、対象節点ごとに1以上の基準節点ペアを抽出する。例えば、節点Aが対象節点に設定された場合には、抽出部14は、節点B,C、節点B,H、節点B,I、節点C,G、節点C,H、節点G,H、節点G,I、節点H,I、節点C,D、節点C,J、節点C,K、節点D,I、節点D,J、節点I,J、節点I,K、節点J,K、節点D,E、節点D,L、節点D,M、節点E,K、節点E,L、節点K,L、節点K,M、節点L,M、節点E,F、節点E,N、節点E,O、節点F,M,節点F,N、節点M,N、節点M,O、および節点N,Oの32通りの節点の組み合わせを、基準節点ペアとして抽出する。
In this embodiment, the
また、上述したように、本実施形態では、中間節点を対象節点とすることもできる。この場合も、抽出部14は、上述したルールに従って基準節点ペアを抽出する。例えば、中間節点Iが対象節点に設定された場合には、抽出部14は、対象節点Iとともに対象節点Iの周囲の要素e1,e2に属する複数の節点A~D、G、H、JおよびKから、複数の基準節点ペアを抽出する。具体的には、抽出部14は、節点A,B、節点A,D、節点A,G、節点A,H,節点A,J、節点A,K、節点B,C、節点B,G、節点B,H、節点C,D、節点C,G、節点C,H、節点C,J、節点C,K、節点D,J、節点D,K、節点G,Hおよび節点J,Kの18通りの節点の組み合わせを、基準節点ペアとして抽出する。
Further, as described above, in the present embodiment, intermediate nodes can also be set as target nodes. Also in this case, the
なお、詳細な説明は省略するが、例えば、評価装置10の負担を軽減するために、上述の基準節点ペアのうちの一部を、基準節点ペアから除外してもよい。例えば、対象節点と同一の要素辺上に位置する中間節点と、対象節点と異なる要素辺上に位置する中間節点との組み合わせを、基準節点ペアから除外してもよい。したがって、図3において節点Aを対象節点とした場合に、対象節点Aと同一の要素辺上に位置する中間節点G(または中間節点I)と、対象節点Aと異なる要素辺上に位置する中間節点Hとの組み合わせを、基準節点ペアから除外してもよい。なお、基準節点ペアから除外する節点の組み合わせは、ユーザが適宜設定することができる。
Although detailed description is omitted, for example, in order to reduce the burden on the
規定部16は、上記のようにして抽出した基準節点ペアごとに、仮想節点を規定する。具体的には、規定部16は、基準節点ペアを構成する2つの節点のうち、一方を第1基準節点、他方を第2基準節点とし、第1基準節点と前記第2基準節点との間に仮想節点を規定する。
The defining
以下、規定部16の機能を具体的に説明する。図4は、規定部16の機能を説明するための図である。なお、図4には、図3に示した要素e1が示されている。また、以下においては、節点Aを対象節点とする基準節点ペアG,Iに対して、仮想節点を規定する場合について説明する。
The functions of the defining
図4を参照して、基準節点ペアG,Iに対して仮想節点を規定する場合、規定部16は、例えば、節点Gを第1基準節点とし、節点Iを第2基準節点とし、第1基準節点Gと第2基準節点Iとを結ぶ線分GI上に、仮想節点Vを規定する。なお、本実施形態では、規定部16は、第1基準節点Gと第2基準節点Iとの距離D0に対する第1基準節点Gと仮想節点Vとの距離D1の比を0~1の範囲で変化させることによって、仮想節点の位置を変更する。これにより、規定部16は、複数の基準節点ペアに対してそれぞれ、複数の仮想節点を規定することができる。
Referring to FIG. 4, when defining virtual nodes for reference node pair G and I,
なお、本明細書において、第1基準節点と第2基準節点との距離に対する第1基準節点と仮想節点との距離の比が0になる場合とは、第1基準節点を仮想節点として規定することを意味し、第1基準節点と第2基準節点との距離に対する第1基準節点と仮想節点との距離の比が1になる場合とは、第2基準節点を仮想節点として規定することを意味する。 In this specification, when the ratio of the distance between the first reference node and the virtual node to the distance between the first reference node and the second reference node is 0, the first reference node is defined as a virtual node. When the ratio of the distance between the first reference node and the virtual node to the distance between the first reference node and the second reference node is 1, it means that the second reference node is defined as the virtual node. means.
なお、第1基準節点と第2基準節点との距離に対する第1基準節点と仮想節点との距離の比については、0~1の範囲で連続的に変化させてもよく、複数の値をユーザが適宜設定してもよい。また、例えば、1つの基準節点ペアに対して規定される仮想節点の数が予め設定されていてもよい。この場合、第1基準節点と第2基準節点との間において、予め設定された数の仮想節点を等間隔で規定してもよい。 Note that the ratio of the distance between the first reference node and the virtual node to the distance between the first reference node and the second reference node may be continuously changed within the range of 0 to 1, and a plurality of values may be set by the user. may be set as appropriate. Also, for example, the number of virtual nodes defined for one reference node pair may be set in advance. In this case, a preset number of virtual nodes may be defined at regular intervals between the first reference node and the second reference node.
荷重負荷部18は、モデル取得部12が取得した解析モデルに荷重を与える。なお、荷重負荷部18の機能は、公知のFEM解析ソフトの機能と同様である。したがって、荷重負荷部18は、公知のFEM解析ソフトと同様に構成することができるので、詳細な説明は省略する。
The
変位取得部20は、荷重負荷部18が荷重を与える前後における解析モデルの表面上の各節点(コーナー節点、中間節点および仮想節点)の変位を取得する。本実施形態では、変位取得部20は、公知のFEM解析ソフトと同様の機能により、コーナー節点および中間節点の変位を取得する。
The
また、本実施形態では、変位取得部20は、第1基準節点および第2基準節点の変位に基づいて、内挿(直接内挿)により仮想節点の変位を取得する。図4に示した例では、変位取得部20は、第1基準節点Gの変位および第2基準節点Iの変位に基づいて、仮想節点Vの変位を算出する。本実施形態では、変位取得部20は、例えば、仮想節点Vの変位を変位ベクトルとして算出する。具体的には、第1基準節点Gの変位をΔG(=(Δx1,Δy1,Δz1))とし、第2基準節点Iの変位をΔI(=(Δx2,Δy2,Δz2))とした場合、仮想節点Vの変位ΔVは、下記式(1)によって算出される。
ΔV=(1-(D1/D0))×ΔG+(D1/D0)× ΔI ・・・(1)
Further, in the present embodiment, the
ΔV=(1−(D 1 /D 0 ))×ΔG+(D 1 /D 0 )×ΔI (1)
算出部22は、変位取得部20が取得した対象節点の変位および仮想節点の変位に基づいて、台車枠100において対象節点に対応する部位の強度に関する値を算出する。本実施形態では、算出部22は、強度に関する値として、ひずみ、応力および安全率を算出する。図4に示した例では、算出部22は、対象節点Aの変位および仮想節点Vの変位に基づいて、台車枠100において対象節点Aに対応する部位のひずみ、応力および安全率を算出する。
The
本実施形態では、算出部22は、まず、対象節点Aの変位および仮想節点Vの変位に基づいて、荷重負荷部18が解析モデルに荷重を与えることによって対象節点Aと仮想節点Vとの間に生じるひずみを算出する。本実施形態では、算出部22は、例えば、下記式(2)によってひずみλを算出する。
λ=(AV1-AV0)/AV0 ・・・(2)
ただし、上記式において、AV0は、荷重負荷部18による荷重負荷前における対象節点Aと仮想節点Vとの距離を示し、AV1は、荷重負荷部18による荷重負荷後における対象節点Aと仮想節点Vとの距離を示す。
In the present embodiment, the
λ=(AV 1 -AV 0 )/AV 0 (2)
However, in the above formula, AV 0 indicates the distance between the target node A and the virtual node V before the load is applied by the
次に、算出部22は、上記のようにして算出したひずみと、解析モデルにおいて予め設定されたヤング率とに基づいて、対象節点Aと仮想節点Vとの間に生じる応力を算出する。
Next, the
最後に、算出部22は、上記のようにして算出した応力、および予め設定された応力と応力限界との関係に基づいて、台車枠100において対象節点Aに対応する部位の安全率を算出する。
Finally, the
なお、本実施形態では、算出部22は、仮想節点の位置ごとに、対象節点に対応する部位の安全率を算出する。算出部22は、例えば、仮想節点の位置ごとに算出した安全率のうち、最も低い安全率を、対象節点に対応する部位の安全率として出力してもよく、算出した全ての安全率を、仮想節点の位置に紐付けて出力してもよい。
Note that, in the present embodiment, the
以下、本実施形態に係る評価装置10の作用効果を説明する。図5は、本実施形態に係る評価装置10の作用効果を説明するための図である。具体的には、図5は、台車枠の解析モデル106の一部を示す斜視図である。解析モデル106は、底部106aおよび底部106aから立ち上がる壁部106bを有している。なお、図5においては、図面が煩雑になることを避けるために、解析モデル106の複数の要素(三角要素)のうち、底部106aおよび壁部106bの内表面の要素(三角要素)のみを示している。
The effects of the
上述したように、本実施形態に係る評価装置10では、対象節点と共通の要素に属する複数の節点から基準節点ペアを抽出する際に、対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない場合には、その節点の組み合わせは基準節点ペアから除外される。
As described above, in the
例えば、図5に示す解析モデル106において、底部106aと壁部106bとの境界に位置する節点aが対象節点に設定されたとする。この場合、対象節点aと共通の要素E1,E2に属する中間節点b、および対象節点と共通の要素E3,E4に属する中間節点cは、対象節点aとともに三角形を形成することができる。しかしながら、中間節点b,cが属する共通の要素は存在しない。したがって、対象節点aの安全率を算出する際には、中間節点b,cの組み合わせは基準節点ペアとして抽出されない。これにより、中間節点b,cを結ぶ線分bc上、すなわち要素が存在していない空間上に仮想節点が配置されることを防止することができる。その結果、対象節点aの周囲に生じるひずみおよび対象節点aに作用する応力を適切に算出することができ、解析対象(台車枠)において対象節点aに対応する部位の安全率を適切に算出することができる。
For example, in the
次に、評価装置10の具体的な構成について説明する。図6は、本発明の一実施形態に係る評価装置の構成を具体的に示すブロック図である。
Next, a specific configuration of the
本実施形態に係る評価装置10においても、上述の実施形態と同様に、モデル取得部12が台車枠のFEM解析用の解析モデルを取得し、抽出部14は、複数の基準節点ペアを抽出し、規定部16は、基準節点ペアごとに仮想節点を規定する。
Also in the
本実施形態では、荷重負荷部18は、異なる複数の荷重条件で解析モデルに荷重を与える。荷重条件の区分としては、例えば、上下方向の荷重、左右方向の荷重、前後方向の荷重、およびねじり方向の荷重等が挙げられる。本実施形態では、例えば、JIS E 4207:2019に示された荷重条件に従って、解析モデルに荷重を与えることができる。
In this embodiment, the
なお、本実施形態では、荷重条件ごとに、平均荷重および変動荷重が予め設定されている。また、本実施形態では、荷重負荷部18は、荷重条件ごとに、平均荷重および変動荷重を加算して得られる荷重を、解析モデルに与える。
In addition, in this embodiment, the average load and the fluctuating load are set in advance for each load condition. Further, in the present embodiment, the
変位取得部20は、複数の荷重条件ごとに、上述の実施形態と同様に、荷重負荷部18が荷重を与える前後における解析モデルの表面上の各節点(コーナー節点、中間節点および仮想節点)の変位を取得する。
For each of a plurality of load conditions, the
算出部22は、第1算出部22a、第2算出部22b、第3算出部22cおよび第4算出部22dを有する。第1算出部22aは、複数の荷重条件ごとに、対象節点と仮想節点との間に生じるひずみを算出する。具体的には、第1算出部22aは、例えば、上述の実施形態と同様に、上記式(2)によって、対象節点と仮想節点との間に生じるひずみを、規定部16が規定した仮想節点の位置ごとに算出する。
The
第2算出部22bは、第1算出部22aが算出したひずみと、解析モデルにおいて予め設定されたヤング率とに基づいて、対象節点と仮想節点との間に生じる応力を荷重条件ごとに算出する。本実施形態では、第2算出部22bは、規定部16が規定した仮想節点の位置ごとに、対象節点と仮想節点との間に生じる応力を算出する。
The
第3算出部22cは、第2算出部22bが荷重条件ごとに算出した応力の平均応力成分および変動応力成分を算出する。なお、上述したように、本実施形態では、荷重負荷部18は、平均荷重および変動荷重を加算して得られる荷重を解析モデルに与えている。すなわち、第2算出部22bによって算出される応力は、平均荷重および変動荷重を加算して得られる荷重に基づいて算出された応力である。
The
したがって、本実施形態では、第3算出部22cは、第2算出部22bが算出した応力と、荷重負荷部18が解析モデルに与えた荷重に対する平均荷重の割合とを乗算することにより、平均応力成分を算出することができる。また、第3算出部22cは、第2算出部22bが算出した応力と、荷重負荷部18が解析モデルに与えた荷重に対する変動荷重の割合とを乗算することにより、変動応力成分を算出することができる。
Therefore, in the present embodiment, the
例えば、荷重負荷部18が解析モデルに与えた荷重が130kNで、そのうち、平均荷重が100kNで、変動荷重が30kNであったとする。この条件で、第2算出部22bによって算出された応力が130MPaであったとする。この場合、第2算出部22bが算出した応力(130MPa)と、荷重負荷部18が解析モデルに与えた荷重(130kN)に対する平均荷重(100kN)の割合(100/130)とが乗算され、平均応力成分が100MPaであると算出される。また、第2算出部22bが算出した応力(130MPa)と、荷重負荷部18が解析モデルに与えた荷重(130kN)に対する変動荷重(30kN)の割合(30/130)とが乗算され、変動応力成分が30MPaであると算出される。
For example, assume that the load applied to the analysis model by the
第3算出部22cは、上記のようにして荷重条件ごとに算出した平均応力成分の和を取って、対象節点と仮想節点との間の複合平均応力を算出する。また、第3算出部22cは、上記のようにして荷重条件ごとに算出した変動応力成分の二乗和平方根を取って、対象節点と仮想節点との間の複合変動応力を算出する。本実施形態では、第3算出部22cは、規定部16が規定した仮想節点の位置ごとに、複合平均応力および複合変動応力を算出する。
The
第4算出部22dは、予め設定された複合平均応力、複合変動応力および応力限界の関係(以下、応力限界関係と記載する。)と、第3算出部22cが算出した複合平均応力および複合変動応力とに基づいて、台車枠100において対象節点に対応する部位の安全率を算出する。本実施形態では、応力限界関係は、図示しないデータベースに予め記憶されていてもよく、評価装置10による解析を実行する際に、ユーザによって入力されてもよい。以下、第4算出部22dによる安全率の算出方法の一例について図面を用いて説明する。
The
図7は、応力限界関係の一例を示す図(応力限界図)である。図7に示す応力限界関係では、平均応力を第1軸(横軸)とし、変動応力を第2軸(縦軸)とするグラフ上に、応力限界線が規定されている。 FIG. 7 is a diagram (stress limit diagram) showing an example of the stress limit relationship. In the stress limit relationship shown in FIG. 7, a stress limit line is defined on a graph having mean stress as the first axis (horizontal axis) and fluctuating stress as the second axis (vertical axis).
図7に示すように、本実施形態では、第4算出部22dは、まず、応力限界線が示されたグラフ上に、第3算出部22cが算出した複合平均応力および複合変動応力を示す点をプロットする。なお、図7には、複合平均応力が「-50MPa」で、複合変動応力が「50MPa」の場合を示している。
As shown in FIG. 7, in the present embodiment, the
次に、第4算出部22dは、プロットした点と原点とを通る直線(以下、判定線と記載する。)を算出する。次に、第4算出部22dは、下記式(3)に基づいて、安全率SFを算出する。
SF=OP1/OP2 ・・・(3)
なお、上記式(3)において、OP1は、判定線および応力限界線の交点と原点とを結ぶ線分の長さであり、OP2は、判定線において、第4算出部22dがプロットした点と原点とを結ぶ線分の長さである。
Next, the
SF=OP1/OP2 (3)
In the above formula (3), OP1 is the length of the line segment connecting the point of intersection of the decision line and the stress limit line and the origin, and OP2 is the point plotted by the
なお、応力限界関係としては、例えば、JIS E 4207:2019に示された応力限界図を用いることができる。 As the stress limit relationship, for example, the stress limit diagram shown in JIS E 4207:2019 can be used.
本実施形態では、第4算出部22dは、規定部16が規定した仮想節点の位置ごとに、上述の(3)式を用いて、台車枠100において対象節点に対応する部位の安全率を算出する。また、本実施形態では、第4算出部22dは、仮想節点の位置ごとに算出した安全率のうち、最も低い安全率を、対象節点に対応する部位の安全率として出力する。
In the present embodiment, the
なお、本実施形態では、例えば、解析モデルに予め形成された全ての節点(仮想節点を除く節点:コーナー節点および中間節点)が対象節点に設定される。したがって、第4算出部22dは、解析モデルに予め設定された全ての節点について、安全率を算出することができる。また、本実施形態では、第4算出部22dは、台車枠100において解析モデルの各節点に対応する部位の安全率を、安全率分布図として出力してもよい。この場合、第4算出部22dは、台車枠100の各部位の安全率を、その大きさに基づいてコンター表示してもよい。
In this embodiment, for example, all nodes (nodes other than virtual nodes: corner nodes and intermediate nodes) previously formed in the analysis model are set as target nodes. Therefore, the
以上のように、本実施形態においては、複合平均応力および複合変動応力を用いて、台車枠100の各部の安全率が算出される。ここで、本実施形態では、複合平均応力および複合変動応力を算出するに際して、要素が存在していない空間上に仮想節点が配置されることが防止される。これにより、各対象節点に作用する複合平均応力および複合変動応力を適切に算出することができ、台車枠100の各部の安全率を適切に算出することができる。
As described above, in the present embodiment, the safety factor of each part of the
[装置動作]
次に、本実施形態に係る評価装置10の動作について図8を用いて説明する。図8は、本発明の一実施形態に係る評価装置10の動作を示すフロー図である。なお、本実施形態に係る評価方法は、評価装置10を動作させることによって実施される。
[Device operation]
Next, the operation of the
図8に示すように、本実施形態では、上述したように、まず、モデル取得部12が、台車枠のFEM解析用の解析モデルを取得する(ステップS1)。次に、上述したように、抽出部14が、解析モデルに予め形成された複数の節点をそれぞれ対象節点に設定するとともに、対象節点ごとに複数の基準節点ペアを抽出する(ステップS2)。
As shown in FIG. 8, in this embodiment, as described above, first, the
次に、上述したように、規定部16が、基準節点ペアごとに仮想節点を規定する(ステップS3)。次に、上述したように、荷重負荷部18が解析モデルに荷重を与える(ステップS4)。
Next, as described above, the
次に、上述したように、変位取得部20が、荷重負荷部18が荷重を与える前後における解析モデルの表面上の各節点(コーナー節点、中間節点および仮想節点)の変位を取得する(ステップS5)。
Next, as described above, the
最後に、上述したように、算出部22が、台車枠100において対象節点に対応する部位の強度に関する値を算出する(ステップS6)。ステップS6では、ステップS2において対象節点に設定された複数の節点についてそれぞれ、強度に関する値(ひずみ、応力および安全率)が算出される。なお、ステップS6では、上述したように、台車枠100において解析モデルの各節点に対応する部位の安全率が、安全率分布図として出力されてもよい。
Finally, as described above, the
[プログラム]
本実施の形態におけるプログラムは、コンピュータに、図8に示すステップS1~S6を実行させるプログラムであればよい。このプログラムをコンピュータにインストールし、実行することによって、本実施の形態における評価装置と評価方法とを実現することができる。この場合、コンピュータのプロセッサは、モデル取得部12、抽出部14、規定部16、荷重負荷部18、変位取得部20、および算出部22(第1算出部22a、第2算出部22b、第3算出部22cおよび第4算出部22d)として機能し、処理を行う。
[program]
The program in this embodiment may be any program that causes a computer to execute steps S1 to S6 shown in FIG. By installing this program in a computer and executing it, the evaluation apparatus and evaluation method according to the present embodiment can be realized. In this case, the processor of the computer includes the
また、本実施形態に係るプログラムは、複数のコンピュータによって構築されたコンピュータシステムによって実行されてもよい。この場合、各コンピュータがそれぞれ、モデル取得部12、抽出部14、規定部16、荷重負荷部18、変位取得部20、および算出部22(第1算出部22a、第2算出部22b、第3算出部22cおよび第4算出部22d)のいずれかとして機能してもよい。
Also, the program according to the present embodiment may be executed by a computer system constructed by a plurality of computers. In this case, each computer has a
[物理構成]
図9は、本発明の一実施形態に係る評価装置を実現するコンピュータの一例を示すブロック図である。コンピュータ110は、本実施形態に係るプログラムを実行することによって、本実施形態に係る評価装置10を実現する。
[Physical configuration]
FIG. 9 is a block diagram showing an example of a computer that implements the evaluation device according to one embodiment of the present invention. The
図9に示すように、コンピュータ110は、CPU111と、メインメモリ112と、記憶装置113と、入力インターフェイス114と、表示コントローラ115と、データリーダ/ライタ116と、通信インターフェイス117とを備える。これらの各部は、バス121を介して、互いにデータ通信可能に接続される。なお、コンピュータ110は、CPU111に加えて、又はCPU111に代えて、GPU(Graphics Processing Unit)、又はFPGA(Field-Programmable Gate Array)を備えていてもよい。
As shown in FIG. 9,
CPU111は、記憶装置113に格納された、本実施の形態におけるプログラム(コード)をメインメモリ112に展開し、これらを所定順序で実行することにより、各種の演算を実施する。メインメモリ112は、典型的には、DRAM(Dynamic Random Access Memory)等の揮発性の記憶装置である。また、本実施の形態におけるプログラムは、コンピュータ読み取り可能な記録媒体120に格納された状態で提供される。なお、本実施の形態におけるプログラムは、通信インターフェイス117を介して接続されたインターネット上で流通するものであってもよい。
The
また、記憶装置113の具体例としては、ハードディスクドライブの他、フラッシュメモリ等の半導体記憶装置が挙げられる。入力インターフェイス114は、CPU111と、キーボードおよびマウスといった入力機器118との間のデータ伝送を仲介する。表示コントローラ115は、ディスプレイ装置119と接続され、ディスプレイ装置119での表示を制御する。
Further, as a specific example of the
データリーダ/ライタ116は、CPU111と記録媒体120との間のデータ伝送を仲介し、記録媒体120からのプログラムの読み出し、およびコンピュータ110における処理結果の記録媒体120への書き込みを実行する。通信インターフェイス117は、CPU111と、他のコンピュータとの間のデータ伝送を仲介する。
Data reader/
また、記録媒体120の具体例としては、CF(Compact Flash(登録商標))およびSD(Secure Digital)等の汎用的な半導体記憶デバイス、フレキシブルディスク(Flexible Disk)等の磁気記録媒体、又はCD-ROM(Compact Disk Read Only Memory)などの光学記録媒体が挙げられる。
Specific examples of the
なお、本実施形態に係る評価装置10は、プログラムがインストールされたコンピュータではなく、各部に対応したハードウェアを用いることによって実現されてもよい、また、評価装置10は、一部がプログラムで実現され、残りの部分がハードウェアで実現されていてもよい。
Note that the
[他の実施形態]
上述の実施形態では、評価装置が、鉄道車両用台車枠の各部の強度に関する値を算出する場合について説明したが、本発明に係る評価装置の評価対象となる構造部材は、鉄道車両用台車枠に限定されない。例えば、評価装置が、他の機械構造部材(自動車用車体等)、建築構造部材および土木構造部材等の種々の構造部材の各部の強度に関する値を算出してもよい。
[Other embodiments]
In the above-described embodiment, the case where the evaluation device calculates the strength value of each part of the railroad vehicle bogie frame has been described. is not limited to For example, the evaluation device may calculate the strength of each part of various structural members such as other mechanical structural members (vehicle bodies for automobiles, etc.), building structural members, and civil engineering structural members.
また、上述の実施形態では、算出部22が強度に関する値として安全率を出力する場合について説明したが、算出部22が強度に関する値としてひずみまたは応力(複合平均応力および複合変動応力を含む。)を出力してもよい。すなわち、算出部22が、ひずみ、応力および安全率のうちの少なくとも一つを出力してもよい。なお、上述したように、算出部22は、仮想節点の位置ごとに、対象節点に対応する部位のひずみおよび応力を算出する。算出部22は、例えば、仮想節点の位置ごとに算出したひずみまたは応力のうち、最も大きいひずみまたは応力を、対象節点に対応する部位のひずみまたは応力として出力してもよく、算出した全てのひずみまたは応力を、仮想節点の位置に紐付けて出力してもよい。
Further, in the above-described embodiment, the case where the
また、上述の実施形態では、算出部22が強度に関する値として、ひずみ、応力および安全率を算出する場合について説明したが、算出部が強度に関する値として、安全率を算出しなくてもよく、安全率および応力を算出しなくてもよい。
Further, in the above-described embodiment, the case where the
本発明によれば、構造部材の強度を適切に評価することができる。 ADVANTAGE OF THE INVENTION According to this invention, the intensity|strength of a structural member can be evaluated appropriately.
10 評価装置
12 モデル取得部
14 抽出部
16 規定部
18 荷重負荷部
20 変位取得部
22 算出部
10
Claims (11)
前記解析モデルの表面上の任意の節点を対象節点として、前記表面上において前記対象節点と共通の要素に属する複数の節点から、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせを基準節点ペアとして抽出する、抽出部と、
前記抽出部に抽出された前記基準節点ペアを構成する2つの節点の間に仮想節点を規定する、規定部と、
前記解析モデルに荷重を与える、荷重負荷部と、
前記荷重負荷部が荷重を与える前後における前記解析モデルの表面上の各節点の変位を取得する、変位取得部と、
前記対象節点および前記仮想節点の変位に基づいて、前記構造部材の前記対象節点に対応する部位の強度に関する値を算出する、算出部と、
を備え、
前記抽出部は、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない場合には、基準節点ペアとして抽出しない、構造部材の評価装置。 a model acquisition unit that acquires an analysis model for FEM analysis of a structural member to be evaluated;
With an arbitrary node on the surface of the analysis model as a target node, a combination of two nodes that can form a triangle together with the target node is selected from a plurality of nodes on the surface that belong to a common element with the target node. an extraction unit for extracting as a reference node pair;
a defining unit that defines a virtual node between two nodes forming the reference node pair extracted by the extracting unit;
a load applying unit that applies a load to the analysis model;
a displacement acquisition unit that acquires the displacement of each node on the surface of the analysis model before and after the load applying unit applies the load;
a calculation unit that calculates a strength value of a portion of the structural member corresponding to the target node based on the displacement of the target node and the virtual node;
with
The extraction unit does not extract a reference node pair even if the combination of two nodes that can form a triangle with the target node does not have a common element to which both nodes belong. evaluation equipment.
前記変位取得部は、前記複数の荷重条件ごとに各節点の変位を取得し、
前記算出部は、
前記複数の荷重条件ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出する、第1算出部と、
前記第1算出部によって算出された前記ひずみに基づいて、前記対象節点と前記仮想節点との間に生じる応力を前記荷重条件ごとに算出する、第2算出部と、
前記第2算出部によって荷重条件ごとに算出された応力に基づいて、それらの平均応力成分の和を取って前記対象節点と前記仮想節点との間の複合平均応力を算出するとともに、それらの変動応力成分の二乗和平方根を取って前記対象節点と前記仮想節点との間の複合変動応力を算出する、第3算出部と、
予め設定された複合平均応力、複合変動応力および応力限界の関係と、前記第3算出部によって算出された前記複合平均応力および前記複合変動応力とに基づいて、前記構造部材の前記対象節点に対応する部位の安全率を算出する、第4算出部と、
を含む、請求項1または2に記載の構造部材の評価装置。 The load applying unit applies a load to the analysis model for each of a plurality of different load conditions,
The displacement acquisition unit acquires the displacement of each node for each of the plurality of load conditions,
The calculation unit
a first calculator that calculates the strain generated between the target node and the virtual node for each of the plurality of load conditions;
a second calculator that calculates the stress generated between the target node and the virtual node for each load condition based on the strain calculated by the first calculator;
Based on the stress calculated for each load condition by the second calculation unit, the sum of the average stress components is calculated to calculate the composite average stress between the target node and the virtual node, and their fluctuations a third calculation unit that calculates a composite fluctuating stress between the target node and the virtual node by taking the square root of the sum of squares of the stress components;
Corresponding to the target node of the structural member based on the preset relationship between the composite mean stress, the composite fluctuating stress, and the stress limit, and the composite mean stress and the composite fluctuating stress calculated by the third calculation unit A fourth calculation unit that calculates the safety factor of the part to
The structural member evaluation device according to claim 1 or 2, comprising:
前記第1算出部は、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出し、
前記第2算出部は、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じる応力を算出し、
前記第3算出部は、前記仮想節点の位置ごとに、前記複合平均応力および前記複合変動応力を算出し、
前記第4算出部は、前記仮想節点の位置ごとに前記安全率を算出し、算出された安全率のうち最も低い安全率を出力する、請求項3に記載の構造部材の評価装置。 The defining unit defines one of two nodes constituting the reference node pair as a first reference node and the other as a second reference node, and defines a first reference node for a distance between the first reference node and the second reference node. changing the position of the virtual node by changing the ratio of the distance between the reference node and the virtual node in the range of 0 to 1;
The first calculation unit calculates a strain generated between the target node and the virtual node for each position of the virtual node,
The second calculation unit calculates a stress generated between the target node and the virtual node for each position of the virtual node,
The third calculation unit calculates the composite average stress and the composite fluctuating stress for each position of the virtual node,
4. The structural member evaluation apparatus according to claim 3, wherein said fourth calculator calculates said safety factor for each position of said virtual node, and outputs the lowest safety factor among the calculated safety factors.
評価対象となる構造部材のFEM解析用の解析モデルを取得する、モデル取得ステップと、
前記解析モデルの表面上の任意の節点を対象節点として、前記表面上において前記対象節点と共通の要素に属する複数の節点から、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせを基準節点ペアとして抽出する、抽出ステップと、
前記抽出ステップで抽出された前記基準節点ペアを構成する2つの節点の間に仮想節点を規定する、規定ステップと、
前記解析モデルに荷重を与える、荷重負荷ステップと、
前記荷重負荷ステップで荷重を与える前後における前記解析モデルの表面上の各節点の変位を取得する、変位取得ステップと、
前記対象節点および前記仮想節点の変位に基づいて、前記構造部材の前記対象節点に対応する部位の強度に関する値を算出する、算出ステップと、
を備え、
前記抽出ステップでは、前記対象節点とともに三角形を形成することができる2つの節点の組み合わせであっても、双方の節点が属する共通の要素が存在しない場合には、基準節点ペアとして抽出しない、構造部材の評価方法。 A computer implemented evaluation method comprising:
a model acquisition step of acquiring an analysis model for FEM analysis of a structural member to be evaluated;
With an arbitrary node on the surface of the analysis model as a target node, a combination of two nodes that can form a triangle together with the target node is selected from a plurality of nodes on the surface that belong to a common element with the target node. an extraction step of extracting as reference node pairs;
a defining step of defining a virtual node between two nodes constituting the reference node pair extracted in the extracting step;
a loading step that applies a load to the analytical model;
a displacement acquisition step of acquiring the displacement of each node on the surface of the analysis model before and after the load is applied in the load application step;
a calculation step of calculating a strength value of a portion of the structural member corresponding to the target node based on the displacements of the target node and the virtual node;
with
In the extraction step, even if a combination of two nodes that can form a triangle with the target node, if there is no common element to which both nodes belong, the structural member is not extracted as a reference node pair. evaluation method.
前記変位取得ステップでは、前記複数の荷重条件ごとに各節点の変位を取得し、
前記算出ステップは、
前記複数の荷重条件ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出する、第1算出ステップと、
前記第1算出ステップで算出された前記ひずみに基づいて、前記対象節点と前記仮想節点との間に生じる応力を前記荷重条件ごとに算出する、第2算出ステップと、
前記第2算出ステップで荷重条件ごとに算出された応力に基づいて、それらの平均応力成分の和を取って前記対象節点と前記仮想節点との間の複合平均応力を算出するとともに、それらの変動応力成分の二乗和平方根を取って前記対象節点と前記仮想節点との間の複合変動応力を算出する、第3算出ステップと、
予め設定された複合平均応力、複合変動応力および応力限界の関係と、前記第3算出ステップで算出された前記複合平均応力および前記複合変動応力とに基づいて、前記構造部材の前記対象節点に対応する部位の安全率を算出する、第4算出ステップと、
を含む、請求項6または7に記載の構造部材の評価方法。 In the load applying step, a load is applied to the analysis model for each of a plurality of different load conditions,
In the displacement acquisition step, the displacement of each node is acquired for each of the plurality of load conditions;
The calculating step includes:
a first calculation step of calculating the strain generated between the target node and the virtual node for each of the plurality of load conditions;
a second calculation step of calculating the stress generated between the target node and the virtual node for each load condition based on the strain calculated in the first calculation step;
Based on the stresses calculated for each load condition in the second calculation step, the sum of the average stress components is calculated to calculate the composite average stress between the target node and the virtual node, and their fluctuations a third calculating step of taking the square root of the sum of squares of the stress components to calculate a composite fluctuating stress between the target node and the virtual node;
Corresponding to the target node of the structural member based on the preset relationship between the composite mean stress, the composite fluctuating stress, and the stress limit, and the composite mean stress and the composite fluctuating stress calculated in the third calculation step A fourth calculation step for calculating the safety factor of the part to be
The method for evaluating a structural member according to claim 6 or 7, comprising:
前記第1算出ステップでは、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じるひずみを算出し、
前記第2算出ステップでは、前記仮想節点の位置ごとに、前記対象節点と前記仮想節点との間に生じる応力を算出し、
前記第3算出ステップでは、前記仮想節点の位置ごとに、前記複合平均応力および前記複合変動応力を算出し、
前記第4算出ステップでは、前記仮想節点の位置ごとに前記安全率を算出し、算出された安全率のうち最も低い安全率を出力する、請求項8に記載の構造部材の評価方法。 In the defining step, one of the two nodes forming the reference node pair is set as a first reference node and the other is set as a second reference node, and the distance between the first reference node and the second reference node is defined as a first reference node. changing the position of the virtual node by changing the ratio of the distance between the reference node and the virtual node in the range of 0 to 1;
In the first calculation step, a strain generated between the target node and the virtual node is calculated for each position of the virtual node;
In the second calculation step, a stress generated between the target node and the virtual node is calculated for each position of the virtual node;
In the third calculation step, the composite mean stress and the composite fluctuating stress are calculated for each position of the virtual node;
9. The structural member evaluation method according to claim 8, wherein in said fourth calculating step, said safety factor is calculated for each position of said virtual node, and the lowest safety factor among the calculated safety factors is output.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019190732A JP7284402B2 (en) | 2019-10-18 | 2019-10-18 | Evaluation device, evaluation method and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019190732A JP7284402B2 (en) | 2019-10-18 | 2019-10-18 | Evaluation device, evaluation method and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021067991A JP2021067991A (en) | 2021-04-30 |
JP7284402B2 true JP7284402B2 (en) | 2023-05-31 |
Family
ID=75637187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019190732A Active JP7284402B2 (en) | 2019-10-18 | 2019-10-18 | Evaluation device, evaluation method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7284402B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114036689B (en) * | 2021-11-26 | 2023-03-14 | 天津大学 | Iteration-based component strength stress optimization method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002073705A (en) | 2000-08-29 | 2002-03-12 | Kawasaki Heavy Ind Ltd | Method for preparing stress analysis model for bogie frame of rolling stock |
JP2005190242A (en) | 2003-12-26 | 2005-07-14 | Sumitomo Metal Ind Ltd | Strength evaluation method for truck frame, and manufacturing method of truck frame |
JP2011173573A (en) | 2010-02-25 | 2011-09-08 | Bridgestone Corp | Simulation method and simulation device |
-
2019
- 2019-10-18 JP JP2019190732A patent/JP7284402B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002073705A (en) | 2000-08-29 | 2002-03-12 | Kawasaki Heavy Ind Ltd | Method for preparing stress analysis model for bogie frame of rolling stock |
JP2005190242A (en) | 2003-12-26 | 2005-07-14 | Sumitomo Metal Ind Ltd | Strength evaluation method for truck frame, and manufacturing method of truck frame |
JP2011173573A (en) | 2010-02-25 | 2011-09-08 | Bridgestone Corp | Simulation method and simulation device |
Also Published As
Publication number | Publication date |
---|---|
JP2021067991A (en) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Balduzzi et al. | Non-prismatic beams: a simple and effective Timoshenko-like model | |
Montenegro et al. | Wheel–rail contact formulation for analyzing the lateral train–structure dynamic interaction | |
Guo et al. | Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profiles | |
Limebeer et al. | Optimal control of a formula one car on a three-dimensional track—Part 2: Optimal control | |
Morfidis | Vibration of Timoshenko beams on three-parameter elastic foundation | |
Auricchio et al. | The dimensional reduction approach for 2D non-prismatic beam modelling: a solution based on Hellinger–Reissner principle | |
JP6453627B2 (en) | Seismic analysis apparatus, method and program | |
Amiri et al. | A procedure for in situ wind load reconstruction from structural response only based on field testing data | |
Ghasemzadeh et al. | Dynamic high order numerical manifold method based on weighted residual method | |
Rodríguez-Tembleque et al. | A FEM–BEM fast methodology for 3D frictional contact problems | |
JP7284402B2 (en) | Evaluation device, evaluation method and program | |
Jang et al. | Remesh-free shape optimization using the wavelet-Galerkin method | |
Sreenath et al. | Beam and shell element model for advanced analysis of steel structural members | |
Dong et al. | On principal shear axes for correction factors in Timoshenko beam theory | |
Ignesti et al. | Development of a model for the prediction of wheel and rail wear in the railway field | |
JP2018200288A (en) | Earthquake response analysis method and earthquake response analysis program | |
JP7469653B2 (en) | Evaluation device, evaluation method, and program | |
Amaral et al. | Computational modeling and constructal design method applied to the geometric evaluation of stiffened thin steel plates considering symmetry boundary condition | |
Blázquez et al. | Effect of numerical artificial corners appearing when using BEM on contact stresses | |
Dreher et al. | Gerrymandering for interfaces: modeling the mechanics of jointed structures | |
Contractor et al. | Design and analysis of ladder frame chassis considering support at contact region of leaf spring and chassis frame | |
JP6086793B2 (en) | Tire wear simulation method and tire wear simulation program | |
Wang et al. | A laminated ring on elastic foundation model and a feedback contact algorithm for tire-road contact analysis | |
Sasikala et al. | Nonsingular boundary element flexural analysis of stiffened plates | |
Prado et al. | An analysis of parametric instability of risers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220606 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230405 |
|
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: 20230418 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230501 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 7284402 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |