JP7292755B2 - 弾塑性解析方法及び弾塑性解析プログラム - Google Patents
弾塑性解析方法及び弾塑性解析プログラム Download PDFInfo
- Publication number
- JP7292755B2 JP7292755B2 JP2021562425A JP2021562425A JP7292755B2 JP 7292755 B2 JP7292755 B2 JP 7292755B2 JP 2021562425 A JP2021562425 A JP 2021562425A JP 2021562425 A JP2021562425 A JP 2021562425A JP 7292755 B2 JP7292755 B2 JP 7292755B2
- Authority
- JP
- Japan
- Prior art keywords
- plastic
- elastic
- tensor
- increment
- strain
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Description
以下、本発明の実施形態1に係る弾塑性体の解析方法について図2を参照して説明する。弾塑性解析の初期値境界値問題の支配方程式は、運動量保存則、角運動量保存則、変位-ひずみ関係式、応力-ひずみ関係式、塑性ひずみに関する方程式、からなる。従来技術においては、塑性ひずみに関する方程式として流れ則が用いられている。
εe=εtotal-εp …式(2)
である。試行弾性ひずみテンソルとは、繰り返し計算のために試行的に与える微小な弾性ひずみテンソルである。次に、試行弾性ひずみテンソルεeに基づいて一意に定義される試行弾性ひずみ指標εeiを決定する(ステップS14)。試行弾性ひずみ指標εeiはスカラー又はベクトルである。
εei-εcri=dεpi …式(3)
である。ただし、εei≦εcriである場合は、dεpi=0とする。即ち、
dεpi=εei-εcri (εei>εcri) …式(3A)
dεpi=0 (εei≦εcri) …式(3B)
とする。式(3A)および式(3B)の条件式の中の不等号は、指標がベクトル値である場合はベクトルに適切な数学的処理を施して得られるスカラー値に基づいて評価される。
実施形態1で説明した、微小変形の場合における一般的な弾塑性解析方法についての具体的な一例である実施形態2について、図面を参照して説明する。静的状態における弾塑性体の微小変形を考える。実施形態2では、実施形態1で説明した弾性限界ひずみ指標εcriの一例である弾性限界せん断ひずみεcrを考慮する。なお、弾性限界せん断ひずみεcrは、弾塑性体の材料試験により得られる数値である。実施形態2では、弾性限界せん断ひずみεcrを正の定数とする。
(ケース1)
(ケース2)
(ケース2)
(ケース3)
次に実施形態3として、弾塑性体の形状の変形が無視できない有限変形の場合について図面を参照して説明する。この場合の数値解析方法は、次の2つの方法をとることができる。
動的あるいは静的平衡状態にある弾塑性体に、荷重増分を負荷し(図7のステップS30)、塑性変形勾配テンソルは変えずに、荷重増分を負荷した後の擬似的平衡状態における全変形勾配テンソルを得る。この全変形勾配テンソルFtotalに、塑性変形勾配テンソルFpの逆テンソル(Fp)-1を適切に作用させたものを、試行弾性変形勾配テンソルFeとする(ステップS32)。即ち、
Fe=(Fp)-1〇Ftotal …式(23)
である。ここで(Fp)-1〇Ftotalは、塑性ひずみ変形勾配テンソルFpの逆テンソル(Fp)-1を全変形勾配テンソルFtotalに適切に作用させる数学的処理を表す。次に、試行弾性変形勾配テンソルFeに基づいて一意に定義される試行弾性変形勾配指標Feiを決定する(ステップS34)。試行弾性変形勾配指標Feiはスカラー又はベクトルである。
Fei/Fcri=dFpi …式(24)
である。ただし、Fei≦Fcriである場合は、dFpi=0とする。即ち、
dFpi=Fei/Fcri (Fei>Fcri) …式(24A)
dFpi=0 (Fei≦Fcri) …式(24B)
とする。式(24A)および式(24B)の条件式の中の不等号は、指標がベクトル値である場合はベクトルに適切な数学的処理を施して得られるスカラー値に基づいて評価される。
動的あるいは静的平衡状態にある弾塑性体に、荷重増分を負荷し(図8のステップS50)、塑性変形勾配テンソルは変えずに、荷重増分を負荷した後の擬似的平衡状態における全変形勾配テンソルを得る。この全変形勾配テンソルFtotalの対数を取って得られる対数全変形勾配テンソルLnFtotalと、塑性変形勾配テンソルFpの対数を取って得られる対数塑性変形勾配テンソルLnFpとの差に基づいて、対数試行弾性変形勾配テンソルLnFeを決定する。例えば、対数全変形勾配テンソルLnFtotalと対数塑性変形勾配テンソルLnFpとの差を対数試行弾性変形勾配テンソルLnFeとする(ステップS52)。即ち、
LnFe=LnFtotal-LnFp …式(25)
である。対数試行弾性変形勾配テンソルLnFeに基づいて一意に定義される対数試行弾性変形勾配指標LnFeiを決定する(ステップS54)。対数試行弾性変形勾配指標LnFeiはスカラー又はベクトルである。
LnFei-LnFcri=dLnFpi …式(26)
である。ただし、LnFei≦LnFcriである場合は、dLnFpi=0とする。即ち、
dLnFpi=LnFei-LnFcri (LnFei>LnFcri) …式(26A)
dLnFpi=0 (LnFei≦LnFcri) …式(26B)
とする。式(26A)および式(26B)の条件式の中の不等号は、指標がベクトル値である場合はベクトルに適切な数学的処理を施して得られるスカラー値に基づいて評価される。
弾塑性解析方法は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求めるステップと、前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定するステップと、前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標を決定するステップと、前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定するステップと、を含む。これにより、バウシンガー効果と、特定の方向の塑性変形履歴に対する弾塑性変形挙動の不変性と、を同時に再現することができる。
前記試行弾性ひずみ指標から前記弾性限界ひずみ指標を差し引いて得られる値を、前記荷重増分に対応する前記塑性ひずみ増分指標として決定し、前記試行弾性ひずみ指標が前記弾性限界ひずみ指標以下である場合は、前記塑性ひずみ増分指標をゼロとすることを趣旨とする。
弾塑性解析方法は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求めるステップと、前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定するステップと、前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定するステップと、前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定するステップと、決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定するステップと、を含む。
前記試行弾性変形勾配指標と前記弾性限界変形勾配指標との比を前記荷重増分に対応する前記塑性変形勾配増分指標とし、前記比が1以下である場合は、前記塑性変形勾配増分指標をゼロとすることを趣旨とする。
弾塑性解析方法は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求めるステップと、前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定するステップと、前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定するステップと、前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定するステップと、決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定するステップと、を含む。
前記対数試行弾性変形勾配指標から前記対数弾性限界変形勾配指標を差し引いて得られる値を前記荷重増分に対応する前記対数塑性変形勾配増分指標とし、前記対数試行弾性変形勾配指標が前記対数弾性限界変形勾配指標以下である場合は、前記対数塑性変形勾配増分指標をゼロとすることを趣旨とする。
弾塑性解析プログラムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求める手順と、前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定する手順と、前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標として決定する手順と、前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムである。
弾塑性解析プログラムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求める手順と、前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定する手順と、前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定する手順と、前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定する手順と、決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムである。
弾塑性解析プログラムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求める手順と、前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定する手順と、前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定する手順と、前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定する手順と、決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムである。
弾塑性解析システムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析システムであって、動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求める試行弾性ひずみテンソル取得部と、前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定する試行弾性ひずみ指標決定部と、前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標を決定する塑性ひずみ増分指標決定部と、前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定する塑性ひずみ増分テンソル決定部と、を含む弾塑性解析システムである。
弾塑性解析システムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析システムであって、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求める試行弾性変形勾配テンソル取得部と、前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定する試行弾性変形勾配指標決定部と、前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定する塑性変形勾配増分指標決定部と、前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定する塑性変形勾配増分テンソル決定部と、決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定する塑性ひずみ増分テンソル決定部と、を含む弾塑性解析システムである。
弾塑性解析システムは、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析システムであって、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求める対数試行弾性変形勾配テンソル取得部と、前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定する対数試行弾性変形勾配指標決定部と、前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定する対数塑性変形勾配増分指標決定部と、前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定する対数塑性変形勾配増分テンソル決定部と、決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定する対数塑性ひずみ増分テンソル決定部と、を含む弾塑性解析システムである。
記憶媒体は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求める手順と、前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定する手順と、前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標として決定する手順と、前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムを記憶させた、コンピュータ読み取り可能な記憶媒体である。
記憶媒体は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求める手順と、前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定する手順と、前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定する手順と、前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定する手順と、決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムを記憶させた、コンピュータ読み取り可能な記憶媒体である。
記憶媒体は、境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、コンピュータに、動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求める手順と、前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定する手順と、前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定する手順と、前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定する手順と、決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定する手順と、を含む手順を実行させるための弾塑性解析プログラムを記憶させた、コンピュータ読み取り可能な記憶媒体である。
S1,S2:外部せん断荷重
E:引張荷重
C:圧縮荷重
Claims (9)
- 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求めるステップと、
前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定するステップと、
前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標を決定するステップと、
前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定するステップと、
を含むプロセスを実行させる弾塑性解析方法。 - 前記コンピュータに、
前記試行弾性ひずみ指標から前記弾性限界ひずみ指標を差し引いて得られる値を、前記荷重増分に対応する前記塑性ひずみ増分指標として決定し、前記試行弾性ひずみ指標が前記弾性限界ひずみ指標以下である場合は、前記塑性ひずみ増分指標をゼロとするプロセスを実行させる、請求項1に記載の弾塑性解析方法。 - 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求めるステップと、
前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定するステップと、
前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定するステップと、
前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定するステップと、
決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定するステップと、
を含むプロセスを実行させる弾塑性解析方法。 - 前記コンピュータに、
前記試行弾性変形勾配指標と前記弾性限界変形勾配指標との比を前記荷重増分に対応する前記塑性変形勾配増分指標とし、前記比が1以下である場合は、前記塑性変形勾配増分指標をゼロとするプロセスを実行させる、請求項3に記載の弾塑性解析方法。 - 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析方法であって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求めるステップと、
前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定するステップと、
前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定するステップと、
前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定するステップと、
決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定するステップと、
を含むプロセスを実行させる弾塑性解析方法。 - コンピュータに、
前記対数試行弾性変形勾配指標から前記対数弾性限界変形勾配指標を差し引いて得られる値を前記荷重増分に対応する前記対数塑性変形勾配増分指標とし、前記対数試行弾性変形勾配指標が前記対数弾性限界変形勾配指標以下である場合は、前記対数塑性変形勾配増分指標をゼロとするプロセスを実行させる、請求項5に記載の弾塑性解析方法。 - 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性ひずみテンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全ひずみテンソルから前記塑性ひずみテンソルを差し引いて得られる試行弾性ひずみテンソルを求める手順と、
前記試行弾性ひずみテンソルに基づいて試行弾性ひずみ指標を決定する手順と、
前記弾塑性体の材料特性に基づいて決定される弾性限界ひずみ指標を、前記試行弾性ひずみ指標から差し引いて得られる値に基づいて、前記荷重増分に対応する塑性ひずみ増分指標として決定する手順と、
前記塑性ひずみ増分指標に基づいて塑性ひずみ増分テンソルを決定する手順と、
を含む手順を実行させるための弾塑性解析プログラム。 - 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における全変形勾配テンソルに、前記塑性変形勾配テンソルの逆テンソルを作用させる数学的処理をほどこして試行弾性変形勾配テンソルを求める手順と、
前記試行弾性変形勾配テンソルに基づいて試行弾性変形勾配指標を決定する手順と、
前記試行弾性変形勾配指標と、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標との比に基づいて、前記荷重増分に対応する塑性変形勾配増分指標を決定する手順と、
前記塑性変形勾配増分指標に基づいて塑性変形勾配増分テンソルを決定する手順と、
決定された前記塑性変形勾配増分テンソルを用いて塑性ひずみ増分テンソルを決定する手順と、
を含む手順を実行させるための弾塑性解析プログラム。 - 境界値問題を数値的に解く数値解析手法を用いて行う弾塑性体の弾塑性解析プログラムであって、
コンピュータに、
動的又は静的平衡状態にある弾塑性体に塑性変形勾配テンソルは変えずに荷重増分を負荷後の擬似的平衡状態における対数全変形勾配テンソルから対数塑性変形勾配テンソルを差し引いて得られる対数試行弾性変形勾配テンソルを求める手順と、
前記対数試行弾性変形勾配テンソルに基づいて対数試行弾性変形勾配指標を決定する手順と、
前記対数試行弾性変形勾配指標から、前記弾塑性体の材料特性に基づいて決定される弾性限界変形勾配指標の対数を取って得られる対数弾性限界変形勾配指標を差し引いて得られる値に基づいて、前記荷重増分に対応する対数塑性変形勾配増分指標を決定する手順と、
前記対数塑性変形勾配増分指標に基づいて対数塑性変形勾配増分テンソルを決定する手順と、
決定された前記対数塑性変形勾配増分テンソルを用いて対数塑性ひずみ増分テンソルを決定する手順と、
を含む手順を実行させるための弾塑性解析プログラム。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2019/047888 WO2021111625A1 (ja) | 2019-12-06 | 2019-12-06 | 弾塑性解析方法及び弾塑性解析プログラム |
Publications (3)
Publication Number | Publication Date |
---|---|
JPWO2021111625A1 JPWO2021111625A1 (ja) | 2021-06-10 |
JPWO2021111625A5 JPWO2021111625A5 (ja) | 2022-09-14 |
JP7292755B2 true JP7292755B2 (ja) | 2023-06-19 |
Family
ID=76221151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021562425A Active JP7292755B2 (ja) | 2019-12-06 | 2019-12-06 | 弾塑性解析方法及び弾塑性解析プログラム |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7292755B2 (ja) |
WO (1) | WO2021111625A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114912314B (zh) * | 2022-04-21 | 2024-04-02 | 中国科学院武汉岩土力学研究所 | 岩土介质弹塑性本构关系隐式自适应应力积分计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001201408A (ja) | 2000-01-18 | 2001-07-27 | Toyota Motor Corp | 応力振動を抑制する動的陽解法有限要素法 |
JP2008142774A (ja) | 2006-11-14 | 2008-06-26 | Jfe Steel Kk | 応力−ひずみ関係シミュレート方法、応力−ひずみ関係シミュレーションシステム、応力−ひずみ関係シミュレーションプログラム、及び当該プログラムを記録した記録媒体 |
US20150213164A1 (en) | 2014-01-27 | 2015-07-30 | GM Global Technology Operations LLC | Product design reliability with consideration of material property changes during service |
-
2019
- 2019-12-06 WO PCT/JP2019/047888 patent/WO2021111625A1/ja active Application Filing
- 2019-12-06 JP JP2021562425A patent/JP7292755B2/ja active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001201408A (ja) | 2000-01-18 | 2001-07-27 | Toyota Motor Corp | 応力振動を抑制する動的陽解法有限要素法 |
JP2008142774A (ja) | 2006-11-14 | 2008-06-26 | Jfe Steel Kk | 応力−ひずみ関係シミュレート方法、応力−ひずみ関係シミュレーションシステム、応力−ひずみ関係シミュレーションプログラム、及び当該プログラムを記録した記録媒体 |
US20150213164A1 (en) | 2014-01-27 | 2015-07-30 | GM Global Technology Operations LLC | Product design reliability with consideration of material property changes during service |
Also Published As
Publication number | Publication date |
---|---|
WO2021111625A1 (ja) | 2021-06-10 |
JPWO2021111625A1 (ja) | 2021-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2014141794A1 (ja) | 応力-ひずみ関係シミュレート方法、スプリングバック量予測方法およびスプリングバック解析装置 | |
Burchitz | Improvement of springback prediction in sheet metal forming | |
Malag et al. | Problems determining of the mechanical properties of metallic materials from the tensile test in the aspect of numerical calculations of the technological processes | |
Song et al. | Springback prediction in T-section beam bending process using neural networks and finite element method | |
CN111539144A (zh) | 一种含裂纹结构件的断裂韧性计算方法及系统 | |
JP3809374B2 (ja) | 応力−ひずみ関係シミュレート方法および除荷過程における降伏点を求める方法 | |
JP7292755B2 (ja) | 弾塑性解析方法及び弾塑性解析プログラム | |
Kos et al. | The development of prediction model for failure force of damaged reinforced-concrete slender columns | |
Voyiadjis et al. | Effects of stress invariants and reverse loading on ductile fracture initiation | |
CN111539142A (zh) | 一种管道裂纹扩展驱动力计算方法及系统 | |
Patel et al. | N-Layer concentric cylinder model (NCYL): an extended micromechanics-based multiscale model for nonlinear composites | |
Meinders et al. | A sensitivity analysis on the springback behavior of the unconstrained bending problem | |
CN111539071A (zh) | 一种差厚板晶体塑性本构模型建立方法、系统及电子设备 | |
US5437190A (en) | Method for determining the effects of stress | |
JP6325865B2 (ja) | 部品の破断判定方法、システム及びプログラム、並びに、理論成形限界線図の作成方法 | |
WO2015037067A1 (ja) | 延性破壊評価システム及び延性破壊評価方法 | |
JP2009259043A (ja) | ゴム材料のシミュレーション方法 | |
Kutt et al. | Non‐linear finite element analysis of springback | |
Yuan et al. | A configurational force based anisotropic damage model for original isotropic materials | |
Wang et al. | Explicit dynamic analysis of sheet metal forming processes using linear prismatic and hexahedral solid-shell elements | |
Rossi | Mechanical behavior of ferritic grade 3Cr12 stainless steel—Part 2: Yield locus and hardening laws | |
Mendoza et al. | Ultimate local strength of a submarine structure considering the influence of localized reduction of thickness | |
Hora et al. | Damage dependent stress limit model for failure prediction in bulk forming processes | |
Altai et al. | Effect of the post-peak behavior on collapse of structural systems | |
Lazzari et al. | Buckling Analysis, Structural Behavior and Strength of Distortional-Global Interaction in Cold-Formed Steel Lipped Channel Columns |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20220620 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220906 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220930 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230509 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230517 |
|
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: 20230530 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230531 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7292755 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |