A. P. S. Selvadurai¹, Wenjun Dong

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 41-54, 2006, DOI:10.3970/cmes.2006.012.041

Abstract A Fourier analysis conducted on both the spatial and the temporal discretizations of the governing partial differential equation shows that the Courant number as well as the time marching scheme have significant influences on the numerical behaviour of a Modified Least Squares (MLS) method for the solution of the advection equation. The variations of the amplification factor and the relative phase velocity with the Courant number and the dimensionless wave number indicate that when Courant number is equal to unity, the MLS method with the specified time-weighting and upwind function gives accurate results. This conclusion… More >

S.Yu. Reutskiy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 27-40, 2006, DOI:10.3970/cmes.2006.012.027

Abstract In this paper a new boundary method for eigenproblems with the Helmholtz equation in simply and multiply connected domains is presented. The solution of an eigenvalue problem is reduced to a sequence of inhomogeneous problems with the differential operator studied. The method shows a high precision in simply and multiply connected domains and does not generate spurious eigenvalues. The results of the numerical experiments justifying the method are presented. More >

A. Dimitrov¹, F.-G. Buchholz², E. Schnack³

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 1-26, 2006, DOI:10.3970/cmes.2006.012.001

Abstract Crack propagation in 3D-structures cannot be reduced (in general) to a series of plane problems along the crack front edge, due to the existence of some "corners'' on the crack front, where the elastic fields are of a real three-dimensional nature. The most important example for such a corner ist the point, where the crack front intersects a free surface of the body. According to the concept of weak and strong singularities, it is possible to obtain the asymptotics for the stress intensity factor (SIF) as well as the strain energy release rate (SERR) in… More >

S. Masuda¹, H. Noguchi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 131-144, 2006, DOI:10.3970/cmes.2006.011.131

Abstract This paper presents a novel and accurate technique for modeling discontinuous derivatives in meshfree methods, which will be used in the analysis of structures with material interfaces. The novelty lies in the formulation of the Moving Least Squares Approximation (MLSA) scheme where an introduced discontinuous derivative basis function replaces the conventional linear basis function. Furthermore, it is easy to implement this novelty into existing meshfree methods, such as the Element Free Galerkin (EFG) method, which are based on the MLSA scheme. The successful analyses of one and two-dimensional structures with material interfaces demonstrate the potential More >

T. Ishizawa¹, M. Iura¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 121-130, 2006, DOI:10.3970/cmes.2006.011.121

Abstract A one-dimensional model is proposed for numerical analysis of partially concrete-filled steel tubular (PCFST) columns subjected to cyclic loadings. The present formulation does not require experimental results nor shell analysis to obtain the constitutive equation of the model. The material properties and dimensions of PCFST columns are required for numerical analysis of the present model. The PCFST columns are assumed to consist of elastic beam and base plastic-hinge region in which steel local buckling is observed. Two parameters are introduced in order to express hardening phenomena of PCFST columns subjected to cyclic loading. Resisting forces More >

H. Kawai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 111-120, 2006, DOI:10.3970/cmes.2006.011.111

Abstract ADVENTURE AutoGL (pronounced as ‘Ote- ga-lu’) is a graphics and GUI library, dedicated for simulation-based research and development. It is designed for the simulation users to develop their own data viewers and editors. Currently, the library is used among many researchers and simulation users, mainly in universities and national research centers. Its functionalities and supported platforms are explained. AutoGL applications of various kinds of simulation methods are demonstrated also. More >

Y. Suetake¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 103-110, 2006, DOI:10.3970/cmes.2006.011.103

Abstract Since Reissner and Mindlin proposed their classical thick plate theories, many authors have presented refined theories including transverse shear deformation. Most of those plate theories have tended to use higher order power series for displacements and stresses along the thickness in order to achieve the higher accuracy. However, they have not carefully noticed lateral load effect. In this paper, we pay attention to constitution of the lateral loads: a body force and upper and lower surface tractions. Especially we formulate a modified theory for plate bending, in which the effect of a body force is More >

T. Fujimoto¹ and T. Nishioka¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 91-102, 2006, DOI:10.3970/cmes.2006.011.091

Abstract In the dynamic fracture of metallic material, some cracks propagate with the incidence of plastic deformation, and distinct plastic strain remains near the post-propagation area. In order to elucidate these dynamic nonlinear fracture processes, the moving finite element method is developed for nonlinear crack propagation. The T* integral is used as the parameter to estimate crack tip condition. First, the effect of material viscosity and crack propagation velocity have been discussed based on the numerical results for fracture under pure mode I high speed loading. Under mixed mode loading, numerical simulations for fracture path prediction More >

K. Yashiro¹, Y. Nakashima¹, Y. Tomita¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 73-80, 2006, DOI:10.3970/cmes.2006.011.073

Abstract A simple back force model is proposed for a dislocation cutting into γ' precipitate, taking the work formaking and recovering an anti-phase boundary (APB) into account. The first dislocation, or a leading partial of a superdislocation, is acted upon by a back force whose magnitude is equal to the APB energy. The second dislocation, or a trailing partial of a superdislocation, is attracted by the APB with a force of the same magnitude. The model is encoded in the 3D discrete dislocation dynamics (DDD) code and applied to the cutting behavior of dislocations at a… More >

Masanori Kikuchi ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 49-60, 2006, DOI:10.3970/cmes.2006.011.049

Abstract Three kinds of fracture specimens are tested under different constraint conditions. By the SEM(Scanning Electron Microscope) observation, it is shown that the roughness of fracture surface is different from each other largely. This is the effect of constraint condition. The dimple fracture process is simulated by the finite element method using Gurson’s constitutive equation, and the crack tip stress fields are obtained. The distributions of stress triaxiality qualitatively agree with the experimental results. The J-R curves obtained also qualitatively agree with those of experiments, and the fracture surface roughness is well simulated. More >