The effective permeability of the underfill flow domain in flip-chip packaging

In order to reach the goals of high electrical performance and dense packaging within the limited space, the flip-chip technology becomes popular in electronics packaging. In the flip-chip assembly, difference between thermal expansion coefficients of the chip and substrate may cause thermal fatigue at solder joints. To avoid this thermal fatigue, epoxy encapsulant is filled into the gap between the substrate and chip by the capillary force. Because of the small space in the flow domain, the underfilling flow can be assumed as a flow in porous medium. Permeability is used to characterize the flow field of the space among the substrate, chip, and solder bumps. In this study, a numerical method is used to determine the effective permeability for the underfilling flow domain. Analysis of the three dimensional flow in a unit cell of the underfill flow domain is performed. The resulting average velocity and pressure gradient are used to calculate the apparent permeability. Comparison with the analytical approximation for the permeability in literature is also performed. The effective permeability calculated using the proposed numerical method gives reasonable prediction of the underfill flow as compared to the experimental result.     

Investigation of the variation of the rubber shearing force during extended stationary contact under constant compressive strain

The effect of extended stationary contact on the shearing force has been investigated in the case of rubbers subjected to constant compressive deformation at high and low temperatures. The dependence of the shearing force on the logarithm of the stationary contact time is linear up to the moment of accumulation of 100% relative residual deformation. It is shown that the time-dependent processes that determine the values of the shearing forces can be intensified by raising the temperature. The results of the investigation can be used in predicting the shearing force in the case of rubbers in extended contact with metal under conditions of constant compressive strain.Laboratory for Problems in Polymer Physics, Lenin Moscow State Pedagogical Institute; Leningrad Branch of the Scientific-Research Institute of the Rubber Industry. Translated from Mekhanika Polimerov, No. 4, pp. 629–633, July–August, 1971.     

Optimum design of compressive residual stress field caused by ultrasonic surface rolling with a mathematical model

Compressive residual stress fields induced by ultrasonic surface rolling process play a key role in determining the fatigue performance of machine parts. The present work is an analytical approach to conducting optimum design of this field to obtain an optimum fatigue resistance. Thus, a mathematical model was presented to predict residual stresses based on circular and elliptical Hertz contact areas. Moreover, to validate the proposed analytical approach, experimental verification was carried out on 18CrNiMo7-6 steel. Analytic solutions were derived from the mathematical model and optimum characteristic parameters were obtained by investigating the characteristic parameters in this field, such as surface residual stress, maximum residual stress and its depths. Results showed that increasing the total force, Hertz contact area and ratio of the radius of tool tip to that of target body could significantly enhance the peak of compressive residual stress and its depth.     

人体动脉瘤生成与破裂的力学分析

在大变形超弹性理论框架下研究了内压、轴向拉伸和扭转联合作用下人体动脉壁的力学响应,应用结构不稳定性理论对动脉瘤生成的可能性进行了解释,应用材料强度理论对动脉瘤破裂的可能性进行了分析．考虑动脉壁中残余应力和平滑肌主动作用的影响,用纤维加强各向异性不可压超弹性复合材料两层厚壁圆筒模型来模拟动脉壁的力学特性．给出了正常和几种非正常状态下动脉壁的变形曲线和应力分布．变形和稳定性分析结果表明该文模型可以模拟正常状态下动脉壁的均匀变形,还可以模拟在动脉壁中弹性蛋白纤维和胶原蛋白纤维强度降低的非正常状态下动脉瘤生成的可能性及动脉瘤的增长．应力和强度分析结果表明该文模型可以模拟当动脉瘤中的最大应力超过管壁的强度时动脉瘤破裂的可能性．     

Novel approach for modelling of nanomachining using a mesh-less method

