JP2007188454A - 設計支援装置及びコンピュータプログラム - Google Patents
設計支援装置及びコンピュータプログラム Download PDFInfo
- Publication number
- JP2007188454A JP2007188454A JP2006008057A JP2006008057A JP2007188454A JP 2007188454 A JP2007188454 A JP 2007188454A JP 2006008057 A JP2006008057 A JP 2006008057A JP 2006008057 A JP2006008057 A JP 2006008057A JP 2007188454 A JP2007188454 A JP 2007188454A
- Authority
- JP
- Japan
- Prior art keywords
- equation
- procedure
- error
- linear combination
- solution
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Abstract
【解決手段】設計支援のためのコンピュータプログラムにおいて、未知関数u、初期条件、境界条件を入力する第1手順と、uを正規直交基底の線形結合で近似する第2手順と、線形結合に含まれる線形結合定数を根とする代数方程式を生成する第3手順と、初期条件、境界条件を正規直交基底で表わし、1以上の拘束条件を得る第4手順と、代数方程式の根を拘束条件の下で求める第5手順と、根が代数方程式の2乗誤差を最小ならしめるか否かを判定する第6手順とを含み、2乗誤差が最小であると判定されれば、根を方程式の解とし、2乗誤差が最小でないと判定されれば、第5手順に戻る。
【選択図】図6
Description
誤差が最小であると判定されれば、根を方程式の解とし、誤差が最小でないと判定されれば、第5手順に戻る。
(t−b1)(t−b2)........(t−bn)
の根がすべて実でなければならない。そこで、nが小さいときは根の判別式を利用することも出来る。ここで、m個の線形結合係数が、拘束条件(初期条件、境界条件)によって関係付けられているため、代数方程式については、(nーm)の実根を求めればよいことになる。代数方程式が解けるということは、元の関数方程式が解けるということと同じであるが、代数方程式は、未知方程式f(x)のONB(x)展開近似に基づくため、解の正確さまでが保障されるものではない。
u(x、y)= X(x)・Y(y)として、
X(x)=(p1ψ1(x)+p2ψ2(x)+p3ψ3(x)) (1)
Y(y)=(k1φ(y)+k2φ(y)) (2)
u(x、y)は、(ψiφj)で表現できる。ここで、τk(x、y)=(ψi(x)・φj(y))(k=1to 6)も正規直交系をなすことは簡単に証明できる。なぜなら、(ψiφj)を積分すると、ψiψjφkφm(i,j=1 to3;k、m=1 to2)の積分となるので、ψi、φiの直交関数性を使えば、τkも正規直交性を有するからである。したがって、u(x、y)は、正規直交基底OCN(τk(x、y))を使って、
u(x、y)=a1τ1(x、y)
+a2τ2(x、y)
+a3τ3(x、y)
+a4τ4(x、y)
+a5τ5(x、y)
+a6τ6(x、y) (3)
ここで 、
τ1=ψ1(x)・φ1(y)
τ2=ψ1(x)・φ2(y)
τ3=ψ2(x)・φ1(y)
τ4=ψ2(x)・φ2(y)
τ5=ψ3(x)・φ1(y)
τ6=ψ3(x)・φ1(y) (4)
dX(x)/dx=q1ψ1(x)+q2ψ2(x)+q3ψ3(X)
dY(y)/dy=r1φ1(y)+r2φ2(y) (5)
同様に、
d2X(x)/dx2=l1ψ1(x)+l2ψ2(x)+l3ψ3(X)
d2Y(y)/dy2=m1φ1(y)+m2φ2(y) (6)
PE=b1τ1(x、y)
+b2τ2(x、y)
+b3τ3(x、y)
+b4τ4(x、y)
+b5τ5(x、y)
+b6τ6(x、y) (7)
以上のように、ポアソン方程式を正規直交基底ONBの線形結合で代数化できた。
ここに、X(x)を正規直交基底の線形結合で近似する。
X(x)≒k1ψ1(x)+k2ψ2(x) (8)
この近似を使って、原点での境界条件を表現すると、
k1ψ1(0)+k2ψ2(0)=4 (9)
となるが、これは(X(0)-4)とディラックのデルタ関数δ(x)の内積に等しいことを下式に示す。
下式は、線素AB上の関数空間の基底ベクトルδABを表す式である。
近似関数u(x、γ)と境界条件式g1(x)との関係は、
u(x、γ)=a1τ1(x、γ)
+a2τ2(x、γ)
+a3τ3(x、γ)
+a4τ4(x、γ)
+a5τ5(x、γ)
+a6τ6(x、γ)=g1(x) (10)
を上述した基底ベクトルδABを使用して積分したものである。
同様に、境界線素BC,CD,AD上からも定数1次式が得られる。
このように、ベクトルAの成分が満たす4つの拘束条件が得られた。したがって、拘束条件が表す領域をビジュアルに表示することにより、領域の谷底や山頂などに、安定解、不安定解などの存在を予測することが出来る。このように、本発明のポアソン方程式の解法は、極めて見とおしのよいものである。
(t−a1)(t−a2)(t−a3)(t−a4)(t−a5)(t−a6)=0 (11)
2実根を求めればよい理由は、6つの線形結合定数のうち4つは境界条件により既に関連付けられているためである。
しかし、方程式の解を確定するためには、ポアソン方程式の誤差又は2乗誤差などを最小化しなければならない。
y(x)=b1ψ1(x)+b2ψ2(x)+b3ψ3(x)+b4ψ4(x)(12)
と近似する。
ψとして、フーリエ級数、ラゲール多項式、ルジャンドル多項式などが境界条件に対応して使用される。
I=(2/3)(b12+b22+b32+b42) (13)
y“は、正規直交系で表現して、y“=d1ψ1+d2ψ2+d3ψ3+d4ψ4としてもよいが、y”を0から1まで積分する都合上、
y‘=c1ψ1(x)+c2ψ2(x)+c3ψ3(x)+c4ψ4(x)
を使用する。
y”の0から1までの定積分は、y’(1)とy’(0)の差分である。
ところで、定数係数c1〜c3は、b1〜b3に基づいて生成されたものであり、b1、b2、b3を成分とする列ベクトルに、3行3列の変換行列(R)を乗算して生成される。したがって、y”の0から1までの定積分I‘は、
I‘=c1(ψ1(1)−ψ1(0))
+c2(ψ2(1)-ψ2(0))
+c3(ψ3(1)-ψ3(0))
+c4(ψ4(1)-ψ4(0)) (14)
微分方程式は、式(13)と式(14)を等置したものである。
b1ψ1(1)+b2ψ2(1)+b3ψ3(1)+b4ψ4(1)=1 (15)
b1ψ1(0)+b2ψ2(0)+b3ψ3(0)+b4ψ4(0)=4 (16)
である。これを0から1まで積分して、下式の境界条件を得る。
(t−b1)(t−b2)(t−b3)(t−b4)=0
境界条件(15)、(16)の下で未知の根は2つに減少する。
b12+b22+b32+b42=定数
となる(パーセヴアル等式近似)。この意味で、線形結合係数ベクトルBの先端は、高次元球面上にある。解Bの先端は、目的関数Lを最小にしており、高次元曲面上で原点を中心とする最も半径の小さな超球と接すると見ても良い。しかし、クランク軸を分析したMCP手法(特開2005−9323号公報)の場合のように、半径SQRT(n)(クランク軸数は(n+1))の球上にはない。MCP手法のときは、球面を固定して拘束面を動かすが、本発明では、拘束面を固定して、球の半径を動かしている。
方程式d2y/dt2+(1-2qcos(2t))y=0
を取り上げる。図6は、実施例3の解法のフローチャートである。
境界条件は、出発時刻t=0として、
y(0)=y(π)
y’(0)=y’(π)
である。
未知関数yの正規直交基底展開は、
y=b1ψ1(t)+b2ψ2(t)+b3ψ3(t)+b4ψ4(t) (17)
y’=c1ψ1(t)+c2ψ2(t)+c3ψ3(t)+c4ψ4(t) (18)
y”=d1ψ1(t)+d2ψ2(t)+d3ψ3(t)+d4ψ4(t) (19)
ここに、b1〜b4を成分とする列ベクトルを(B)、変換行列を(K1)とすると、
c1〜c4を成分とする列ベクトル(C)=(K1)(B)であり、
d1〜d4を成分とする列ベクトル(D)=(K2)(B)である。
式(15)及び式(17)を用いると、Mathieu方程式が代数化される。
(ψT(0))(B)=(ψT(π))(B) (20)
(ψT(0))(K1)(B)=(ψT(π))(K1)(B) (21)
これを0からπまで積分したものが境界条件である。
(t-b1)(t-b2)(t−b3)(t−b4)
となる。この4次代数方程式が、拘束条件(20)、(21)の下で解かれる。したがって、2つの実根を求めればよい。
目的関数Lは、ベクトル形式で、
L=∫(y”+(1-2qcos2t)y)2dt
=(B)T(K1)T・(∫((ψ(t))・(ψ(t))Tdt)・(K2)(B)
+2(B)T(K2)T・(∫((ψ(t))・(1-2qcos2t)・(ψ(t))Tdt)・(B)+(B)T・(∫((ψ(t))・(1−2qcos2t)・(ψ(t))Tdt)・(B)
となる。積分は、0からπまでである。目的関数Lは、パラメトリック発振であろうが、時間tの関数であろうが、(B)の2次曲面形式になる。
y(t+π)=exp(iπν)y(t)
で与えられることがあり(Floquetの解)、この場合の境界条件のベクトル表現は、
(ψ(π))T・(B)=exp(iπν)・(ψ(0))T・(B)
(ψ(π))T・(K1)・(B)=exp(iπν)・(ψ(0))T・(B)
となる。b1〜b4の間に成立する2つの1次形式の係数は一般に複素数となるが、方程式(18)の根がすべて実根であることに変わりはない。
y1=b1φ1(t)+b2φ2(t) (22)
y2=c1φ1(t)+c2φ2(t) (23)
y1(T)=y1(0)
y2(T)=y2(0)
dy1/dt(0)=dy2/dt(0)
dy1/dt(T)=dy2/dt(T) (24)
式(21)、式(22)を、4つの境界条件(24)に代入すれば、境界条件が代数化され、それを0からTまで積分して境界条件を得る。
(dy1/dt(T)=dy1/dt(0)、dy2dt(T)=dy2/dt(0))
を課せば、基底数をそれだけ増やしやすくなる。
y=b1ψ1(x)+b2ψ2(x)+b3ψ3(x) (25)
b1ψ1(0)+b2ψ2(0)+b3ψ3(0)=C (26)
これが境界条件である。
代数方程式は、
(t-b1)(t−b2)(t−b3)
=t3−(b1+b2+b3)t2+(b1b2+b2b3+b3b1)t+b1b2b3
(27)
この代数方程式に、拘束条件式(26)から、例えば、
b3=(C−b1ψ1(0)−b2ψ2(0))/ψ3
を代入すれば、求める根はb1、b2だけとなる。したがって、2次元実空間である許容領域の中で、2つの実根を探せばよい。
S56では、目的関数Lは、
2乗誤差(b1ψ1(x)+b2ψ2(x)+b3ψ3(x)−f(x))2
を1〜αまで積分したものである。αは有限であるが、必要に応じて充分大なる値とする。
T(t)=a1φ1(t)+a2φ2(t)
X(x)=c1ζ1(x)+c2ζ2(x)
すると、u=a1c1φ1(t)ζ1(x)
+a1c2φ1(t)ζ2(x)
+a2c1φ2(t)ζ1(x)
+a2c2φ2(t)ζ2(x)
=b1τ1(t、x)
+b2τ1(t、x)
+b3τ1(t、x)
+b4τ4(t、x) (28)
これは、ポアソン方程式の場合と同様の手順である。
2 近似次数入力部
3 許容領域表示部
4 解出力部
5 CPU
Claims (11)
- 未知関数uの方程式の解を求める設計支援装置において、
前記方程式、及び、その初期条件及び/又は境界条件、若しくは前記方程式に伴う拘束条件を入力する方程式入力部と、
前記未知関数を、有限個数nの規格直交基底の線形結合で近似するときの前記nを入力する近似次数入力部と、
前記近似により生じた前記方程式の誤差であって、maxノルム、1−ノルム、2−ノルム、ミンコフスキーノルム、ヘルダーノルム、又は距離の公理を満足するノルムで表わされる誤差を計算するCPUと、
前記初期条件及び/又は境界条件によって関係付けられた前記線形結合の結合係数の数をm個とするとき、(前記n−前記m)次元実空間の領域を許容領域として表示する許容領域表示部と、
前記方程式の解を出力する解出力部とを備え、
前記許容領域内の任意の点を候補点として、前記候補点における前記誤差を計算し、
さらに、別の前記候補点における前記誤差を計算し、
前記誤差を最小ならしめる前記許容領域の点を前記方程式の解として求めることを特徴とする設計支援装置。 - 請求項1において、
前記方程式は、微分作用素、積分作用要素、漸化式の写像作用素、又は非線形方程式における非線形作用素をD、定数をC、未知関数を前記uとするとき、Du=Cであり、前記uを前記線形結合で近似し且つ前記方程式全体を代数化した後のDu=Cを使って、前記誤差である前記ノルム||Du−C||を計算することを特徴とする設計支援装置。 - 請求項1において、
前記許容領域表示部は、
前記許容領域が2次元以下のときは2次元以下のグラフィック表示を行い、
前記許容領域が3次元以上のときは3次元のグラフィック表示を行い、又は、2次元領域に階層化された複数のグラフとして表示し、
前記許容領域の曲面が3次元以上のときは、3次元のグラフィック表示を行い、又は、2次元領域に階層化された立体グラフ又は等高線図として表示することを特徴とする設計支援装置。 - 請求項1において、
前記CPUは、前記線形結合の係数bi(iは1から前記nまで)を根とする代数方程式において、
前記根のすべてを実数に制限するか、
又は、前記根が1対をなす共役根であるときは複素平面内の一部の領域に存在するものに限定することを特徴とする設計支援装置。 - コンピュータを、
前記方程式、及び、その初期条件及び/又は境界条件を入力する方程式入力部、
前記未知関数を、有限個数nの規格直交基底の線形結合で近似するときの前記nを入力する近似次数入力部、
前記近似により生じた前記方程式の誤差であって、maxノルム、1−ノルム、2−ノルム、ミンコフスキーノルム、ヘルダーノルム、又は距離の公理を満足するノルムで表わされる誤差を計算するCPU、
前記初期条件及び/又は境界条件によって関係付けられた前記線形結合の結合係数の数をm個とするとき、(前記n−前記m)次元実空間の領域を許容領域として表示する許容領域表示部、
前記方程式の解を出力する解出力部として機能させるコンピュータプログラムであって、
前記許容領域内の任意の点を候補点として、前記候補点における前記誤差を計算し、
さらに、別の前記候補点における前記誤差を計算し、
前記誤差を最小ならしめる前記許容領域の点を前記方程式の解として求めることを特徴とするコンピュータプログラム。 - 請求項5において、
前記方程式は、微分作用素、積分作用要素、漸化式の写像作用素、又は非線形方程式における非線形作用素をD、定数をC、未知関数を前記uとするとき、Du=Cであり、前記uを前記線形結合で近似し且つ前記方程式全体を代数化した後のDu=Cを使って、前記誤差である前記ノルム||Du−C||を計算することを特徴とするコンピュータプログラム。 - 請求項5において、
前記許容領域表示部は、
前記許容領域が2次元以下のときは2次元以下のグラフィック表示を行い、
前記許容領域が3次元以上のときは3次元のグラフィック表示を行い、又は、2次元領域に階層化された複数のグラフとして表示し、
前記許容領域の曲面が3次元以上のときは、3次元のグラフィック表示を行い、又は、2次元領域に階層化された立体グラフ又は等高線図として表示することを特徴とするコンピュータプログラム。 - 請求項5において、
前記CPUは、前記線形結合の係数bi(iは1から前記nまで)を根とする代数方程式において、
前記根のすべてを実数に制限するか、
又は、前記根が1対をなす共役根であるときは複素平面内の一部の領域に存在するものに限定することを特徴とするコンピュータプログラム。 - 請求項5において、
未知関数u、初期条件、境界条件を入力する第1手順と、
前記uを変数分離し、各変数の関数を、各々異なる正規直交基底の線形結合で近似する第2手順と、
前記方程式を前記正規直交基底で表わし、その表現に含まれる線形結合定数を根とする代数方程式を生成する第3手順と、
前記初期条件、前記境界条件を前記正規直交基底で表わし、1以上の拘束条件を得る第4手順と、
前記代数方程式の根を前記拘束条件の下で求める第5手順と、
前記根が前記代数方程式の前記誤差を最小ならしめるか否かを判定する第6手順とを含み、
前記誤差が最小であると判定されれば、前記根を前記方程式の解とし、
前記誤差が最小でないと判定されれば、第5手順に戻ることを特徴とするコンピュータプログラム。 - 請求項5において、
未知関数u、初期条件、境界条件を入力する第1手順と、
前記uを正規直交基底の線形結合で近似する第2手順と、
前記方程式を定積分し、前記定積分に含まれる線形結合定数を根とする代数方程式を生成する第3手順と、
前記初期条件、前記境界条件を前記正規直交基底で表わし、1以上の拘束条件を得る第4手順と、
前記代数方程式の根を前記拘束条件の下で求める第5手順と、
前記根が前記代数方程式の前記誤差を最小ならしめるか否かを判定する第6手順とを含み、
前記誤差が最小であると判定されれば、前記根を前記方程式の解とし、
前記誤差が最小でないと判定されれば、第5手順に戻ることを特徴とするコンピュータプログラム。 - 請求項5において、
未知関数u、初期条件、境界条件を入力する第1手順と、
前記uを正規直交基底の線形結合で近似する第2手順と、
前記線形結合に含まれる線形結合定数を根とする代数方程式を生成する第3手順と、
前記初期条件、前記境界条件を前記正規直交基底で表わし、1以上の拘束条件を得る第4手順と、
前記代数方程式の根を前記拘束条件の下で求める第5手順と、
前記根が前記代数方程式の前記誤差を最小ならしめるか否かを判定する第6手順とを含み、
前記誤差が最小であると判定されれば、前記根を前記方程式の解とし、
前記誤差が最小でないと判定されれば、第5手順に戻ることを特徴とするコンピュータプログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006008057A JP4599301B2 (ja) | 2006-01-16 | 2006-01-16 | 設計支援装置及びコンピュータプログラム |
EP06027001.4A EP1808775A3 (en) | 2006-01-16 | 2006-12-28 | Design aiding apparatus and computer program |
US11/646,285 US20070179762A1 (en) | 2006-01-16 | 2006-12-28 | Design aiding apparatus and computer program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006008057A JP4599301B2 (ja) | 2006-01-16 | 2006-01-16 | 設計支援装置及びコンピュータプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007188454A true JP2007188454A (ja) | 2007-07-26 |
JP4599301B2 JP4599301B2 (ja) | 2010-12-15 |
Family
ID=37909778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006008057A Expired - Fee Related JP4599301B2 (ja) | 2006-01-16 | 2006-01-16 | 設計支援装置及びコンピュータプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US20070179762A1 (ja) |
EP (1) | EP1808775A3 (ja) |
JP (1) | JP4599301B2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009020619A (ja) * | 2007-07-10 | 2009-01-29 | Nippon Telegr & Teleph Corp <Ntt> | 多項式間距離算出装置、方法、プログラムおよび記録媒体、一変数最近実多項式算出装置、方法、プログラムおよび記録媒体 |
JP2010102680A (ja) * | 2008-09-25 | 2010-05-06 | Fujitsu Ltd | レイアウト評価装置、レイアウト評価プログラム、ダミールール生成装置及びダミールール生成プログラム |
JP2011233062A (ja) * | 2010-04-30 | 2011-11-17 | Shuzo Maeda | 3次元形状を3次元形状へ変換する方法 |
JP2012083958A (ja) * | 2010-10-12 | 2012-04-26 | Sumitomo Rubber Ind Ltd | 構造格子を用いたシミュレーション方法 |
JP2022506325A (ja) * | 2018-11-13 | 2022-01-17 | スーチョウ レインメド メディカル テクノロジー カンパニー リミテッド | 造影画像に基づき血管評定パラメータを取得する方法、装置及びシステム |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8238675B2 (en) * | 2008-03-24 | 2012-08-07 | Microsoft Corporation | Spectral information recovery for compressed image restoration with nonlinear partial differential equation regularization |
CN103615739B (zh) * | 2013-11-27 | 2015-11-18 | 广东电网公司电力科学研究院 | 燃烧锅炉运行控制方法与系统 |
CN103646143B (zh) * | 2013-12-17 | 2016-04-06 | 辽宁石油化工大学 | 一种火灾环境下lpg储罐涂层防护系统性能分析的方法 |
US9256709B2 (en) * | 2014-02-13 | 2016-02-09 | Taiwan Semiconductor Manufacturing Company, Ltd. | Method for integrated circuit mask patterning |
CN108491646A (zh) * | 2018-03-28 | 2018-09-04 | 青岛理工大学 | 一种隧道近距穿越重要建筑物爆破引起振动损伤的鉴定方法 |
CN112036039B (zh) * | 2020-09-01 | 2022-12-23 | 内蒙古科技大学 | 一种共轭传热材料热性能的高精度数值分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01145723A (ja) * | 1987-08-28 | 1989-06-07 | Hitachi Ltd | プログラム生成方法 |
JPH06148038A (ja) * | 1992-11-09 | 1994-05-27 | Hitachi Ltd | 振動解析装置 |
JP2005009323A (ja) * | 2003-06-16 | 2005-01-13 | Mitsubishi Heavy Ind Ltd | 多気筒内燃機関、そのクランク軸のクランク配置決定方法 |
JP2005202634A (ja) * | 2004-01-15 | 2005-07-28 | Keio Gijuku | 六面体有限要素数値解析方法 |
JP2005267214A (ja) * | 2004-03-18 | 2005-09-29 | Canon Inc | 有限要素法による制御プログラムおよび記憶媒体 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2070847A1 (en) * | 1992-06-10 | 1993-12-11 | Jyrki Kauppinen | Multicomponent analysis of ft-ir spectra |
US5966524A (en) * | 1997-07-24 | 1999-10-12 | Lucent Technologies Inc. | 3-D electromagnetic infinite element |
US6963824B1 (en) * | 1999-02-19 | 2005-11-08 | Davidson Joseph K | Method and apparatus for geometric variations to integrate parametric computer aided design with tolerance analyses and optimization |
ATE377209T1 (de) * | 2001-05-25 | 2007-11-15 | Parametric Optimization Soluti | Verbesserte prozesssteuerung |
JP3861012B2 (ja) * | 2002-01-30 | 2006-12-20 | 三菱重工業株式会社 | 多気筒内燃機関 |
US7283140B2 (en) * | 2005-06-21 | 2007-10-16 | Microsoft Corporation | Texture montage |
-
2006
- 2006-01-16 JP JP2006008057A patent/JP4599301B2/ja not_active Expired - Fee Related
- 2006-12-28 EP EP06027001.4A patent/EP1808775A3/en not_active Withdrawn
- 2006-12-28 US US11/646,285 patent/US20070179762A1/en not_active Abandoned
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01145723A (ja) * | 1987-08-28 | 1989-06-07 | Hitachi Ltd | プログラム生成方法 |
JPH06148038A (ja) * | 1992-11-09 | 1994-05-27 | Hitachi Ltd | 振動解析装置 |
JP2005009323A (ja) * | 2003-06-16 | 2005-01-13 | Mitsubishi Heavy Ind Ltd | 多気筒内燃機関、そのクランク軸のクランク配置決定方法 |
JP2005202634A (ja) * | 2004-01-15 | 2005-07-28 | Keio Gijuku | 六面体有限要素数値解析方法 |
JP2005267214A (ja) * | 2004-03-18 | 2005-09-29 | Canon Inc | 有限要素法による制御プログラムおよび記憶媒体 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009020619A (ja) * | 2007-07-10 | 2009-01-29 | Nippon Telegr & Teleph Corp <Ntt> | 多項式間距離算出装置、方法、プログラムおよび記録媒体、一変数最近実多項式算出装置、方法、プログラムおよび記録媒体 |
JP2010102680A (ja) * | 2008-09-25 | 2010-05-06 | Fujitsu Ltd | レイアウト評価装置、レイアウト評価プログラム、ダミールール生成装置及びダミールール生成プログラム |
JP2011233062A (ja) * | 2010-04-30 | 2011-11-17 | Shuzo Maeda | 3次元形状を3次元形状へ変換する方法 |
JP2012083958A (ja) * | 2010-10-12 | 2012-04-26 | Sumitomo Rubber Ind Ltd | 構造格子を用いたシミュレーション方法 |
JP2022506325A (ja) * | 2018-11-13 | 2022-01-17 | スーチョウ レインメド メディカル テクノロジー カンパニー リミテッド | 造影画像に基づき血管評定パラメータを取得する方法、装置及びシステム |
JP7162934B2 (ja) | 2018-11-13 | 2022-10-31 | スーチョウ レインメド メディカル テクノロジー カンパニー リミテッド | 造影画像に基づき血管評定パラメータを取得する方法、装置及びシステム |
Also Published As
Publication number | Publication date |
---|---|
EP1808775A3 (en) | 2013-05-01 |
EP1808775A2 (en) | 2007-07-18 |
US20070179762A1 (en) | 2007-08-02 |
JP4599301B2 (ja) | 2010-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4599301B2 (ja) | 設計支援装置及びコンピュータプログラム | |
US7647211B2 (en) | Constraint-based method of designing a route for a transport element | |
Boukharfane et al. | A combined ghost-point-forcing/direct-forcing immersed boundary method (IBM) for compressible flow simulations | |
US7668700B2 (en) | Adaptive distance field constraint for designing a route for a transport element | |
Masters et al. | Review of aerofoil parameterisation methods for aerodynamic shape optimisation | |
Colonna et al. | Numerical simulation of dense gas flows on unstructured grids with an implicit high resolution upwind Euler solver | |
Kim et al. | Analysis of hydroelasticity of floating shiplike structure in time domain using a fully coupled hybrid BEM-FEM | |
Firrone et al. | Passive control of vibration of thin-walled gears: advanced modelling of ring dampers | |
Gong et al. | Hybrid nearly singular integration for isogeometric boundary element analysis of coatings and other thin 2D structures | |
Shershnev et al. | HyCFS, a high-resolution shock capturing code for numerical simulation on hybrid computational clusters | |
Brahimi et al. | Blade dynamical response based on aeroelastic analysis of fluid structure interaction in turbomachinery | |
US9361413B1 (en) | Systems and methods for simulating contact between physical objects | |
Barrett et al. | Integrated free‐form method for aerostructural optimization of wind turbine blades | |
Scholten et al. | Uncoupled method for static aeroelastic analysis | |
Deng et al. | Reduced-order multiscale modeling of plastic deformations in 3D alloys with spatially varying porosity by deflated clustering analysis | |
Singh et al. | Mixed-variable optimization strategy employing multifidelity simulation and surrogate models | |
Chen et al. | Towards high-accuracy deep learning inference of compressible turbulent flows over aerofoils | |
Zhang et al. | Surrogate-based structural optimization design of large-scale rectangular pressure vessel using radial point interpolation method | |
US10878147B1 (en) | Systems and methods for simulating contact between physical objects | |
Tissera et al. | Computational fluid dynamics methods for hypersonic flow around blunted-cone-cylinder-flare | |
Fenrich et al. | Reliable multidisciplinary design of a supersonic nozzle using multifidelity surrogates | |
Vandevoorde et al. | Comparison of algorithms for unsteady flow calculations in inlet and exhaust systems of IC engines | |
Yuanqiu et al. | Identification of crystal orientation for turbine blades with anisotropy materials | |
Saldía et al. | Numerical simulations using tvd schemes of two-dimensional supersonic flow in chemical equilibrium | |
Jiang et al. | Hybrid central-upwind finite volume schemes for solving the Euler and Navier–Stokes equations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080819 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090417 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090626 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090824 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100115 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20100201 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100316 |
|
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: 20100910 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100927 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131001 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |