CN1034699C - 形状模拟方法 - Google Patents
形状模拟方法 Download PDFInfo
- Publication number
- CN1034699C CN1034699C CN91109020A CN91109020A CN1034699C CN 1034699 C CN1034699 C CN 1034699C CN 91109020 A CN91109020 A CN 91109020A CN 91109020 A CN91109020 A CN 91109020A CN 1034699 C CN1034699 C CN 1034699C
- Authority
- CN
- China
- Prior art keywords
- unit
- volume fraction
- volume
- instantaneous
- rate
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000000463 material Substances 0.000 claims description 47
- 239000007787 solid Substances 0.000 claims description 9
- 230000004941 influx Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 abstract description 10
- 239000000126 substance Substances 0.000 abstract description 7
- 210000004027 cell Anatomy 0.000 description 13
- 238000001556 precipitation Methods 0.000 description 12
- 210000000540 fraction c Anatomy 0.000 description 10
- 239000002245 particle Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000001681 protective effect Effects 0.000 description 8
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 6
- 229910052782 aluminium Inorganic materials 0.000 description 6
- 239000004411 aluminium Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000003475 lamination Methods 0.000 description 6
- 238000007747 plating Methods 0.000 description 5
- 238000005507 spraying Methods 0.000 description 5
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 4
- 235000017491 Bambusa tulda Nutrition 0.000 description 4
- 241001330002 Bambuseae Species 0.000 description 4
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 4
- 239000011425 bamboo Substances 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 239000002244 precipitate Substances 0.000 description 3
- 239000000758 substrate Substances 0.000 description 3
- 229910004298 SiO 2 Inorganic materials 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000839 emulsion Substances 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 241000721047 Danaus plexippus Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- AZDRQVAHHNSJOQ-UHFFFAOYSA-N alumane Chemical group [AlH3] AZDRQVAHHNSJOQ-UHFFFAOYSA-N 0.000 description 1
- 238000005267 amalgamation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000012456 homogeneous solution Substances 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01L—SEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
- H01L21/00—Processes or apparatus adapted for the manufacture or treatment of semiconductor or solid state devices or of parts thereof
- H01L21/02—Manufacture or treatment of semiconductor devices or of parts thereof
- H01L21/04—Manufacture or treatment of semiconductor devices or of parts thereof the devices having potential barriers, e.g. a PN junction, depletion layer or carrier concentration layer
- H01L21/18—Manufacture or treatment of semiconductor devices or of parts thereof the devices having potential barriers, e.g. a PN junction, depletion layer or carrier concentration layer the devices having semiconductor bodies comprising elements of Group IV of the Periodic Table or AIIIBV compounds with or without impurities, e.g. doping materials
- H01L21/30—Treatment of semiconductor bodies using processes or apparatus not provided for in groups H01L21/20 - H01L21/26
- H01L21/302—Treatment of semiconductor bodies using processes or apparatus not provided for in groups H01L21/20 - H01L21/26 to change their surface-physical characteristics or shape, e.g. etching, polishing, cutting
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Manufacturing & Machinery (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Power Engineering (AREA)
- Drying Of Semiconductors (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明为一种形状模拟方法,其特征在于:将解析领域分成若干个单元;定义每个单元中物质最初的体积率;对于各个单元,算出经过微小时间后,物质粒子的流入量及流出量;以算出流入量及流出量为基础,算出在微小时间内关于各单元的物质的体积率;根据所定的体积率的值,以等体积率面为基准模拟物质的形状。
Description
本发明涉及形状的模拟方法,特别是关于LSI等半导体装置加工形状的模拟方法。
图11A~11C是“Line-Profile ResistDeVelopment Simulation Techniques”:R.E.JEWETT etc.;以及Polymer Engineering andScience,June 1977,Vol.17,No.6,P.381~384中所示的现有的形状表示模式的概念图,第11A图是弦线模式,笫11B图是单元分割模式,第11C图显示的是射线描记模式。
在第11A图的弦线模型中,物质的形状用短线(112)互相连接各接线点(111)的方式来表示,在微小时间内,使各线段移动。以表示物质形状的时间变化。一般情况下,使接点朝看相互连接的两根短线所形成的角度的二等分线方向移动,因而线段也就移动了。
在笫11B图的单元分割模型中,用小的正方形单元(113)表示物质的形状,以单元的添加和消除来表示形状的时间变化。举个例子,在这个模型中,根据(i、j、k)的单元o(真空),1(硅)及2(酸化物)内所代表内容的记忆方法,可以简单而且明了地记忆前道工序的物质形状。
此外,在第11C图的射线描记模型中,同喷流模型一样,用在相互连接点(114)上连接短线的模式来表示物质的形状,接线点(114)沿着光线前进的方向移动。
然而,喷流及射线描记模型中,一旦产生第12A图那样的过长线段,就应该象第12B图那样分割成两个线段(122)及(123);当出现如第图12C那样的过短线段(124)或者是第12E图那样的环形(1.25)时,就必须笫12D图或笫12F图那样进行消除特异点的处理。而且这种处理每经过一个微小的时间就必须进行。在计算立体形状的情况下,近行这种处理是很复杂的时间也很长,这就存在一定的问题。此外,在上述这些模型由于使用了一维的近似,因此要想正确地模拟具有粒子量变化的腐刻或沉积是极其困难的。
另外,在这些射线自动描记模型中,也存在看形状的记忆方法很难的问题。例如,象第13A图显示的那样,使硅(131)的凹部(132)沉积酸化物随着堆积的进行,就成为图13B所示的酸化物(133)表面的线段(134)和(135)交错那样的情形这样就形成了环(136),要用前面所提到的特异点处理进行消除,最终留下了必要的线段,离实际的形状就相差太远了。
另一方面,在单元分割模式中,为了不形成倾斜表面,形状是呈阶段形成的;特别是为了高精度地进行立体形状的模拟,单元数就需要很多,而且还存在计算时间增大的问题。
本发明的目的就是为了解决这样一些问题,目的是提供一种能够进行高精度而且快速模似立体形状的模拟方法。
在本发明的形状的模拟方法中,把解析领域分割成多个单元,定义出各单元物质的最初的体积率;每经过极短的时间就算出各单元中物质粒子的流入量和流出量;每经过极短的时间根据算出的流入量和流出量各单元物质的体积率;根据具有既定的体积率的值的等体积率面模拟物质的形状。
在本发明中,从每经过极短的时间的各单元的物质粒子的流入量和流出量计算出单元的物质的体积率,根据等体积率面模拟物质的形状。
图1是关于这个发明的一个形状模拟方法实例的流程图。图2是沉淀计算的流程图。图3A~3C是真对图2的流程图。具体地计算物质的体积率的实例。图4A~4D显示了解析领域中几种物质混一起其边界的体积率的关系。图5是腐刻计算的流程图。图6A~6C是真对图5的流程图。具体地计算物质体积率的实例。图7A和图77B分别显示了各单元的体积率,及关于图7A的表所表现的物质的形状。图8A~8F是立体角的计算方法。图9A~9E是根据这个发明所进行的形状模拟的连续工序,图10A是图9所用的靶和极板位置关系。图10B和图10C分别显示了图10A极板中央的槽及偏离中央的槽的铝层堆积情况的模拟结果。图11A~11C是最初的形状模拟方法的概念图。图12A~12F是关于弦线模式的特异点处理。图13A~13C显示的是采用弦线模式进行处理模拟时存在的问题。
此外,各图中有相当的部分都采用同一符号,同一注脚。
以下是这个发明的实施例,并根据附图加以说明。
图1是表示本发明的一个实施例的形状模拟方法的流程图。首先,在步骤S1中输入处理器参数,然后,在步骤S2中,将解析领域分割成若干个正方体单元,设各单元的x方向为i,y方向为j,z方向为k,用(i、j、k)来表示。
接下来,在步骤S3对处理工序进行判断,如果是沉淀工序就进入步骤S4,进行沉淀方面的计算,如果是腐刻工序就进入步骤S5进行腐刻方面的计算,经过步骤S4和S5。在沉淀工序或腐刻工序终了时,各单元物质的体积率就被算出来了。在这之后,在步骤S6存贮各单元体积率,根据需要在步骤S7根据既定体积率值所确定的等体积率而来表示物质的形状。步骤S8是在一连串的程序全部终了时,判断是否返回执行S3~S7程序。
在这里,图2显示的是关于步骤S4中的沉淀计算的流程图。下面,参照图2具体说明一下沉淀的计算方法。图3A是用x-z二维断面表示的沉淀工序前的状态,其解析领域被分割成若干个正方体单元(i、j、k)。各单元内的数字表示的是物质的体积率Ct(i、j、k)。首先,在图2中的步骤11导出表面单元,这个表面单元如单元(i、j、k)周围的单元的体积率Ct(i+1j,k),Ct(i,j+1,k),Ct(i,j,k+1)内的不足0.5时,那么,这个单元(i、j、k)就作为表面单元。但是,体积率不足0.5的单元不作为表面单元,图3A中标着斜线的单元为表面单元。
下一步。在步骤S12进行沉淀种类的判别。如果是喷镀沉淀,则执行步骤S13,算出中间电极和晶片的位置关系、及从立体角算出沉淀速度。如果是均质沉淀,就执行步骤S14,由立体角计算出沉淀速度。然后,根据沉淀速度求出单位时间内由被堆积物质x穿过表示单元(i,j,k)的而形成的该单元的流入体积率Ri、j、k。
此外,执行步骤S15,根据下面的公式计算出在微小时间Δt秒之后单元(i、j、k)中的全部物质的体积率Ct+Δt
Ct+Δt(i、j、k)=Ct(i,j,k)+Rijk·Δt各单元Δt秒后的体积率表示在图3B中。
在这之后,执行步骤S16,进行调整各单元体积率的处理,也就是说,在图3B中画圈的单元(i=1,k=2)及(i=2,k=3)中全物质的体积率超过1的情况下,设这些单元周围的单元(i+1,j,k)、(i,j+1,k)、(i,j,k+1)的体积率都在0.5以下的单元为α,β,γ……,这些单元与单元(i,j,k)邻接的面的面积为Sβ,Sβ,Sγ……,这样,
Ct+Δt(i、j、k)←1
Ct+Δt(α)←Ct+Δt(α)+ηSα/(Sα+Sβ+Sγ+……) ……)
Ct+Δt(γ)←Ct+Δt(γ)+ηSγ/(Sα+Sβ+Sγ+……)。
其中η=Ct+Δt(i,j,k)-1。
同样,物质X的情况, Sγ+……) +Sγ+……) +Sγ+……)。
经过这样的调整处理,全物质的体积率就变成图3C的样子。在步骤S17和S18对沉淀时间t0进行判定。如果经过微小时间Δt后仍少于时间t0,则以上一连串的步骤S11~S16被被重复执行。
另外,例如在解析领域里,有三个物质X1、X2及X3混合的情况下,图4A~4C分别表示的各物质X1、X2及X3的体积率(CX1(i,j,k)CX2(i,j,k),CX3(i,j,k)进行加法运算,图4D显示了全物质的体积率C(i,j,k)。
关于图1的步骤5腐刻计算的流程图表示在图5中。下面,参照图5具体说明腐刻的计算方法。图6A显示的是腐刻工序前x-z二维断面的状态,解析领域被分割成许多的正方体单元(i,j,k)。首先,在图5的步骤S21导出表面单元。这个表面单元如单元(i,j,k)的体积率Ct(i,j,k)>0且周围单元的体积率Ct(i+1,j,k)、Ct(i,j+1,k)、Ct(i,j,k+1)为0时,这个单元就被定为表面单元。图6A中画有斜线的单元为表面单元。
下一步,步骤S22进行处理种类的判别,在印刷的情况下,执行步骤S23·进行光强度及感光剂分布的计算;接看执行步骤S2S24,计算显象速度;之后,执行步骤S25。根据表面积,计算出单位时间的流出重。另外,在进行等方性腐刻的情况下,由步骤S22直接进入步骤S25计算流出量。在进行异方性腐刻的情况下执行步骤S26,从腐刻的角度分散及立体角计算出单位时间的流出量。然后,根据步骤S25或S26算出的流出量,求出单位时间内通过被腐刻物质X的表面单元(i,j,k)的面从这个面流出的体积率Rijk。
然后,执行步骤S27,根据下面的公式计算出微小时间Δt秒后的单元(i,j,k)的物质X的体积率Cx t+Δt(i,j,k)
图6B中表示了各单元Δt秒后的体积率。
在这之后,步骤S28,进行各单元体积率的调整处理,也就是说,如图6B中画圈的单元(i=2,k=2)那样,在
(i,j,k<0的情况下,设该单元(i,j,k)周围的单元(i+1,j,k)、(i,j+1,k)及(i,j,k+1)内的物质X的体积率在0.5以上的单元α,β,γ,……与单元(i,j,k)邻接的面的面积Sα、Sβ,Sγ……,那么, +Sγ+……) +Sγ+……) +Sγ+……)。
其中 。
根据上面的式子,单元α、β、γ……的体积率为负数时,该单元的体积率就被认为是0 。经过这样的调整处理,全物质的体积率就变成图6C的样子。
在解析领域中许多物质X1、X2、X3……混在一起的情况下,根据各物质进行以上的处理,对于单元(i,j,k)中各物质X1、X2、X3……的体积率CX1(i,j,k)CX2(ij,k)、CX3(i,j,k)……相加后算出的值就被定为全物质的体积率Ct+Δt(i,j,k)
以上一连串步骤S21~S28以微小的时间间隔重复执行,直到步骤S29和S30判断出已超过腐刻时间t0为止。
假如上述所做的沉淀计算及腐刻计算中的微小时间为Δt,单位时间流入或流出的体积率Rijk的最大值为Rmax,最好
Δt=1/2·Rmax。
但是,沉淀时间或腐刻时间t0的最后的时间步骤的时间间隔Δtend必须定为:
Δtend=t0-η·Δt。
按照以上的作法,标出沉淀工序或腐刻工序终了后的体积率,在图1的步骤S6把它存入附图中未示出的文件,以及在步骤S7根据体积率表示物质的形状。此时,对于各单元(i,j,k)体积率C(i,j,k)采用线性内插的方法进行近似运算。根据(比方说)C(i,j,k)=0.5的等体积率而来表现物质的形状。图7A和图7B中示出了表示形状的一个实例。对图7A中的各单元的体积率进行内插近似以后,就得到了图7B中所示的物质形状。
接下来,说明图2的步骤S13,S14及图5的步骤S26所采用的立体角的计算方法的一个实例。首先,单元(i、j、k)有六个面,有必要计算各面入射的物质量、即粒子数。计算一个面的入射粒子数时,要象图8A和图8B显示的那样把这个面的中心点O作为起始点确定矢量V=(Vx,Vy,Vz)。而,而,Vx=sinθ·cos
Vy=sinθ·sin
Vz=-cosθ
0≤θ≤π/2
0≤≤2π。
此时,微小立体角ΔΩ给定为:
ΔΩ =sinθdde。
O方向的分割数为Nθ、方向的分割数为N时,
O=π(I-1/2)/2Nθ
=2π(J-1/2)/N。而1≤I≤Nθ1≤J≤N。
这个矢量V通过被物质来填充的其他的单元时,由其方向粒子是不能入射的。这里讲叙一种能发现矢量V从哪个单元通过、计算时间又比较短的方法。
首先,如图8C所示,把x y平面分割成四个象限R1~R4,判别出矢量V通过的单元在哪个象限。在这里,说明一下单元位于象限R1的情况,即由于Vx>0,Vy>0,Vz<0,假设在单元(i,j,k)的下面通过的单元定为(u、v、w),
则u>i,v>j,w≤k成立。
而且,粒子的起始点定为γ0=(x0,y0,z0),用粒子之速度来判断矢量V,t秒后的粒子位置γ=(x,y,z)为:
γ=V·t+γ0
x=Vx·t+X0
y=Vy·t+y0
z=Vz·t +z0
另外,如图8D和图8E所示,相邻单元的边界的坐标用BX(i),BY(j),BZ(k)表示,粒子通过单元(i,j,k)的边界BX(i+1),BY(j+1),BZ(k)的时刻分别为:
tx(i+1)=〔BX(i+1)-X0〕/Vx,
ty(j+1)=〔BY(j+1)-y0〕/Vy,
tz(k)=〔BZ(k)-z0〕/Vz,
在图8D和图8E所示的例子中,由于
tz(k)<ty(j+1)<tx(i+1),
0<t<tx(i+1)时刻的情况下,x方向单元的代号为i,0<t<ty(j+1)时刻的情况下,y方向单元的代号为j,在这种情况下,到时刻t变为ty(j+1)为止,z方向的单元边界坐标以BZ(k-1)、BZ(k-2)、……变动,如果时刻t超过了ty(j+1),下一次y方向的单元边界坐标一个个地移动,并与tx(i+1)和ty(j+2)进行比较。
通过各单元边界处的时刻tx、ty和tz如图8F的直线所示。根据图8F,矢量V通过的单元为
(i,j,k-1)、(i,j,k-2)、(i,j+1,k-2)、
(i,j+1,k-3)、(i,j+2,k-3)、(ij+2,k-4)、
(i,j+3,k-4、(i+1,j+3,k-4)、(j+1,j+3,k-5)、……
同时地写出位于x、y、z方向上的单元的代号的话当然很好。但是,此时,考虑了tx、ty和tz各周期的大小,首先tz开始移动,然后是ty,最后是tx移动,这对于研究同时置位的单元代号是很好的。
另外,判别位于最高位置的填满物质的单元、即正坐标最小的单元,与这个单元相比较位于上面的单元(一正方向),就不考虑有关立体角的判断。由此,计算时间更加缩短。
图9A~图9E显示的是采用本发明的一个形状模拟的连续工序的例子。这个模拟过程是根据(1)保护层显影,(2)等方性腐刻(3)异方性腐刻,(4)保护层去除及(5)喷镀沉淀的连续进行来形成触孔和铝层的。
(1)保护层显影
在君9A中边的厚度d1=1μm计算使用一个1μm×1μm的矩形屏蔽、边界厚度为d1=1μm的保护层(91)表面的光强度分布,下一步计算保护层(91)中的感光剂浓度,最后根据感光计的浓度算出显象的速度。据此,可以决定每个单元的显象速度(腐刻规度)。在这个保护层显象的模拟中,由于从厚度d2=1μm的SiO2基板(92)来的反射光干涉效果,可以进行明确的模拟。
(2)等方性腐刻
在图9B所示的等方性腐刻中,只有SiO2基板本身被腐刻,这是根据表面单元的面积计算出从单元流出的SiO2的量。
(3)异方性腐刻
在图9C所示的异方性腐刻中,腐刻剂从上方垂直流向下向,来模拟SiO2基板的腐刻。
(4)除去保护层
在图9D中,等方性腐刻的保护层全部除去。
(5)喷镀沉淀
在图9E中,中间电极的中央有接触孔存在,由于铝的表面扩散效果,形成铝层。
图10A显示的是上述的喷镀模拟中所使用的靶(94)与极板(95)之间的位置关系。靶是一个直径为300mm的圆板,距它的中心d3=50mm处形成一个宽度为d4=90mm的侵蚀孔(95a),在靶(94)的下方d5=85mm处放置了一个直径为150mm的极板(95)。
上述的喷镀沉淀模拟中,触孔位于板的中部,而图10A所示的位于板中央的槽(95a)和距中央d6=50mm处的槽(95b)堆积铝层后的模拟结果分别显示在图10B和图10C中,这两个槽宽2μm,深2μm。板中央的槽(95a)形成左右对称的铝层,而偏离板中心的槽形成非对称的铝层,这是很明显的模拟结果。
在上述的实例中采用了正交的网格式分割解析领域也可以分割成三角形或多角形。
正如以上所说明的,这个发明涉及形状的模拟方法,解析领域被分成若干个单元,各单元内物质初期的体积率的定义,关于各单元经过微小时间后物质粒子的流入量和流出量的计算,根据所算出的流入量和流出量算出经过微小时间后物质的体积率。根据所定体积率值的等体积率而来模拟物质的形状,使高精度且高速度地模拟立体形状成为可能。
Claims (6)
1.一种描述经受改变表面形状的处理的物体表面的方法,其特征在于:
将一个解析体积分成许多单元,所述解析体积中包括一个具有表面的物体以及该解析体积的体积填充部分,所述的每个单元均具有一个单元体积,
定义每个单元最初的体积率,每个单元最初的体积率等于该单元中的物体体积与该单元的体积之比,
在处理期间t0内每隔时间周期Δt计算在改变物体表面形状的处理期间每个单元的物质流入速率及流出速率,这里t0远大于Δt;
每隔时间周期Δt以最初体积率和各个单元的物质流入速率及流出速率计算出每个单元的体积率,
在经过时间to之后,将体积率为预定值的单元连接起来以便描述物体的表面,以及
以图形方式将所描述的表面显示出来。
2.如权利要求1所述的方法,其特征在于:
将每个所具有的体积率至少为0.5并且具有一个体积率小于0.5的邻接单元的单元视为表面单元,
在进行向物体的表面加入物质的处理期间,每经过一个时间周期Δt即计算每个表面单元中的瞬时体积率,
将瞬时体积率超过1的每个表面单元的体积率置为1,并将超过1的体积率部分从瞬时体积率分配至体积率不大于0.5内邻接单元中。
3.如权利要求1所述的方法,其特征在于:
将体积率至少为0并且具有一个体积率为0的邻接单元的每个单元视为表面单元,
在进行从物体表面去除物质的处理期间,每经过时间周期Δt即计算每个表面单元的瞬时体积率,以及
将瞬时体积率为已变成负值的每个表面单元的体积率置为0,将将这些负瞬时体积率分配至体积率至少为0.5的邻接单元。
4.如权利要求1所述的方法,其特征在于:在检索出每个单元的体积率之后,通过内插方法对体积率进行近似计算。
5.如权利要求1所述的方法,其特征在于:每个单元的物质流入率和流出率通过使用立体角来计算。
6.如权利要求5所述的方法,其特征在于:处于最高位置的单元从体积率为1的单元中确定;当确定上述立体角时,位于最高位置的单元以上的单元不在考虑之列。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP254113/90 | 1990-09-26 | ||
JP2254113A JP2507696B2 (ja) | 1990-09-26 | 1990-09-26 | 形状シミュレ―ション方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1060167A CN1060167A (zh) | 1992-04-08 |
CN1034699C true CN1034699C (zh) | 1997-04-23 |
Family
ID=17260408
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN91109020A Expired - Fee Related CN1034699C (zh) | 1990-09-26 | 1991-09-19 | 形状模拟方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US5293557A (zh) |
JP (1) | JP2507696B2 (zh) |
KR (1) | KR940010503B1 (zh) |
CN (1) | CN1034699C (zh) |
DE (1) | DE4132102A1 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2851447B2 (ja) * | 1991-03-08 | 1999-01-27 | 三菱電機株式会社 | 形状シミュレーション方法 |
JP3526909B2 (ja) * | 1993-09-21 | 2004-05-17 | 株式会社ルネサステクノロジ | 形状シミュレーション方法 |
US5812435A (en) * | 1993-09-21 | 1998-09-22 | Mitsubishi Denki Kabushiki Kaisha | Shape simulation method allowing simulation of processed shape during steps of manufacturing a semiconductor device in a short period of time |
US6138052A (en) * | 1997-02-10 | 2000-10-24 | Betalaser Mike, Inc. | Product forming apparatus having computer-based device for non-contact gauging of a product size |
JPH1187197A (ja) * | 1997-09-14 | 1999-03-30 | Toshiba Corp | シミュレーション方法及びこの方法を実施するシミュレーション装置 |
US6091845A (en) | 1998-02-24 | 2000-07-18 | Micron Technology, Inc. | Inspection technique of photomask |
KR100280555B1 (ko) * | 1999-01-27 | 2001-01-15 | 김영환 | 모디파이드 셀 모델을 이용한 식각 및 증착 모사방법 |
SE0101984D0 (sv) | 2001-05-31 | 2001-05-31 | Skf Ab | A device, computer program product and method for indicating a function deviation of one or more details of manufacturing equipment using frequency component analyses |
US7419418B2 (en) * | 2003-08-26 | 2008-09-02 | Ormond, Llc | CNC abrasive fluid-jet milling |
EP1860582B1 (en) | 2005-03-17 | 2011-09-14 | Fujitsu Ltd. | Simulation apparatus, simulation method, simulation program, and computer readable recording medium in which that program has been recorded |
JP6089590B2 (ja) * | 2012-10-30 | 2017-03-08 | 大日本印刷株式会社 | 形状変化シミュレーション装置および方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4821214A (en) * | 1986-04-17 | 1989-04-11 | Brigham Young University | Computer graphics method for changing the shape of a geometric model using free-form deformation |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
UST944007I4 (zh) * | 1973-08-27 | 1976-03-02 | ||
US4534003A (en) * | 1981-08-24 | 1985-08-06 | At&T Bell Laboratories | Optimized reaction injection molding |
US4463380A (en) * | 1981-09-25 | 1984-07-31 | Vought Corporation | Image processing system |
US4882763A (en) * | 1984-12-31 | 1989-11-21 | The Standard Oil Company | Method of making a rock-pore micromodel involving generation of an accurate and reliable template image of an actual reservoir rock pore system |
US4785399A (en) * | 1987-03-03 | 1988-11-15 | International Business Machines Corporation | Shaping geometric objects by cumulative translational sweeps |
JPH0616475B2 (ja) * | 1987-04-03 | 1994-03-02 | 三菱電機株式会社 | 物品の製造システム及び物品の製造方法 |
DE3830570A1 (de) * | 1987-09-08 | 1989-03-16 | Toshiba Machine Co Ltd | Berechnungsverfahren fuer die stroemungsanalyse beim spritzgiessen |
JPH01141021A (ja) * | 1987-11-27 | 1989-06-02 | Toshiba Mach Co Ltd | 溶融材料の金型成形における流動解析結果の表示方法 |
US5031127A (en) * | 1987-11-27 | 1991-07-09 | Toshiba Machine Co., Ltd. | Molten injection-molding method |
US5070469A (en) * | 1988-11-29 | 1991-12-03 | Mitsubishi Denki Kabushiki Kaisha | Topography simulation method |
US5067101A (en) * | 1988-11-29 | 1991-11-19 | Mitsubishi Denki Kabushiki Kaisha | Topography simulation method |
US5071597A (en) * | 1989-06-02 | 1991-12-10 | American Bank Note Holographics, Inc. | Plastic molding of articles including a hologram or other microstructure |
-
1990
- 1990-09-26 JP JP2254113A patent/JP2507696B2/ja not_active Expired - Fee Related
-
1991
- 1991-04-16 US US07/685,655 patent/US5293557A/en not_active Expired - Lifetime
- 1991-09-19 CN CN91109020A patent/CN1034699C/zh not_active Expired - Fee Related
- 1991-09-20 KR KR1019910016481A patent/KR940010503B1/ko not_active IP Right Cessation
- 1991-09-26 DE DE4132102A patent/DE4132102A1/de not_active Ceased
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4821214A (en) * | 1986-04-17 | 1989-04-11 | Brigham Young University | Computer graphics method for changing the shape of a geometric model using free-form deformation |
Non-Patent Citations (1)
Title |
---|
IEE,CAT.NO.87,0189-1 1987.10.1 Y.HIRAI ET ALI SYMPOSUM ON VLSI TECHNOLOGY * |
Also Published As
Publication number | Publication date |
---|---|
JPH04133326A (ja) | 1992-05-07 |
CN1060167A (zh) | 1992-04-08 |
DE4132102A1 (de) | 1992-04-09 |
KR940010503B1 (ko) | 1994-10-24 |
US5293557A (en) | 1994-03-08 |
JP2507696B2 (ja) | 1996-06-12 |
KR920007107A (ko) | 1992-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1034699C (zh) | 形状模拟方法 | |
Aref et al. | Vortex dynamics of the two-dimensional turbulent shear layer | |
Mason et al. | Trailing vortices in the wakes of surface-mounted obstacles | |
CN106503154B (zh) | 批量河道断面形态数据的自动化提取方法 | |
CN104626585A (zh) | 一种用于sla3d打印机的平面分割方法及装置 | |
CN106777917B (zh) | 水工建筑物对海事交管雷达遮蔽区域计算和影响评价方法 | |
CN108180918A (zh) | 一种点云测地路径正向跟踪生成方法及装置 | |
CN106840116A (zh) | 三维放样方法 | |
CN114330501B (zh) | 一种基于动态时间规整的轨迹样式识别方法及设备 | |
McDonald et al. | Kinematic simulation of the uncut chip thickness and surface finish using a reduced set of 3D grinding wheel measurements | |
CN1702824A (zh) | 形貌仿真系统、形貌仿真方法以及计算机产品 | |
CN1255751C (zh) | 用于沉积的信息处理方法和设备 | |
Pakkanen | The Temple of Zeus at Stratos: New observations on the building design | |
CN102930143A (zh) | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 | |
CN109812260A (zh) | 用于海上勘探的定向井轨迹设计方法和系统、电子设备 | |
CN1065625C (zh) | 图形测定装置 | |
CN106056124A (zh) | 一种基于三维曲面循环引导的浮雕提取方法 | |
CN206311103U (zh) | 三维放样杆件 | |
CN206311021U (zh) | 三维放样平台 | |
JPH11279759A (ja) | スパッタ形状シミュレーション方法 | |
JP2006091907A (ja) | 地形データ処理プログラム | |
JP4203601B2 (ja) | 地形データ処理プログラム | |
JP2006120167A (ja) | 地形データ処理プログラム | |
Bork | Geometry and scenography in the late Gothic choir of Metz Cathedral | |
Krakovsky | A model restoration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C15 | Extension of patent right duration from 15 to 20 years for appl. with date before 31.12.1992 and still valid on 11.12.2001 (patent law change 1993) | ||
OR01 | Other related matters | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |