Experimental and Numerical Analysis of Thermal Deformation of Electronic Packages,QFP and MCM

0Introduction Anelectronicpackageisgenerallyconstructedwithanactivesiliconchip,mountisland,gold wires,leadframesandsoldersasshowninFig.1(a).Toprotectfromtheenvironment,thesilicon chipisusuallyencapsulatedinresin.SincethesematerialshavedifferentCTE(coefficientofthermal expansion),stressandstrainaregeneratedbyinternalheatinginserviceoperationorthermalloading fromtheenvironment[1].Insomecases,thismaycausefailureofthewireandcrackinganddebonding insidethepackage.Thereliabilityofelectronicpackages…     

Parametric studies on relationships between flutter derivatives of slender bridge （Ⅰ）

Relationships between flutter derivatives of slender bridge are investigated based on our previously proposed semi-analytical flutter derivatives of flexible structure. The intrinsic relations are validated with test data of flutter derivatives of two bridges. Changes in flutter derivatives with the aerodynamic center, rotation speed, and angle variation are also studied by using a parametric method. The results show correctness of the proposed expressions of flutter derivatives given by authors in Ref. [1], and indicate that certain relations exist between these derivatives. It is also shown that semi-analytical flutter derivatives are applicable to bridges with a streamlined cross-section.     

Numerical simulation of solid tumor angiogenesis with Endostatin treatment： a combined analysis of inhibiting effect of anti-angiogenic factor and micro mechanical environment of extracellular matrix

To investigate the influence of anti-angiogenesis drug Endostatin on solid tumor angiogenesis, a mathematical model of tumor angiogenesis was developed with combined influences of local extra-cellular matrix mechanical environment, and the inhibiting effects of Angiostatin and Endostatin. Simulation results show that Angiostatin and Endostatin can effectively inhibit the process of tumor angiogenesis, and decrease the number of blood vessels in the tumor. The present model could be used as a valid theoretical method in the investigation of anti-angiogenic therapy of tumors.     

Carrier fringe method of moiré interferometry for tiny strain measurements in micro-field

In this paper, we demonstrate a new optical method for tiny strain measurements based on the principle of carrier fringes of moire interferometry. A cross-line grating with frequency of 1200 lp/mm is replicated on the specimen surface, and the strain can be deduced from the changes in carrier fringes before and after the deformation of an object. Four coherent laser beams are used to obtain the carrier fringe patterns of field U and V. Both theoretical analysis and numerical simulation indicate that the ideal accuracy of strain can be controlled within a range of ±1με. Case study of a plane extension experiment shows that the measurement accuracy of strain can be controlled within the range of ±10με. The average strain values of every row of field U and every column of field V can be obtained by using this method, and approximated strain of every pixel in the whole-field can be further acquired, and thus it is possible to measure tiny strains occurred in a micro-field. The technology in this paper can provide comprehensive information for analyzing related mechanical content in the field of MEMS.     

基于变量修正分布的混凝土结构长期变形预测Prediction of long-term deformation for concrete structure based on revised distribution of variables

由于随机因素的影响,混凝土结构长期变形通常表现出很强的离散性,Bayesian理论提供了改善这种离散性的方法。基于Bayesian理论的混凝土结构长期变形预测方法的核心是引入短期变形构造似然函数,通过修正先验概率得到长期变形的后验概率分布。但是对于实际结构而言,在施工之前短期变形及其标准差无法获取,这就使得这种方法在开展时机和实际应用方面存在一定的限制性。为了改善这种限制性,在随机变量修正分布的基础上,结合拉丁超立方抽样技术,采用CEB‐FIP（MC90）模型建立了钢筋混凝土梁长期变形的随机分析模型。采用该模型进行混凝土梁长期变形随机分析,得到基于变量修正分布的混凝土梁长期变形预测结果,并分别与先验预测结果和Bayesian预测结果进行比较。研究结果表明,基于变量修正分布的预测结果与Bayesian预测结果十分接近,比先验预测结果不确定性降低50％左右,与试验结果吻合良好。     

Global analysis of Ivlev＇s type predator-prey dynamic systems

Consider a class of Ivlev＇s type predator-prey dynamic systems with prey and predator both having linear density restricts. By using the qualitative methods of ODE, the global stability of positive equilibrium and existence and uniqueness of non-small amplitude stable limit cycle are obtained. Especially under certain conditions, it shows that existence and uniqueness of non-small amplitude stable limit cycle is equivalent to the local un-stability of positive equilibrium and the local stability of positive equilibrium implies its global stability. That is to say, the global dynamic of the system is entirely determined by the local stability of the positive equilibrium.     

《Stochastic Dynamics of Structures》评介

<正>现代科学理论、科学试验和计算技术的迅速发展极大地提高了人类认识自然、改造自然并与自然和谐相处的能力.在这一进程中,人们日益认识到:必须兼顾系统与环境中客观存在的各种随机要素,才能更准确地反映工程结构的动力响应和性态.这一客观需求,促进了结构随机动力学的深入研究与持续发展.     

The application of three-dimensional hydroelastic analysis of ship structures in Pekeris hydro-acoustic waveguide environment简