In this paper, a mesh-less method, called Smoothed Particle Hydrodynamics (SPH) is used to simulate the nanomachining operation in order to assist with the understanding of the fundamental mechanisms of nano scale material deformation and the characteristics of the post machined surface. An elasto-plastic nano-machining analysis is used to form a nano-groove using a conical tool on a copper specimen. The SPH solutions are validated against nano scale machining experiments conducted using a nanoindenter. The simulated results showed that the normal force is greater than the cutting force in this nano scale machining operation, which is consistent with the experimental results. Both the ploughing and cutting mechanisms were observed in these machining conditions and increased with the increase of the depth of cut. Moreover, the results reveal that the larger negative rake angle reduced the ploughing mechanism and caused higher residual stress and strain along the machined surface. Therefore, the effect of machining parameters on the nano deformation mechanism and the quality of the machined surface can be rapidly assessed using SPH.     

液滴对弹性梯度薄基变形的影响

液滴润湿现象在细胞的变形和软器件的设计和制作中具有潜在的指导意义.该文在考虑三相接触线处线张力的情况下,研究了液滴引起的梯度薄基变形问题.首先,利用积分变换法求解了基底变形的本构方程,给出了变形的法向位移表达式.其次,讨论了基底弹性模量非梯度、指数型梯度和幂型梯度变化时基底的变形情况.最后,给出了液滴大小、弹性模量、线张力及梯度指数变化时位移的变化情况.数值结果表明弹性模量逐渐减小和梯度指数逐渐增大时,湿润脊逐渐变高,变形也越大;线张力和特征深度越小,位移的峰值越高,变形也越大;液滴半径较小时,湿润脊的对称性会变得更好.     

Mathematical modelling of chip thickness in micro-end- milling: A Fourier modelling

This study developed a model of undeformed chip thickness in micro-end-milling for the use in estimating cutting constants based on measured cutting forces. The proposed estimation method is based upon the invertibility of the average milling force model. In this paper, chip thickness in micro-end-milling was estimated by summing the thicknesses of the conventional chip component and the additional chip component. Thickness was then expressed in terms of Fourier series. The analyses showed that the fast convergence of Fourier series gives the Fourier chip thickness model sufficient accuracy when using only five terms of the truncated Fourier series for common micro-end-milling processes. The Fourier coefficients can be expressed in terms of the ratio of feed per tooth to cutter radius for different numbers of cutter teeth. The accuracy and conciseness of the chip thickness model enables the modelling of average cutting force in a closed form, which can be applied to identify the cutting constants. Cutting force experiments verify that the model prediction agrees very well with the experimental results.     

Plastic deformation of rough surfaces

Friction forces are only transferred within the the real area of contact A_real, which is usually smaller than the apparent area of contact A_o. The maximum friction stress τ_fric is therefore determined by the shear limit τ_max in the area of real contact and the fraction of the real area of contact (τ_fric = τ_max (A_real/A_o)). For rough surfaces the size of A_real is governed o by the plastic deformation of the surface roughness. We present a fully elasto-plastic halfspace contact formulation based on the work of Jacq et al. [1]. Linear elastic-plastic material behavior is modeled based on v. Mises plasticity with isotropic hardening. The algorithm gives the residual stress as well as the full plastic deformation field due to a frictionless normal contact. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Oldroyd-B流体在孔隙壁面残余油上的作用力分析

为了分析计算粘弹性流体驱替残余油的微尺度力,从水动力学角度探索非牛顿流体的流变特性,选取Oldroyd-B本构方程来模拟粘弹性流体,并结合连续性方程和运动方程得到了粘弹性流体在微孔道中的流动方程,利用边界条件计算得到流动的流场,结合应力张量理论,计算出粘弹性流体作用在残余油上的法向偏应力和水平应力差,计算结果表明:沿流动方向,粘弹性流体的弹性越大,法向偏应力越大;垂直于流动方向,法向偏应力近似对称分布;随着粘弹性流体的弹性变化,水平应力差的变化趋势发生了变化,威森伯格数We越大,残余油所受的水平应力差先逐渐增加,达到峰值后降低,这种趋势更有利于残余油的变形,为下一步分析残余油的变形,并从主体上分离奠定了基础.     

Moving contact problems involving a rigid punch and a functionally graded coating

Frictional contact mechanics analysis for a rigid moving punch of an arbitrary profile and a functionally graded coating/homogeneous substrate system is carried out. The rigid punch slides over the coating at a constant subsonic speed. Smooth variation of the shear modulus of the graded coating is defined by an exponential function and the variation of the Poisson's ratio is assumed negligible. Coulomb's friction law is adopted. Hence, tangential force is proportional to the normal applied force through the coefficient of friction. An analytical method is developed utilizing the singular integral equation approach. Governing partial differential equations are derived in accordance with the theory of elastodynamics. The mixed boundary value problem is reduced to a singular integral equation of the second kind, which is solved numerically by an expansion-collocation technique. Presented results illustrate the effects of punch speed, coefficient of friction, material inhomogeneity and coating thickness on contact stress distributions and stress intensity factors. Comparisons indicate that the difference between elastodynamic and elastostatic solutions tends to be quite larger especially at higher punch speeds. It is shown that use of the elastodynamic theory provides more realistic results in contact problems involving a moving punch.     

Pull-in instability and free vibration of electrically actuated poly-SiGe graded micro-beams with a curved ground electrode

This paper investigates the pull-in instability and free vibration of functionally graded poly-SiGe micro-beams under combined electrostatic force, intermolecular force and axial residual stress, with an emphasis on the effects of ground electrode shape, position-dependent material composition, and geometrically nonlinear deformation of the micro-beam. The differential quadrature (DQ) method is employed to solve the nonlinear differential governing equations to obtain the pull-in voltage and vibration frequencies of the clamped poly-SiGe micro-beams. The present analysis is validated through direct comparisons with published experimental results and excellent agreement has been achieved. A parametric study is conducted to investigate the effects of material composition, ground electrode shape, axial residual stress and geometrical nonlinearity on the pull-in voltage and frequency characteristics.     

Modeling of large deformation frictional contact in powder compaction processes

In this paper, the computational aspects of large deformation frictional contact are presented in powder forming processes. The influence of powder–tool friction on the mechanical properties of the final product is investigated in pressing metal powders. A general formulation of continuum model is developed for frictional contact and the computational algorithm is presented for analyzing the phenomena. It is particularly concerned with the numerical modeling of frictional contact between a rigid tool and a deformable material. The finite element approach adopted is characterized by the use of penalty approach in which a plasticity theory of friction is incorporated to simulate sliding resistance at the powder–tool interface. The constitutive relations for friction are derived from a Coulomb friction law. The frictional contact formulation is performed within the framework of large FE deformation in order to predict the non-uniform relative density distribution during large deformation of powder die pressing. A double-surface cap plasticity model is employed together with the nonlinear contact friction behavior in numerical simulation of powder material. Finally, the numerical schemes are examined for efficiency and accuracy in modeling of several powder compaction processes.     

Investigation of the relation between the depth of penetration and indentation diameter in polymer hardness measurements

It has been experimentally demonstrated that for polymeric materials the area of the indentation should not be calculated from the depth of penetration of the spherical indenter, since this gives hardness values that are too low. It is shown that when the surface of a polymeric specimen is indented by a rigid ball, the effect of the spherical stress tensor on the yield point is important; the average stress on the contact area at the beginning for forced high-elastic deformation is close to the compressive yield stress of the polymeric material and considerably higher than its tensile yield stress.Institute of Polymer Mechanics, Academy of Sciences of the Latvian SSR, Riga. Translated from Mekhanika Polimerov, No. 6, pp. 995–998, November–December, 1971.     

Simulation of micro-cutting considering finite deformation crystal plasticity

Micro-machining processes on metalic microstructures are influenced by the crystal structure, i. e. the grain orientation. Furthermore, the chip formation underlies large deformations. To perform finite element simulations of micro-cutting processes, a large deformation material model is necessary in order to model the hyperelastic and finite plastic material behaviour. In the case of cp-titanium material with hcp-crystal structure the anisotropic behaviour must be considered by an appropriate set of slip planes and slip directions. In the present work the impact of the grain orientation on the plastic deformation is demonstrated by means of finite element simulations of a finite deformation single slip crystal plasticity model. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Determination of a stochastic friction coefficient depending on a randomly rough surface

The goal of this contribution is to calculate the the friction coefficient for a scanned surface of a worn brake pad. The data shows that the asperities can be approximated by paraboloids which allows to calculate the contact force and area with the Hertz contact model if the deformation is elastic. The friction force is calculated with the Bowden-Tabor approach which suggests that the friction force is the force to shear apart contacting asperities. This is considered to be the dominant friction cause in dry contact. To generate many surfaces with similar peak statistics the spectral decomposition is used. The friction coefficient and it's stochastic properties is calculated for these surfaces. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling the Sedimentation of Red Blood Cells in Flow under Strong External Magnetic Body Force Using a Lattice Boltzmann Fictitious Domain Method

Experimental observations show that a strong magnetic field has a dramatic influence on the sedimentation of RBCs, which motivates us to model the sedimentation of red blood cell (RBC) under strong external magnetic body force. To model the sedimentation of a RBC in a square duct and a circular pipe, a recently developed technique derived from the lattice Boltzmann and the distributed Lagrange multiplier/fictitious domain methods (LBM-DLM/FD) is extended to employ the mesoscopic network model for simulations of the sedimentation of a RBC in flow. The flow is simulated by the LBM with a strong magnetic body force, while the network model is used for modeling RBC deformation. The fluid-RBC interactions are enforced by the Lagrange multiplier. The sedimentation of RBC in a square duct and a circular pipe is simulated, which demonstrates the developed method's capability to model the sedimentation of RBCs in various flows. Numerical results illustrate that the terminal settling velocity increases incrementally with the exerted body force. The deformation of RBC has a significant effect on the terminal settling velocity due to the change in the frontal area. The larger the exerted force, the smaller the frontal area and the larger the RBC deformation become. Additionally, the wall effect on the motion and deformation of RBC is also investigated.     

The axisymmetric contact problem for an elastic layer loaded with a deformable cover plate

The axisymmetric problem of the contact interaction of an elastic cover plate with an elastic layer, loaded at infinity with a uniform stretching force, directed parallel to the boundaries of the layer, is considered. The cover plate resists stretching but does not resist bending. The contact shearing stress under the cover plate, the displacement of the points of the cover plate and the deformation distortion coefficient of the elastic layer are determined.     

Nonlinear analysis of residual stresses in a rail manufacturing process by FEM

The aim of this paper is to study the residual stresses in an UIC-60 rail and their reduction by means of roller straightening. Both experimental and numerical investigations have been carried out in the past to reveal the formation of dominant longitudinal residual stresses. However, the agreement between both investigations was not particularly good. The finite element method (FEM) has also been used to simulate one, two and three-dimensional analyses of a rail during roller straightening processes. The present model considers the longitudinal movement of a rail through the straightening machine, contact conditions between rail and rollers and kinematic hardening so as to take into account the plastic behaviour of the rail material (steel). These results were compared with the experimental investigations and good agreement was observed. In this respect, this paper presents a novel, more realistic numerical simulation by FEM for the roller straightening process. Finally, an improvement of the straightening process in order to obtain smaller residual stress in the rail section is proposed.     

Numerical analysis of chip formation during machining for different value of failure strain

Grinding is a very complicated processing. To increase quality of product and minimize the cost of abrasive machining, we should know physical phenomena which exist during the process. The first step to solution of this problem is analysis of machining process with a single abrasive grain. In the papers [1, 2] the thermo-mechanical models of this process are presented, but in this work attention is concentrated on chip formation and his separation from object. The influence of failure strain ε_f on states of strain and stress in surface layer during machining is explained. The phenomena on a typical incremental step were described using step-by-step incremental procedure, with updated Lagrangian formulation. Then, the Finite Element Method (FEM) and Dynamic Explicit Method (DEM) were used to obtain the solution. Application was developed in the ANSYS system, which makes possible a complex time analysis of the physical phenomena: states of displacements, strains and stress. Numerical computations of the strain have been conducted with the use of two methodologies. The first one requires an introduction of boundary conditions for displacements in the contact area determined in modeling investigation, while the second – a proper definition of the contact zone through the introduction of finite elements of TARGET and CONTACT types, without the necessity to introduce boundary conditions. This model includes variational equations of the object's motion and deformation. Examples of calculations for the displacement, strain and stress field in the surface layer zones were presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)