A suitable computational strategy for the parametric analysis of problems with multiple contact

The aim of the present work is to develop an application of the LArge Time INcrement (LATIN) approach for the parametric analysis of static problems with multiple contacts. The methodology adopted was originally introduced to solve viscoplastic and large‐transformation problems. Here, the applications concern elastic, quasi‐static structural assemblies with local non‐linearities such as unilateral contact with friction. Our approach is based on a decomposition of the assembly into substructures and interfaces. The interfaces play the vital role of enabling the local non‐linearities, such as contact and friction, to be modelled easily and accurately. The problem on each substructure is solved by the finite element method and an iterative scheme based on the LATIN method is used for the global resolution. More specifically, the objective is to calculate a large number of design configurations. Each design configuration corresponds to a set of values of all the variable parameters (friction coefficients, prestress) which are introduced into the mechanical analysis. A full computation is needed for each set of parameters. Here we propose, as an alternative to carrying out these full computations, to use the capability of the LATIN method to re‐use the solution to a given problem (for one set of parameters) in order to solve similar problems (for the other sets of parameters). Copyright © 2003 John Wiley & Sons, Ltd.     

A three-dimensional FE analysis of large deformations for impact loadings using tetrahedral elements

 A three-dimensional dynamic program for the anaysis of large deformations in contact-penetration problems is developed using the finite element Lagrangian method with explicit time integration. By incorporating a tetrahedral element, which allows a single-point integration without a special hourglass control scheme, this program can be more effective to the present problem. The position code algorithm is used to search contact surface. Eroding surfaces are also considered. The defense node algorithm was slightly modified for the calculation of contact forces. A study of obliquity effects on metallic plate perforation and ricochet processes in thin plates impacted by a sphere was conducted. It is well simulated that on separation of two parts of the sphere, the portion still within the crater tends to perforate, while the portion in contact with the plate surface ricochets. This deformation pattern is observed in experiments, especially at high obliquities. A long rod that impacts an oblique steel plate at high impact velocity was also simulated in order to study the dynamics of the rod caused by the three dimensional asymmetric contact. The agreement between simulated and experimental results is quite good. Fracture phenomena occuring at high obliquity deserves further investigations. Received: 20 February 2002 / Accepted: 20 September 2002     

任意取向多个介质旋转椭球的散射场

本文提出了一种用于计算任意取向多个介质椭球电磁散射的混合方法.利用扩展边界条件(EBC)导出具有旋转轴对称单个椭球以球面波矢量函数为基底的T矩阵,结合球矢量波函数旋转定理和递推T矩阵方法,给出多个任意取向介质椭球的散射场.数值例子验证了该混合方法的计算精度和有效性.     

Ca掺杂钨镁酸铅陶瓷材料的制备及其介电性能

研究了Ca掺杂钨镁酸铅(PMW)陶瓷材料的合成、结构、烧结以及介电性能。结果发现：在Ca^2 摩尔分数小于15％时，能形成单相的PCMW钙铁矿相，结构由原来的斜方相向立方相转变。用二步合成法制备的样品容易致密烧结，气孔率比一步法制备的样品小。Ca的加入降低了材料的介电损耗，在频率为1MHz时，介质损耗达到了10^-4。当Ca^2 摩尔分数大于10％时，材料的Curie峰宽化显著，介电常数温度系数降低。     

Development of geometrically exact new shell elements based on general curvilinear co‐ordinates

In the present study first‐order shear deformable shell finite elements based on general curvilinear co‐ordinates are proposed. For the development of the present shell elements, a partial mixed variational functional with independently assumed strains is provided in order to avoid the severe locking troubles known as transverse shear and membrane lockings. Bubble functions are included in the shape function of displacement to improve the performance of the developed element. The proposed assumed strain four‐ and nine‐node elements based on the general tensor shell theory provide an efficient linkage framework for shell surface modelling and finite element analysis. In the several benchmark problems, the present shell elements with exact geometric representations demonstrate their performance compared to previously reported results. Copyright © 2002 John Wiley & Sons, Ltd.     

Range of the Applicability of the Wilkins Method to the Investigation of the Dynamic Stress-Strain State of Anisotropic Elastic Axisymmetric Shells

The exactness and stability of the Wilkins method in the investigation of the stress-strain state of axisymmetric anisotropic elastic shells, including thick-walled ones, have been studied for different values of artificial viscosity, pulse load rise and fall time, and different deformation process durations. A comparison of numerical calculations, performed by means of an application software package developed by V. A. Romashchenko on the basis of the Wilkins method, with experimental results, known numerical data, and with calculations using the computation kernel LS-DYNA 3D has been carried out.     

防止干燥塔塔顶及锥节产生焦粉的方法

描述及分析了干燥塔在生产乳粉过程中出现的塔顶及锥节焦粉现象，并提出了改进方法。     

平行线束法的仿射研究

用 3 个仿射图的投射方向构成 1 个仿射坐标系,同时也构成 1 个三棱锥,在此仿射坐标系中用球面三角重新研究由两个仿射图求第 3 个仿射图的投影作图理论 ——平行线束法。进一步还给出了由两个仿射图得到正投影的条件、由两个正投影得到第 3 个正投影的条件、以及由 V,H 两投影得到指定方向上的正投影的条件,简化了研究过程,得到了更简单的公式。     

一类连续模糊动态系统控制方法在Furuta摆中的应用

为了研究Furuta摆在倒立点附近的稳定控制问题，采用基于T-S模糊模型的一类连续模糊动态系统控制方法，首次进行Furuta摆的模糊状态空间建模，并采用模糊切换的思想，设计了Furuta摆在倒立点附近的稳定控制器，进而分析了系统稳定性。与已有的方法相比，该方法的优点是稳定区域大、调节时间短。数字仿真结果验证了方法的有效性。     

NUMERICAL INVESTIGATION ON MUD BED-GENERATION AND IMPACTING A WALL ALONG A SLOPE

In this paper mud was treated as the Bingham fluid. In staggered grids, two-dimensional incompressible Navier-Stokes equations for non-Newtonian fluid was solved by the MAC method. Numerical simulations were conducted on mud bed-generating phenomenon and mud impacting a wall along a slope. The distributions of free surface, pressure and velocity of mud and water were obtained. The results indicate that the computed layer thickness of mud bed almost equals the theoretically predicted value. Because of the differences in constitutive relationship, the distributions of free surface are different for water and mud. The distortion of water free surface is much more complicated.