N.S. Mera, L. Elliott, D.B. Ingham, D. Lesnic¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 101-106, 2000, DOI:10.3970/cmes.2000.001.403

Abstract In this paper the iterative algorithm proposed by [Kozlov and Maz'ya (1990)] for the backward heat conduction problem is extended in order to solve the Cauchy steady state heat conduction problem and the accuracy, convergence and stability of the numerical algorithm are investigated. The numerical results which are obtained confirm that this new iterative BEM procedure is accurate, convergent and stable with respect to increasing the number of boundary elements and decreasing the amount of noise which is added into the input data. More >

Y.C. Shiah¹, C.L. Tan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 91-99, 2000, DOI:10.3970/cmes.2000.001.393

Abstract In the direct formulation of the boundary element method (BEM), a volume integral arises in the resulting integral equation if thermal effects are present. The steps to transform this volume integral into boundary ones in an exact analytical manner are reviewed in this paper for two- dimensional anisotropic thermoelasticity. The general applicability of the BEM algorithm for fracture mechanics applications is demonstrated by three crack problems with slanted cracks. The numerical results of the stress intensity factors are presented and compared with those obtained using superposition. More >

F. Simonetti¹, R. M. Ardito Marretta²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 81-90, 2000, DOI:10.3970/cmes.2000.001.383

Abstract Advanced propellers are being developed to improve the performance and fuel economy of future transport aircraft. To study them, various aerodynamic prediction models and systems (from theory to experiment) have been developed via several approaches (Free Wake Analysis, helicoidal source methods, scale model tests). This study focuses on the development of an efficient numerical method to predict the behaviour of rotor or propeller in forward flight. Based on a variational approach, the present numerical technique allows a significant reduction of computer resources used in the calculation of instantaneous velocities to determine the wake geometry and More >

B. Pichler¹, H.A. Mang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 43-55, 2000, DOI:10.3970/cmes.2000.001.345

Abstract In order to avoid a fully nonlinear analysis to obtain stability limits on nonlinear load-displacement paths, linear eigenvalue problems may be used to compute estimates of such limits. In this paper an asymptotic approach for assessment of the errors resulting from such estimates is presented. Based on the consistent linearization of the geometrically nonlinear static stability criterion – the so-called consistently linearized eigenvalue problem – higher-order estimation functions can be calculated. They are obtained from a scalar post-calculation performed after the solution of the eigenproblem. Different extensions of these higher-order estimation functions are presented. An ab More >

F. M. A. El-Saeidy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 33-42, 2000, DOI:10.3970/cmes.2000.001.335

Abstract In rotating radial ball bearings supported on elastic casings with the bearing outer ring lightly fitted into the housing, the force due to the ball elastic contact is indeed a rotating load rolling over the housing. For accurate estimation of the dynamic deformations of the casing annulus (hole), which in turn affect the bearing tolerances and hence the magnitudes of the generated forces, effect of the load rotation (motion) should be considered. Considering the integral casing and the outer ring to be a plate, an isoparametric plane stress finite-element (FE) based analytical procedure is presented… More >

J. Lufrano¹, P. Sofronis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 119-132, 2000, DOI:10.3970/cmes.2000.001.279

Abstract Transient hydrogen diffusion and hydride formation coupled with material deformation are studied in Zr-2.5Nb alloys used in the pressure tubes of CANDU nuclear generating stations. The energetics of the hydride formation is revisited and the terminal solid solubility of hydrogen in solution is defined on the basis of the total elastoplastic work done on the system by the forming hydride and the external loads. Probabilistic precipitation of hydride is modeled in the neighborhood of a crack tip under mode I plane strain loading and a uniform initial hydrogen concentration below the stress free terminal solid More >

Igor Patlashenko¹, Dan Givoli²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 61-70, 2000, DOI:10.3970/cmes.2000.001.221

Abstract The method of Absorbing Boundary Conditions (ABCs) is considered for the numerical solution of a class of nonlinear exterior wave scattering problems. Recently, a scheme based on the exact nonlocal Dirichlet-to-Neumann (DtN) ABC has been proposed for such problems. Although this method is very accurate, it is also highly expensive computationally. In this paper, the nonlocal ABC is replaced by a low-order local ABC, which is obtained by localizing the DtN condition in a certain "optimal'' way. The performance of the new local scheme is compared to that of the nonlocal scheme via numerical experiments More >

Selden B. Crary¹, Peter Cousseau², David Armstrong¹, David M. Woodcock³, Eva H. Mok¹, Olivier Dubochet⁴, Philippe Lerch⁴, Philippe Renaud²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.1, pp. 127-140, 2000, DOI:10.3970/cmes.2000.001.127

Abstract We present a new and unique software capability for finding statistical optimal designs of deterministic experiments on continuous cuboidal regions. The objective function for the design optimization is the minimization of the expected integrated mean squared error of prediction of the metamodel that will be found, subsequent to the running of the computer simulations, using the best linear unbiased predictor (BLUP). The assumed response-model function includes an unknown, stochastic term, Z. We prove that this criterion, which we name I_Z-optimality, is equivalent to I-optimality for non-deterministic experiments, in the limit of zero correlations among the Z's for More >

Wenjing Ye¹, Subrata Mukherjee²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.1, pp. 111-120, 2000, DOI:10.3970/cmes.2000.001.111

Abstract Polynomial driving-force comb drives are designed using numerical simulation. The electrode shapes are obtained using the indirect boundary element method. Variable gap comb drives that produce combinations of linear, quadratic, and cubic driving-force profiles are synthesized. This inverse problem is solved by an optimization procedure. Sensitivity analysis is carried out by the direct differentiation approach (DDA) in order to compute design sensitivity coefficients (DSCs) of force profiles with respect to parameters that define the shapes of the fingers of a comb drive. The DSCs are then used to drive iterative optimization procedures. Designs of variable More >

Toru Ikeda¹, Isao Arase, Yuya Ueno, Noriyuki Miyazaki

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.1, pp. 91-98, 2000, DOI:10.3970/cmes.2000.001.091

Abstract In electronic devices, the corners of joined dissimilar materials exist between plastic resin and a die pad or a chip. Failure of the plastic resin is often caused from these corners during the assembly process or the operation of products. The strength evaluation of the corner is important to protect the failure of plastic packages. To evaluate the singular stress field around a corner, we utilize the stress intensity factors of the asymptotic solution for a corner of joined dissimilar materials. We show that the accurate stress intensity factor can be analyzed by the displacement More >