Homogenization of microsystem interconnects based on micropolar theory and discontinuous kinematics

In the present paper, the homogenized mechanical response of an interface in a microsystem interconnection is established on the basis of micropolar theory. The interface is treated as a finite RVE (representative volume element), across which macroscopic discontinuities occur as expressed in terms of the regularized discontinuous displacement and rotation fields. For the microstructure within the interfacial RVE, the micro-macro kinematical coupling is introduced as a second-order Taylor series expansion, along with a fluctuation term representing the microscopic displacement variation. In the second-order term of the expansion a restriction for the curvature is made, which motivates the adopted micropolar kinematics. Explicit expressions for the homogenized traction vector and the couple stress traction, associated with displacement and rotational jumps across the interface surface, are derived. A planar elastic interface is subjected to three basic deformation modes, i.e. the standard modes I, II and a non-conventional rotation mode, which are considered in the numerical examples representing a typical interconnect. A comparison to the results from the Taylor assumption is made, which shows that the Taylor assumption method produces an overstiffening of the interface.     

准脆性材料断裂过程区中的圆饼微裂纹损伤计算

准脆性工程材料及结构在外力作用下,不仅引起内部缺陷变化和微裂纹的出现及发展,且使得其结构承载能力降低或性能劣化．在其材料失效过程中常存在裂缝与断裂损伤过程区．为研究材料细观缺陷或微裂纹与宏观破坏的规律,通过细观力学方法,对于代表性体积单元RVE中的圆饼型微裂纹的尺寸与密度变化,探讨其宏观断裂过程区力学参量与损伤之间的量化关系．借助宏观断裂过程区的黏聚裂纹模型,将损伤单元RVE嵌入到宏观裂缝端部的断裂过程区中,对其进行联接细观损伤到宏观破坏的力学多尺度研究．文中也通过实验数据,对其理论计算结果进行了算例的讨论与分析．     

Homogenization-based continuum plasticity-damage model for ductile failure of materials containing heterogeneities

This paper develops an accurate and computationally efficient homogenization-based continuum plasticity-damage (HCPD) model for macroscopic analysis of ductile failure in porous ductile materials containing brittle inclusions. Example of these materials are cast alloys such as aluminum and metal matrix composites. The overall framework of the HCPD model follows the structure of the anisotropic Gurson-Tvergaard-Needleman (GTN) type elasto-plasticity model for porous ductile materials. The HCPD model is assumed to be orthotropic in an evolving material principal coordinate system throughout the deformation history. The GTN model parameters are calibrated from homogenization of evolving variables in representative volume elements (RVE) of the microstructure containing inclusions and voids. Micromechanical analyses for this purpose are conducted by the locally enriched Voronoi cell finite element model (LE-VCFEM) [Hu, C., Ghosh, S., 2008. Locally enhanced Voronoi cell finite element model (LE-VCFEM) for simulating evolving fracture in ductile microstructures containing inclusions. Int. J. Numer. Methods Eng. 76(12), 1955-1992]. The model also introduces a novel void nucleation criterion from micromechanical damage evolution due to combined inclusion and matrix cracking. The paper discusses methods for estimating RVE length scales in microstructures with non-uniform dispersions, as well as macroscopic characteristic length scales for non-local constitutive models. Comparison of results from the anisotropic HCPD model with homogenized micromechanics shows excellent agreement. The HCPD model has a huge efficiency advantage over micromechanics models. Hence, it is a very effective tool in predicting macroscopic damage in structures with direct reference to microstructural composition.     

不连续岩体变形破坏过程中AE的空间演化与意义

通过一些典型节理化岩体模型变形破坏过程中的AE事件的定位分析和破裂面形态的研究,探讨了不连续岩体变形破坏过程的空间演化特征。结果表明：单节理岩体变形破坏的空间动态特征受节理角控制,随着节理角从小变大,岩体的变形破坏方式从AE在局部集中的突发式失稳向AE弥漫于整个节理层面的类似渐进式的破坏形式发展；失稳破裂面从简单的整体状岩-膏界面变成形态复杂的层内复合破裂面；复杂节理岩体的变形破坏受一条（组）主导构造控制,强AE事件集中发生在不连续构造交汇部位,在包含多条主构造的岩体模型中最终失稳破裂面只与最后阶段AE的位置一致。     

角铺设复合材料层合板内部脱层破坏的研究

在面内剪切外载作用下，角铺设复合材料层板板最终的宏观破坏模式是脱层，然而从细观角度来看，宏观的脱层破坏可以对应不同的细观损伤过程，以破坏面的形貌为例，有些破坏面的形貌为例，有些破坏面主要由裸露的纤维和纤维迹组成，而有些破坏面则主要由矩齿形基体材料组成。不同的过程对应着不同的力学性能、诸如脱层强度、韧性等、本文从细观角度研究了脱层破坏过程，并就铺设角、界面强度、基体开裂强度对该过程的影响进行了讨论。     

Thermomechanical Multiscale Constitutive Modeling: Accounting for Microstructural Thermal Effects

In this work we present a thermomechanical multiscale constitutive model for materials with microstructure. In these materials thermal effects at microscale have an impact on the effective macroscopic stress. As a result, it turns out that the homogenized stress depends upon the macroscopic temperature and its gradient. In order to allow this interplay to be thermodynamically valid, we resort to a macroscopic extended thermodynamics whose elements are derived from the microscopic behavior using homogenization concepts. Hence, the thermodynamics implications of this new class of multiscale models are discussed. A variational approach based on the Hill–Mandel Principle of Macro-homogeneity, and which makes use of the volume averaging concept over a local representative volume element (RVE), is employed to derive the thermal and mechanical equilibrium problems at the RVE level and the corresponding homogenization expressions for the effective heat flux and stress. The material behavior at the RVE level is described through standard phenomenological constitutive models. To sum up, the novel contribution of the model presented here is that it allows to include the microscopic temperature fluctuation field, obtained from the multiscale thermal analysis, in the micro-mechanical problem at the RVE level while keeping thermodynamic consistency.     

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution is proposed. The basic assumption is that strain strength of geomaterial comply with a certain distribution law in space. Failure of tensile fracture plane and shear fracture plane in representative volume element (RVE) with iso-strain are discussed, and generalized failure constitutive model of fracture plane in RVE is established considering combined effect of tension and shear. Fracture plane consists of elastic microplanes and fractured microplanes. Elastic microplanes are intact parts of the fracture plane, and fractured microplanes are the rest parts of the fracture plane whose strain have ever exceeded their strain strength. Interaction mode on elastic microplanes maintains linear elasticity, while on fractured microplanes it turns into contact and complies with Coulomb’s friction law. Intact factor and fracture factor are defined to describe damage state of the fracture plane which can be easily expressed with cumulative integration of distribution density function of strain strength. Strong nonlinear macroscopic behavior such as yielding and strain softening can be naturally obtained through statistical microstructural damage of fracture plane due to distribution of strain strength. Elastic–brittle fracture model and ideal elastoplastic model are special cases of this model when upper and lower limit of distribution interval are equal.     

基于介观力学信息的颗粒材料损伤——愈合与塑性宏观表征

本文在二阶计算均匀化框架下提出了颗粒材料损伤--愈合与塑性的多尺度表征方法. 颗粒材料结构在宏观尺度模型化为梯度Cosserat连续体,在其有限元网格的每个积分点处定义具有离散颗粒介观结构的表征元. 建立了表征元离散颗粒系统的非线性增量本构关系. 表征元周边介质作用于表征元边界颗粒的增量力与增量力偶矩以表征元边界颗粒的增量线位移与增量转动角位移、当前变形状态下表征元离散介观结构弹性刚度、以及凝聚到表征元边界颗粒的增量耗散摩擦力表示. 基于平均场理论与Hill定理,导出了基于介观力学信息的梯度Cosserat连续体增量非线性本构关系. 在等温热动力学框架下定义了表征颗粒材料各向异性损伤--愈合和塑性的损伤、愈合张量因子与综合损伤、愈合效应的净损伤张量因子和塑性应变. 此外,定义了损伤和塑性耗散能密度与愈合能密度,以定量比较材料损伤、愈合、塑性对材料失效的效应. 应变局部化数值例题结果显示了所建议的颗粒材料损伤--愈合--塑性表征方法的有效性.      

核电690合金传热管微动疲劳试验及损伤机理分析

