プレス成形で用いる金型に存在するドロービードを解析上省略し、その代替としてドロービード位置で被加工材にドロービード張力を付与した等価ドロービードが設定された、コンピュータを用いたシェル要素による成形解析方法であって、実験又は計算に基づき予め求められた、前記金型のダイ肩R部における前記被加工材の通過抵抗力を記憶媒体に記憶する通過抵抗力記憶工程と、成形条件を設定する成形条件設定工程と、前記成形条件に基づいて成形解析演算する演算工程と、を備え、前記成形条件設定工程において、前記記憶媒体に記憶された前記通過抵抗力を前記成形条件に反映させることを特徴とする成形解析方法。前記通過抵抗力は、実験又は計算により求められた実測通過抵抗力から、前記通過抵抗力を前記成形条件に反映しないで成形解析した場合に生じる前記ダイ肩R部における仮想通過抵抗力を減算した値である請求項1に記載の成形解析方法。前記成形条件設定工程では、前記ドロービード張力に前記通過抵抗力を加える請求項2に記載の成形解析方法。前記成形条件設定工程では、前記成形条件における前記ダイ肩R部位置に対して前記通過抵抗力を加える請求項2に記載の成形解析方法。前記実測通過抵抗力は、前記ダイ肩R部に生じる摩擦力と、前記被加工材が前記ダイ肩R部を通過する際の曲げ-曲げ戻し変形で生じる応力との和である請求項2~4の何れか一項に記載の成形解析方法。プレス成形で用いる金型に存在するドロービードを解析上省略し、その代替としてドロービード位置で被加工材にドロービード張力を付与した等価ドロービードが設定され、シェル要素により成形解析を行う、コンピュータを備えた成形解析装置であって、実験又は計算に基づき予め求められた、前記金型のダイ肩R部における前記被加工材の通過抵抗力を記憶する通過抵抗力記憶部と、前記通過抵抗力を反映した成形条件を記憶する成形条件記憶部と、前記成形条件に基づいて成形解析演算する演算部と、を備えることを特徴とする成形解析装置。請求項1~5の何れか一項に記載の方法を実施するためにコンピュータにより実行されるプログラムであって、前記通過抵抗力記憶工程では、ユーザーの操作に応じて、前記通過抵抗力を前記記憶媒体に記憶し、前記成形条件設定工程では、ユーザーの操作に応じて、前記成形条件を設定することを特徴とするプログラム。請求項7に記載のプログラムをコンピュータ読み取り可能に記録した記録媒体。しかしながら、上記のような従来の成形解析では、金型のダイ肩R部におけるR形状が正確に再現されておらず、その分、解析精度が低下していた。つまり、計算時間を短縮するためにメッシュを大きくとるため、ダイ肩R部の形状は、大まか例えば、直角又は緩やかな曲線の両極端に表される。換言すると、ダイ肩R部の形状は、メッシュが粗いため、そのままの曲率では表現できず、例えばR部を拡大して表している。例えば、自動車の大物プレス部品を、一般に普及している性能のコンピュータを用いて数時間程度の実用的な計算時間で解析しようとする場合、3mm以上のメッシュを用いなければならないため、R1曲率半径が1mmやR2曲率半径が2mmのダイ肩R部形状は精度良く再現されない。これによれば、ダイ肩R部で生じる抵抗力は、解析上と実際上とで大きく異なる可能性があり、解析精度向上の面で問題があった。本発明は、このような事情に鑑みて為されたものであり、ダイ肩R部に生じる抵抗力に着目し、解析時間の増加を抑えつつ、解析精度を向上させることができる成形解析方法及び成形解析装置を提供することを目的とする。本発明は、プレス成形で用いる金型に存在するドロービードを解析上省略し、その代替としてドロービード位置で被加工材にドロービード張力を付与した等価ドロービードが設定された、コンピュータを用いたシェル要素による成形解析方法であって、実験又は計算に基づき予め求められた、前記金型のダイ肩R部における前記被加工材の通過抵抗力を記憶媒体に記憶する通過抵抗力記憶工程と、成形条件を設定する成形条件設定工程と、前記成形条件に基づいて成形解析演算する演算工程と、を備え、前記成形条件設定工程において、前記記憶媒体に記憶された前記通過抵抗力を前記成形条件に反映させることを特徴とする。本方法によれば、解析上、ダイ肩R部の形状が大まかに表されていたとしても、実験又は計算に基づいて予め求められた通過抵抗力が解析結果に反映される。メッシュを細かくしていないため計算時間の増加はほぼない。しかも、通過抵抗力を反映するため被加工材の流入量計算を実際の流入量に近づけることができ、結果として成形解析精度は向上する。ここで、前記通過抵抗力は、実験又は計算により求められた実測通過抵抗力から、前記通過抵抗力を考慮しないで成形解析した場合に生じる前記ダイ肩R部における仮想通過抵抗力を減算した値であることが好ましい。これによれば、成形条件をより正確に修正することができる。ここで、前記成形条件設定工程では、前記ドロービード張力に前記通過抵抗力を加えるようにしてもよい。これによれば、成形条件設定工程において、元々設定されていたドロービード張力を修正するだけでよく、ソフトウェア設計が容易となる。一方、前記成形条件設定工程では、前記成形条件における前記ダイ肩R部位置に対して前記通過抵抗力を加えるようにしてもよい。換言すると、本方法では、通過抵抗力を、ドロービード張力同様、成形解析における演算処理の条件としてダイ肩R部位置に加える。これによれば、通過抵抗力を加える位置が実際に抵抗力が生じる位置と同等になり、より正確な解析が可能となる。ここで、前記実測通過抵抗力は、前記ダイ肩R部に生じる摩擦力と、前記被加工材が前記ダイ肩R部を通過する際の曲げ-曲げ戻し変形で生じる応力との和であることが好ましい。これにより、実状に即した成形解析が可能となる。また、本方法は、コンピュータを備える成形解析装置としても記載することができる。すなわち、本発明は、プレス成形で用いる金型に存在するドロービードを解析上省略し、その代替としてドロービード位置で被加工材にドロービード張力を付与した等価ドロービードが設定され、シェル要素により成形解析を行う、コンピュータを備えた成形解析装置であって、実験又は計算に基づき予め求められた、前記金型のダイ肩R部における前記被加工材の通過抵抗力を記憶する通過抵抗力記憶部と、前記通過抵抗力を反映した成形条件を記憶する成形条件記憶部と、前記成形条件に基づいて成形解析演算する演算部と、を備えることを特徴とする。この構成によれば、本方法と同様の効果が発揮される。上記いずれの方法においても同様に装置として記載できる。また、本発明は、上記方法を実施するためにコンピュータにより実行されるプログラムとしても記載できる。この場合、前記通過抵抗力記憶工程では、ユーザーの操作に応じて、前記通過抵抗力を前記記憶媒体に記憶し、前記成形条件設定工程では、ユーザーの操作に応じて、成形条件を設定する。また、本発明は、このプログラムをコンピュータ読み取り可能に記録した記録媒体としても記載できる。本発明によれば、解析時間の増加を抑えつつ、解析精度を向上させることができる。次に、好ましい実施形態を挙げ、図1~図6を参照して本発明をより詳しく説明する。本実施形態の成形解析方法は、CPUや記憶媒体メモリ、ハードディスク等を備えたコンピュータ成形解析装置を用いて行うものであり、主に、通過抵抗力記憶工程と、成形条件設定工程と、演算工程と、を備えている。まず、通過抵抗力記憶工程では、実験又は計算に基づき予め求められた、プレス成形における金型ダイのダイ肩R部9における被加工材Wの通過抵抗力を、コンピュータのメモリ又はハードディスク等の記憶媒体に記憶する。通過抵抗力は、図1に示すように、プレス時に流入しようとする被加工材Wに対し、流入を阻止しようとする力図1の矢印参照である。ここでの通過抵抗力とは、実験又は計算により求められた実測通過抵抗力から、通過抵抗力を考慮しないで成形解析した場合に生じるダイ肩R部における仮想通過抵抗力を減算した値である。実測通過抵抗力は、被加工材の流入に対するダイ肩R部の摩擦力と、被加工材がダイ肩R部9を通過する際の曲げ-曲げ戻し変形で生じる応力との和である。仮想通過抵抗力は、従来の成形解析のように、通過抵抗力を反映させないで成形解析した場合のダイ肩R部9の抵抗力である。仮想通過抵抗力は、粗いメッシュで且つ拡大されたダイ肩R部9に生じる抵抗力、すなわち正確でない抵抗力の値である。メッシュを細かくすることで再現度が高まるが、計算時間が大きくなってしまう。仮想通過抵抗力は設定メッシュ粗さ等に応じて計算過程で生じる値であるため、当該値を実測通過抵抗力から減算する差分をとることでより正確な値通過抵抗力となる。なお、金型及び被加工材Wの材料に応じた通過抵抗力のデータベースをコンピュータに記憶してもよい。実験計算については後述する。成形条件設定工程は、ユーザー操作により、成形解析演算の基準となる成形条件を設定するいわゆる前工程の1つである。成形条件設定工程では、ユーザーにより、例えば、摩擦係数、しわ押え力、工具移動速度、移動量、移動方向、各種変位拘束、及び数値ビード力ドロービード張力などの境界条件がコンピュータに設定される。また、成形条件設定工程では、例えば、被加工材の縦弾性係数、密度、降伏応力、引張強さ、及びr値塑性ひずみ比、ランクフォード値などの材料特性値も設定される。なお、成形解析における他の前工程は、公知の方法と同様である。つまり、前工程は、成形条件設定工程のほかに、例えば、成形工具形状モデリング工程、工具メッシュデータ作成工程、被加工材メッシュデータ作成工程、及びブランクホールド工程などを含んでいる。成形工具形状モデリング工程は、製品形状に対し余肉、しわ押え面、ビードなどの成形に必要な形状を設計する工程である。工具メッシュデータ作成工程は、工具サーフェスデータから金型のメッシュデータを作成する工程である。被加工材メッシュデータ作成工程は、材料サイズ定義や材料メッシュデータを作成する工程である。ブランクホールド工程は、ブランクホルダで板材を挟んだ際の状態を解析する工程である。ブランクホールド工程の解析結果は、成形条件設定工程に反映される。このように、実際の演算処理後処理の前には、前処理工程として各種条件等の設定が行われる。演算工程は、コンピュータが、成形条件に基づいて、公知の静的陰解法等により成形解析演算を行う工程である。これにより、プレス成形後の被加工材の形状や応力分布などがシミュレーションされる。本実施形態では、演算工程において、通過抵抗力を考慮した成形解析演算が為される。具体的に、コンピュータは、演算工程の前に、メモリ等に記憶された通過抵抗力を成形条件に加える処理を行う適用工程。その後、コンピュータは、通過抵抗力が成形条件に反映された状態で、成形解析演算を実行する演算工程。なお、適用工程は、メモリ中からユーザーが通過抵抗力を選択して実行するものでもよい。また、順番として、前工程が為された後に、通過抵抗力をメモリ等に記憶し、通過抵抗力を成形条件に反映させてもよい。ここで、通過抵抗力を求めるための実験予備実験の一例について説明する。まず、被加工材と同じ材料の板材を準備する。図2に示すように、板材W1の中央部分を、成形で用いる金型のダイ肩R部と同形状の金型91を支点に曲げる。そして、板材W1の一方端部Yに一定荷重を与えた状態で、板材W1の他方端部Zを一定距離ストロークさせる図2矢印参照。その際、当該他方端部Zにかかる荷重を測定する。これにより、金型91ダイ肩R部9に生じる摩擦力と、板材W1が金型91ダイ肩R部9を通過する際の曲げ-曲げ戻し変形で生じる応力との和である実測通過抵抗力実測した通過抵抗力が得られる。予備実験において一方端部Yに与える一定荷重についてはいくつかの条件荷重値で行われ、コンピュータ演算時には、一方端部Yに相当する部分の荷重を参照して通過抵抗力が選択される。なお、コンピュータは、一方端部Yに相当する部分の荷重が実験条件同士間の荷重である場合は、内挿により通過抵抗力を決定する。また、通過抵抗力を求めるための計算の一例は、ダイ肩R部9についての実験図2の実験を再現した成形解析演算シミュレーションによって求める方法である。例えば、ソリッド要素を用いて、ダイ肩R部9についてのみシミュレーションを行う。ユーザーは摩擦係数や板材の材料モデル応力とひずみの関係を特定し、コンピュータがソリッド要素で計算する。これによっても実測通過抵抗力を求めることができる。ここで、処理フローについて簡単に説明する。図3に示すように、まず、予備実験又は計算により実測通過抵抗力を求めるS101。続いて、通過抵抗力を加えない状態で成形解析した場合に生じる解析上の通過抵抗力仮想通過抵抗力を求め、実測通過抵抗力から仮想通過抵抗力を引き減算し、通過抵抗力を求めるS102。続いて、求めた通過抵抗力をコンピュータの記憶部例えばメモリに記憶させるS103。続いて、前工程を行い、各種成形条件を設定するS104。続いて、設定した成形条件に対し、記憶部に記憶された通過抵抗力を反映させるS105。具体的に、ここでは、ドロービード張力に通過抵抗力を加算する。続いて、通過抵抗力が反映された成形条件に基づき、コンピュータが動的陽解法を用いて成形解析演算を行うS106。演算終了後、ディスプレイ等の表示装置に解析結果を表示するS107。通過抵抗力は、同じ材料で成形解析する場合、新たに予備実験等をすることなく適用することができる。実施例ここで、本方法にて実際に成形解析した結果本発明と、通過抵抗力を考慮しない従来方法で成形解析した結果従来とを比較する。図4に示すようなプレス品に対して、フランジ部分の長さを比較することで、材料の流入量を比較する。本方法にあっては、B側伸び側のダイ肩R部がR1曲率半径が1mmでメッシュが1.5mmであるため、R部形状をR3曲率半径が3mmに拡大している。そして、R1とR3の通過抵抗力の差をB側のドロービード張力に加算して成形解析した。図5に示すように、本発明によれば、A側縮み側の中央及び端において、従来よりも実際にプレス成形した実験値に近似した値が算出された。また、B側の中央及び端においても、従来よりも実験値に近似した値が算出された。今回の解析において、本発明はB側に通過抵抗力を適用しており、B側の精度が従来よりも遥かに向上している。B側の精度が向上することで、A側の精度も向上し、全体として解析精度が向上した。このように、本実施形態によれば、メッシュ粗さを細かくせずとも解析精度を向上させることができる。つまり、本実施形態によれば、解析時間の増加を抑えつつ、解析精度を向上させることができる。なお、上記実施例では、ドロービード張力に通過抵抗力を加えていたが、ダイ肩R部9位置に加えるように設計してもよい。この方法は、例えば、等価ドロービード設定においてビード位置にドロービード張力を加えるように、ダイ肩R部9位置に通過抵抗力を加える設定を行うことで実現することができる。これによれば、実際に抵抗力が生じるダイ肩R部9位置に対して通過抵抗力が設定されるため、より高精度な成形解析が可能となる。この方法は、例えば、ビード部8とダイ肩R部9の位置が大きく離れている場合に特に有効である。また、本実施形態は、図6に示すように、コンピュータを備えた成形解析装置1を用いて実行している。成形解析装置1は、主に、データ等を記憶する記憶部2と、演算処理を実行する演算部3と、を備えている。記憶部2は、メモリ又はハードディスク等の記憶媒体であって、通過抵抗力、成形条件、及びプログラム等を記憶している。ユーザーの選択又は演算部3のプログラム処理により、通過抵抗力は成形条件に組み込まれる。演算部3は、CPU等であり、解析プログラムにそって演算処理を実行する。具体的に、演算部3は、記憶部2に記憶され且つ通過抵抗力が反映された成形条件に基づき、解析プログラムに従って成形解析演算をする。これにより、本方法が実行され、上記効果が発揮される。本発明の方法は、それを実施するために各工程が規定されたプログラムにより、コンピュータにより実行できる。このプログラムは、CD-ROMやDVD-ROM等の記憶媒体に記憶されていてもよい。なお、本実施形態では、成形解析装置1は、ディスプレイ等の表示装置4と、キーボード等の入力装置5と、をさらに備えている。また、通過抵抗力やプログラム等は、それぞれCD-ROMやDVD-ROM等の記憶媒体記憶部2に相当するに記憶されていてもよい。通過抵抗力と成形条件とは別々の記憶媒体に保存されていてもよい。