The hydroelastic analysis and sonoelastic analysis methods are incorporated with the Green＇s function of the Pekeris ocean hydro-acoustic waveguide model to produce a three-dimensional sonoelastic analysis method for ships in the ocean hydro-acoustic environment. The seabed condition is represented by a penetrable boundary of prescribed density and sound speed. This method is employed in this paper to predict the vibration and acoustic radiation of a 1 500 t Small Water Area Twin Hull （SWATH） ship in shallow sea acoustic environment. The wet resonant frequencies and radiation sound source levels are predicted and compared with the measured results of the ship in trial.     

A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis简

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.     

Experimental analysis of sand particles＇ lift-off and incident velocities in wind-blown sand flux

The probability distributions of sand particles' lift-off and incident velocities in a wind-blown sand flux play very important roles in the simulation of the wind-blown sand movement. In this paper, the vertical and the horizontal speeds of sand particles located at 1.0 mm above a sand-bed in a wind-blown sand flux are observed with the aid of Phase Doppler Anemometry (PDA) in a wind tunnel. Based on the experimental data, the probability distributions of not only the vertical lift-off speed but also the lift-off velocity as well as its horizontal component and the incident velocity as well as its vertical and horizontal components can be obtained by the equal distance histogram method. It is found, according to the results of the χ ²-test for these probability distributions, that the probability density functions (pdf's) of the sand particles' lift-off and incident velocities as well as their vertical components are described by the Gamma density function with different peak values and shapes and the downwind incident and lift-off horizontal speeds, respectively, can be described by the lognormal and the Gamma density functions. These pdf's depend on not only the sand particle diameter but also the wind speed.     

P-FEMOL ANALYSIS OF PLATE BENDING PROBLEMS

In this paper,the p-version of the finite element method of lines(FEMOL) for the analysis of the Mindlin-Reissner plate bending problems is presentedand a class of p-FEMOL elements with polynomial degrees as high as nine isdeveloped.Numerical examples given in this paper show tremendous performance ofthe present method;namely,rapid convergence rate,high accuracy for bothdisplacements and stress resultants,removal of shear-locking trouble,capability ofdealing with difficult problems such as the boundary layer behavior near a free edgeand stress concentration around a hole.     

Three-dimensional mixed-mode stress intensity factors for cracks in functionally graded materials using enriched finite elements

Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading and constraint conditions. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previous developments applied so far on isotropic homogeneous and isotropic interface cracks. The spatial variation of FGM material properties is taken into account at the level of element integration points. To validate the developed method, two- and three-dimensional mixed-mode fracture problems are selected from the literature for comparison. Two-dimensional cases include: inclined central crack in a large FGM medium under uniform tensile strain and stress loadings, a slanted crack in a finite-size FGM plate under exponentially varying tensile stress loading and an edge crack in a finite-size plate under shear traction load. The three-dimensional example models a deflected surface crack in a finite-size FGM plate under uniform tensile stress loading. Comparisons between current results and those from analytical and other numerical methods yield good agreement. Thus, it is concluded that the developed three-dimensional enriched finite elements are capable of accurately computing mixed-mode fracture parameters for cracks in FGMs.     

具有待定形函数的有限元线法

有限元线法是一种优良的半解析法，但其解函数存在解析方向和离散方向上的精度不相称的问题。本文提出了一种基于有限元线法的双向解析法——具有待定形函数的有限元线法。该方法设形函数为未知待定并借助能量原理对形函数进行解析性求解，使解函数在两个方向上的精度达到相称，从而大幅度提升了解答的整体质量和精度。文中，以二维Poisson方程问题为例，具体给出了待定形函数的有限元线法的算法和算例，用以表明本法的可行性和有效性。     

基于DIC技术的爆炸应力波过异质界面应变场演化规律实验研究

采用氯仿粘结聚碳酸酯（polycarbonate, PC）板和聚甲基丙烯酸甲酯（polymethylmethacrylate, PMMA）板模拟含异质界面模型;在PC介质中布置柱状炮孔并与界面呈一定角度,根据炮孔端部与界面相对位置,分别于柱状炮孔两个端部设置起爆点,起爆点远离界面端部时定义为孔口起爆,靠近界面端部时定义为孔底起爆;借助数字图像相关实验系统,研究爆炸应力波通过异质界面后PMMA介质应变场演化过程及炮孔底部区域拉、压应变变化规律。实验结果表明,异质界面改变了爆炸应力波过界面后的传播形态。孔口起爆时,异质界面受爆破荷载作用后易形成应力集中区,界面处产生开裂,横向拉伸波作用是造成异质界面开裂的主要原因。起爆方式对过界面后介质PMMA的横/纵向拉、压应变场作用贡献不同,主要体现在应变场强度、拉/压应变场位置分布2个方面。在炮孔底部区域,起爆方式对应变场时程特性的影响主要体现在作用时效长短和应变强度2个方面。孔口起爆时,横/纵向应变体现出短时效、高强度的变化特征。就应变强度而言,起爆方式对横向压应变的影响显著强于其对纵向拉应变的影响。对空间分布特性影响主要体现在衰减程度,起爆方式对纵向应变衰减程度影响较大。无论采用何种起爆方式,爆炸应变场在PC介质中衰减速度较快,进入PMMA介质后衰减速度显著降低。     

