CN109558610A - 膜层刻蚀区域等效力学参数的计算方法和设备 - Google Patents
膜层刻蚀区域等效力学参数的计算方法和设备 Download PDFInfo
- Publication number
- CN109558610A CN109558610A CN201710883448.6A CN201710883448A CN109558610A CN 109558610 A CN109558610 A CN 109558610A CN 201710883448 A CN201710883448 A CN 201710883448A CN 109558610 A CN109558610 A CN 109558610A
- Authority
- CN
- China
- Prior art keywords
- stress
- mechanical parameter
- equivalent mechanical
- etch areas
- film layer
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 34
- 238000004088 simulation Methods 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 48
- 239000000463 material Substances 0.000 claims abstract description 36
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 4
- 238000009826 distribution Methods 0.000 description 9
- 230000001788 irregular Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 238000005530 etching Methods 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000264877 Hippospongia communis Species 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
-
- 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]
-
- 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/70—Manufacture or treatment of devices consisting of a plurality of solid state components formed in or on a common substrate or of parts thereof; Manufacture of integrated circuit devices or of parts thereof
- H01L21/71—Manufacture of specific parts of devices defined in group H01L21/70
- H01L21/768—Applying interconnections to be used for carrying current between separate components within a device comprising conductors and dielectrics
- H01L21/76801—Applying interconnections to be used for carrying current between separate components within a device comprising conductors and dielectrics characterised by the formation and the after-treatment of the dielectrics, e.g. smoothing
- H01L21/76829—Applying interconnections to be used for carrying current between separate components within a device comprising conductors and dielectrics characterised by the formation and the after-treatment of the dielectrics, e.g. smoothing characterised by the formation of thin functional dielectric layers, e.g. dielectric etch-stop, barrier, capping or liner layers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
-
- 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/70—Manufacture or treatment of devices consisting of a plurality of solid state components formed in or on a common substrate or of parts thereof; Manufacture of integrated circuit devices or of parts thereof
- H01L21/71—Manufacture of specific parts of devices defined in group H01L21/70
- H01L21/76—Making of isolation regions between components
- H01L21/762—Dielectric regions, e.g. EPIC dielectric isolation, LOCOS; Trench refilling techniques, SOI technology, use of channel stoppers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Power Engineering (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Micromachines (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种膜层刻蚀区域等效力学参数的计算方法和设备,属于刻蚀区域设计技术领域,其可至少部分解决现有的获取刻蚀区域等效力学参数的方法难度大或只能用于规则区域的问题。本发明的膜层刻蚀区域等效力学参数的计算方法包括:选取膜层刻蚀区域的至少一部分为分析区域;建立与分析区域对应的平面模型;以第一密度对平面模型进行网格划分;根据膜层材料的实际力学参数和第一密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第一模拟应力;计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
Description
技术领域
本发明属于刻蚀区域设计技术领域,具体涉及一种膜层刻蚀区域等效力学参数的计算方法和设备。
背景技术
显示基板的许多膜层都要刻蚀出开口(图案),形成一定刻蚀区域(如微小刻蚀区域)。不同刻蚀区域有不同的应力分布和变形,若变形不合适可能造成不良。因此,对设计出的刻蚀区域先要计算其应力分布和变形,若不合适则要对开口的个数、尺寸、形状、位置等进行调整。
以上应力分布和变形的一种计算方式为建立刻蚀区域的三维立体模型,并根据膜层材料的实际力学参数计算应力分布和变形,这种方式的建模、网格划分等非常复杂,耗时长,较难实用。另一种常用的计算方式为将刻蚀区域视为各向异性的平面盲板(即各方向力学性能不同的、无厚度的、无开口的完整二维板)并建立模型,再根据等效力学参数通过有限元法计算应力分布和变形,这种方式运算简单。其中,等效力学参数是指将刻蚀区域视为各方向具有不同性能的平面盲板时,该平面盲板在各方向上所具有的等效泊松比、等效弹性模量等。由于刻蚀区域并非平面盲板,故等效力学参数与材料的实际力学参数必然不同;而只有根据准确的等效力学参数,才能准确算出刻蚀区域的应力分布和变形。
获得等效力学参数的方法包括实验法和理论法。实验法是指制造标准试样产品并对其进行测试,但由于工艺限制,标准试样经常较难得到,且会带来成本、时间的增加。理论法是通过理论推算获得等效力学参数,但其只能用于规则的刻蚀区域,对大量不规则的刻蚀区域则不可用。
发明内容
本发明至少部分解决现有的获取刻蚀区域等效力学参数的方法难度大或只能用于规则区域的问题,提供一种实现简单且适用于不规则区域的膜层刻蚀区域等效力学参数的计算方法和设备。
解决本发明技术问题所采用的技术方案是一种膜层刻蚀区域等效力学参数的计算方法,其包括:
S01、选取膜层刻蚀区域的至少一部分为分析区域;
S02、建立与所述分析区域对应的平面模型;
S03、以第一密度对所述平面模型进行网格划分;
S04、根据膜层材料的实际力学参数和第一密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第一模拟应力;
S05、计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
优选的是,所述分析区域为长方形区域,其具有平行于X方向的第一边和第三边,以及平行于Y方向的第二边和第四边,所述第一边和第三边的长度为L1,所述第二边和第四边的长度为L2。
进一步优选的是,所述L1在0.001~0.5mm;
所述L2在0.001~0.5mm。
进一步优选的是,所述模拟边界条件包括第一边界条件和第二边界条件,其中,
所述第一边界条件为第一边、第三边、第四边固定,第二边沿X方向远离第四边发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
所述第二边界条件为第一边、第二边、第四边固定,第三边沿X方向远离第一边的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2。
进一步优选的是,所述等效力学参数包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX,其通过以下公式计算:
进一步优选的是,所述C1在0.0001~0.005;
所述C2在0.0001~0.005。
优选的是,所述网格为正方形,其边长在0.003~0.02mm。
优选的是,所述膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
优选的是,在所述步骤S04和S05之间,还包括:
以第二密度对平面模型进行网格划分,所述第二密度为第一密度的二倍;
根据膜层材料的实际力学参数和第二密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第二模拟应力;
比较所述第二模拟应力与第一模拟应力的差别是否超过阈值,若是则返回S03步骤,并重新选取更大的第一密度;若否则进入下一步。
进一步优选的是,所述第二模拟应力与第一模拟应力的差别为:第二模拟应力与第一模拟应力差值的绝对值与第二模拟应力的比值;
所述阈值为3%。
解决本发明技术问题所采用的技术方案是一种膜层刻蚀区域等效力学参数的计算设备,其包括:
区域选取单元,用于选取膜层刻蚀区域的至少一部分为分析区域;
模型建立单元,用于建立与所述分析区域对应的平面模型;
网格划分单元,用于对所述平面模型进行网格划分;
模拟应力计算单元,用于根据膜层材料的实际力学参数和当前的网格划分,通过有限元法分析平面模型在模拟边界条件下的模拟应力;
等效力学参数计算单元,用于计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
优选的是,所述分析区域为长方形区域,其具有平行于X方向的第一边和第三边,以及平行于Y方向的第二边和第四边,所述第一边和第三边的长度为L1,所述第二边和第四边的长度为L2。
进一步优选的是,所述模拟边界条件包括第一边界条件和第二边界条件,其中,
所述第一边界条件为第一边、第三边、第四边固定,第二边沿X方向远离第四边发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
所述第二边界条件为第一边、第二边、第四边固定,第三边沿X方向远离第一边的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2;
所述等效力学参数包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX,其通过以下公式计算:
优选的是,所述膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
优选的是,所述膜层刻蚀区域等效力学参数的计算设备还包括:
比较单元,用于比较第二模拟应力与第一模拟应力的差别是否超过阈值;其中,第一模拟应力为在第一密度的网格划分下计算得到的模拟应力,第二模拟应力为在第二密度的网格划分下计算得到的模拟应力,且所述第二密度为第一密度的二倍。
本发明的膜层刻蚀区域等效力学参数的计算方法中,先将刻蚀区域视为各向异性的平面盲板并用材料的实际力学参数(三维宏观材料实际具有的力学参数)计算其第一模拟应力,再用该第一模拟应力反推等效力学参数;因此该方法不需要制备标准试样产品,易实现,成本和时间消耗均小;而该第一模拟应力是通过有限元法计算的,故对不规则的刻蚀区域也适用;同时,经过试验证明,该方法的计算过程简单,且得到的等效力学参数与通过现有方法得到的等效力学参数差别很小,这表明该方法的计算结果准确,用其所得的等效力学参数可准确计算出刻蚀区域的实际应力分布和变形,并最终设计出合理的刻蚀区域。
附图说明
图1为本发明的实施例的一种膜层刻蚀区域等效力学参数的计算方法的流程图;
图2为本发明的实施例的方法中的一种从分析区域得到平面盲板的过程示意图;
图3为本发明的实施例的方法中的一种分析区域的示意图;
图4为本发明的实施例的方法中第一边界条件下平面盲板的状态示意图;
图5为本发明的实施例的方法中第二边界条件下平面盲板的状态示意图;
图6为本发明的实施例的方法中的一种刻蚀区域的示意图;
图7为本发明的实施例的方法中的另一种刻蚀区域的示意图;
图8为通过两种方式得到的图7的刻蚀区域的等效弹性模量的对比图;
其中,附图标记为:B1、第一边;B2、第二边;B3、第三边;B4、第四边。
具体实施方式
为使本领域技术人员更好地理解本发明的技术方案,下面结合附图和具体实施方式对本发明作进一步详细描述。
其中,X方向的等效泊松比VXY是指:当材料受到沿X方向的应力时,其在Y方向(垂直于X方向)产生的应变与在X方向产生的应变的比值。当然,Y方向的等效泊松比VYX的意义与此相对,在此不再详细描述。
实施例1:
本实施例提供一种膜层刻蚀区域等效力学参数的计算方法,其包括:
S11、选取膜层刻蚀区域的至少一部分为分析区域;
S12、建立与分析区域对应的平面模型;
S13、以第一密度对平面模型进行网格划分;
S14、根据膜层材料的实际力学参数和第一密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第一模拟应力;
S15、计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
本实施例的膜层刻蚀区域等效力学参数的计算方法中,先将刻蚀区域视为各向异性的平面盲板并用材料的实际力学参数(三维宏观材料实际具有的力学参数)计算其第一模拟应力,再用该第一模拟应力反推等效力学参数;因此该方法不需要制备标准试样产品,易实现,成本和时间消耗均小;而该第一模拟应力是通过有限元法计算的,故对不规则的刻蚀区域也适用;同时,经过试验证明,该方法的计算过程简单,且得到的等效力学参数与通过现有方法得到的等效力学参数差别很小,这表明该方法的计算结果准确,用其所得的等效力学参数可准确计算出刻蚀区域的实际应力分布和变形,并最终设计出合理的刻蚀区域。
实施例2:
如图1至图8所示,本实施例提供一种膜层刻蚀区域等效力学参数的计算方法。
其中,刻蚀区域是指:基板(如显示基板)的膜层中具有一定刻蚀开口(图案)的区域。
等效力学参数是指:若将刻蚀区域视为各向异性的平面盲板(即各方向力学性能不同的、无厚度的、无开口的完整二维板)时,能使该平面盲板表现出与实际刻蚀区域相同的力学行为的力学参数。具体的,等效力学参数可包括等效弹性模量和等效泊松比等。
本实施例的膜层刻蚀区域等效力学参数的计算方法具体包括以下步骤:
S21、选取膜层刻蚀区域的至少一部分为分析区域。
也就是说,在膜层具有开口的刻蚀区域(如微小刻蚀区域)中选取部分或全部,作为后续步骤的分析对象。
其中,若将整个刻蚀区域作为一个分析区域,则可得到对应整个刻蚀区域的等效力学参数。
若将刻蚀区域分为多个分析区域,且每个分析区域有不同的开口状况,则每个分析区域可分别计算得到不同的等效力学参数,从而刻蚀区域可视为由多个不同的平面盲板拼接形成。
若刻蚀区域如图6所示,是由多个周期性分布的如图2所示的分析区域组合而成的,则通过图2所示的分析区域计算得到的等效力学参数,通常与以图6整体为分析区域计算得到的等效力学参数相等。
优选的,分析区域为长方形区域,其具有平行于X方向的第一边B1和第三边B3,以及平行于Y方向的第二边B2和第四边B4,其中,第一边B1和第三边B3的长度为L1,第二边B2和第四边B4的长度为L2。
也就是说,优选如图2所示,选择长方形的区域为分析区域(故后续得到的平面盲板也是长方形的),这样的分析区域整体形状规则,便于计算。如图3所示,由于分析区域为长方形,故以上X方向和Y方向必然相互垂直,且L1与L2必然不相等(图中以L2更长为例)。
当然,以上描述只是限定分析区域整体形状为长方形,而分析区域中开口的形状、位置、个数、尺寸等都可以是不规则的,由此,本实施例的方法适用于不规则的刻蚀区域。
更优选的,L1在0.001~0.5mm;L2在0.001~0.5mm。
如果分析区域过大,则其划分出的网格也多,会给有限元计算带来困难,而若分析区域太小,则往往需要将刻蚀区域划分为成很多分析区域,同样计算麻烦,故其尺寸优选在以上范围内。
S22、建立与分析区域对应的平面模型。
如图2所示,将分析区域转换为平面模型,即转变为具有相同边界尺寸,在相同位置有相同开口,但却无厚度的二维板状模型(如Plane 42模型),作为有限元分析的对象。
S23、以第一密度对平面模型进行网格划分。
如图2所示,对以上建立的平面模型,以特定的第一密度对其进行网格划分;显然,此时只会对平面模型的实体部分进行网格划分,而不会其开口进行网格划分(开口处实际没有平面模型)。
优选的,网格为正方形,其边长在0.003~0.02mm。
有限元的网格通常可为正方形,其密度越大(尺寸约小)则计算结果越准确,但所需的计算量也就越大。经研究发现,以上尺寸范围内的网格,可在相对较低的计算量下达到较高的准确性。
S24、根据膜层材料的实际力学参数和第一密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第一模拟应力。
本步骤为根据膜材所用的材料的实际力学参数(三维宏观材料实际具有的力学参数,其通常已知),通过有限元法计算以上形式的平面模型在特定的边界条件(约束条件)下所受的应力。也就是说,对于与平面模型形状、尺寸相同的无厚度的平面板,假设其发生特定的变形,则通过有限元法计算其在该变形下所产生的应力,作为第一模拟应力。
优选的,膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
也就是说,当要计算在不同温度下的等效力学参数时,则本步骤中计算第一模拟应力所用的材料的实际力学参数,也应为材料在相应温度下的实际力学参数。由此,本实施例的方法可考虑温度的影响,从而得到更准确的等效力学参数。
优选的,第一边界条件为第一边B1、第三边B3、第四边B4固定,第二边B2沿X方向远离第四边B4发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
第二边界条件为第一边B1、第二边B2、第四边B4固定,第三边B3沿X方向远离第一边B1的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2。
也就是说,可设定两种不同的边界条件。其中,第一边界条为三条边固定,一条长度为L2的边沿X方向位移,此时通过有限元法可计算得到平面模型在X方向和Y方向的应力分别为σX1和σY1。而第二边界条件为三条边固定,一条长度为L1的边沿Y方向位移,此时通过有限元计算得到的平面模型在X方向和Y方向的应力分别为σX2和σY2。
其中,以上C1和C2即为平面模型在相应方向上的应变,其可根据需要人为设定。但优选的,C1在0.0001~0.005,C2在0.0001~0.005。这是因为在以上较小的应变范围内,绝大多数的材料都是弹性变形,有限元计算的结果与实际更吻合。
S25、以第二密度对平面模型进行网格划分,第二密度为第一密度的二倍;根据膜层材料的实际力学参数和第二密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第二模拟应力;比较第二模拟应力与第一模拟应力的差别是否超过阈值,若是则返回S23步骤,并重新选取更大的第一密度;若否则进入下一步。
如前,网格划分会对计算得结果产生影响,故为保证网格划分是合适的,可在用第一密度的网格划分计算得到一定的模拟应力(第一模拟应力)后,再用增大一倍的第二密度(同样面积内的网格数是第一网格密度的二倍)的网格划分计算模拟应力(第二模拟应力),并对两个模拟应力进行比较,若二者的差别超过预定值,则认为第一密度过小导致计算结果不够准确,应返回S23步骤选取更大的第一密度重新计算。
优选的,第二模拟应力与第一模拟应力的差别为:第二模拟应力与第一模拟应力差值的绝对值与第二模拟应力的比值;阈值为3%。
也就是说,以第一模拟应力为100%,若第二模拟应力超过103%或低于97%,则认为二者差别过大,应返回S23步骤。
S26、计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
本步骤中,在观念上(不用进行实际操作)将平面模型看成一个各向异性的、没有开口的平面板,并且计算该平面盲板应当具有怎样的力学参数(等效力学参数)才能使其在相同的边界条件(模拟边界条件)下受到相同的应力(第一模拟应力);或者说,认为该平面盲板在模拟边界条件下正好受到第一模拟应力,并据此反推其具有的力学参数,作为等效力学参数。
具体的,待计算的等效力学参数可包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX。
如前,第一边界条件下,平面模型沿X方向和Y方向的应变分别为C1和0,而沿X方向和Y方向的应力分别为σX1和σY1。
由此,可如图4所示,认为有一个各向异性的平面盲板,其在X方向和Y方向的应力分别为σX1和σY1,而应变分别为C1和0,由此可得其等效力学参数应满足以下方程:
类似的,根据第二边界条件,可得:
联立以上式1至式4,可得方程:
求解以上式5,可知等效力学参数可通过以下公式计算:
可见,通过以上式6,即可求出分析区域(为刻蚀区域或刻蚀区域的一部分)的等效力学参数。在该计算方法中,设置的模拟边界条件十分简单,故其所需的运算量和时间少,容易实现;而且,其计算过程又同时考虑到了泊松比和弹性模量两个参数的综合作用,因此其计算结果准确,误差小。
本实施例的验证:
验证1:
制备如图6所示的刻蚀区域的实际产品,其中,刻蚀区域其两边长分别为L1=100mm、L2=60mm的矩形,膜层厚度为0.03mm,膜层材料为SUS304,其中有多个排成矩阵的正方形开口,每个开口边长为0.05mm,相邻开口在X方向和Y方向的间距均为0.05mm mm。
对刻蚀区域的三条边进行约束,使一条长度为L1的边向X方向发生不同位移,实际测定沿X方向的应力;同时,在模拟边界条件下,按照本实施例的方法计算出沿X方向的第一模拟应力(其中网格为正方形,边长为0.025mm),二者具体结果如下表:
可见,在不同程度的变形下,实测应力与第一模拟应力之间的差别均很小(小于4%),这表明本实施例的方法算出的第一模拟应力与实际应力是很接近的,据此算出的等效力学参数准确性很高。
验证2:
对以上验证1的刻蚀区域,分别在不同的网格边长(网格形状均为正方形)下求得其等效力学参数,并用等效力学参数通过有限元法计算刻蚀区域在沿X方向发生0.01的应变时的应力,具体结果如下表:
网格边长(mm) | 计算的应力(MPa) |
0.02 | 5.2527 |
0.01 | 5.0891 |
0.005 | 4.9357 |
0.004 | 4.9127 |
0.003 | 4.8954 |
可见,在网格边长增大进10倍(0.003mm至0.02mm)的情况下,计算得到的应力的差别不足10%。这表明,本实施例的方法适用的网格密度范围很大,可选用较低的网格密度,从而在保证准确性的前提下明显缩短计算时间,提高效率。
验证3:
选取如图7所示的L=5mm的正方形刻蚀区域,其正中有正方形开口,开口与正方形边缘间的距离为t/2,膜层材料的实际弹性模量取194GPa,通过“Wang A J,McDowell D L,In-plane stiffness and yield strength of periodic metal honeycombs,ASME.J.Eng.Mater.Tech”提供的理论法计算其在不同的t值下的沿X方向(因是正方形,故X方向和Y方向等价)的等效弹性模量;同时,通过本实施例的方法计算出沿X方向的等效弹性模量(其中网格为正方形,其边长为0.01mm),具体结果如下表和图8所示:
可见,在不同程度的t值下,通过理论法和本实施例的方法计算得到的等效弹性模量的差均很小(小于7%),这表明本实施例的方法计算得到的等效力学参数准确性很高。
实施例3:
本实施例提供一种膜层刻蚀区域等效力学参数的计算设备,其包括:
区域选取单元,用于选取膜层刻蚀区域的至少一部分为分析区域;
模型建立单元,用于建立与分析区域对应的平面模型;
网格划分单元,用于对平面模型进行网格划分;
模拟应力计算单元,用于根据膜层材料的实际力学参数和当前的网格划分,通过有限元法分析平面模型在模拟边界条件下的模拟应力;
等效力学参数计算单元,用于计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
优选的,分析区域为长方形区域,其具有平行于X方向的第一边和第三边,以及平行于Y方向的第二边和第四边,第一边和第三边的长度为L1,第二边和第四边的长度为L2。
进一步优选的,模拟边界条件包括第一边界条件和第二边界条件,其中,
第一边界条件为第一边、第三边、第四边固定,第二边沿X方向远离第四边发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
第二边界条件为第一边、第二边、第四边固定,第三边沿X方向远离第一边的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2;
等效力学参数包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX,其通过以下公式计算:
优选的,膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
优选的,膜层刻蚀区域等效力学参数的计算设备还包括:
比较单元,用于比较第二模拟应力与第一模拟应力的差别是否超过阈值;其中,第一模拟应力为在第一密度的网格划分下计算得到的模拟应力,第二模拟应力为在第二密度的网格划分下计算得到的模拟应力,且第二密度为第一密度的二倍。
本实施例的膜层刻蚀区域等效力学参数的计算设备可执行以上实施例的方法,故可简单的获取准确的等效力学参数,且适用于不规则的刻蚀区域。
可以理解的是,以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式,然而本发明并不局限于此。对于本领域内的普通技术人员而言,在不脱离本发明的精神和实质的情况下,可以做出各种变型和改进,这些变型和改进也视为本发明的保护范围。
Claims (15)
1.一种膜层刻蚀区域等效力学参数的计算方法,其特征在于,包括:
S01、选取膜层刻蚀区域的至少一部分为分析区域;
S02、建立与所述分析区域对应的平面模型;
S03、以第一密度对所述平面模型进行网格划分;
S04、根据膜层材料的实际力学参数和第一密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第一模拟应力;
S05、计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
2.根据权利要求1所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述分析区域为长方形区域,其具有平行于X方向的第一边和第三边,以及平行于Y方向的第二边和第四边,所述第一边和第三边的长度为L1,所述第二边和第四边的长度为L2。
3.根据权利要求2所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述L1在0.001~0.5mm。
所述L2在0.001~0.5mm。
4.根据权利要求3所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,所述模拟边界条件包括第一边界条件和第二边界条件,其中,
所述第一边界条件为第一边、第三边、第四边固定,第二边沿X方向远离第四边发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
所述第二边界条件为第一边、第二边、第四边固定,第三边沿X方向远离第一边的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2。
5.根据权利要求4所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,所述等效力学参数包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX,其通过以下公式计算:
6.根据权利要求4所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述C1在0.0001~0.005;
所述C2在0.0001~0.005。
7.根据权利要求1所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述网格为正方形,其边长在0.003~0.02mm。
8.根据权利要求1所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
9.根据权利要求1所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,在所述步骤S04和S05之间,还包括:
以第二密度对平面模型进行网格划分,所述第二密度为第一密度的二倍;
根据膜层材料的实际力学参数和第二密度的网格划分,通过有限元法分析平面模型在模拟边界条件下的第二模拟应力;
比较所述第二模拟应力与第一模拟应力的差别是否超过阈值,若是则返回S03步骤,并重新选取更大的第一密度;若否则进入下一步。
10.根据权利要求9所述的膜层刻蚀区域等效力学参数的计算方法,其特征在于,
所述第二模拟应力与第一模拟应力的差别为:第二模拟应力与第一模拟应力差值的绝对值与第二模拟应力的比值;
所述阈值为3%。
11.一种膜层刻蚀区域等效力学参数的计算设备,其特征在于,包括:
区域选取单元,用于选取膜层刻蚀区域的至少一部分为分析区域;
模型建立单元,用于建立与所述分析区域对应的平面模型;
网格划分单元,用于对所述平面模型进行网格划分;
模拟应力计算单元,用于根据膜层材料的实际力学参数和当前的网格划分,通过有限元法分析平面模型在模拟边界条件下的模拟应力;
等效力学参数计算单元,用于计算等效力学参数,其中,在模拟边界条件下,等效力学参数能使与平面模型边界尺寸相同的各向异性的平面盲板达到第一模拟应力。
12.根据权利要求11所述的膜层刻蚀区域等效力学参数的计算设备,其特征在于,
所述分析区域为长方形区域,其具有平行于X方向的第一边和第三边,以及平行于Y方向的第二边和第四边,所述第一边和第三边的长度为L1,所述第二边和第四边的长度为L2。
13.根据权利要求12所述的膜层刻蚀区域等效力学参数的计算设备,其特征在于,
所述模拟边界条件包括第一边界条件和第二边界条件,其中,
所述第一边界条件为第一边、第三边、第四边固定,第二边沿X方向远离第四边发生C1L1的位移,在第一边界条件下沿X方向的应力为σX1,而沿Y方向的应力为σY1;
所述第二边界条件为第一边、第二边、第四边固定,第三边沿X方向远离第一边的方向发生C2L2的位移,在第二边界条件下沿X方向的应力为σX2,而沿Y方向的应力为σY2;
所述等效力学参数包括X方向的等效弹性模量EX和等效泊松比VXY,以及Y方向的等效弹性模量EY和等效泊松比VYX,其通过以下公式计算:
14.根据权利要求11所述的膜层刻蚀区域等效力学参数的计算设备,其特征在于,
所述膜层材料的实际力学参数为膜层材料在指定温度下的实际力学参数。
15.根据权利要求11所述的膜层刻蚀区域等效力学参数的计算设备,其特征在于,还包括:
比较单元,用于比较第二模拟应力与第一模拟应力的差别是否超过阈值;其中,第一模拟应力为在第一密度的网格划分下计算得到的模拟应力,第二模拟应力为在第二密度的网格划分下计算得到的模拟应力,且所述第二密度为第一密度的二倍。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710883448.6A CN109558610B (zh) | 2017-09-26 | 2017-09-26 | 膜层刻蚀区域等效力学参数的计算方法和设备 |
PCT/CN2018/107412 WO2019062727A1 (zh) | 2017-09-26 | 2018-09-25 | 膜层刻蚀区域等效力学参数的计算方法和设备 |
US16/479,182 US11250185B2 (en) | 2017-09-26 | 2018-09-25 | Method and apparatus for calculating equivalent mechanical parameters of film layer etching region |
EP18860934.1A EP3690681B1 (en) | 2017-09-26 | 2018-09-25 | Method and device for calculating equivalent mechanical parameters of film layer etching region |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710883448.6A CN109558610B (zh) | 2017-09-26 | 2017-09-26 | 膜层刻蚀区域等效力学参数的计算方法和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109558610A true CN109558610A (zh) | 2019-04-02 |
CN109558610B CN109558610B (zh) | 2021-01-29 |
Family
ID=65863285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710883448.6A Active CN109558610B (zh) | 2017-09-26 | 2017-09-26 | 膜层刻蚀区域等效力学参数的计算方法和设备 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11250185B2 (zh) |
EP (1) | EP3690681B1 (zh) |
CN (1) | CN109558610B (zh) |
WO (1) | WO2019062727A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112581868A (zh) * | 2020-12-09 | 2021-03-30 | 惠州市华星光电技术有限公司 | 柔性显示面板及其制备方法 |
WO2022082936A1 (zh) * | 2020-10-23 | 2022-04-28 | 东方晶源微电子科技(北京)有限公司深圳分公司 | 一种负显影光刻工艺的全芯片快速仿真方法、负显影光刻胶模型、opc模型及电子设备 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111625982B (zh) * | 2020-05-29 | 2023-07-25 | 成都赫尔墨斯科技股份有限公司 | 一种基于Ansys的便携导航设备的冲击与坠撞分析方法 |
CN112036068A (zh) * | 2020-09-07 | 2020-12-04 | 中国十七冶集团有限公司 | 一种基于bim技术的深基坑减震抗浮锚杆施工优化方法 |
CN113591353B (zh) * | 2021-08-03 | 2023-10-27 | 中国科学院上海应用物理研究所 | 一种基于ANSYS Workbench的多层膜热分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793833A (zh) * | 2005-12-02 | 2006-06-28 | 北京大学 | 薄膜应力测试方法及测试结构 |
CN103076116A (zh) * | 2012-12-29 | 2013-05-01 | 上海捷锦电力新材料有限公司 | 预制式电缆附件界面压力检测方法 |
CN105443102A (zh) * | 2014-07-14 | 2016-03-30 | 中国石油化工股份有限公司 | 一种水平井层内爆炸裂缝体建模方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7026249B2 (en) * | 2003-05-30 | 2006-04-11 | International Business Machines Corporation | SiGe lattice engineering using a combination of oxidation, thinning and epitaxial regrowth |
CN101479730B (zh) * | 2006-06-27 | 2011-06-08 | 日本电气株式会社 | 基板或电子部件的翘曲分析方法、基板或电子部件的翘曲分析系统及基板或电子部件的翘曲分析程序 |
CN104280155B (zh) * | 2013-07-14 | 2019-03-12 | 付康 | 一种确定多层膜薄膜应力的系统与方法 |
CN105512360B (zh) * | 2015-11-24 | 2018-11-23 | 广州兴森快捷电路科技有限公司 | 预测封装基板蚀刻后翘曲的方法 |
US10002222B2 (en) * | 2016-07-14 | 2018-06-19 | Arm Limited | System and method for perforating redundant metal in self-aligned multiple patterning |
US10339250B2 (en) * | 2016-11-29 | 2019-07-02 | Taiwan Semiconductor Manufacturing Company, Ltd. | Method of generating engineering change order (ECO) layout of base cell and computer-readable medium comprising executable instructions for carrying out said method |
-
2017
- 2017-09-26 CN CN201710883448.6A patent/CN109558610B/zh active Active
-
2018
- 2018-09-25 WO PCT/CN2018/107412 patent/WO2019062727A1/zh unknown
- 2018-09-25 EP EP18860934.1A patent/EP3690681B1/en active Active
- 2018-09-25 US US16/479,182 patent/US11250185B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793833A (zh) * | 2005-12-02 | 2006-06-28 | 北京大学 | 薄膜应力测试方法及测试结构 |
CN103076116A (zh) * | 2012-12-29 | 2013-05-01 | 上海捷锦电力新材料有限公司 | 预制式电缆附件界面压力检测方法 |
CN105443102A (zh) * | 2014-07-14 | 2016-03-30 | 中国石油化工股份有限公司 | 一种水平井层内爆炸裂缝体建模方法 |
Non-Patent Citations (3)
Title |
---|
FRIEDER LUCKLUM 等: "Highly insulating, fully porous silicon substrates for high temperature micro-hotplates", 《SENSORS AND ACTUATORS A: PHYSICAL》 * |
成广庆: "微尺度多层膜力学响应的有限元分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技I辑》 * |
王经伟: "多层膜材料力学性能的有限元模拟与分析", 《沈阳建筑大学硕士学位论文》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022082936A1 (zh) * | 2020-10-23 | 2022-04-28 | 东方晶源微电子科技(北京)有限公司深圳分公司 | 一种负显影光刻工艺的全芯片快速仿真方法、负显影光刻胶模型、opc模型及电子设备 |
CN112581868A (zh) * | 2020-12-09 | 2021-03-30 | 惠州市华星光电技术有限公司 | 柔性显示面板及其制备方法 |
CN112581868B (zh) * | 2020-12-09 | 2021-11-02 | 惠州市华星光电技术有限公司 | 柔性显示面板及其制备方法 |
Also Published As
Publication number | Publication date |
---|---|
US20190370423A1 (en) | 2019-12-05 |
EP3690681B1 (en) | 2024-01-17 |
US11250185B2 (en) | 2022-02-15 |
EP3690681A4 (en) | 2021-05-26 |
CN109558610B (zh) | 2021-01-29 |
EP3690681A1 (en) | 2020-08-05 |
WO2019062727A1 (zh) | 2019-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109558610A (zh) | 膜层刻蚀区域等效力学参数的计算方法和设备 | |
Yang et al. | A three-node triangular element with continuous nodal stress | |
CN103245437B (zh) | 一种确定非线性薄膜应力的系统与方法 | |
CN101629859B (zh) | 基于变形测量与数值反求确定薄膜应力的系统与方法 | |
Zheng et al. | A novel twice-interpolation finite element method for solid mechanics problems | |
CN104280155A (zh) | 一种确定多层膜薄膜应力的系统与方法 | |
CN104236768A (zh) | 一种基于有限元校正系数的钻孔法测量残余应力的方法 | |
CN107729605A (zh) | 一种基于板壳理论的剥层法残余应力测量值的修正方法 | |
Abbasi et al. | FLAC3D mesh and zone quality | |
Van der Kelen et al. | Identification of the full anisotropic flow resistivity tensor for multiple glass wool and melamine foam samples | |
CN104149340A (zh) | 基于3d打印技术获得金属零件性能的方法 | |
CN109323791A (zh) | 基于增量切割法的复合材料板材的残余应力分布测量方法 | |
CN107704657A (zh) | 一种由残余应力释放引起的板类零件翘曲变形的计算方法 | |
Negrea et al. | The elastic contact of a sphere with an elastic half-space, a comparison between analytical and finite elements solutions | |
Tang et al. | A novel four-node quadrilateral element with continuous nodal stress | |
CN104036085B (zh) | 用于cae仿真分析的复杂型面汽车声学部件的厚度统计方法 | |
CN107808021B (zh) | 基于cfd的流体设备的阻力计算方法 | |
CN105488263A (zh) | 预测封装基板阻焊后翘曲的方法 | |
JP3994884B2 (ja) | 薄膜の膜応力評価方法と弾性率・線膨張率同定方法 | |
CN103033467A (zh) | 二维多孔材料共面动态缓冲性能的测定方法 | |
Litvinenko et al. | Numerical methods for uncertainty quantification and bayesian update in aerodynamics | |
CN105093933B (zh) | 一种确定lpv变增益控制器的方法 | |
Liu et al. | Physics-informed sparse identification of bistable structures | |
Kianian et al. | Buckling analysis of functionally graded plates based on two-variable refined plate theory using the bubble finite strip method | |
CN106709167A (zh) | 一种基于相似度理论的装配有效性评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |