JPH06209101A - Device/process simulating method - Google Patents

Device/process simulating method

Info

Publication number
JPH06209101A
JPH06209101A JP268193A JP268193A JPH06209101A JP H06209101 A JPH06209101 A JP H06209101A JP 268193 A JP268193 A JP 268193A JP 268193 A JP268193 A JP 268193A JP H06209101 A JPH06209101 A JP H06209101A
Authority
JP
Japan
Prior art keywords
discretized
physical quantity
contour lines
regions
region
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.)
Withdrawn
Application number
JP268193A
Other languages
Japanese (ja)
Inventor
Kenichiro Sonoda
賢一郎 園田
Kyohiko Kotani
教彦 小谷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP268193A priority Critical patent/JPH06209101A/en
Publication of JPH06209101A publication Critical patent/JPH06209101A/en
Withdrawn legal-status Critical Current

Links

Abstract

PURPOSE:To enable a device/process simulating method to be applied even in the rapid change in the physical quantity by a method wherein discrete lattice points and regions are formed as it filling up the space between adjacent contour lines or surfaces in arbitrary physical quantity in analytic region with polygons or polyhedrons to be analyzed. CONSTITUTION:Discrete lattice points 10 are determined along contour lines 9 while triangular discrete regions 10a are formed as it filling up the spaces between adjacent contour lines 9. That is, since the concentration distribution of the contour lines 9 is to be densified in the parts where the distribution rapidly changes, respective discrete regions 10a are made inevitably smaller so as to increase the density thereof. Accordingly, the device/process simulating method faithfully following after the change in physical quantity with high precision can be applied even within the analytic region including the parts in the rapid change in the physical quantity.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は、半導体デバイスの物
理的特性や製造プロセスに対して任意の物理量の変数が
与える影響を解析するためのデバイス/プロセスシミュ
レーション方法に関し、特に、解析領域内に離散化格子
点と離散化領域を発生させて解析するシミュレーション
方法の精度の改善に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a device / process simulation method for analyzing the influence of a variable of an arbitrary physical quantity on a physical characteristic of a semiconductor device or a manufacturing process, and more particularly to a device / process simulation method. The present invention relates to improvement in accuracy of a simulation method for generating and analyzing a lattice point and a discretized region.

【0002】[0002]

【従来の技術】一般に、デバイス/プロセスシミュレー
ション方法においては、方程式を離散化して解くため
に、解析領域内に離散化格子点と離散化領域を発生させ
ることが必要である。
2. Description of the Related Art Generally, in a device / process simulation method, it is necessary to generate discretized grid points and discretized regions in an analysis region in order to discretize and solve equations.

【0003】図4は、従来のデバイス/プロセスシミュ
レーション方法において発生させられた離散化格子点と
離散化領域を有するMOSFETの概略的な断面図であ
る。この図において、太い実線1−2−3で囲まれた領
域と太い実線4−5−6で囲まれた領域がソース/ドレ
イン領域を表わしている。太い実線2−3−6−5−8
−7で囲まれた領域は半導体基板領域である。細い実線
の交点10は離散化格子点であり、お互いに近接する離
散化格子点間を結ぶ細い実線によって囲まれた領域10
aが離散化領域である。この場合、離散化領域10aは
三角形の形状を有している。
FIG. 4 is a schematic sectional view of a MOSFET having discretized lattice points and discretized regions generated in a conventional device / process simulation method. In this figure, a region surrounded by a thick solid line 1-2-3 and a region surrounded by a thick solid line 4-5-6 represent source / drain regions. Thick solid line 2-3-6-5-8
The region surrounded by -7 is the semiconductor substrate region. The intersections 10 of the thin solid lines are the discretized grid points, and the region 10 surrounded by the thin solid lines connecting the discretized grid points adjacent to each other.
a is a discretized area. In this case, the discretized area 10a has a triangular shape.

【0004】[0004]

