S.K. Pal¹, D.A. Linkens²

CMC-Computers, Materials & Continua, Vol.2, No.2, pp. 113-118, 2005, DOI:10.3970/cmc.2005.002.113

Abstract In hot rolling, the final thickness of the strip is achieved through plastic deformation of the original stock by a series of counter-rotating rollers. In this work, static recrystallization kinetics in between two stages of steel rolling has been modelled, and simulation studies have also been performed to find out the effect of entry temperature on the recrystallization kinetics. A viable bond graph model has been developed to study the kinetics of the process. Low-carbon steel has been considered for this purpose. More >

B. Amaziane¹, A. Naji², D. Ouazar³, A. H.-D. Cheng⁴

CMC-Computers, Materials & Continua, Vol.2, No.2, pp. 85-96, 2005, DOI:10.3970/cmc.2005.002.085

Abstract In this paper the optimization of groundwater pumping in coastal aquifers under the threat of saltwater intrusion is investigated. The aquifer is inhomogeneous and contains several hydraulic conductivities zones. The aquifer data such as the hydraulic conductivities are uncertain, but with their expected mean and standard deviation values given. A stochastic boundary element method based on the perturbation technique is employed as the simulation tool. The stochastic optimization is handled by the chance-constrained programming. Genetic algorithm is selected as the optimization tool. Numerical examples of deterministic and stochastic problems are provided to demonstrate the feasibility More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 193-210, 2005, DOI:10.3970/cmes.2005.009.193

Abstract Recently, the MLPG (Meshless Local Petrov-Galerkin Method) method has been successfully extended to simulating nonlinear water waves [Ma, (2005)]. In that paper, the author employed the Heaviside step function as the test function to formulate the weak form over local sub-domains, acquiring an expression in terms of pressure gradient. In this paper, the solution for Rankine sources is taken as the test function and the local weak form is expressed in term of pressure rather than pressure gradient. Apart from not including pressure gradient, velocity gradient is also eliminated from the weak form. In addition, More >

Chein-Shan Liu, Yu-Ling Ku

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 151-178, 2005, DOI:10.3970/cmes.2005.009.151

Abstract In this paper we are concerned with the integration of a semi-discretized version of the Landau-Lifshitz equation, which is fundamental to describe the magnetization dynamics in micro/nano-scale magnetic systems. The resulting ordinary differential equations at the interior grid points are numerically integrated by a combination of the group preserving scheme derived by Liu (2004a) and the fourth-order Runge-Kutta method, abbreviated as GPS-RK4. The new method not only conserves the magnetization magnitude and has the fourth-order accuracy, but also preserves the Lyapunov property of the Landau-Lifshitz equation, namely the free energy is decreasing with time. In More >

U. Andreaus^1,3, R.C. Batra², M. Porfiri^{2, 3}

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 111-132, 2005, DOI:10.3970/cmes.2005.009.111

Abstract Structural health monitoring techniques based on vibration data have received increasing attention in recent years. Since the measured modal characteristics and the transient motion of a beam exhibit low sensitivity to damage, numerical techniques for accurately computing vibration characteristics are needed. Here we use a Meshless Local Petrov-Galerkin (MLPG) method to analyze vibrations of a beam with multiple cracks. The trial and the test functions are constructed using the Generalized Moving Least Squares (GMLS) approximation. The smoothness of the GMLS basis functions requires the use of special techniques to account for the slope discontinuities at More >

Wenzhong Tang¹, Suresh G. Advani¹

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 165-176, 2005, DOI:10.3970/cmes.2005.008.165

Abstract An optimization method is proposed to simulate the motion of long flexible fibers in shear flow. The fiber is modeled as spheres connected by massless rigid rods and ball-socket joints. The optimization method is mathematically justified and used to obtain the position of a fiber at the next time step from its current position. Results for a single fiber in simple shear flow agree well with those reported in the literature. The usefulness of the method is demonstrated by simulating the motion of two interactive fibers subjected to shear flow field, and by studying the More >

D. Soares Jr.¹, W.J. Mansur^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 153-164, 2005, DOI:10.3970/cmes.2005.008.153

Abstract A coupling procedure is described to perform time-domain numerical analyses of dynamic fluid-structure interaction. The fluid sub-domains, where acoustic waves propagate, are modeled by the Boundary Element Method (BEM), which is quite suitable to deal with linear homogeneous unbounded domain problems. The Finite Element Method (FEM), on the other hand, models the structure sub-domains, adopting a time marching scheme based on implicit Green's functions. The BEM/FEM coupling algorithm here developed is very efficient, eliminating the drawbacks of standard and iterative coupling procedures. Stability and accuracy features are improved by the adoption of different time steps More >

J. Ma¹, H. Lu¹, B. Wang¹, S. Roy¹, R. Hornung², A. Wissink², R. Komanduri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 135-152, 2005, DOI:10.3970/cmes.2005.008.135

Abstract In the simulation of a wide range of mechanics problems including impact/contact/penetration and fracture, the material point method (MPM), Sulsky, Zhou and Shreyer (1995), demonstrated its computational capabilities. To resolve alternating stress sign and instability problems associated with conventional MPM, Bardenhagen and Kober (2004) introduced recently the generalized interpolation material point (GIMP) method and implemented for one-dimensional simulations. In this paper we have extended GIMP to 2D and applied to simulate simple tension and indentation problems. For simulations spanning multiple length scales, based on the continuum mechanics approach, we present a parallel GIMP computational method… More >

G. Himpel, E. Kuhl, A. Menzel, P. Steinmann¹

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 119-134, 2005, DOI:10.3970/cmes.2005.008.119

Abstract The changing mass of biomaterials can either be modelled at the constitutive level or at the kinematic level. This contribution attends on the description of growth at the kinematic level. The deformation gradient will be multiplicatively split into a growth part and an elastic part. Hence, in addition to the material and the spatial configuration, we consider an intermediate configuration or grown configuration without any elastic deformations. With such an ansatz at hand, contrary to the modelling of mass changes at the constitutive level, both a change in density and a change in volume can More >

B. Chrasekhar¹, S. M. Rao²

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 105-118, 2005, DOI:10.3970/cmes.2005.008.105

Abstract In this work, a simple, robust, and an efficient numerical algorithm to calculate the scattered acoustic fields from complex shaped objects such as aircrafts and missiles, subjected to a plane wave incidence is presented. The work is based on the recently proposed method of moments (MoM) and the potential theory, unlike the standard Helmholtz integral equation (HIE) solution method. For the numerical solution, the scattering structure is approximated by planar triangular patches. For the MoM solution of complex bodies involving open/closed/intersecting surfaces, a unified set of basis functions to approximate the source distribution is defined. More >