STUDY ON THE SURFACE CRACK GROWTH BEHAVIOR IN 14MnNbq BRIDGE STEEL

Three-dimensional crack closure correction methods are investigated in this paper.The fatigue crack growth tests of surface cracks in 14MnNbq steel for bridge plate subjected to tensile and bending loadings are systematically conducted.The experimentally measured fatigue crack growth rates of surface cracks are compared with those of through-thickness cracks in detail.It is found that the crack growth rates of surface cracks are lower than those of through-thickness cracks.In order to correct their differences in fatigue crack growth rates, a dimensionless crack closure correction model is proposed.Although this correction model is determined only by the experimental data of surface cracks under tensile loading with a constant ratio R=0.05, it can correlate the surface crack growth rates with reasonable accuracy under tensile and bending loadings with various stress ratios ranging from 0 to 0.5.Furthermore, predictions of fatigue life and crack aspect ratio for surface cracks are discussed, and the predicted results are also compared with those obtained from other prediction approaches.Comparison results show that the proposed crack closure correction model gives better prediction of fatigue life than other models.     

A new analytical-numerical method for calculating interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses

Based on the elementary solutions and new integral equations, a new analytical-numerical method is proposed to calculate the interacting stresses of multiple circular holes in an infinite elastic plate under both remote stresses and arbitrarily distributed stresses applied to the circular boundaries. The validity of this new analytical-numerical method is verified by the analytical solution of the bi-harmonic stress function method, the numerical solution of the finite element method, and the analytical-numerical solutions of the series expansion and Laurent series methods. Some numerical examples are presented to investigate the effects of the hole geometry parameters (radii and relative positions) and loading conditions (remote stresses and surface stresses) on the interacting tangential stresses and interacting stress concentration factors (SCFs). The results show that whether the interference effect is shielding (k <1) or amplifying (k> 1) depends on the relative orientation of holes (α) and remote stresses (σ_x^∞, σ_y^∞). When the maximum principal stress is aligned with the connecting line of two-hole centers and σ_y^∞ <0.5σ_x^∞, the plate containing two circular holes has greater stability than that containing one circular hole, and the smaller circular hole has greater stability than the bigger one. This new method not only has a simple formulation and high accuracy, but also has an advantage of wide applications over common analytical methods and analytical-numerical methods in calculating the interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses.     

COMPUTATION OF STRESS INTENSITY FACTORS BY THE SUB-REGION MIXED FINITE ELEMENT METHOD OF LINES

Based on the sub-region generalized variational principle,a sub-region mixed ver- sion of the newly-developed semi-analytical‘finite element method of lines’(FEMOL)is pro- posed in this paper for accurate and efficient computation of stress intensity factors(SIFs)of two-dimensional notches/cracks.The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used,with the sought SIFs being among the unknown coefficients.The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements.A mixed system of ordinary differential equations(ODEs) and al- gebraic equations is derived via the sub-region generalized variational principle.A singularity removal technique that eliminates the stress parameters from the mixed equation system even- tually yields a standard FEMOL ODE system,the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver.A number of numerical examples,including bi-material notches/cracks in anti-plane and plane elasticity,are given to show the generally excellent performance of the proposed method.     

Elastoplastic solution for an eccentric crack loaded by two pairs of point tensile forces

The strain energy density theory and the near crack line analysis method are applied to investigate an eccentric crack loaded by two pairs of tensile point forces in a finite plate. The minimum values of SED in the vicinity of the crack tip are determined, the initial growth orientation of crack are determined. Obtained is the elastic-plastic solution near the crack line of an eccentric crack loaded by two pairs of point tensile forces under large scale yielding condition. More specifically, the near field solution contains the unit normal vector of the elastic-plastic boundary and the elastic-plastic stress field. The length of the plastic zone along the crack line is found to vary with the external load and the bearing capacity of a finite plate with an eccentric crack loaded by two pairs of tensile point forces. Compared with small scale yielding condition, the normalized load obtained is higher than those under small scale yielding condition when the length of the plastic zone is the same.     

Symmetrical contact problem of a thick transversely isotropic plate

The symmetrical problem of the contact between a spherical indenter and a thick transversely isotropic plate is solved using the techniques of the Hankel transform. Solutions are written as the sums of the associated half-space solution and plate thickness effect terms. The normalized normal contact stress and the surface radial stress are obtained and calculated numerically for both composite materials and metallic substances. The example materials give both real and complex characteristic roots. A general method of calculation is described to determine the maximum tensile radial stress and the maximum compressive normal contact stress.The plate thickness effects on the contact stress and on the tensile radial stress are shown to be increasingly important with decreasing thickness. The effects that the material anisotropy has on the magnification of the contact stress and the maximum tensile radial stress are clearly revealed in the numerical results presented. The effects of material anisotropy for the composite are compared to those for the metallic substances.