Timon Rabczuk¹, Pedro Areias²

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 115-130, 2006, DOI:10.3970/cmes.2006.016.115

Abstract This paper proposes a meshfree method for arbitrary evolving cracks in thin shells. The approach is an improvement of the method proposed by Rabczuk T., Areias P.M.A., Belytschko T. (A meshfree thin shell for large deformation, finite strain and arbitrary evolving cracks, International Journal for Numerical Methods in Engineering). In the above cited paper, a shell was developed based on an intrinsic basis of third order completeness. Third order completeness was necessary to remove membrane locking. This resulted in the use of very large domains of influence that made the method computationally expensive. If the More >

K.H. Chen¹, J.T. Chen², J.H. Kao³

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

Abstract In this paper, we employ the regularized meshless method (RMM) to search for eigenfrequency of two-dimension acoustics with multiply-connected domain. The solution is represented by using the double layer potentials. The source points can be located on the physical boundary not alike method of fundamental solutions (MFS) after using the proposed technique to regularize the singularity and hypersingularity of the kernel functions. The troublesome singularity in the MFS methods is desingularized and the diagonal terms of influence matrices are determined by employing the subtracting and adding-back technique. Spurious eigenvalues are filtered out by using singular More >

S.-Y. Leu¹, J.T. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 129-140, 2006, DOI:10.3970/cmes.2006.014.129

Abstract Plastic limit angular velocity of rotating hollow cylinders made of the von Mises materials with nonlinear isotropic hardening is investigated numerically and analytically in the paper. The paper applies sequential limit analysis to deal with the rotating problems involving hardening material property and weakening behavior resulted from the widening deformation. By sequential limit analysis, the paper treats the plasticity problems as a sequence of limit analysis problems stated in the upper bound formulation. Rigorous upper bounds are acquired iteratively through a computational optimization procedure with the angular velocity factor as the objective function. Especially, rigorous More >

J. Ma², H. Lu², B. Wang², R. Hornung³, A. Wissink³, R. Komanduri^2,*

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 101-118, 2006, DOI:10.3970/cmes.2006.014.101

Abstract A new method for multiscale simulation bridging two scales, namely, the continuum scale using the generalized interpolation material point (GIMP) method and the atomistic scale using the molecular dynamics (MD), is presented and verified in 2D. The atomistic strain from the molecular dynamics simulation is determined through interpolation of the displacement field into an Eulerian background grid using the same generalized interpolation functions as that in the GIMP method. The atomistic strain is consistent with that determined from the virial theorem for interior points but provides more accurate values at the boundary of the MD… More >

Taeyoung Ha¹, Sangwon Seo², Dongwoo Sheen³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.1, pp. 57-76, 2006, DOI:10.3970/cmes.2006.014.057

Abstract We present nonoverlapping domain decomposition methods for the approximation of both electromagnetic fields in a three-dimensional bounded domain satisfying absorbing boundary conditions. A Seidel-type domain decomposition iterative method is introduced based on a hybridization of a nonconforming mixed finite element method. Convergence results for the numerical procedure are proved by introducing a suitable pseudo-energy. The spectral radius of the iterative procedure is estimated and a method for choosing an optimal matching parameter is given. A red-black Seidel-type method which is readily parallelizable is also introduced and analyzed. Numerical experiments confirm that the presented algorithms are faster than More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.2, pp. 83-108, 2006, DOI:10.3970/cmes.2006.012.083

Abstract The system we consider consists of two parts: a purely algebraic system describing the manifold of constraints and a differential part describing the dynamics on this manifold. For the constrained dynamical problem in its engineering application, it is utmost important to developing numerical methods that can preserve the constraints. We embed the nonlinear dynamical system with dimensions n and with k constraints into a mathematically equivalent n + k-dimensional nonlinear system, which including k integrating factors. Each subsystem of the k independent sets constitutes a Lie type system of X^˙_i = A_iX_i with A_i ∈ so(n_i,1) and n₁ +···+n_k = n.… More >

T. Hasebe¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 145-156, 2006, DOI:10.3970/cmes.2006.011.145

Abstract This paper presents recent achievements in field theoretical approach toward substantial linkage among key hieratical scales dominating polycrystalline plasticity of metals and alloys. Major ingredients of the theory are briefly shown first, which is followed by several overwhelming results and some implications including key factors for dislocation cell structure evolution, key features of polycrystalline plasticity and their rational modeling in crystal plasticity-based constitutive equation. 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 >

Jong Youn Park¹, Young-Sam Cho², Sung Youb Kim¹, Sukky Jun³, Seyoung Im¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 61-72, 2006, DOI:10.3970/cmes.2006.011.061

Abstract We present a coarse-graining computation for deformations of CNTs (carbon nanotubes) via QC (quasicontinuum), particularly targeting analysis of multi-walled carbon nanotubes. Higher order triangular elements are utilized for proper interpolation of atom positions of the CNT on the basis of QC approach. The computing scheme enables one to differentiate between the fully atomistic zone and the coarse-grained zone in the framework of the multiscale computing. Several numerical examples demonstrate the effectiveness and accuracy of the present methodology. 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 >