海水中的直流水平电偶极子在界面处产生的电磁场

舰船底部附近的电化学电流会产生恒定电磁场,研究这种电磁场分布对水下侦测舰船非常重要.文中将该场源理想化为一个直流水平电偶极子,这一直流水平电偶板子是低频电偶极子在频率趋于0时的极限情况.利用贝赛尔函数的性质,将这些积分表达式简化为一系列简单的解析表达式,并在该表达式的基础上进行分析讨论算出电磁场的空间分布.     

内压作用下发动机反向封头屈曲参数化研究

内压作用下发动机反向封头的稳定性是设计中的考虑重点.采用轴对称有限元方法对不同设计参数的等厚度和变厚度两大类反向封头进行稳定性研究,模型考虑封头开口加强环及封头与圆柱段连接处的局部加强影响.参数化研究椭球比、壁厚对封头临界屈曲载荷的影响.研究结果表明,等厚度反向封头的外形和变厚度的壁厚对结构临界屈曲载荷的影响十分敏感,在最小厚度一定的情况下,椭球比为1.413的变厚度封头是最小质量设计.另外文中给出可供工程实际设计参考的图表及曲线.     

非间断软件回归测试模型研究

软件测试是开发过程的重要环节。回归错误在开发过程中严重影响着开发的效率,其隐匿特性也使得测试比较困难。因此,本文提出一种非间断测试的模型,用于加强对回归错误的管理,节约不必要的无谓工作时间,提高软件开发的效率。这种模型对用户开发行为进行模拟和推测,在经过大量的数据训练后,寻找测试行为、错误发生等开发动作的捕获点,从而找到最佳的测试策略。在不影响开发人员工作的前提下,可对当前最新的代码进行测试,最大程度上避免回归错误长期隐藏在前期代码中,保持整个开发过程的稳步前进。     

仪器安装板振动环境虚拟试验和验证技术研究

基于商用有限元软件MSC.Patran和Nastran,研究了导弹仪器安装板结构的振动环境预示技术。建立了仪器安装板的三维实体模型,单机仪器使用质点模拟,并通过MPC法与安装板连接;设置了多种阻尼特性,通过模态叠加法预示了该仪器安装板在振动条件下的结构动力学环境,通过振动试验验证了预示结果。研究发现,使用大阻尼能较好地预示该仪器安装板结构的动力学环境。     

弹性抗滑桩内力计算的有限差分法

在传统抗滑桩内力计算的悬臂桩法的基础上,将滑动面以上的部分视为定向铰支的悬臂梁,以使滑动面上下位移符合连续性条件,滑动面以下采用地基系数"m"法计算桩身内力,并在此基础上推导了有限差分计算公式。用MATLAB编制了弹性桩全桩的内力计算程序,可用于滑坡推力和桩前剩余抗滑力为梯形、三角形和矩形的情况。     

Study on the fluid/structure interaction at different inlet pressures in molded packaging

This paper presents the computational study of fluid/structure interaction (FSI) analysis in the molding process using the Mesh-based parallel Code Coupling Interface (MpCCI) method with finite volume coding (FLUENT 6.3) and finite element coding (ABAQUS 6.9). The FSI analysis is implemented on the molded package during the encapsulation process with different inlet pressures. Real-time flow visualization, deformation and stress of the silicon die during the encapsulation process are presented in this paper. A fluctuation phenomenon of the silicon die is found in the encapsulation process when the inlet pressure increases. The maximum deformation during the process is determined at different locations on the silicon die, calculated during the final stage of the filling process. The deformation and stress of the die is exponentially increased with increasing inlet pressure. The maximum stress on the solder bump is concentrated near to the inlet gate. Thus, the present FSI analysis approach is expected to be a guideline or reference and provides better understanding of the encapsulation process for package design in the microelectronic industry.     

On the proper characterization of tooling motions and initial conditions in powder die compaction modeling

