Jihan Kim¹, Yong Hyup Kim¹, Sung Won Lee²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 263-276, 2004, DOI:10.3970/cmes.2004.006.263

Abstract The Koiter's asymptotic method is combined with the assumed strain solid shell element formulation for postbuckling analysis of composite and sandwich structures. The assumed strain solid shell element is free of locking and the small angle assumption, and it allows multiple plies through the element thickness. While laminated composite structures are modeled with single element through the thickness, sandwich structures are modeled with three elements stacked through the thickness to model the face sheets and the core independently. The Koiter's method is used to trace initial postbuckling path. Subsequently, the Koiter's method is switched to More >

B. Jin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 253-262, 2004, DOI:10.3970/cmes.2004.006.253

Abstract In this paper, we propose a new numerical scheme for the solution of the Laplace and biharmonic equations subjected to noisy boundary data. The equations are discretized by the method of fundamental solutions. Since the resulting matrix equation is highly ill-conditioned, a regularized solution is obtained using the truncated singular value decomposition, with the regularization parameter given by the L-curve method. Numerical experiments show that the method is stable with respect to the noise in the data, highly accurate and computationally very efficient. More >

L. G. Zhou¹, L. M. Hsiung², Hanchen Huang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 245-252, 2004, DOI:10.3970/cmes.2004.006.245

Abstract Using molecular dynamics simulations, we study the deformation of polysynthetically twinned (PST) TiAl at room temperature. The simulation cell is pre-strained and thermodynamically relaxed to zero stress, so that no dislocations pre-exist in γ−α₂ interfaces. A uniaxial compression is then applied along one 1/6<112] direction. Our results show that interfacial dislocation pairs nucleate at the γ−α₂ interface under the compression. The glide and agglomeration of these dislocations lead to the nucleation of deformation twins from the interface. Based on our studies, twins may nucleate without pre-existing interfacial dislocations. Further we have monitored the propagation of the More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 227-244, 2004, DOI:10.3970/cmes.2004.006.227

Abstract Defining the crack path in brittle and non-brittle crack is not easy, due to several unknowns. If the direction of crack propagation can be computed by means of one of the existing criteria, it is not known whether this direction will remain constant during crack propagation. A crack initiation leads to an enhanced stress field at crack tip. During propagation, the enhanced tip stress field propagates into the solid, locally interacting with the pre-existing stress field. This interaction can lead to modifications of the propagation direction, depending on the domain and crack geometry. Moreover, trajectory… More >

F. Auricchio¹, L. Petrini², R. Pietrabissa³, E. Sacco⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 365-380, 2003, DOI:10.3970/cmes.2003.004.365

Abstract Since 80’s many devices were developed to exploit the unique blend of mechanical and biocompatibility properties of shape memory alloys in orthodontic applications. It results in a high clinical effectiveness, but also in a spreading of technical knowledge on the properties of the single appliances. The goal of the present contribution is to contrast this sense of bewilderness and to prepare the basis for a simulationtool able to support the orthodontist choice. In particular a finite-element beam with a one-dimensional constitutive law, able to describe the SMA super elasticity and shape memory effect, is presented: More >

P. Vena¹, R. Contro

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 497-506, 2003, DOI:10.3970/cmes.2003.004.497

Abstract The adoption of artificial ligaments in current surgery is still characterised by a low success rate due to the fact that mechanical properties of the biomedical devices are such that a biomechanical compatibility is not fully satisfied. A durable artificial ligament should exhibit stiffness as well as strength properties which are such that a full articulation functionality is guaranteed. To this purpose, reliable numerical methods able to predict the mechanical behaviour of such devices both in the elastic and in inelastic range until complete rupture, could be used for designing of devices with tailored mechanical More >

M.C. Ho Ba Tho¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 489-496, 2003, DOI:10.3970/cmes.2003.004.489

Abstract The objective of the paper is to address the methodology developed to model bone and joints with individualised geometric and material properties from medical image data. An atlas of mechanical properties of human bone has been investigated demonstrating individual differences. From these data, predictive relationships have been established between mechanical properties and quantitative data derived from measurements on medical images. Subsequently, geometric and numerical models of bones with individualised geometrical and mechanical properties have been developed from the same source of image data. The advantages of this modelling technique are its ability to study the More >

M. Thiriet¹, S. Naili², C. Ribreau²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 473-488, 2003, DOI:10.3970/cmes.2003.004.473

Abstract The laminar steady flow of incompressible Newtonian fluid is studied in rigid pipes with cross configuration of a collapsed tube to determine both the entry length and the wall shear stress (WSS). The cross section shapes have been defined from the collapse of an infinitely long elastic tube subjected to an uniform transmural pressure. Five characteristic collapsed configurations, from the unstressed down to the point-contact states, with a finite and infinite curvature radius at the contact point, are investigated, although the wall contact is not necessary observed in veins. Such collapsed shapes induce cross gradient More >

A.R.Skovoroda¹, R.V.Goldstein²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 457-472, 2003, DOI:10.3970/cmes.2003.004.457

Abstract We discuss the inverse problem of the recovery of the distribution of the elastic moduli of a stratified layer, based on measurements of the surface displacement under localized surface loads. A general parametric solution and a numerical procedure for computing the parameters are presented. Examples of numerical results are given. The problem and its solution are related to the monitoring of elastic properties of living tissues. More >

I. J. Rao¹, J.D. Humphrey², K.R. Rajagopal³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 439-456, 2003, DOI:10.3970/cmes.2003.004.439

Abstract Recent discoveries in molecular and cell biology reveal that many cell types sense and respond (via altered gene expression) to changes in their mechanical environment. Such mechanotransduction mechanisms are responsible for many changes in structure and function, including the growth and remodeling process. To understand better, and ultimately to use (e.g., in tissue engineering), biological growth and remodeling, there is a need for mathematical models that have predictive and not just descriptive capability. In contrast to prior models based on reaction-diffusion equations or the concept of volumetric growth, we examine here a newly proposed constrained More >