Ru-Min Chao, Yen-Ji Chen, F.C. Lin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 73-86, 2001, DOI:10.3970/cmes.2001.002.073

Abstract The two-dimensional elasticity problem of an isotropic material, containing a centered-crack with unknown boundary traction is studied by the inverse procedure. The dual boundary integral equations are used to analyze the problem. While solving the ill-posed inverse problem, both of the conjugate gradient method and the regularization method are used. A scaling factor depending upon the material constant μ is introduced into the sensitivity matrix in order to keep the order of magnitude the same throughout the formulation. The result by using the displacement measurement will be compared with those by stress measurement, and an extensive More >

M. Kögl, L. Gaul²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.4, pp. 27-44, 2000, DOI:10.3970/cmes.2000.001.479

Abstract A Boundary Element formulation is presented for the solution of three-dimensional problems of anisotropic elastodynamics. Due to the complexity of the dynamic fundamental solutions for anisotropic materials and the resulting high computational costs, the approach at hand uses the fundamental solution of the static operator. This leads to a domain integral in the representation formula which contains the inertia term. The domain integral can be transformed to the boundary using the Dual Reciprocity Method. This results in a system of ordinary differential equations in time with time-independent matrices. Several general questions concerning the anisotropic solutions, More >

José A. González, Ramón Abascal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 141-150, 2000, DOI:10.3970/cmes.2000.001.443

Abstract In this work an approach to the two-dimensional steady-state rolling contact problem, with and without force transmission, is presented. The problem is solved by the combination of the Boundary Element Method with a formulation of the variational inequalities that govern the problem in the contact area, producing finally a mathematical programming problem. This formulation avoids the direct use of the contact constrains, but it drives to the minimisation of a non-differentiable function, being necessary the use of an specific numerical tool as the modified Newton's method. More >

C.P . Providakis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 123-130, 2000, DOI:10.3970/cmes.2000.001.425

Abstract Boundary and Finite Element methodologies for the determination of the inelastic response of thick plates resting on Winkler-type elastic foundations are compared and critically discussed. For comparison reasons the domain/boundary element and the finite element methodology use isoparametric elements of the same accuracy level. After a discretizaton of the integral equations of motion in both methodologies an efficient step-by-step time integration algorithm is used to solve the resulting matrix equations. Comparison studies are shown for impacted elastoplastic thick plates with smooth boundaries and supported on different Winkler-type foundations. The numerical results reveal that boundary element More >

R. Mustata¹, S. D. Harris², L. Elliott¹, D. Lesnic¹, D. B. Ingham¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 107-116, 2000, DOI:10.3970/cmes.2000.001.409

Abstract An inverse boundary element method is developed to characterise the components of the hydraulic conductivity tensor K of anisotropic materials. Surface measurements at exposed boundaries serve as additional input to a Genetic Algorithm (GA) using a modified least squares functional that minimises the difference between observed and BEM-predicted boundary pressure and/or hydraulic flux measurements under current hydraulic conductivity tensor component estimates. More >

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 >

J.A.F. Santiago¹, L.C. Wrobel²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 73-80, 2000, DOI:10.3970/cmes.2000.001.375

Abstract This work presents a boundary element formulation for two-dimensional acoustic wave propagation in shallow water. It is assumed that the velocity of sound in water is constant, the free surface is horizontal, and the seabed is irregular. The boundary conditions of the problem are that the sea bottom is rigid and the free surface pressure is atmospheric.
For regions of constant depth, fundamental solutions in the form of infinite series can be employed in order to avoid the discretisation of both the free surface and bottom boundaries. When the seabed topography is irregular, it is More >

A.J. Yiotis¹, J.T. Katsikadelis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 95-104, 2000, DOI:10.3970/cmes.2000.001.255

Abstract The Analog Equation Method is applied to the static and dynamic analysis of thin cylindrical shell panels. The Fl\"{u}gge theory is adopted. The three displacement components are established by solving two membrane and one plate bending problems under the same boundary conditions subjected to "appropriate'' (equivalent) fictitious loads. Numerical results are presented which illustrate the efficiency and the accuracy of the proposed method. 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 >