Numerical Study of Finite Element Method Based Solutions for Propagation of Wetting Fronts in Unsaturated Soil

The accurate prediction of the propagation of a wetting front in an unsaturated soil subjected to surficial infiltration is of practical importance to many geotechnical and geoenvironmental problems. The finite element method is the most common solution technique as the hydraulic soil properties are highly nonlinear. Two important issues are often found to create difficulties in such analyses. First, numerical oscillations are usually observed in the calculated pore pressures at the wetting front. Second, when a reasonable mesh size and time step are used, the elevation of the wetting front may be seriously overpredicted. This paper is focused on the second issue. The under-relaxation (UR) technique used in the iterative process within each time step is found to have a serious impact on rate of convergence with refinement in mesh size and time step. Two different techniques are typically used; the first evaluates the hydraulic conductivity using an average of heads calculated from the preceding time node and the most recent iteration of the current time node (UR1), and the second evaluates the hydraulic conductivity using the average of heads calculated from the two most recent iterations of the current time nodes (UR2). The study shows that UR1, which is adopted in programs such as SEEP/W, ensures that the solution converges rapidly to a stable solution within a time step, but may converge to the wrong wetting front at a given elapsed time unless a sufficiently refined mesh is used. UR2 converges much more slowly within a time step, but the error in the wetting front is smaller than that generated by UR1.     

Numerical Simulation of Waveguides of Arbitrary Cross-Section

Finite difference equations are derived for the simulation of dielectric waveguides using an H_z-E_z formulation defined on a nonuniform triangular grid. The resulting equations may be solved as a banded eigenproblem for waveguide structures of arbitrary shape composed of regions of piecewise constant isotropic dielectric, and all transverse fields then computed from the solutions. Benchmark comparisons are presented for problems with analytic solutions, as well as a sample calculation of the propagation loss of a hollow Bragg fiber.     

Finite element analysis of femoral components paper III – hip joints

Finite element analysis of the hip implant is conducted in this paper for representative femoral cross-section geometries and development of stress in the presence of bone cement is elucidated. Differences in cement stresses generated by varying implant cross-sections were compared with conventional features derived from representative implants. The analysis was performed under idealized implant assemblies by constraining the implant movement in the assembly. The cross-sections and implant geometries used are generic and intended to be representative of available geometries. This paper is a part of student research projects, prepared from a series of activities in biomedical implant research currently underway at Ohio Northern and Arkansas Tech Universities.     

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     

Nonreactive atom/molecule-surface scattering within the finite basis wave packet method

A quantum wave packet code for studying nonreactive scattering of closed-shell atoms or diatomic molecules from a rigid surface is described. The time evolution relies on the Chebychev propagator. Up to 5 collider degrees of freedom, 3 in translation and 2 in rotation, are treated in a pseudospectral way with the momentum or finite basis representation as the primary space. Potential matrix elements are efficiently evaluated by means of sequential 1D transformations between momentum and coordinate spaces. Fast Fourier transforms are performed for the translational and azimuthal coordinates whereas a Gauss-Legendre transform is used for the polar coordinate. This pseudospectral strategy minimizes memory requirements because no off-diagonal Hamiltonian matrix elements need to be stored. In addition, a wide variety of physical systems can be studied since no particular functional form is imposed for the interaction potential.     

射孔孔眼中散射超声脉冲波的FEM模拟

采用有限元素法 ,分别数值计算了平底孔道、台阶式孔道和尖劈状孔道等 3种孔眼模型中不同中心频率 (10 0 k Hz、5 0 0 k Hz、1MHz)的散射超声脉冲波。计算结果表明 ,无论选用哪一种孔底模型 ,均可以得到较为明显的对应于孔底的声波信号。但是 ,孔道越不规则 ,使用的超声信号的主频越高 ,则孔底信号越差。为了兼顾孔底信号质量和信号采集、传输的代价 ,在未来的检测仪器中应采用尽可能低的测量频率 (低于 5 0 0 k Hz)     

���������ز��׹��ƻ�������Ԫ����

文中建立了套管、水泥环和地层互作用的弹塑性的有限元力学模型，用MSC.Marc软件的高度非线性功能对本文的复杂力学系统进行了一系列的计算机仿真模拟研究，为复杂地层套管破坏机理研究提供了简便，可行的方法。本方法可以反推地层载荷，为复杂地层套管破坏的预防措施提供了理论依据。     

A continuum-based modeling of MEMS devices for estimating their resonant frequencies

A continuum-based modeling of coupled electrostatics-structure interactions is presented for the frequency computations of MEMS devices. The present general formulation of electrostatics accounting for free space is validated first by specializing it to one-dimensional uniform motion of conducting surfaces and comparing the resulting electrostatics to conventional lumped models. The general coupled electrostatics-structure interactions are then applied for the prediction of resonant frequencies of MEMS devices due to bias-voltage changes and temperature variations. Comparisons of predicted resonant frequencies obtained by the present coupled electrostatics-structure interaction models with experimental results available in the literature demonstrate that the proposed continuum-based interaction modeling yields high-confidence predictions of resonant frequencies of MEMS devices.     

Model Validation for Bridge-Road-Vehicle Dynamic Interaction System

The dynamic response of highway bridges subjected to moving truckloads has been observed to be dependent on (1) dynamic characteristics of the bridge; (2) truck configuration, speed, and lane position on the bridge; and (3) road surface roughness profile of the bridge and its approach. Historically, truckloads were measured to determine the load spectra for girder bridges. However, truckload measurements are either made for a short period of time [for example, weigh-in-motion (WIM) data] or are statistically biased (for example, weigh stations) and cost prohibitive. The objective of this paper is to present results of a 3D computer-based model for the simulation of multiple trucks on girder bridges. The model is based on the grillage approach and is applied to four steel girder bridges tested under normal truck traffic. Actual truckload data collected using a discrete bridge WIM system are used in the model. The data include axle loads, truck gross weight, axle configuration, and statistical data on multiple presence (side by side or following). The results are presented as a function of the static and dynamic stresses in each girder and compared with code provisions for dynamic load factor. The study provides an alternate method for the development of live-load models for bridge design and evaluation.