An axisymmetrical quasi-Kirchhoff-type shell element for large plastic deformations

Summary An axisymmetrical shell element for large plastic strains is developed. The theory is based on the multiplicative decomposition of the material deformation gradient into an elastic and a plastic part. For quasi-Kirchhoff-type axisymmetric shells this leads to a product of the elastic and plastic stretches. By introduction of logarithmic strains the decomposition becomes additive. Plastic incompressibility is fulfilled in an exact manner.Support for this work was provided by the Deutsche Forschungsgemeinschaft (DFG) under contract Wr 19/7-1. The financial aid for the first and second author is gratefully acknowledged.     

芯片封装焊球连接疲劳寿命预测分析--能量法和有效应变法之比较

研究了目前常用的两种芯片封装焊球连接的寿命预测方法一能量法和有效应变法。对焊球连接的材料本构关系进行了对比分析，并利用FORTRAN语言编制了相应的材料本构模型子程序，将其与有限元仿真工具ANSYS耦合，实现了将焊球连接的材料本构模型用户子程序导入到ANSYS的材料库。在此基础上模拟了三维芯片叠层球栅阵列尺寸封装的焊球结构，在热循环条件下(-40～ 125℃的工作状态，并分别利用能量法和有效应变法对焊球连接的寿命预测进行了比较分析，最后对两种方法作出了评价。     

A Low Order Nonconforming Mixed Finite Element Method for Non-Stationary Incompressible Magnetohydrodynamics System

A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by $H^1(\Omega)$-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of $L^2(H^1)$-norm for the velocity field, $L^2(L^2)$-norm for the pressure and the broken $L^2(H^1)$-norm for the magnetic field are derived.     

有限元/边界元耦合方法在海洋生物目标强度特性研究的应用*

目标强度特性是海洋生物声学识别与资源量评估的重要依据,其中,基于近似几何体和声阻抗特性的理论模型法是研究海洋生物目标强度的重要手段。由于对几何形态近似处理以及数值求解方法的限制,传统理论模型对声波频率、入射方位以及目标声阻抗、形态尺寸等均有各自不同的适用范围,单一模型难以满足不同种类或同一种类但不同尺寸海洋生物的目标强度求解。本文尝试将逐渐见诸应用的有限元/边界元耦合方法用于海洋生物目标强度特性研究,分别以球形生物、纺锤形鱼类尾明角灯鱼(Ceratoscopelus warmingii)和细长形浮游动物南极大磷虾(Euphausia superba)为例进行仿真计算,并与相适应的经典理论模型进行对比分析。结果表明,对于球形生物,有限元/边界元耦合方法与解析模型的目标强度频响曲线完全吻合;对于纺锤形鱼类,有限元/边界元耦合方法可有效弥补基于模态级数解的形变圆柱体模型在中低频和两端入射时的准确性问题;对于细长形浮游动物,有限元/边界元耦合方法与畸变波玻恩近似模型高度吻合。综上,有限元/边界元耦合方法对多种海洋生物目标强度求解均具有较好的适用性,未来有待进一步结合实验测定进行验证。     

液固对流相变问题的焓式有限元法

本文基于焓的对流扩散方程,提出一种求解具有糊状区液固相变对流传热问题的有限元方法,避免了通常采用温度对流扩散方程进行数值求解所带来的许多困难,对问题的解决更加自然、简便。算例表明,本文提出的方法是有效且可行的。     

一种四面体网格局部自适应方法在有限元中的应用

通过提出一种针对三维四面体网格的区域整形技术，给出了一种简便实用的四面体网格局部自适应加密方法，并将其与流量修正有限元法结合，应用于三维高速流的计算     

平板网架结构分析的超级有限元法

本文给出一种二维矩形超级单元，用以分析平板网架结构，其做法是将大型平板网架结构离散成一系列矩形超级单元，考虑弯曲、剪切、挤压、拉压、翘曲等多种非经典变形效应，采用构件自由度向超级元自由度的转换把大型多构件问题的求解变为二维问题的求解，可大大减少未知量，又能保证精度，并可方便地与一般有限元软件连接。     

Experimental and FEM study of windshield subjected to high speed bird impact

Both experimental and finite element model (FEM) results are presented for the dynamic strength behavior of windshield subjected to bird impact. The experimental data taken from a series of high speed photographs are compared with the numerical results predicted by using FEM in which the windshield was modeled entirely with solid elements and the bird body was approximately simulated by an elastic-plastic material with failure element behavior. Effective plastic strain and element pressure were adopted as the failure criteria and once the pressure or the effective plastic strain of an element reached the critical value, the element would lose the tensile resistance capability completely. The deflection and stress distribution in the windshield were obtained. It is shown that the result from the finite element analysis agrees with those from the full-scale bird impact test. The project supported by the National Natural Science Foundation of China (10272011)     

用均匀化理论分析蜂窝结构的等效弹性参数

在线弹性范围内,根据均匀化理论,并结合有限元方法推导出适用于二维周期性结构的均匀化的有限元格式（Homo FEM）,计算出不同相对密度下的规则蜂窝结构的等效弹性模量Ee和泊松比νe.同时,利用蜂窝结构的代表胞元模型,用常规的有限元方法（FEM）模拟计算出相应的等效弹性参数.最后将两种数值计算结果与己有的理论公式进行了比较和分析讨论.结果表明:在考察的相对密度全场范围内（0～0.4）,HOmO FEM得到的蜂窝结构的 Ee和νe 与 FEM使用平面实体单元模拟计算得到的结果一致吻合,反映出 Homo FEM数值方法的客观准确性和可行性.而 Gibson公式和 W-K得到的等效弹性模量值 Ee只是在较小相对密度的情况下（小于0.15）与数值计算结果吻合.当结构的相对密度较大时,必须考虑胞棱附近区域由应力集中导致的复杂的应力和应变分布的影响.     

动态J积分及Rice J积分一致性比较与分析

首先回顾了动态J积分的发展进程,然后给出了两种常用动态J积分和Rice J积分的数学计算式,并对它们的一致性加以讨论,由此提出一种广义J积分的定义。然后通过算例,基于ANSYS软件,建立平面应变情形下纯I型中心裂纹模型,并利用编写的APDL程序计算了各个J积分,通过结果比较,进一步验证了在动态加载条件下,动态J积分比Rice J积分具有更好的路径守恒性,但仍可以看出它们的取值具有一定的一致性。最后,通过J积分间接法和裂尖位移外推法分别计算裂纹动态应力强度因子,所得结果吻合很好,验证了编写程序的正确性。