T. Nishioka¹, J. Furutuka¹, S. Tchouikov¹, T. Fujimoto¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.1, pp. 129-146, 2002, DOI:10.3970/cmes.2002.003.129

Abstract The governing condition of dynamic crack bifurcation phenomena had not been fully elucidated until our recent experimental studies. We found from the experimental results that the energy flux per unit time into a propagating crack tip or into a fracture process zone governs the crack bifurcation. Regarding the numerical simulation of dynamic crack bifurcation, to the authors' knowledge, no accurate simulations have been carried out, due to several unresolved difficulties. In order to overcome the difficulties, for the analysis of dynamic crack bifurcation, we developed a moving finite element method based on Delaunay automatic triangulation. More >

Chyou-Chi Chien, Tong-Yue Wu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 213-226, 2001, DOI:10.3970/cmes.2001.002.213

Abstract This study presents a novel computational method for implementing the time finite element formulation for the equations of linear structural dynamics. The proposed method adopts the time-discontinuous Galerkin method, in which both the displacement and velocity variables are represented independently by second-order interpolation functions in the time domain. The solution algorithm derived utilizes a predictor/multi-corrector technique that can effectively obtain the solutions for the resulting system of coupled equations. The numerical implementation of the time-discontinuous Galerkin finite element method is verified through several benchmark problems. Numerical results are compared with exact and accepted solutions from More >

S. Rugonyi, K. J. Bathe¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 195-212, 2001, DOI:10.3970/cmes.2001.002.195

Abstract The solution of fluid flows, modeled using the Navier-Stokes or Euler equations, fully coupled with structures/solids is considered. Simultaneous and partitioned solution procedures, used in the solution of the coupled equations, are briefly discussed, and advantages and disadvantages of their use are mentioned. In addition, a simplified stability analysis of the interface equations is presented, and unconditional stability for certain choices of time integration schemes is shown. Furthermore, the long-term dynamic stability of fluid-structure interaction systems is assessed by the use of Lyapunov characteristic exponents, which allow differentiating between a chaotic and a regular system 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 >

N. D. Ngo, K. K. Tamma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 57-72, 2000, DOI:10.3970/cmes.2000.001.359

Abstract In the process modeling via Resin Transfer Molding (RTM) for thick composite sections, multi-layer preforms with varying thermophysical characteristics across the different layers, or for geometrically complex mold geometries with varying thicknesses, the assumption of a thin shell-like geometry is no longer valid. The flow in the through thickness direction is no longer negligible and current practices of treating the continuously moving flow front as two-dimensional and the temperature and cure as three-dimensional are not representative of the underlying physics. In view of these considerations, in the present study, the focus is on the non-isothermal… More >

H. Okumura¹, M. Kawahara¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 71-78, 2000, DOI:10.3970/cmes.2000.001.231

Abstract This paper presents a new approach to a shape optimization problem of a body located in the unsteady incompressible viscous flow field based on an optimal control theory. The optimal state is defined by the reduction of drag and lift forces subjected to the body. The state equation used is the transient incompressible Navier--Stokes equations. The shape optimization problem can be formulated to find out geometrical coordinates of the body to minimize the performance function that is defined to evaluate forces subjected to the body. The fractional step method with the implicit temporal integration and More >

K. Kaminishi¹, M. Iino², H. Bessho², M. Taneda³

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

Abstract An FEA (finite element analysis) program employing a new scheme for crack growth analysis is developed and a prediction method for crack growth life is proposed. The FEA program consists of the subroutines for the automatic element re-generation using the Delaunay Triangulation technique, the element configuration in the near-tip region being provided by a super-element, elasto-inelastic stress analyses, prediction of crack extension path and calculation of fatigue life. The FEA results show that crack extension rate and path are controlled by a maximum opening stress range, Δσ_θmax, at a small radial distance of r = d, where More >

N. Miyazaki¹, T. Tamura², K. Yamamoto¹

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

Abstract Quantitative estimation of the failure of a gadolinium orthosilicate (Gd₂SiO₅, hereafter abbreviated as GSO) single crystal induced by thermal stress was investigated. A GSO cylindrical test specimen was heated in a silicone oil bath, then subjected to large thermal stress by room temperature silicone oil. Cracking occurred during cooling. The transient heat conduction analysis was performed to obtain temperature distribution in the test specimen at the time of cracking, using the surface temperatures measured in the test. Then the thermal stress was calculated using the temperature profile of the test specimen obtained from the heat conduction… More >

Philippe Matagne¹, Jean-Pierre Leburton², Jacques Destine, Guy Cantraine³

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

Abstract We investigate the quantum mechanical properties and single-electron charging effects in vertical semiconductor quantum dots by solving the Schrödinger and Poisson (SP) equations, self-consistently. We use the finite element method (FEM), specifically the Bubnov-Galerkin technique to discretize the SP equations. Owing to the cylindrical symmetry of the structure, the mesh is generated from hexahedral volume elements. The fine details of the electron spectrum and wavefunctions in the quantum dot are obtained as a function of macroscopic parameters such as the gate voltage, device geometry and doping level. The simulations provide comprehensive data for the analysis More >