【発明が解決しようとする課題】従来のデバイス/プロ
セスシミュレーション方法においては、離散化格子点1
0と離散化領域10aは物理量の変化を考慮することな
く、それぞれの領域1−2−3,4−5−6,2−3−
6−5−8−7の各々を複数の離散化領域に分割するの
に都合のよいように、それらの領域の境界の形状に応じ
て離散化格子点および離散化領域を発生させていた。し
かし、物理量が大きく変化する部分において離散化格子
点の密度が小さくて各離散化領域の面積が大きい場合に
は、精度の高いシミュレーションが期待できない。した
がって、従来のデバイス/シミュレーション方法におい
ては、物理量が大きく変化する部分においては、一旦発
生させられた格子点と離散化領域を修正する必要があっ
た。
DISCLOSURE OF THE INVENTION In the conventional device / process simulation method, the discretized lattice points 1
0 and the discretized region 10a are the respective regions 1-2-3, 4-5, 2-3-, without considering the change of the physical quantity.
In order to conveniently divide each of 6-5-8-7 into a plurality of discretized regions, the discretized grid points and the discretized regions are generated according to the shape of the boundary between those regions. However, if the density of the discretized lattice points is small and the area of each discretized region is large in the portion where the physical quantity changes greatly, a highly accurate simulation cannot be expected. Therefore, in the conventional device / simulation method, it is necessary to correct the once generated grid points and the discretized area in the portion where the physical quantity changes greatly.

【0005】本発明は、このような先行技術の課題を解
消するためになされたもので、解析領域内の物理量の変
化の度合を考慮して離散化格子点と離散化領域を発生す
ることによって、精度の高いデバイス/プロセスシミュ
レーション方法を提供することを目的としている。
The present invention has been made in order to solve the above-mentioned problems of the prior art, and generates a discretized grid point and a discretized region in consideration of the degree of change of the physical quantity in the analysis region. It is intended to provide a highly accurate device / process simulation method.

【0006】[0006]

【課題を解決するための手段】半導体デバイスの物理的
特性や製造プロセスに対して任意の物理量の変数が与え
る影響を解析するための本発明によるデバイス/プロセ
スシミュレーション方法は、解析領域内の任意の物理量
の隣り合う等高線または等高面の間を多角形または多面
体で埋めつくすように離散化格子点と離散化領域を発生
させて解析することを特徴としている。
A device / process simulation method according to the present invention for analyzing the influence of a variable of an arbitrary physical quantity on a physical characteristic of a semiconductor device or a manufacturing process is performed by an arbitrary method within an analysis area. It is characterized by generating and analyzing discretized grid points and discretized regions so that adjacent contour lines or contour surfaces of physical quantities are filled with polygons or polyhedra.

【0007】[0007]

【作用】本発明によるデバイス/プロセスシミュレーシ
ョン方法においては、物理量が急激に変化する領域では
その物理量の等高線または等高面の間隔が小さくなる。
したがって、これらの等高線または等高面の間を埋める
ように離散化格子点と離散化領域を発生させるので、物
理量の変化の急激な領域をも含めて精度の高いシミュレ
ーション結果を得ることができる。
In the device / process simulation method according to the present invention, the interval between the contour lines or contour surfaces of the physical quantity becomes small in the region where the physical quantity changes rapidly.
Therefore, since the discretized grid points and the discretized regions are generated so as to fill the space between these contour lines or contour surfaces, it is possible to obtain a highly accurate simulation result including a region where the physical quantity changes rapidly.

【0008】[0008]

【実施例】図1は、本発明の一実施例によるデバイス/
プロセスシミュレーション方法を実行するためのシミュ
レーション装置を概略的に示すブロック図である。この
シミュレーション装置は、読取装置11,中央演算装置
12,出力装置13,記憶装置14を含んでいる。読取
装置11は、解析領域の形状やその領域中の物理量の変
化などを読取る。中央演算装置12は、読取装置によっ
て読取られた入力データに基づき、その入力データの入
力形式規則をチェックし、物理量の等高線を求め、離散
化格子点を発生し、離散化領域を発生する。出力装置1
3は、中央演算装置によって演算されたシミュレーショ
ン結果を出力する。記憶装置14は、入力データ,物理
量データ,等高線を表わすデータ,発生された離散化格
子点を表わすデータ,および発生された離散化領域を表
わすデータなどをデータ記憶部に記憶する。
FIG. 1 shows a device / device according to an embodiment of the present invention.
It is a block diagram which shows schematically the simulation apparatus for performing a process simulation method. This simulation device includes a reading device 11, a central processing unit 12, an output device 13, and a storage device 14. The reading device 11 reads the shape of the analysis area and changes in physical quantities in the area. Based on the input data read by the reading device, the central processing unit 12 checks the input format rule of the input data, obtains contour lines of physical quantities, generates discretized grid points, and generates discretized regions. Output device 1
3 outputs the simulation result calculated by the central processing unit. The storage device 14 stores input data, physical quantity data, data representing contour lines, data representing generated discretized grid points, data representing generated discretized areas, and the like in a data storage unit.

