Francesca Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.3&4, pp. 269-282, 2011, DOI:10.3970/cmes.2011.074.269

Abstract In industrial applications, microstructure inhomogeneities can derive from the manufacturing process and the final mechanical properties of the material depend on the resulting, complex, structural pattern of the constituents. In this paper, Cell Method plane models in the elastic and plastic fields are presented and applied to predict the behaviour of four sintered alloys, where the spatial arrangement of voids within the base material contributes to determine the mechanical behaviour. Unlike the Finite Elements and other methods, the Cell Method is a numerical method based on a direct discrete formulation of equilibrium equations, so that More >

G. R. Liu^{1, 2}

The International Conference on Computational & Experimental Engineering and Sciences, Vol.12, No.4, pp. 113-114, 2009, DOI:10.3970/icces.2009.012.113

Abstract This paper introduces first a weakened weakform (W2) using a generalized gradient smoothing technique for an unified formulation of a wide class of compatible and incompatible displacement methods including settings of the finite element methods (FEM) and meshfree methods of special properties including the upper bound properties. A G space is first defined to include discontinuous functions allowing the use of much more types of methods/techniques to create shape functions for numerical models; Properties and a set of important inequalities for G spaces are then proven in theory and analyzed in detail. We prove that… More >

X.Y. Cui^1,2, G. R. Liu^2,3, G. Y. Li^1,4, G. Zheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 217-230, 2008, DOI:10.3970/cmes.2008.038.217

Abstract In this paper, a gradient smoothed formulation is proposed to deal with a fourth-order differential equation of Bernoulli-Euler beam problems for static and dynamic analysis. Through the smoothing operation, the C¹ continuity requirement for fourth-order boundary value and initial value problems can be easily relaxed, and C⁰ interpolating function can be employed to solve C¹ problems. In present thin beam problems, linear shape functions are employed to approximate the displacement field, and smoothing domains are further formed for computing the smoothed curvature and bending moment field. Numerical examples indicate that very accurate results can be yielded when More >

P. Bettini, E. Brusa, M. Munteanu, R. Specogna, F. Trevisan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.3, pp. 215-242, 2008, DOI:10.3970/cmes.2008.033.215

Abstract Electrostatic microactuator is a paradigm of MEMS. Cantilever and double clamped microbeams are often used in microswitches, microresonators and varactors. An efficient numerical prediction of their mechanical behaviour is affected by the nonlinearity of the electromechanical coupling. Sometimes an additional nonlinearity is due to the large displacement or to the axial-flexural coupling exhibited in bending. To overcome the computational limits of the available numerical methods two new formulations are here proposed and compared. Modifying the classical beam element in the Finite Element Method to allow the implementation of a \emph {Non incremental sequential approach} is… More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 187-196, 2008, DOI:10.3970/cmes.2008.027.187

Abstract The Cell Method is a recent numerical method that can be applied in several fields of physics and engineering. In this paper, the elastodynamics formulation is extended to include system internal damping, highlighting some interesting characteristics of the method. The developed formulation leads to an explicit solving system. The mass matrix is diagonal (without lumping) and in the most general case a time-dependent damping coefficient can be defined for each node. \newline Accuracy and convergence rate have been tested with reference to the classical problem of a particle free vibration with viscous damping.
An application More >

M-A Lavoie¹, A. Gakwaya¹, M. Nejad Ensan², D.G. Zimcik²

The International Conference on Computational & Experimental Engineering and Sciences, Vol.2, No.4, pp. 111-118, 2007, DOI:10.3970/icces.2007.002.111

Abstract This paper reviews numerical methods that are currently available to simulate bird strike as well as the theory of the event. It also summarizes important parameters and provides guidelines as to how to set up the analysis and how to evaluate a model. The information provided is based on physical properties and available results regarding a bird and its behaviour upon impact. The simulations have been performed with LS-DYNA 970 but can be done in similar dynamic finite elements analysis codes. 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 >

V. Cristini¹, Y. Renardy²

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.2, pp. 77-94, 2006, DOI:10.3970/fdmp.2006.002.077

Abstract We review studies of a drop of viscous liquid, suspended in another liquid, and undergoing breakup in an impulsively started shear flow. Stokes flow conditions as well as the effects of inertia are reported. They reveal a universal scaling for the fragments, which allows one to use sheared emulsions to produce monodispersity as an alternative to microfluidic devices. More >

R.L. Huston¹, C.-Q. Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.2, pp. 143-152, 2005, DOI:10.3970/cmes.2005.010.143

Abstract This paper presents a summary of recent developments in computational methods for multibody dynamics analyses. The developments are presented within the context of an automated numerical analysis. The intent of the paper is to provide a basis for the easy development of computational algorithms. The principal concepts discussed are: differentiation algorithms, partial velocities and partial angular velocities, generalized speeds, Euler parameters, Kane's equations, orthogonal complement arrays, lower body arrays and accuracy testing functions. More >

J. A. Nairn¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 649-664, 2003, DOI:10.3970/cmes.2003.004.649

Abstract A new algorithm is described which extends the material point method (MPM) to allow explicit cracks within the model material. Conventional MPM enforces velocity and displacement continuity through its background grid. This approach is incompatible with cracks which are displacement and velocity discontinuities. By allowing multiple velocity fields at special nodes near cracks, the new method (called CRAMP) can model cracks. The results provide an ``exact'' MPM analysis for cracks. Comparison to finite element analysis and to experiments show it gets good results for crack problems. The intersection of crack surfaces is prevented by implementing More >