核电蒸汽发生器传热管在微幅磨损与交变载荷的作用下形成微动疲劳,导致其表面裂纹萌生和扩展乃至破裂,从而影响反应堆的安全. 为研究径向载荷以及轴向交变应力对690合金管微动疲劳寿命的影响规律,开展690合金管管材的微动疲劳试验,获得690合金管管材的微动疲劳寿命曲线,并与相关研究数据进行对比分析,以便探讨材料在微动疲劳下的寿命模型. 对不同载荷下的690合金管试样的磨痕表面进行三维形貌和扫描电镜观测,分析磨损表面的损伤机理;对不同载荷下的690合金管试样断口的宏观与微观形貌进行表征,分析裂纹萌生、起裂过程及其失效机理. 结果表明690合金管与403不锈钢(SS)抗振条间的磨损机理为剥层及磨粒磨损;690合金管在径向载荷作用下于微动磨损处产生裂纹源,裂纹在轴向交变应力的作用下不断向内部扩展,最终导致断裂;其断裂形式为解理疲劳断裂.      

Determination of the size of the representative volume element for random quasi-brittle composites

A representative volume element (RVE) is related to the domain size of a microstructure providing a “good” statistical representation of typical material properties. The size of an RVE for the class of quasi-brittle random heterogeneous materials under dynamic loading is one of the major questions to be answered in this paper. A new statistical strategy is thus proposed for the RVE size determination. The microstructure illustrating the methodology of the RVE size determination is a metal matrix composite with randomly distributed aligned brittle inclusions: the hydrided Zircaloy constituting nuclear claddings. For a given volume fraction of inclusions, the periodic RVE size is found in the case of overall elastic properties and of overall fracture energy. In the latter case, the term “representative” is discussed since the fracture tends to localize. A correlation factor between the “elastic” RVE and the “fracture” RVE is discussed.     

Configurational forces and couples in fracture mechanics accounting for microstructures and dissipation

Configurational forces and couples acting on a dynamically evolving fracture process region as well as their balance are studied with special emphasis to microstructure and dissipation. To be able to investigate fracture process regions preceding cracks of mode I, II and III we choose as underlying continuum model the polar and micropolar, respectively, continuum with dislocation motion on the microlevel. As point of departure balance of macroforces, balance of couples and balance of microforces acting on dislocations are postulated. Taking into account results of the second law of thermodynamics the stress power principle for dissipative processes is derived.Applying this principle to a fracture process region evolving dynamically in the reference configuration with variable rotational and crystallographic structure, the configurational forces and couples are derived generalizing the well-known Eshelby tensor. It is shown that the balance law of configurational forces and couples reflects the structure of the postulated balance laws on macro- and microlevel: the balance law of configurational forces and configurational couples are coupled by field variable, while the balance laws of configurational macro- and microforces are coupled only by the form of the free energy. They can be decoupled by corresponding constitutive assumption.Finally, it is shown that the second law of thermodynamics leads to the result that the generalized Eshelby tensor for micropolar continua with dislocation motion consists of a non-dissipative part, derivable from free and kinetic energy, and a dissipative part, derivable from a dissipation pseudo-potential.     

沿指定应力路径加载的有限元RVE模型研究

在细观尺度上建立能反映材料微观组织结构又能反映统计意义上宏观力学性能的代表性体积单元（Representative Volume Element, RVE）,对其进行复杂加载下的数值研究,是目前预测材料宏观力学性能较有效的方法。本文从理论上分析并提供了对正六面体RVE在任意应力状态及任意应力路径下加载及宏观应力、应变计算的方法,用有限元软件ABAQUS实现了数值计算过程,并用此方法对循环加载下缺口圆棒颈部中心和边缘位置进行了RVE分析。结果表明：（1）此方法能准确的控制并实现正六面体RVE在任意应力状态及应力状态路径下加载;（2）通过RVE分析,可用于复杂加载下试样局部细观结构变化的研究。     

近场动力学方法及其应用

近场动力学(peridynamics,PD)是一种新兴的基于非局部作用思想建立模型并通过求解空间积分方程描述物质力学行为的方法.它兼有分子动力学方法和无网格方法的优点,避免了基于连续性假设建模和求解空间微分方程的传统宏观方法在面临不连续问题时的奇异性,又突破了经典分子动力学方法在计算尺度上的局限,在宏/微观不连续力学问题分析中均表现出很高的求解精度和效率.首先概述了PD方法的理论基础、建模思路和计算体系;进而介绍了PD方法在不同尺度不连续力学问题中的应用,包括均匀与非均匀材料和结构的大变形、损伤、断裂、冲击、穿透和失稳问题,结晶相变动力学问题以及纳米材料和结构的破坏问题;最后讨论了PD方法在理论、计算和应用等方面值得进一步研究的问题.     

Influences of initial porosity,stress triaxiality and Lode parameter on plastic deformation and ductile fracture

Local mechanical properties in aluminum cast components are inhomogeneous as a consequence of spatial distribution of microstructure,e.g.,porosity,inclusions,grain size and arm spacing of secondary dendrites.In this work,the effect of porosity is investigated.Cast components contain voids with different sizes,forms,orientations and distributions.This is approximated by a porosity distribution in the following.The aim of this paper is to investigate the influence of initial porosity,stress triaxiality and Lode parameter on plastic deformation and ductile fracture.A micromechanical model with a spherical void located at the center of the matrix material,called the representative volume element(RVE),is developed.Fully periodic boundary conditions are applied to the RVE and the values of stress triaxiality and Lode parameter are kept constant during the entire course of loading.For this purpose,a multi-point constraint(MPC)user subroutine is developed to prescribe the loading.The results of the RVE model are used to establish the constitutive equations and to further investigate the influences of initial porosity,stress triaxiality and Lode parameter on elastic constant,plastic deformation and ductile fracture of an aluminum die casting alloy.     

考虑晶体滑移面分解正应力的细观损伤模型

材料内部的解理、滑移面剥离等细观损伤是引起宏观失效的根源,从细观尺度研究损伤的发生和发展有助于深入认识材料的变形和失效过程.本文基于晶体塑性理论,从滑移系的受力和变形出发研究材料的细观损伤,建立了考虑滑移面分解正应力的细观损伤模型,为晶体材料解理断裂的分析提供了新方法.首先,在晶体弹塑性变形构型的基础上引入损伤变形梯度...     

基于均匀化理论韧性复合材料塑性极限分析

运用细观力学中的均匀化方法分析了韧性复合材料的塑性极限承载能力.从反映复合材料细观结构的代表性胞元入手,将均匀化理论运用到塑性极限分析中,计算由理想刚塑性、Mises组分材料构成的复合材料的极限承载能力.运用机动极限方法和有限元技术,最终将上述问题归结为求解一组带等式约束的非线性数学规划问题,并采用一种无搜索直接迭代算法求解.为复合材料的强度分析提供了一个有效手段.     

Identification of damage mechanism and validation of a fracture model based on mesoscale approach in spalling of titanium alloy

The subject of this paper is identification of the physical mechanisms of spalling at low impact velocities for Ti–6Al–4V alloy and determination of the macroscopic stress of spalling via meso-macro approach. Spalling is a specific mode of fracture which depends on the loading history. The aspects of the initial microstructure and its evolution during plastic deformation are very important. In order to identify the spalling physical mechanisms in titanium alloy, numerous pictures by the optical microscopy of the spall surfaces created by plate impact technique have been taken. The scenario of failure observed is in complete agreement with known physical micro-mechanisms: namely nucleation, propagation and coalescence by adiabatic shearing of micro-voids. The most interesting point in spall fracture of Ti–6Al–4V alloy is the nucleation of micro-voids. A significant amount of small micro-voids in the region of the expected spall plane has been observed. It appears that microstructural effects are important due to dual α–β phase microstructure, called Widmanstätten structure. The orientation of microstructure has a direct influence on nucleation mechanism by means of distribution of nucleation sites and decohesion between the softer particles (α-phase lamellae) and the harder lattice (β-phase). According to these observations, a fracture model has been developed. This model is based on the numerous post-mortem microscopic observations of spall specimens. The goal is to determine the macroscopic stress of spalling in function of loading time and damage level via a meso-macro approach.     

On the relationship between impact energy and delamination area

Thin glass/epoxy plates fabricated from 3M prepreg tape were subjected to low-velocity impact. Delamination in the impacted composite plates was measured by edge replication and was identified as the major damage mode. The effects of fiber orientation, thickness, and lamination on the delamination resistance were investigated. Experimental results verified three previous findings: (1) a linear relationship holds between delamination area and impact energy, (2) the mismatch of bending stiffness between adjacent laminae can be correlated with the delamination area on the interface, and (3) the behavior of a thin composite plate under low-velocity impact is very similar to that caused by global bending. In addition, based on the calculation of impact energy per unit delamination area, the dynamic fracture energy and the energy of dissipation in the thin glass/epoxy plates can be examined.