【0009】図3は、本発明の一実施例によるデバイス
/プロセスシミュレーション方法において解析されるM
OSFETの概略的な断面図を示している。太い実線で
囲まれた領域1−2−3および領域4−5−6はソース
/ドレイン領域である。太い実線2−3−6−5−8−
7で囲まれた領域は基板領域である。破線で表わされた
曲線9は、不純物濃度の等高線を表わしている。この図
からわかるように、離散化格子点10は等高線9に沿っ
て発生させられ、三角形の離散化領域10aは、隣り合
う等高線9の間を埋めつくすように形成される。
FIG. 3 shows the M analyzed in the device / process simulation method according to one embodiment of the present invention.
Figure 2 shows a schematic cross section of an OSFET. Regions 1-2-3 and 4-5-6 surrounded by thick solid lines are source / drain regions. Thick solid line 2-3-6-5-8-
The area surrounded by 7 is the substrate area. The curve 9 represented by a broken line represents a contour line of the impurity concentration. As can be seen from this figure, the discretized grid points 10 are generated along the contour lines 9, and the triangular discretized regions 10 a are formed so as to fill up between the adjacent contour lines 9.

【0010】すなわち、濃度勾配が急激に変化する部分
においては等高線が密に分布しているので、必然的に離
散化領域10aの各々が小さくなって、離散化領域10
aの密度が上昇することになる。その結果、物理量の変
化の度合が急峻な部分を含む解析領域においても、その
物理量の変化に忠実に追随した精度の高いデバイス/プ
ロセスシミュレーション方法が可能となる。
That is, since the contour lines are densely distributed in the portion where the concentration gradient changes abruptly, each of the discretized regions 10a inevitably becomes smaller and the discretized region 10a becomes smaller.
The density of a will increase. As a result, a highly accurate device / process simulation method that faithfully follows the change in the physical quantity is possible even in an analysis region including a portion where the degree of change in the physical quantity is steep.

【0011】図2は、図3に示されているような離散化
格子点と離散化領域を発生させる手順の一例を示すフロ
ーチャートである。このフローチャートにおける手順が
スタートすれば、まずステップS21において、入力デ
ータが入力形式規則を満たしているか否かが判断され
る。入力データが入力形式規則を満たしていなければス
テップS21aに進み、入力形式規則を満たしていない
旨を出力装置に表示してシミュレーション装置を停止さ
せる。入力データが入力形式規則を満たしていれば、ス
テップS22に進み、適当な離散化格子点を用いて不純
物濃度の等高線を表わす。このとき、一例としての不純
物濃度の等高線のみならず、電位ポテンシャルの等高線
などの他の等高線に基づく解析も可能であることはいう
までもない。また、描かれた等高線に対してスプライン
補間などの補間を行なってもよい。
FIG. 2 is a flow chart showing an example of a procedure for generating the discretized grid points and the discretized regions as shown in FIG. When the procedure in this flowchart starts, first, in step S21, it is determined whether or not the input data satisfies the input format rule. If the input data does not satisfy the input format rule, the process proceeds to step S21a, the fact that the input format rule is not satisfied is displayed on the output device, and the simulation device is stopped. If the input data satisfies the input format rule, the process proceeds to step S22, and the contour lines of the impurity concentration are expressed using appropriate discretized lattice points. At this time, it is needless to say that analysis based not only on the contour line of the impurity concentration as an example but also on the contour line of the potential potential is possible. Further, interpolation such as spline interpolation may be performed on the drawn contour lines.

【0012】等高線が定められれば、ステップS23に
進み、それらの等高線の曲率や隣の等高線との距離など
を考慮して分割し、その分割点を離散化格子点として決
定する。すべての等高線が分割されてすべての離散化格
子点が決定されれば、ステップS24へ進む。
When the contour lines are determined, the process proceeds to step S23, and the contour points are divided in consideration of the curvature of the contour lines and the distance to the adjacent contour line, and the division points are determined as discretized grid points. If all the contour lines are divided and all the discretized grid points are determined, the process proceeds to step S24.

【0013】ステップS24においては、すべての離散
化格子点に関して、或る離散化格子点からその離散化格
子点が属しない他の等高線上の離散化格子点のうちで最
も近い離散化格子点までを線分で結び、それによって離
散化領域10aを形成する。すべての離散化格子点から
の線分が形成されればステップS25へ進む。
In step S24, for all the discretized lattice points, from a certain discretized lattice point to the closest discretized lattice point among other discretized lattice points on other contour lines to which the discretized lattice point does not belong. Are connected by a line segment, thereby forming the discretized region 10a. If line segments from all the discretized grid points are formed, the process proceeds to step S25.

【0014】ステップS25においては、離散化領域1
0aのうち三角形以外の領域を探し、もし三角形以外の
領域が存在すればステップS26へ進む。ステップS2
6においては、三角形以外の離散化領域を最も短い対角
線で分割し、三角形の離散化領域を生成する。すなわ
ち、三角形以外の離散化領域がなくなってすべての離散
化領域10aが三角形になるまで、ステップS25とS
26が繰返される。そして、解析領域のすべてが三角形
の離散化領域10aによって分割されれば、離散化領域
の生成手順は終了する。
In step S25, the discretized area 1
A region other than a triangle is searched for in 0a, and if a region other than a triangle exists, the process proceeds to step S26. Step S2
In 6, the discretized area other than the triangle is divided by the shortest diagonal to generate the triangular discretized area. That is, steps S25 and S are performed until there are no discretized areas other than triangles and all the discretized areas 10a are triangular.
26 is repeated. Then, when all of the analysis regions are divided by the triangular discretized region 10a, the discretized region generation procedure ends.

【0015】なお、以上の実施例では、二次元の解析領
域を三角形の離散化領域10aに分割する方法を説明し
たが、三次元の領域を分割する場合には、以上の説明の
「等高線」を「等高面」に置換え、「三角形の離散化領
域」を「多面体の離散化領域」に変更すればよいことが
理解されよう。そのとき、多面体の離散化領域10aと
して、たとえば四面体や三角柱などを用いることができ
る。
In the above embodiment, the method of dividing the two-dimensional analysis region into the triangular discretized regions 10a has been described. However, when dividing the three-dimensional region, the "contour line" described above is used. It will be understood that can be replaced by "contour surface" and "triangular discretized region" changed to "polyhedral discretized region". At this time, for example, a tetrahedron or a triangular prism can be used as the polyhedral discretized region 10a.

【0016】さらに、以上の実施例の説明において、等
高線を規定する物理量として不純物濃度を例にして説明
したが、電気的ポテンシャル,磁気的ポテンシャル,温
度,圧力などの他の物理量を表わす等高線に関して解析
することも可能であり、それらの等高線の間隔は必ずし
も一様である必要はないことも理解されよう。
Further, in the above description of the embodiment, the impurity concentration is taken as an example of the physical quantity defining the contour line, but the contour line representing other physical quantities such as electric potential, magnetic potential, temperature and pressure is analyzed. It will also be appreciated that the contours may not necessarily be evenly spaced.

【0017】[0017]

【発明の効果】以上のように、本発明によれば、任意の
物理量の等高線または等高面を利用して離散化格子点を
発生させ、それに応じて離散化領域を形成するので、物
理量の急峻な変化にも追随し得る精度の高いデバイス/
プロセスシミュレーション方法を提供することができ
る。
As described above, according to the present invention, the discretized grid points are generated by utilizing contour lines or contour surfaces of arbitrary physical quantities, and the discretized areas are formed accordingly. Highly accurate device that can follow sharp changes /
A process simulation method can be provided.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の一実施例によるデバイス/プロセスシ
ミュレーション方法を実行し得る装置の一例を示すブロ
ック図である。
FIG. 1 is a block diagram showing an example of an apparatus that can execute a device / process simulation method according to an embodiment of the present invention.

【図2】本発明の一実施例によるデバイス/プロセスシ
ミュレーション方法の処理手順の一例を示すフローチャ
ートである。
FIG. 2 is a flowchart showing an example of a processing procedure of a device / process simulation method according to an embodiment of the present invention.

【図3】本発明のデバイス/プロセスシミュレーション
方法によって発生された離散化格子点と離散化領域の一
例を示すMOSFETの概略的な断面図である。
FIG. 3 is a schematic cross-sectional view of a MOSFET showing an example of discretized lattice points and discretized regions generated by the device / process simulation method of the present invention.

【図4】従来のデバイス/プロセスシミュレーション方
法によって発生された離散化格子点と離散化領域の一例
を示すMOSFETの概略的な断面図である。
FIG. 4 is a schematic cross-sectional view of a MOSFET showing an example of discretized lattice points and discretized regions generated by a conventional device / process simulation method.

【符号の説明】[Explanation of symbols]