The present paper is concerned with the finite element modeling of Powder Metallurgy (P/M) cold die compaction process. Rather than on material constitutive theories or on numerical algorithmic issues, attention is confined exclusively on an scarcely addressed issue in the P/M modeling literature: the proper characterization of the boundary (tooling motions) and initial conditions of the problem. A case study of the compaction of an axially symmetric multilevel adapter in an advanced CNC press machine is used to convey the relevance of the accurate representation of these input data in the quality of model predictions. It is shown that unawareness or deliberate simplification of apparently insignificant details in this respect may cause errors far overshadowing those introduced by deficiencies in either the constitutive model or in the corresponding algorithmic solution procedure. The discussion of this case study serves also to provide useful modeling guidelines; illustrate frequent difficulties, as the unavailability of some information when guessing starting conditions; and reveal subtle, yet relevant for modeling purposes, technical details of advanced CNC press machines.     

The orthogonal meshless finite volume method for solving Euler–Bernoulli beam and thin plate problems

In this paper a new method untitled “orthogonal meshless finite volume method” (OMFVM) is developed for solving elastostatic problems in Euler–Bernoulli beam and thin plate. In this method, the weak formulation of a conservation law is discretized by restricting it to a discrete set of test functions. In contrast to the usual finite volume approach, the test functions are not taken as characteristic functions of the control volumes in a spatial grid, but are chosen from a Heaviside step function. The present approach eliminates the expensive process of directly differentiating the OMLS interpolations in the entire domain. This method was evaluated by applying the formulation to a variety of patch test and thin beam problems. The formulation successfully reproduced exact solutions. Numerical examples demonstrate the advantages of the present methods: (i) lower-order polynomial basis can be used in the OMLS interpolations; (ii) smaller support sizes can be used in the OMFVM approach; and (iii) higher accuracies and computational efficiencies are obtained.     

Bending of Delaminated Composite Conoidal Shells under Uniformly Distributed Load

Conoidal shells are very popular roofing units owing to their aesthetic elegance and stiffness. Many parts of the globe, which were earlier assumed to be seismologically stable, are now being considered as earthquake prone. Hence the necessity to build light structures using composites has become very important. In this paper an eight-noded isoparametric shell element is applied for analyzing the bending behavior of delaminated composite conoidal shells under a uniformly distributed load with different practical boundary conditions. To ensure compatibility of deformation and equilibrium of forces and moments at the delamination crack front, a multipoint constraint algorithm is incorporated, which leads to an unsymmetrical stiffness matrix. This formulation is validated through the solution of benchmark problems. Lamination, curvature, and extent of delamination area are varied to compare the performances of delaminated conoidal shells against those with no damage. The results are carefully observed, and a set of conclusions is presented at the end of the paper.     

A note on the divergence problem arising in calculation of Green’s function for a 2D piezoelectric thin film strip containing straight line defects and free boundaries: Fourier approach for bodies of unrestricted anisotropy

In this paper we derive a set of novel formulas for computation of the Green’s function and the coupled electro-elastic fields in a 2D piezoelectric strip with free boundaries and containing a distribution of straight line defects. The strip is assumed to be of unrestricted anisotropy, but allowing piezoelectricity, and in this sense situation is more general than in the available literature where only cubic symmetry was investigated. We employ a set of already known analytic formulas for the Fourier amplitude of the Green’s function and the corresponding electro-elastic fields. The key novelty of this paper is solution for the divergence problem occurring during integration of the Fourier amplitude. This problem is caused by poles at k = 0 in various matrix components of the amplitude. From purely mathematical point of view such poles lead to quantities which do not tend to zero at infinity, and this situation is clearly unphysical. To resolve this issue it is demonstrated by means of rigorous analysis that when some additional physical conditions are imposed, physical fields exhibit regular behavior at infinity - the poles do not contribute. Nevertheless, they lead to irremovable numerical ∞ − ∞ uncertainties spreading over the whole domain of integration. This motivates us to compute exact formulas for all these poles to enable engineering calculations involving the system in question.