JP2022102135A - プレス機の剛性予測方法 - Google Patents
プレス機の剛性予測方法 Download PDFInfo
- Publication number
- JP2022102135A JP2022102135A JP2020216687A JP2020216687A JP2022102135A JP 2022102135 A JP2022102135 A JP 2022102135A JP 2020216687 A JP2020216687 A JP 2020216687A JP 2020216687 A JP2020216687 A JP 2020216687A JP 2022102135 A JP2022102135 A JP 2022102135A
- Authority
- JP
- Japan
- Prior art keywords
- rigidity
- press machine
- press
- side direction
- long side
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005259 measurement Methods 0.000 claims abstract description 45
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000009795 derivation Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 6
- 238000005452 bending Methods 0.000 abstract 1
- 238000005094 computer simulation Methods 0.000 abstract 1
- 239000002184 metal Substances 0.000 abstract 1
- 230000006870 function Effects 0.000 description 38
- 238000009826 distribution Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000000465 moulding Methods 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
【課題】新規なプレス機剛性の予測方法を提供する。【解決手段】プレス機の剛性予測方法は、代表金型を介して前記プレス機に荷重を載荷して、荷重と前記プレス機における複数の測定点での実変位との関係を実測する実測ステップS1と、プレス機の長辺長さ及び短辺長さを反映したプレス機モデルを用意する解析モデル作成ステップS2と、プレス機モデルのたわみと、荷重との関係を、プレス機の剛性を示す剛性係数を用いた関数で表す理論計算モデル作成ステップS3と、実変位及び荷重に基づいて、剛性係数を導出する剛性係数導出ステップS4と、を含む。【選択図】図4
Description
本発明は、プレス機の剛性予測方法に関する。
従来、ワーク成形時におけるプレス金型に見込まれる変形をプレス金型の設計に反映させるため、ワーク成形時におけるプレス機に装着されたプレス金型の局所の変形を実測し、プレス機とプレス金型とから構成されたプレス成形装置の剛性を予測する方法があった。
しかしながら、従来の方法は、特定のプレス金型を使用した場合の、その特定のプレス金型の局所の変形を実測する方法であるので、プレス機が同じであっても、形状の異なるプレス金型を使用する場合、プレス機の剛性を予測するには、その異なるプレス金型毎にプレス金型の変形を実測する必要があった。
また、従来、プレス機の面方向に沿って連続的な剛性分布を予測する方法はなかった。
しかしながら、従来の方法は、特定のプレス金型を使用した場合の、その特定のプレス金型の局所の変形を実測する方法であるので、プレス機が同じであっても、形状の異なるプレス金型を使用する場合、プレス機の剛性を予測するには、その異なるプレス金型毎にプレス金型の変形を実測する必要があった。
また、従来、プレス機の面方向に沿って連続的な剛性分布を予測する方法はなかった。
本発明は、上記背景技術の問題点に鑑み、新規なプレス機の剛性予測方法を提供することを課題とする。
本発明の要旨は以下の通りである。
(1)本発明の一態様に係るプレス機の剛性予測方法は、代表金型を介して前記プレス機に荷重を載荷して、前記荷重と前記プレス機における複数の測定点での実変位との関係を実測する実測ステップと、前記プレス機の長辺長さ及び短辺長さを反映したプレス機モデルを用意する解析モデル作成ステップと、前記プレス機モデルのたわみと、前記荷重との関係を、前記プレス機の剛性を示す剛性係数を用いた関数で表す理論計算モデル作成ステップと、前記実変位及び前記荷重に基づいて、前記剛性係数を導出する剛性係数導出ステップと、を含む。
(2)上記(1)において、前記代表金型の載荷面積と、指定のプレス金型の載荷面積との比に基づいて、前記剛性係数を補正し、前記指定のプレス金型を備えた前記プレス機の補正剛性係数を算出する補正ステップを含んでよい。
(3)上記(1)又は(2)において、前記関数は、二次元フーリエ級数展開されたものであってよい。
(4)上記(3)において、前記関数は、正弦関数によって二次元フーリエ級数展開されたものであってよい。
(5)上記(3)において、前記関数は、ルジャンドル多項式によって二次元フーリエ級数展開されたものであってよい。
(6)上記(3)において、前記関数は、ラゲール多項式によって二次元フーリエ級数展開されたものであってよい。
(7)上記(3)において、前記関数は、エルミート多項式によって二次元フーリエ級数展開されたものであってよい。
(8)上記(4)において、前記関数は、式(A)で表されてよい。
・・・式(A)
ここで、F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=1,3,...)
rn b:短辺方向スケーリングパラメータ(n=1,3,...)
(9)上記(5)において、前記関数は、式(B)で表されてよい。
・・・式(B)
ここで、Pm:m次のルジャンドル多項式
Pn:n次のルジャンドル多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(10)上記(6)において、前記関数は、式(C)で表されてよい。
・・・式(C)
ここで、Lm:m次のラゲール多項式
Ln:n次のラゲール多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amnは、剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(11)上記(7)において、前記関数は、式(D)で表されてよい。
・・・式(D)
ここで、Hm:m次のエルミート多項式
Hn:n次のエルミート多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(12)上記(3)から(11)のいずれかにおいて、前記二次元フーリエ級数展開の級数は、有限であってよい。
(2)上記(1)において、前記代表金型の載荷面積と、指定のプレス金型の載荷面積との比に基づいて、前記剛性係数を補正し、前記指定のプレス金型を備えた前記プレス機の補正剛性係数を算出する補正ステップを含んでよい。
(3)上記(1)又は(2)において、前記関数は、二次元フーリエ級数展開されたものであってよい。
(4)上記(3)において、前記関数は、正弦関数によって二次元フーリエ級数展開されたものであってよい。
(5)上記(3)において、前記関数は、ルジャンドル多項式によって二次元フーリエ級数展開されたものであってよい。
(6)上記(3)において、前記関数は、ラゲール多項式によって二次元フーリエ級数展開されたものであってよい。
(7)上記(3)において、前記関数は、エルミート多項式によって二次元フーリエ級数展開されたものであってよい。
(8)上記(4)において、前記関数は、式(A)で表されてよい。
ここで、F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=1,3,...)
rn b:短辺方向スケーリングパラメータ(n=1,3,...)
(9)上記(5)において、前記関数は、式(B)で表されてよい。
ここで、Pm:m次のルジャンドル多項式
Pn:n次のルジャンドル多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(10)上記(6)において、前記関数は、式(C)で表されてよい。
ここで、Lm:m次のラゲール多項式
Ln:n次のラゲール多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amnは、剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(11)上記(7)において、前記関数は、式(D)で表されてよい。
ここで、Hm:m次のエルミート多項式
Hn:n次のエルミート多項式
F:荷重
x:測定点における長辺方向の座標
y:測定点における短辺方向の座標
Amn:剛性係数
a:長辺長さ
b:短辺長さ
rm a:長辺方向スケーリングパラメータ(m=0,1,...)
rn b:短辺方向スケーリングパラメータ(n=0,1,...)
(12)上記(3)から(11)のいずれかにおいて、前記二次元フーリエ級数展開の級数は、有限であってよい。
本発明によれば、新規なプレス機の剛性予測方法を提供できる。
図1に示すように、不図示のブランク(例えば、金属材料)に強い圧力を加え、ブランクをプレス金型の形状に変形させてプレス成形品を製造する機械である一般的なプレス成形装置1は、基本的には、プレス機10と、プレス機10に取り付けられたプレス金型20とで構成されている。
プレス機10は、下金型21が上部に取り付けられるボルスタ11と、上金型22が下部に取り付けられるスライド12と、ボルスタ11とスライド12とを近づけたり遠ざけたりするように駆動する駆動装置13とを含んでいる。
プレス金型20は、下金型21及び上金型22を含んでいる。
プレス機10の構造は、そのプレス機10のプレス能力、仕様、設計コンセプト等によって異なる。プレス能力、仕様、設計コンセプト等が同じであっても、プレス機10の構造には、製造誤差等に起因する、機体差がある。また、プレス機10に取り付けられるプレス金型20の構造は、製造するプレス成形品の仕様により異なる。したがって、プレス機10の構造、プレス金型20の構造の相違により、プレス荷重が同じであっても、プレス機10の剛性が異なるので、プレス機10の変形の状態は異なる。
ところで、プレス成形によってプレス成形品を精度良く製造するため、プレス成形品のブランクに接する成形面を含むプレス金型20の形状に、あらかじめプレス機10のたわみω(変位の面的な分布を意味するが、以下、変位又は変形という場合もある。)を見込んでおくことが必要となる場合がある。この場合は、荷重F(プレス荷重)とプレス機10のたわみωとの関係であるプレス機10の剛性を把握しておくことが好ましい。
しかしながら、プレス機10の剛性を把握するため、プレス機10のたわみωを精度良く実測するには、測定点が多くなるので、相応の設備等のコストを要する。また、プレス機10の剛性を、FEM等により計算しても、実際の剛性との乖離があったり、計算の負荷が過大となったりする。
しかしながら、プレス機10の剛性を把握するため、プレス機10のたわみωを精度良く実測するには、測定点が多くなるので、相応の設備等のコストを要する。また、プレス機10の剛性を、FEM等により計算しても、実際の剛性との乖離があったり、計算の負荷が過大となったりする。
そこで、本発明は、代表金型を介してプレス機に荷重を載荷して、荷重とプレス機における複数の測定点での実変位との関係を実測する実測ステップと、プレス機の長辺長さ及び短辺長さを反映したプレス機モデルを用意する解析モデル作成ステップと、プレス機モデルのたわみと、荷重との関係を、プレス機の剛性を示す剛性係数を用いた関数で表す理論計算モデル作成ステップと、実変位及び荷重に基づいて、剛性係数を導出する剛性係数導出ステップと、を含むプレス機の剛性予測方法である。これにより、プレス機剛性を、小さな計算負荷で、実際の剛性との乖離を小さく予測できる。よって、新規なプレス機剛性の予測方法を提供できる。
以下、本発明の実施形態を詳細に説明する。
以下、本発明の実施形態を詳細に説明する。
(実施形態)
図1は、プレス機10とプレス金型20とを備えたプレス成形装置1の説明図である。図2は、プレス機モデル100の説明図である。図3は、図1におけるA矢視断面図であり、平面視におけるプレス機10及びプレス金型20の外形状の説明図である。図4は、プレス機剛性の予測方法のフロー図である。
以下、本実施形態に係るプレス機剛性の予測方法を、図4に示すフロー図に沿って説明する。
図1は、プレス機10とプレス金型20とを備えたプレス成形装置1の説明図である。図2は、プレス機モデル100の説明図である。図3は、図1におけるA矢視断面図であり、平面視におけるプレス機10及びプレス金型20の外形状の説明図である。図4は、プレス機剛性の予測方法のフロー図である。
以下、本実施形態に係るプレス機剛性の予測方法を、図4に示すフロー図に沿って説明する。
(プレス機剛性の予測方法)
図4に示すように、本実施形態に係るプレス機剛性の予測方法は、実測ステップS1と、解析モデル作成ステップS2と、理論計算モデル作成ステップS3と、剛性係数導出ステップS4と、を含んでいる。
図4に示すように、本実施形態に係るプレス機剛性の予測方法は、実測ステップS1と、解析モデル作成ステップS2と、理論計算モデル作成ステップS3と、剛性係数導出ステップS4と、を含んでいる。
前提として、まず、プレス機10を用意する。プレス機10は、上述のような一般的なプレス成形装置1を構成するものであり、スライド12及びボルスタ11を備えている。スライド12又はボルスタ11は、長辺長さa及び短辺長さbを有する平面視において矩形状の版状体である。長辺長さaと短辺長さbとは、同じであっても、異なっていてもよい。すなわち、スライド12又はボルスタ11は、平面視において、長方形又は正方形の外形を有していてよい。なお、スライド12とボルスタ11の長辺長さ又は短辺長さは、一致しなくてもよい。
スライド12は、上金型22を着脱自在な実質的に平坦な下面を有している。ボルスタ11は、下金型21を着脱自在な実質的に平坦な上面を有している。すなわち、スライド12の下面は、平面視において長辺長さa及び短辺長さbを有する矩形状である。
同様に、ボルスタ11の上面は、平面視において長辺長さa及び短辺長さbを有する矩形状である。
スライド12は、上金型22を着脱自在な実質的に平坦な下面を有している。ボルスタ11は、下金型21を着脱自在な実質的に平坦な上面を有している。すなわち、スライド12の下面は、平面視において長辺長さa及び短辺長さbを有する矩形状である。
同様に、ボルスタ11の上面は、平面視において長辺長さa及び短辺長さbを有する矩形状である。
また、プレス機10に取り付ける代表金型を用意する。代表金型は、プレス機10の変形を測定するためのプレス金型20である。
代表金型は、スライド12に取り付けられる上金型22と、ボルスタ11に取り付けられる下金型21とを備えている。
代表金型は、スライド12に取り付けられる上金型22と、ボルスタ11に取り付けられる下金型21とを備えている。
図3に示すように、上金型22及び下金型21は、平面視において金型長辺長さα及び金型短辺長さβを有する外形である。上金型22は、スライド12の下面に対向して面接触し、金型長辺長さα及び金型短辺長さβを有する取付面を有している。下金型21は、ボルスタ11の上面に対向して面接触し、金型長辺長さα及び金型短辺長さβを有する取付面を有している。なお、上金型22と下金型21の長辺長さ又は短辺長さは、一致しなくてもよい。
平面視において、プレス機10の外形と、代表金型(プレス金型20)の外形とは、相似であってよい。すなわち、長辺長さaと短辺長さbの比は、金型長辺長さαと金型短辺長さβの比と等しくてよい。なお、例えば、a/αは、3/1であり、b/βは、3/1である。
続いて、代表金型をプレス機10に取り付ける。具体的には、プレス機10のスライド12の下面に、代表金型の上金型22の取付面を対向させて面接触させた状態で、上金型22を取り付ける。プレス機10のボルスタ11の上面に、代表金型の下金型21の取付面を対向させて面接触させた状態で、下金型21を取り付ける。
(1)そして、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点の座標(xi,yj)での実変位ωijとの関係を実測する(実測ステップS1)。以下、荷重Fのときのプレス機10における複数の座標(xi,yj)での実変位ωijを、ω(xi,yj,F)という場合がある。また、以下、荷重Fのときのプレス機モデル100における複数の座標(xi,yj)での仮想変位ωωijを、ωω(xi,yj,F)という場合がある。
なお、通常、プレス機10のスライド12とボルスタ11とは、板厚、補強構造等が異なり、剛性係数Amn(剛性の分布)が異なるので、この実測は、プレス機10のスライド12と、ボルスタ11とで、個別に実施する。以下の説明においても、特に説明のない限り、「プレス機」は、スライド12又はボルスタ11のいずれかを意味する場合がある。
荷重Fは、鉛直方向に作用する力である。
複数の測定点の座標(xi,yj)は、長辺方向(x方向)にi箇所、短辺方向(y方向)にj箇所設けられる。例えば、長辺方向に2つのラインに沿ってそれぞれ10箇所と、短辺方向に2つのラインに沿ってそれぞれ10箇所であれば、合計40箇所の測定点である。測定点は、図心を含む図心近傍の1点と、図心を通る長辺方向に沿う2点と、図心を通る短辺方向に沿う2点と、を含む、少なくとも5点以上である。現実的には、測定点は、長辺と短辺の近傍付近と、代表金型と対向するプレス機10(スライド12又はボルスタ11)の面を含むように長辺方向及び短辺方向に沿って複数箇所のように、比較的大きな変形が見込まれる箇所を含み、プレス機10の面の変形がわかるように、プレス機10の面上で偏在することなくバランスよく対称的に設ける。
実変位ωijは、測定点の座標(xi,yj)において、鉛直方向の変位を測定可能な変位計で測定する。変位計は、例えば、ばねを利用した接触型又はレーザを利用した非接触型の変位計であってよく、撮像装置で撮像した撮像画像の変化をコンピュータプログラムで解析することによって変位を算出するものであってもよい。
なお、通常、プレス機10のスライド12とボルスタ11とは、板厚、補強構造等が異なり、剛性係数Amn(剛性の分布)が異なるので、この実測は、プレス機10のスライド12と、ボルスタ11とで、個別に実施する。以下の説明においても、特に説明のない限り、「プレス機」は、スライド12又はボルスタ11のいずれかを意味する場合がある。
荷重Fは、鉛直方向に作用する力である。
複数の測定点の座標(xi,yj)は、長辺方向(x方向)にi箇所、短辺方向(y方向)にj箇所設けられる。例えば、長辺方向に2つのラインに沿ってそれぞれ10箇所と、短辺方向に2つのラインに沿ってそれぞれ10箇所であれば、合計40箇所の測定点である。測定点は、図心を含む図心近傍の1点と、図心を通る長辺方向に沿う2点と、図心を通る短辺方向に沿う2点と、を含む、少なくとも5点以上である。現実的には、測定点は、長辺と短辺の近傍付近と、代表金型と対向するプレス機10(スライド12又はボルスタ11)の面を含むように長辺方向及び短辺方向に沿って複数箇所のように、比較的大きな変形が見込まれる箇所を含み、プレス機10の面の変形がわかるように、プレス機10の面上で偏在することなくバランスよく対称的に設ける。
実変位ωijは、測定点の座標(xi,yj)において、鉛直方向の変位を測定可能な変位計で測定する。変位計は、例えば、ばねを利用した接触型又はレーザを利用した非接触型の変位計であってよく、撮像装置で撮像した撮像画像の変化をコンピュータプログラムで解析することによって変位を算出するものであってもよい。
(2)また、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する(解析モデル作成ステップS2)。図2に示すように、プレス機モデル100は、プレス機10(スライド12又はボルスタ11)を模擬した仮想の解析モデルである。プレス機モデル100は、四辺(二つの長辺及び二つの短辺、すなわち、周縁部)を支持された、長辺長さa及び短辺長さbを有する、一様断面(厚さ一定)の、中実な直方体である。プレス機モデル100には、金型長辺長さα及び金型短辺長さβを有するプレス金型モデル200からの反力が作用する。
(3)次に、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す(理論計算モデル作成ステップS3)。
ここで、プレス機モデル100は、上述のように長辺長さa及び短辺長さbを有する、一様断面(厚さ一定)の、中実な直方体であって、四辺(二つの長辺及び二つの短辺)を支持されたものである。ここで、プレス金型におけるプレス機10と対向する取付面は、実際、比較的単純な放物面形状であることから、プレス金型からプレス機10への荷重の伝達にサンヴナンの原理を適用し、プレス機モデル100にその原理を反映して計算したほうがより実際に近いたわみωを計算できると考えられる。プレス機モデル100のような矩形状の版状体のたわみは、長辺方向及び短辺方向の直交する2方向への対称性がある。したがって、長辺方向及び短辺方向のそれぞれに対して、長辺長さa又は短辺長さbを奇数で除した値が半周期となる波を重ね合わせたような形状になることが明らかであることから、式(1)のような、二次元フーリエ級数展開された関数である、一様断面四辺支持平板たわみの解析解により、たわみωを計算することができる。
・・・式(1)
ここで、t:プレス機の版厚
E:プレス機のヤング係数
ν:ポワソン比
ここで、プレス機モデル100は、上述のように長辺長さa及び短辺長さbを有する、一様断面(厚さ一定)の、中実な直方体であって、四辺(二つの長辺及び二つの短辺)を支持されたものである。ここで、プレス金型におけるプレス機10と対向する取付面は、実際、比較的単純な放物面形状であることから、プレス金型からプレス機10への荷重の伝達にサンヴナンの原理を適用し、プレス機モデル100にその原理を反映して計算したほうがより実際に近いたわみωを計算できると考えられる。プレス機モデル100のような矩形状の版状体のたわみは、長辺方向及び短辺方向の直交する2方向への対称性がある。したがって、長辺方向及び短辺方向のそれぞれに対して、長辺長さa又は短辺長さbを奇数で除した値が半周期となる波を重ね合わせたような形状になることが明らかであることから、式(1)のような、二次元フーリエ級数展開された関数である、一様断面四辺支持平板たわみの解析解により、たわみωを計算することができる。
ここで、t:プレス機の版厚
E:プレス機のヤング係数
ν:ポワソン比
さらに、式(2)においてm及びnが大きくなるにつれてAmnは無視できるほど小さくなることから、式(2)を有限級数化すると、下式のようになる。
ここで、長辺長さaに応じて設定する補正値である長辺方向スケーリングパラメータrm
aと、短辺長さbに応じて設定する補正値である短辺方向スケーリングパラメータrn
bとを考慮して、次のように置き換える。
x/a → rm a・x/a,
y/b → rn b・y/b
すると、式(2)は、式(A)で表すことができる。
・・・式(A)
x/a → rm a・x/a,
y/b → rn b・y/b
すると、式(2)は、式(A)で表すことができる。
なお、実際のプレス機10のたわみω又は剛性係数Amnを計算する場合、たわみωの幾何学的に対称な形状を表すため、M及びNは奇数とする。M及びNは、計算負荷と計算精度のバランスをみて、それぞれ、15以下で十分であり、7以下であってよい。
このように、たわみωは、正弦関数によって二次元フーリエ級数展開された関数で表すことができる。これにより、解析解を基本とし、スケーリングパラメータが一定の場合、代数的解法で計算できるので、短い解析時間で、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
このように、たわみωは、正弦関数によって二次元フーリエ級数展開された関数で表すことができる。これにより、解析解を基本とし、スケーリングパラメータが一定の場合、代数的解法で計算できるので、短い解析時間で、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
(4)そして、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する(剛性係数導出ステップS4)。
例えば、スケーリングパラメータが一定の場合、式(A)の代わりに、各座標(xi,yj)におけるたわみωを表す行列を、剛性係数、x方向の変形モード及びy方向の変形モードの各行列の積と荷重Fとの積の形式で表現し、特異値分解を用いて剛性係数Amnを導出できる。
例えば、スケーリングパラメータが一定の場合、式(A)の代わりに、各座標(xi,yj)におけるたわみωを表す行列を、剛性係数、x方向の変形モード及びy方向の変形モードの各行列の積と荷重Fとの積の形式で表現し、特異値分解を用いて剛性係数Amnを導出できる。
一般的には、図5に示すように、次のようなフローに従って、収束計算により、剛性係数Amnを導出する。図5は、剛性係数導出ステップS4のフロー図である。
(4-1)まず、式(4)に示すような仮剛性係数Amnt、仮長辺方向スケーリングパラメータrmt a及び仮短辺方向スケーリングパラメータrnt bを設定する。そして、初期値として、仮剛性係数Amntの初期パラメータAmn0と、初期長辺方向スケーリングパラメータrm0 aと、初期短辺方向スケーリングパラメータrn0 bと、を設定する。
・・・式(4)
(4-1)まず、式(4)に示すような仮剛性係数Amnt、仮長辺方向スケーリングパラメータrmt a及び仮短辺方向スケーリングパラメータrnt bを設定する。そして、初期値として、仮剛性係数Amntの初期パラメータAmn0と、初期長辺方向スケーリングパラメータrm0 aと、初期短辺方向スケーリングパラメータrn0 bと、を設定する。
(4-2)次に、剛性係数Amnの代わりに仮剛性係数Amntを代入した式(A)に、測定点の座標(x,y)を代入する。これを複数の測定点の座標(xi,yj)ごとに行う。すなわち、式(5)のように、プレス機モデル100の測定点の座標(xi,yj)における変位である仮想変位ωωij(=ωω(xi,yj,F))を計算する。
・・・式(5)
(4-3)次に、実変位ωijと仮想変位ωωijとの残差eを評価する。
残差eは、式(6)の第1項のように、実変位ωijと仮想変位ωωijとの差の二乗である二乗誤差として求められる。また、式(5)算出される仮想変位ωωijと荷重Fとの関係は、比較的強い非線形性を有するため、残差eは、計算の再現性及び安定性を高めるため、式(6)の第2項のような、係数ペナルティ項を含んでもよい。
・・・式(6)
ここで、λは、ペナルティ係数である。
残差eは、式(6)の第1項のように、実変位ωijと仮想変位ωωijとの差の二乗である二乗誤差として求められる。また、式(5)算出される仮想変位ωωijと荷重Fとの関係は、比較的強い非線形性を有するため、残差eは、計算の再現性及び安定性を高めるため、式(6)の第2項のような、係数ペナルティ項を含んでもよい。
ここで、λは、ペナルティ係数である。
(4-4)ここで、式(6)で計算された残差eが、所定値tol以上となる場合、前回設定した仮剛性係数Amntの近傍探索を行う。すなわち、前回設定した仮剛性係数Amntと比べて少しだけ(例えば、0.1%)異ならせた新たな仮剛性係数Amntを設定し、再び、式(4)から式(6)を計算する。これを、式(6)で計算された残差eが、所定値tol未満となるまで繰り返し計算する(近傍探索の繰り返し計算ループ)。
一方、式(6)で計算された残差eが、所定値tol未満となる場合、計算を終了する。このときの仮剛性係数Amntを剛性係数Amnとする。
このようにして、剛性係数Amnを導出できる。
一方、式(6)で計算された残差eが、所定値tol未満となる場合、計算を終了する。このときの仮剛性係数Amntを剛性係数Amnとする。
このようにして、剛性係数Amnを導出できる。
(4-5)また、式(5)算出される仮想変位ωωijと荷重Fとの関係は、比較的強い非線形性を有するため、適切でない局所最適解に陥ることを避けることを目的として、(4-4)で述べた近傍探索の繰り返し計算ループの外側に、大域的最適解に近づけることを目的として、大変動の繰り返し計算ループを更に設けることが好ましい(例えば、Basinhopping法)。すなわち、前回設定した仮剛性係数Amntと比べて大きく(例えば、5%)異ならせた新たな仮剛性係数Amntを設定し、再び、式(4)から式(6)を計算する。
(5)そして、適宜、導出した剛性係数Amnを用いた式(A)により、実変位ωijを計算する。
本実施形態によれば、上述のような各ステップ(実測ステップ、解析モデル作成ステップ、理論計算モデル作成ステップ、剛性係数導出ステップ)を含むので、プレス機10の代表点(格子点)における荷重と変位の関係のみならず、プレス機10の面(x-y平面)に沿って連続的な剛性分布を得ることができる。また、プレス機10の剛性(剛性係数Amn)を小さな計算負荷で、実際の剛性との乖離を小さく予測できる。また、このような各ステップを経て計算された剛性係数Amnを用いることにより、ワークに接する成形面を含むプレス金型の形状の設計に、プレス機10のたわみ(変形)をあらかじめ精度良く見込んでおくことができる。よって、プレス成形装置によってワークを精度良く製造できる。
ところで、代表金型と指定のプレス金型の形状が異なる場合、すなわち、代表金型の載荷面積α*βと、指定のプレス金型の載荷面積α2*β2とが異なる場合、例えば、代表金型の載荷面積α*βと、指定のプレス金型の載荷面積α2*β2との比に基づいて、剛性係数Amnを補正し、指定のプレス金型を備えたプレス機10の補正剛性係数を算出してもよい(補正ステップ)。これにより、代表金型とは異なる形状のプレス金型取り付けたプレス機10の剛性を示す剛性係数Amnを、適切に簡単に補正することで計算できる。よって、プレス金型が変わっても、プレス機10の剛性を適切に予測できる。
(ルジャンドル多項式による二次元フーリエ級数展開された関数)
ところで、上述の式(A)は、正弦関数によって二次元フーリエ級数展開された関数であったが、これに換えて、式(B)のような、ルジャンドル多項式によって二次元フーリエ級数展開された関数を用いてもよい。
・・・式(B)
これにより、放物線を示す項の多項式でたわみωを表すことができるので、比較的少ないパラメータで、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
ところで、上述の式(A)は、正弦関数によって二次元フーリエ級数展開された関数であったが、これに換えて、式(B)のような、ルジャンドル多項式によって二次元フーリエ級数展開された関数を用いてもよい。
これにより、放物線を示す項の多項式でたわみωを表すことができるので、比較的少ないパラメータで、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
(ラゲール多項式による二次元フーリエ級数展開された関数)
ところで、上述の式(A)又は式(B)の関数に換えて、式(C)のような、ラゲール多項式によって二次元フーリエ級数展開された関数を用いてもよい。
・・・式(C)
これにより、指数曲線を示す項の多項式でたわみωを表すことができるので、プレス機10の中心から周縁部に向かって離れるに連れて指数関数的にたわみ難く、変形し難くなるような、プレス機10の剛性の分布を算出できる。よって、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
ところで、上述の式(A)又は式(B)の関数に換えて、式(C)のような、ラゲール多項式によって二次元フーリエ級数展開された関数を用いてもよい。
これにより、指数曲線を示す項の多項式でたわみωを表すことができるので、プレス機10の中心から周縁部に向かって離れるに連れて指数関数的にたわみ難く、変形し難くなるような、プレス機10の剛性の分布を算出できる。よって、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
(エルミート多項式による二次元フーリエ級数展開された関数)
ところで、上述の式(A)、式(B)又は式(C)の関数に換えて、式(D)のような、エルミート多項式によって二次元フーリエ級数展開された関数を用いてもよい。
・・・式(D)
これにより、二次の指数曲線を示す項の多項式でたわみωを表すことができるので、プレス機10の中心から周縁部に向かって離れるに連れて指数関数的にたわみ難く、変形し難くなるような、プレス機10の剛性の分布を算出できる。よって、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
ところで、上述の式(A)、式(B)又は式(C)の関数に換えて、式(D)のような、エルミート多項式によって二次元フーリエ級数展開された関数を用いてもよい。
これにより、二次の指数曲線を示す項の多項式でたわみωを表すことができるので、プレス機10の中心から周縁部に向かって離れるに連れて指数関数的にたわみ難く、変形し難くなるような、プレス機10の剛性の分布を算出できる。よって、実際に近いたわみωを算出でき、実際に近い剛性係数Amnを予測できる。
(プレス機剛性の予測プログラム)
次に、上述のプレス機剛性の予測方法を、コンピュータに実行させるプレス機剛性の予測プログラムについて説明する。
プレス機剛性の予測プログラムは、コンピュータに、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点での実変位ωijとの関係を入力する実測値入力ステップと、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する解析モデル作成ステップと、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す理論計算モデル作成ステップと、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する剛性係数導出ステップと、剛性係数Amnを出力する剛性係数出力ステップと、を実行させるものである。
ここで、プレス機剛性の予測プログラムによって出力された剛性係数Amnは、コンピュータと連携するメモリに保存される。メモリに保存された剛性係数Amnは、モニタに表示してもよく、プリンタに印刷することでもよい。剛性係数Amnは、ユーザが視覚的に捉えられるように、例えば、荷重Fの増減に応じたプレス機10のたわみ(変形)の変化を三次元的に表した図を、モニタに表示できる。
このように、プレス機剛性の予測プログラムは、剛性係数Amnを出力できるので、CAEに利用でき、プレス機10の変形を見込んでプレス金型20を設計できる。よって、プレス機10の変形を見込んで設計されたプレス金型20により、プレス成形品を精度良く製造できる。
次に、上述のプレス機剛性の予測方法を、コンピュータに実行させるプレス機剛性の予測プログラムについて説明する。
プレス機剛性の予測プログラムは、コンピュータに、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点での実変位ωijとの関係を入力する実測値入力ステップと、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する解析モデル作成ステップと、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す理論計算モデル作成ステップと、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する剛性係数導出ステップと、剛性係数Amnを出力する剛性係数出力ステップと、を実行させるものである。
ここで、プレス機剛性の予測プログラムによって出力された剛性係数Amnは、コンピュータと連携するメモリに保存される。メモリに保存された剛性係数Amnは、モニタに表示してもよく、プリンタに印刷することでもよい。剛性係数Amnは、ユーザが視覚的に捉えられるように、例えば、荷重Fの増減に応じたプレス機10のたわみ(変形)の変化を三次元的に表した図を、モニタに表示できる。
このように、プレス機剛性の予測プログラムは、剛性係数Amnを出力できるので、CAEに利用でき、プレス機10の変形を見込んでプレス金型20を設計できる。よって、プレス機10の変形を見込んで設計されたプレス金型20により、プレス成形品を精度良く製造できる。
(プレス機剛性の予測装置)
次に、上述のプレス機剛性の予測プログラムを搭載した、プレス機剛性の予測装置について説明する。
プレス機剛性の予測装置は、一般的なコンピュータと同様に、ストレージ、プロセッサ、メモリ、インタフェース及び出力装置を備えている。
詳細には、プレス機剛性の予測装置は、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点での実変位ωijとの関係を入力する実測値入力部と、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する解析モデル作成部と、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す理論計算モデル作成部と、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する剛性係数導出部と、剛性係数Amnを出力する剛性係数出力部と、を備えている。
ここで、剛性係数出力部は、剛性係数Amnを保存するメモリ、剛性係数Amnを表示するモニタ及び剛性係数Amnを印刷するプリンタを含む。
このように、プレス機剛性の予測装置は、剛性係数Amnを出力する剛性係数出力部を備えているので、CAEに利用でき、プレス機10の変形を見込んでプレス金型20を設計できる。よって、プレス機10の変形を見込んで設計されたプレス金型20により、プレス成形品を精度良く製造できる。
次に、上述のプレス機剛性の予測プログラムを搭載した、プレス機剛性の予測装置について説明する。
プレス機剛性の予測装置は、一般的なコンピュータと同様に、ストレージ、プロセッサ、メモリ、インタフェース及び出力装置を備えている。
詳細には、プレス機剛性の予測装置は、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点での実変位ωijとの関係を入力する実測値入力部と、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する解析モデル作成部と、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す理論計算モデル作成部と、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する剛性係数導出部と、剛性係数Amnを出力する剛性係数出力部と、を備えている。
ここで、剛性係数出力部は、剛性係数Amnを保存するメモリ、剛性係数Amnを表示するモニタ及び剛性係数Amnを印刷するプリンタを含む。
このように、プレス機剛性の予測装置は、剛性係数Amnを出力する剛性係数出力部を備えているので、CAEに利用でき、プレス機10の変形を見込んでプレス金型20を設計できる。よって、プレス機10の変形を見込んで設計されたプレス金型20により、プレス成形品を精度良く製造できる。
(その他の実施形態)
以上、図面を参照して実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
例えば、上述の二次元フーリエ級数展開された関数(式(A)、式(B)、式(C)及び式(D))は、長辺方向の基底関数と短辺方向の基底関数とで同じ種類となる多項式(例えば、式(B)はルジャンドル関数とルジャンドル関数とを組み合わせた多項式)であるが、これに限らない。例えば、長辺方向の基底関数をルジャンドル関数(式(B)参照)として短辺方向の基底関数をエルミート関数(式(D)参照)とする等、長辺方向の基底関数と短辺方向の基底関数とで異なる種類の基底関数を組み合わせた多項式であってよい。
以上、図面を参照して実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
例えば、上述の二次元フーリエ級数展開された関数(式(A)、式(B)、式(C)及び式(D))は、長辺方向の基底関数と短辺方向の基底関数とで同じ種類となる多項式(例えば、式(B)はルジャンドル関数とルジャンドル関数とを組み合わせた多項式)であるが、これに限らない。例えば、長辺方向の基底関数をルジャンドル関数(式(B)参照)として短辺方向の基底関数をエルミート関数(式(D)参照)とする等、長辺方向の基底関数と短辺方向の基底関数とで異なる種類の基底関数を組み合わせた多項式であってよい。
本実施形態に係るプレス機10の剛性予測方法は、代表金型を介してプレス機10に荷重Fを載荷して、荷重Fとプレス機10における複数の測定点での実変位ωijとの関係を実測する実測ステップを含んでいる。また、プレス機10の長辺長さa及び短辺長さbを反映したプレス機モデル100を用意する解析モデル作成ステップを含んでいる。また、プレス機モデル100のたわみωωと、荷重Fとの関係を、プレス機10の剛性を示す剛性係数Amnを用いた関数で表す理論計算モデル作成ステップを含んでいる。また、実変位ωij及び荷重Fに基づいて、剛性係数Amnを導出する剛性係数導出ステップと、を含んでいる。本実施形態に係るプレス機10の剛性予測方法は、このような各ステップを経るので、プレス機剛性を、小さな計算負荷で、実際の剛性との乖離を小さく予測できる。また、プレス機10における実測したいくつかの代表点における剛性のみならず、プレス機10の面方向に沿って連続的な(面的な)剛性分布を予測できる、新規なプレス機剛性の予測方法を提供できる。
1 プレス成形装置
10 プレス機
11 ボルスタ
12 スライド
13 駆動装置
20 プレス金型
21 下金型
22 上金型
100 プレス機モデル
200 プレス金型モデル
a 長辺長さ
b 短辺長さ
Amn 剛性係数
Amn0 初期パラメータ
Amnt 仮剛性係数
e 残差
F 荷重
p0 圧力
S1 実測ステップ
S2 解析モデル作成ステップ
S3 理論計算モデル作成ステップ
S4 剛性係数導出ステップ
tol 所定値
(x,y)、(xi,yj) 座標
α 金型長辺長さ
β 金型短辺長さ
ω (プレス機の)たわみ
ωω (プレス機モデルの)たわみ
ωij 実変位
ωωij 仮想変位
rm a 長辺方向スケーリングパラメータ
rn b 短辺方向スケーリングパラメータ
rm0 a 初期長辺方向スケーリングパラメータ
rn0 b 初期短辺方向スケーリングパラメータ
10 プレス機
11 ボルスタ
12 スライド
13 駆動装置
20 プレス金型
21 下金型
22 上金型
100 プレス機モデル
200 プレス金型モデル
a 長辺長さ
b 短辺長さ
Amn 剛性係数
Amn0 初期パラメータ
Amnt 仮剛性係数
e 残差
F 荷重
p0 圧力
S1 実測ステップ
S2 解析モデル作成ステップ
S3 理論計算モデル作成ステップ
S4 剛性係数導出ステップ
tol 所定値
(x,y)、(xi,yj) 座標
α 金型長辺長さ
β 金型短辺長さ
ω (プレス機の)たわみ
ωω (プレス機モデルの)たわみ
ωij 実変位
ωωij 仮想変位
rm a 長辺方向スケーリングパラメータ
rn b 短辺方向スケーリングパラメータ
rm0 a 初期長辺方向スケーリングパラメータ
rn0 b 初期短辺方向スケーリングパラメータ
Claims (12)
- プレス機の剛性予測方法であって、
代表金型を介して前記プレス機に荷重を載荷して、前記荷重と前記プレス機における複数の測定点での実変位との関係を実測する実測ステップと、
前記プレス機の長辺長さ及び短辺長さを反映したプレス機モデルを用意する解析モデル作成ステップと、
前記プレス機モデルのたわみと、前記荷重との関係を、前記プレス機の剛性を示す剛性係数を用いた関数で表す理論計算モデル作成ステップと、
前記実変位及び前記荷重に基づいて、前記剛性係数を導出する剛性係数導出ステップと、を含む
ことを特徴とするプレス機の剛性予測方法。 - 前記代表金型の載荷面積と、指定のプレス金型の載荷面積との比に基づいて、前記剛性係数を補正し、前記指定のプレス金型を備えた前記プレス機の補正剛性係数を算出する補正ステップを含む
ことを特徴とする請求項1に記載のプレス機の剛性予測方法。 - 前記関数は、二次元フーリエ級数展開されたものである
ことを特徴とする請求項1又は請求項2に記載のプレス機の剛性予測方法。 - 前記関数は、正弦関数によって二次元フーリエ級数展開されたものである
ことを特徴とする請求項3に記載のプレス機の剛性予測方法。 - 前記関数は、ルジャンドル多項式によって二次元フーリエ級数展開されたものである
ことを特徴とする請求項3に記載のプレス機の剛性予測方法。 - 前記関数は、ラゲール多項式によって二次元フーリエ級数展開されたものである
ことを特徴とする請求項3に記載のプレス機の剛性予測方法。 - 前記関数は、エルミート多項式によって二次元フーリエ級数展開されたものである
ことを特徴とする請求項3に記載のプレス機の剛性予測方法。 - 前記二次元フーリエ級数展開の級数は、有限である
ことを特徴とする請求項3から請求項11のいずれか1項に記載のプレス機の剛性予測方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020216687A JP2022102135A (ja) | 2020-12-25 | 2020-12-25 | プレス機の剛性予測方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020216687A JP2022102135A (ja) | 2020-12-25 | 2020-12-25 | プレス機の剛性予測方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022102135A true JP2022102135A (ja) | 2022-07-07 |
Family
ID=82272984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020216687A Pending JP2022102135A (ja) | 2020-12-25 | 2020-12-25 | プレス機の剛性予測方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2022102135A (ja) |
-
2020
- 2020-12-25 JP JP2020216687A patent/JP2022102135A/ja active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7194388B2 (en) | Method for determining a die profile for forming a metal part having a desired shape and associated methods | |
US5156053A (en) | Measuring system using a robot | |
US20110246150A1 (en) | Method, device, program, and recording medium of analyzing cause of springback | |
WO2010073756A1 (ja) | スプリングバック発生原因分析方法、スプリングバック発生原因分析装置、スプリングバック発生原因分析プログラム及び記録媒体 | |
WO2010122769A1 (ja) | 成形シミュレーション方法、成形シミュレーション装置、及び成形シミュレーションプログラム、並びにその記録媒体 | |
JPWO2005093612A1 (ja) | 溶接変形算出方法、溶接変形算出装置、及びコンピュータプログラム | |
JP4884839B2 (ja) | プレス金型修正形状データの作成方法 | |
US6539763B1 (en) | Precision press brake | |
JP4884840B2 (ja) | プレス金型修正形状データの作成方法 | |
JP3978377B2 (ja) | 成形シミュレーション解析方法 | |
JPH11179433A (ja) | 曲げ加工方法およびこの曲げ加工システム | |
CN113950650A (zh) | 调整量推定装置、调整量推定方法及调整量推定程序 | |
JP5969886B2 (ja) | 対応点算出システム及びプログラム、並びに金型形状生成システム及びプログラム | |
JP2022102135A (ja) | プレス機の剛性予測方法 | |
JP2007229784A (ja) | プレス金型修正形状データの作成方法 | |
JP2016198997A (ja) | コンピュータ支援の樹脂挙動解析装置 | |
JP4313890B2 (ja) | スプリングバック量予測方法 | |
JP6044606B2 (ja) | 見込み金型形状作成方法及び装置 | |
JP3809527B2 (ja) | 固有変形算出装置、コンピュータプログラム、記録媒体及び溶接時の固有変形の算出方法 | |
JP2005138120A (ja) | プレス成形中の金型たわみ分布のシミュレーション方法 | |
JP2002126834A (ja) | 金型およびその設計装置、設計方法ないし製造装置、並びに、成形品 | |
JPH10154169A (ja) | 成形品の変形評価方法 | |
Maier et al. | Data-driven compensation for bulk formed parts based on material point tracking | |
JP7541241B2 (ja) | プレス成形シミュレーション方法 | |
JP2003344041A (ja) | ワーク形状評価装置およびその方法ならびにそのプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230821 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240612 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240702 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240809 |