1−2−3,4−5−6 ソース/ドレイン領域 2−3−6−5−8−7 半導体基板領域 9 物理量の等高線 10 離散化格子点 10a 離散化領域 11 読取装置 12 中央演算装置 13 出力装置 14 記憶装置 1-2-3, 4-5-6 Source / drain region 2-3-3-6-5-8-7 Semiconductor substrate region 9 Contour line of physical quantity 10 Discretized lattice point 10a Discretized region 11 Reader 12 Central processing unit 13 Output device 14 Storage device

─────────────────────────────────────────────────────
─────────────────────────────────────────────────── ───

【手続補正書】[Procedure amendment]

【提出日】平成5年1月14日[Submission date] January 14, 1993

【手続補正1】[Procedure Amendment 1]

【補正対象書類名】図面[Document name to be corrected] Drawing

【補正対象項目名】図1[Name of item to be corrected] Figure 1

【補正方法】変更[Correction method] Change

【補正内容】[Correction content]

【図1】 [Figure 1]

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】 半導体デバイスの物理的特性や製造プロ
セスに対して任意の物理量の変数が与える影響を解析す
るためのデバイス/プロセスシミュレーション方法にお
いて、 解析領域内の任意の物理量の隣り合う等高線または等高
面の間を多角形または多面体で埋めつくすように離散化
格子点と離散化領域を発生させて解析することを特徴と
するデバイス/プロセスシミュレーション方法。
1. A device / process simulation method for analyzing an influence of a variable of an arbitrary physical quantity on a physical characteristic of a semiconductor device or a manufacturing process, wherein adjacent contour lines or contours of arbitrary physical quantities in an analysis area are analyzed. A device / process simulation method characterized by generating and analyzing discretized lattice points and discretized regions so that the spaces between high surfaces are filled with polygons or polyhedra.
JP268193A 1993-01-11 1993-01-11 Device/process simulating method Withdrawn JPH06209101A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP268193A JPH06209101A (en) 1993-01-11 1993-01-11 Device/process simulating method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP268193A JPH06209101A (en) 1993-01-11 1993-01-11 Device/process simulating method

Publications (1)

Publication Number Publication Date
JPH06209101A true JPH06209101A (en) 1994-07-26

Family

ID=11536045

Family Applications (1)

Application Number Title Priority Date Filing Date
JP268193A Withdrawn JPH06209101A (en) 1993-01-11 1993-01-11 Device/process simulating method

Country Status (1)

Country Link
JP (1) JPH06209101A (en)

Similar Documents

Publication Publication Date Title
KR0160367B1 (en) Mesh generator and generating method
Bott A positive definite advection scheme obtained by nonlinear renormalization of the advective fluxes
KR960013371B1 (en) Shape simulation method
US4969116A (en) Method of generating discretization grid for finite-difference simulation
US5581489A (en) Model generator for constructing and method of generating a model of an object for finite element analysis
US4941114A (en) Triangular mesh generation method
JPH05138265A (en) Method and device for simulating molding
Owen et al. Surface mesh sizing control
KR20060034202A (en) Cad system and cad program
CN103970928A (en) Simulation Program, Simulation Method, And Simulation Device
JP2004038630A (en) Numerical analysis device using mixed lattice type solution-adaptive lattice approach
US7516054B2 (en) Method and apparatus for numerical calculation related to fluids, program for implementing the method, and storage medium storing the program
JPH06209101A (en) Device/process simulating method
JPH09223248A (en) Three-dimensional processing method for two-dimensional graphic data
Ng et al. Incremental tessellation of trimmed parametric surfaces
JPH1064834A (en) Diffusion simulating method
JP2926823B2 (en) Three-dimensional three-dimensional mesh forming method
JP4707216B2 (en) Structure analysis device
EP3511149A1 (en) Curvature deformation prevention design method for resin molded article, program, recording medium, and curvature deformation prevention design apparatus for resin molded article
JP3006374B2 (en) Forging die direction calculator
JPS62165271A (en) Fluid analysis supporting method
JP3805120B2 (en) Volume model generation system, volume model generation method, and computer-readable recording medium
JP2005056288A (en) Numerical analysis device using mixed lattice type solution-adaptive lattice approach
JPH0644350A (en) Method for generating solid-state element of finite element method mainly composed of hexahedron due to plane display
JPH0527901B2 (en)

Legal Events

Date Code Title Description
A300 Withdrawal of application because of no request for examination

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20000404