CN101515307A - 微纳电子器件结构的多物理场界面多尺度设计方法 - Google Patents
微纳电子器件结构的多物理场界面多尺度设计方法 Download PDFInfo
- Publication number
- CN101515307A CN101515307A CNA2009100247179A CN200910024717A CN101515307A CN 101515307 A CN101515307 A CN 101515307A CN A2009100247179 A CNA2009100247179 A CN A2009100247179A CN 200910024717 A CN200910024717 A CN 200910024717A CN 101515307 A CN101515307 A CN 101515307A
- Authority
- CN
- China
- Prior art keywords
- micro
- interface
- physical field
- electronic device
- nano electronic
- 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.)
- Pending
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种微纳电子器件结构的多物理场界面多尺度设计方法,先描述耦合机制,用分子动力学描述界面组织演化和材料缺陷特征;再用有限元为分子动力学模型提供热循环、弯曲、拉伸行为的边界条件;最后进行界面多物理场特征的数学模拟。本发明是针对微纳电子器件结构多物理场界面特别是封装界面特征研究多尺度模型的构架、针对微观界面尺度材料组织的演变、热缺陷的产生机制、材料缺陷如孔洞和裂纹等对微结构物理特性的影响形成,兼顾计算效率和科学准确性;实现微纳电子器件结构多物理场界面全性能特征的宏微纳观模拟,解决了微纳电子器件结构多物理场界面设计难题。
Description
技术领域
本发明涉及微纳观模拟,特指一种用于研究电子器件结构多物理场界面的多尺度耦合设计方法。
背景技术
多层结构和多个界面是微纳电子器件本身以及器件互连和封装中普遍存在的现象,界面分层失效成为产品性能和可靠性方面关心的重要问题。国外研究者通过大量实验发现界面是微系统制造和运行中的关键部分,很多破坏和缺陷都发生在界面附近。但微观材料界面规律的研究刚刚起步,以连续介质力学为基础的宏观理论已不再适用。微观材料界面物理特征不仅与微结构的几何与材料分布形态有关,还与结构的宏观边界条件,各种载荷等多种因素有关。
通过微观途径,可以建立起对材料行为的基本认识,它正逐渐成为发展新材料和高性能器件的不可或缺的重要手段。目前,分子动力学是微观建模的主要手段,许多在实验中无法获得的微观细节,在分子动力学模拟中都可以方便地观察到,分子动力学在原子尺度上比任何其他方法具有更高的时间和空间求解能力,因此,无法用连续介质分析方法求解的微观物理现象都可以用分子动力学进行有效地研究,然而分子动力学计算所需的庞大的时间是一个非常难以解决的问题,其计算时间随着原子数的增加而急剧增加,为了减少计算时间,分子动力学仿真所用的原子数量一般比较少,即便是超级计算机也只能模拟109个原子,也就是不到1平方微米的材料,这样的尺度范围对如裂纹扩展、能量冲击等的仿真显然是不够的。
发明内容
本发明的目的是为进一步促使微纳电子器件结构特别是封装设计向高效及低成本的方向发展而提供一种微纳电子器件结构的多物理场界面多尺度设计方法。
本发明采用的技术方案为采用如下步骤:先描述耦合机制,用分子动力学描述界面组织演化和材料缺陷特征;再用有限元为分子动力学模型提供热循环、弯曲、拉伸行为的边界条件;最后进行界面多物理场特征的数学模拟。
本发明是针对微纳电子器件结构多物理场界面特别是封装界面特征研究多尺度模型的构架、针对微观界面尺度材料组织的演变、热缺陷的产生机制、材料缺陷如孔洞、裂纹等对微结构物理特性的影响,形成了适用于微纳电子器件结构多物理场界面物理特性多尺度设计方法;该设计方法兼顾计算效率和科学准确性;实现微纳电子器件结构多物理场界面全性能特征的宏微纳观模拟,解决了微纳电子器件结构多物理场界面设计难题。
附图说明
图1为不同尺度的耦合机制示意图;
图2为多尺度算法流程图;
图3为实施例中按求解过程编程计算出在150000时间步时多尺度区域的温度分布图;
图4为不同热源温度下分子动力学模型和多尺度模型的热阻对照图。
具体实施方式
本发明用Green-Lindsay/Lord-Shulman(格林-林德赛/劳德-舒尔曼)方程描述耦合机制;用分子动力学描述界面组织演化和材料缺陷特征;包括不同形状、尺寸和数量的纳观缺陷和材料组织演化。用有限元(FE)为分子动力学(MD)模型提供热循环、弯曲、拉伸等行为的边界条件;进行界面多物理场特征的数学模拟。
在多尺度方法中,不同尺度的耦合是非常重要的问题。为了把离散的原子和连续的介质联系起来,应该把原子模型和连续模型用一种恰当的方式混合起来。如图1所示,采用一种耦合区域把两个尺度联系起来。在耦合区域中,原子和有限元网格重叠在一起。在耦合区域,通过有限元的形状函数可以得到原子温度。在多尺度耦合中,有限元区域需要一个独立的能量方程,从而使分子动力学区域产生的能量(比如热能)可以传播到有限元区域中去。采用映射算子技术,根据分子动力学模拟提供的信息,如原子速度和质量,可以获得连续的温度方程。
如图2所示,整个求解的过程可以描述为:
1)求解有限元方程[Kf]{uf}={Qf},f代表单元,K代表刚度,u代表位移,Q代表外加载。
2)通过有限元形状函数 将有限元节点温度转换为原子温度,a和e分别指原子区域和耦合区域的原子,[N]是有限元形状函数矩阵。
3)根据上一步得到的原子温度求解耦合区域原子的新速度和新位置{ra e}。
4)求解原子区域所有原子的新位置。
5)计算耦合区域原子的热流{Qa e}。
6)计算有限元区域的节点热流 其中Ve是包含嵌入原子的单元体积。
7)根据以上的计算结果,得到新的有限元解,然后继续重复第二步。
实施例
本实施例是对界面传热特性的模拟。首先用Green-Lindsay/Lord-Shulman(格林-林德赛/劳德-舒尔曼)方程描述传热耦合机制,再用分子动力学将界面组织结构划分成原子构造集成,然后用有限元为将热加载设置成边界条件,最后按下列步骤编程计算:
1)求解有限元方程[Kf]{uf}={Qf},f代表单元,K代表刚度,u代表位移,Q代表外加载,本实施例中的Q为热加载。
2)通过有限元形状函数 将有限元节点温度转换为原子温度,a和e分别指原子区域和耦合区域的原子,[N]是有限元形状函数矩阵。
3)根据上一步得到的原子温度求解耦合区域原子的新速度{va e}和新位置{ra e}。
4)求解原子区域所有原子的新位置。
5)计算耦合区域原子的热流{Qa e}。
6)计算有限元区域的节点热流 其中Ve是包含嵌入原子的单元体积。
7)根据以上的计算结果,得到新的有限元解,然后继续重复第2)步。
下面是根据上述实施例计算所得的不同特性:
1.界面传热的温度分布
多尺度模拟的一个关键问题是不同尺度的相互作用,如图3所示,可以看出原子区域描述的界面相互作用对有限元区域的影响,并体现出合理的温度梯度,说明了有限元区域和分子动力学区域的相互作用。
2.界面热阻结果与分子动力学模型的对比
为实现与分子动力学模型的量化对比,对不同热源温度下分子动力学模型和多尺度模型的热阻进行了计算,多尺度模型的热阻仍在分子动力学区域范围内进行计算。
如图4所示,多尺度模型计算得到的热阻更大,这与分子动力学区域和有限元区域的耦合机制有关。而在热源温度增加后,两种模型的计算结果有更接近的趋势,两种模型的结果较为接近,本发明提出的多尺度模型具有和分子动力学模型相当精确度,但计算时间更少。
Claims (3)
1、一种微纳电子器件结构的多物理场界面多尺度设计方法,其特征在于采用如下步骤:
1)描述耦合机制;
2)用分子动力学描述界面组织演化和材料缺陷特征;
3)用有限元为分子动力学模型提供热循环、弯曲、拉伸行为的边界条件;
4)进行界面多物理场特征的数学模拟。
2、根据权利要求1所述的微纳电子器件结构的多物理场界面多尺度设计方法,其特征在于:步骤2)中包括不同形状、尺寸和数量的纳观缺陷和材料组织演化。
3、根据权利要求1所述的微纳电子器件结构的多物理场界面多尺度设计方法,其特征在于描述耦合机制的步骤为:
(1)求解有限元方程;
(2)通过有限元形状函数将有限元节点温度转换为原子温度;
(3)根据得到的原子温度求解耦合区域原子的新速度和新位置;
(4)求解原子区域所有原子的新位置;
(5)计算耦合区域原子的热流;
(6)计算有限元区域的节点热流;
(7)根据计算结果,得到新的有限元解,然后继续重复第(2)步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100247179A CN101515307A (zh) | 2009-02-11 | 2009-02-11 | 微纳电子器件结构的多物理场界面多尺度设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100247179A CN101515307A (zh) | 2009-02-11 | 2009-02-11 | 微纳电子器件结构的多物理场界面多尺度设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101515307A true CN101515307A (zh) | 2009-08-26 |
Family
ID=41039759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2009100247179A Pending CN101515307A (zh) | 2009-02-11 | 2009-02-11 | 微纳电子器件结构的多物理场界面多尺度设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101515307A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033995A (zh) * | 2010-12-02 | 2011-04-27 | 中国人民解放军空军工程大学 | 一种耦合计算液压金属管路压力、流量和温度的方法 |
CN102156762A (zh) * | 2010-12-02 | 2011-08-17 | 中国人民解放军空军工程大学 | 一种耦合计算液压舵机压力、流量和温度的方法 |
CN102222142A (zh) * | 2011-06-20 | 2011-10-19 | 江苏大学 | 微/纳/光电子器件界面结构的跨尺度设计方法 |
CN103745028A (zh) * | 2013-12-02 | 2014-04-23 | 温州大学 | 一种微纳器件材料界面热疲劳性能模拟测试方法 |
CN104721152A (zh) * | 2015-01-20 | 2015-06-24 | 哈尔滨医科大学 | 一种羟基红花黄色素a冻干粉针剂的制备方法及其产品 |
CN110491451A (zh) * | 2019-07-09 | 2019-11-22 | 中山大学 | 一种构建拉伸-剪切耦合流场的耗散粒子动力学方法 |
WO2020058110A1 (en) | 2018-09-19 | 2020-03-26 | Signify Holding B.V. | Luminaire mounting |
-
2009
- 2009-02-11 CN CNA2009100247179A patent/CN101515307A/zh active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033995A (zh) * | 2010-12-02 | 2011-04-27 | 中国人民解放军空军工程大学 | 一种耦合计算液压金属管路压力、流量和温度的方法 |
CN102156762A (zh) * | 2010-12-02 | 2011-08-17 | 中国人民解放军空军工程大学 | 一种耦合计算液压舵机压力、流量和温度的方法 |
CN102033995B (zh) * | 2010-12-02 | 2012-07-25 | 中国人民解放军空军工程大学 | 一种耦合计算液压金属管路压力、流量和温度的方法 |
CN102222142A (zh) * | 2011-06-20 | 2011-10-19 | 江苏大学 | 微/纳/光电子器件界面结构的跨尺度设计方法 |
CN103745028A (zh) * | 2013-12-02 | 2014-04-23 | 温州大学 | 一种微纳器件材料界面热疲劳性能模拟测试方法 |
CN104721152A (zh) * | 2015-01-20 | 2015-06-24 | 哈尔滨医科大学 | 一种羟基红花黄色素a冻干粉针剂的制备方法及其产品 |
WO2020058110A1 (en) | 2018-09-19 | 2020-03-26 | Signify Holding B.V. | Luminaire mounting |
CN110491451A (zh) * | 2019-07-09 | 2019-11-22 | 中山大学 | 一种构建拉伸-剪切耦合流场的耗散粒子动力学方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101515307A (zh) | 微纳电子器件结构的多物理场界面多尺度设计方法 | |
Hsieh et al. | The mechanical response of cellular materials with spinodal topologies | |
Niiranen et al. | Isogeometric analysis for sixth-order boundary value problems of gradient-elastic Kirchhoff plates | |
Ghannadpour et al. | Bending, buckling and vibration problems of nonlocal Euler beams using Ritz method | |
Wu | Molecular dynamics study on mechanics of metal nanowire | |
Carrera et al. | One-dimensional finite element formulation with node-dependent kinematics | |
Setoodeh et al. | Application of TW-DQ method to nonlinear free vibration analysis of FG carbon nanotube-reinforced composite quadrilateral plates | |
Olsson et al. | On the importance of surface elastic contributions to the flexural rigidity of nanowires | |
Afrooz et al. | Effects of the carbon nanotube distribution on the macroscopic stiffness of composite materials | |
CN104182585B (zh) | 硅胶弹性体手感分析及疲劳寿命预测的有限元模拟方法 | |
Wang et al. | Effect of CNT volume fractions on nonlinear vibrations of PMMA/CNT composite plates: A multiscale simulation | |
Lv et al. | A novel metamaterial with individually adjustable and sign-switchable Poisson's ratio | |
Viet et al. | Mechanical attributes and wave propagation characteristics of TPMS lattice structures | |
Shariati et al. | A review on stress-driven nonlocal elasticity theory | |
Ramteke et al. | Multidirectional grading influence on static/dynamic deflection and stress responses of porous FG panel structure: a micromechanical approach | |
Chen et al. | A particle-based multiscale simulation procedure within the material point method framework | |
Mobarakian et al. | Developing a thermomechanical and thermochemical model for investigating the cooling rate effects on the distortion of unsymmetrical viscoelastic polymeric composite laminates | |
Sadeghi et al. | TPMS metamaterial structures based on shape memory polymers: Mechanical, thermal and thermomechanical assessment | |
Shen et al. | Multiscale simulation based on the meshless local Petrov-Galerkin (MLPG) method | |
Li et al. | Numerical cross-scale optimization of homogenized composite laminates under impact loading | |
Zhang et al. | Mechanics of stretchy elastomer lattices | |
Xu et al. | Compression behavior of 4D printed metamaterials with various Poisson's ratios | |
Zhang et al. | Structural design and tunable mechanical properties of novel corrugated 3D lattice metamaterials by geometric tailoring | |
Patni et al. | Efficient modelling of beam-like structures with general non-prismatic, curved geometry | |
Thirunavukkarasu et al. | Mechanically robust 3D printed elastomeric lattices inspired by strong and tough hierarchical structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20090826 |