P.C. Chen^∗,†, C.W. Colwell^†, D.D. D’Lima^†,‡

Molecular & Cellular Biomechanics, Vol.8, No.2, pp. 135-148, 2011, DOI:10.3970/mcb.2011.008.135

Abstract A nonlinear viscoelastic finite element model of ultra-high molecular weight polyethylene (UHMWPE) was developed in this study. Eight cylindrical specimens were machined from ram extruded UHMWPE bar stock (GUR 1020) and tested under constant compression at 7% strain for 100 sec. The stress strain data during the initial ramp up to 7% strain was utilized to model the "instantaneous" stress-strain response using a Mooney-Rivlin material model. The viscoelastic behavior was modeled using the time-dependent relaxation in stress seen after the initial maximum stress was achieved using a stored energy formulation. A cylindrical model of similar dimensions was created using a… More >

K. Y. Volokh^*, E. Y. S. Chao^†, M. Armand^‡

Molecular & Cellular Biomechanics, Vol.4, No.2, pp. 67-74, 2007, DOI:10.3970/mcb.2007.004.067

Abstract Information about the stress distribution on contact surfaces of adjacent bones is indispensable for analysis of arthritis, bone fracture and remodeling. Numerical solution of the contact problem based on the classical approaches of solid mechanics is sophisticated and time-consuming. However, the solution can be essentially simplified on the following physical grounds. The bone contact surfaces are covered with a layer of articular cartilage, which is a soft tissue as compared to the hard bone. The latter allows ignoring the bone compliance in analysis of the contact problem, i.e. rigid bones are considered to interact through a compliant cartilage. Moreover, cartilage… More >

Zhao Xilong¹, Zhang Jianxun¹, Chen Hongyuan^1,2, Song Xu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.92, No.6, pp. 557-571, 2013, DOI:10.3970/cmes.2013.092.557

Abstract In order to study failure perspective of Ti-6Al-4V titanium alloy welded joint, fracture behavior of the heat affected zone is researched. Microhardness, tensile test and microstructure are used to study the mechanical properties of the titanium alloy laser welded joint. The tensile and microhardness results show that heat affected zone near the base metal is the weakest in welded joints. This paper is based on the results from in-situ tension test to observe the process of crack expansion in heat affected zone. And then, in-situ test is simulated via the finite element method on ABAQUS Standard; therefore the fracture criteria… More >

Y.J. Ning^1,2,3, Z.Y. Zhao³, J.P. Sun³, W.F. Yuan^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.89, No.3, pp. 221-262, 2012, DOI:10.3970/cmes.2012.089.221

Abstract In this paper, wave propagations in jointed rock masses are modeled by the discontinuous deformation analysis (DDA) method. The selection of the numerical control parameters in the DDA for wave propagation modeling is discussed in detail, and the effects of the joint stiffness, the seismic loading frequency, the joint strength, and the incident angle on the propagations of stress waves in a jointed rock mass are modeled and analyzed. Two nonreflecting boundary conditions including the viscous boundary condition (VBC) and the superposition boundary condition (SBC) are coupled into the DDA. The applicability of the two nonreflecting boundary conditions for simple… More >

Anant Diwakar¹, D. N.Srinath¹, Sanjay Mittal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.1, pp. 61-90, 2010, DOI:10.3970/cmes.2010.069.061

Abstract Aerodynamic shape optimization of airfoils is carried out for two values of Reynolds numbers: 10³ and 10⁴, for an angle of attack of 5^o. The objective functions used are (a) maximization of lift (b) minimization of drag and (c) minimization of drag to lift ratio. The surface of the airfoil is parametrized by a 4^th order non-uniform rational B-Spline (NURBS) curve with 61 control points. Unlike the efforts in the past, the relatively large number of control points used in this study offer a rich design shape with the possibility of local bumps and valleys on the airfoil surface. The… More >

D.N. Srinath¹, Sanjay Mittal¹, Veera Manek²

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.2, pp. 169-190, 2009, DOI:10.3970/cmes.2009.051.169

Abstract A continuous adjoint method is formulated and implemented for the multi-point shape optimization of airfoils at low Re. The airfoil shape is parametrized with a non-uniform rational B-Spline (NURBS). Optimization studies are carried out for two different objective functions. The first involves an inverse function on the lift coefficient over a range of Re. The objective is to determine a shape that results in a lift coefficient of 0.4 at three values of Re: 10, 100 and 500. The second objective involves a direct function on the lift coefficient over a range of angles of attack,a. The lift coefficient is… More >

G. Cocchetti², G. Maier², X. P. Shen³

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.3, pp. 279-298, 2002, DOI:10.3970/cmes.2002.003.279

Abstract Interface models mean here relationships between displacement jumps and tractions across a locus of displacement discontinuities. Frictional contact and quasi-brittle fracture interpreted by cohesive crack models are typical mechanical situations concerned by the present unifying approach. Plastic-softening multidissipative interface models are studied in piecewise linear formats, i.e. assuming linearity for yield functions, plastic potentials and relationships between static and kinematic internal variables. The properties and the pros and cons of such simplified models in a variety of formulations (fully non-holonomic in rates, holonomic and in finite steps), all mathematically described as linear complementarity problems, are comparatively investigated in view of… More >

K.I. Tserpes¹, A.S. Koumpias¹

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.2, pp. 169-182, 2012, DOI:10.32604/cmes.2012.083.169

Abstract In this work, a comparison between a cohesive zone model and a continuum damage model in predicting the mode-I fracture behavior of adhesively bonded joints is performed on the basis of reliability and applicability. The cohesive zone model (CZM) is based on an exponential traction law characterizing the behavior of the interface elements. The continuum damage model (CDM) is based on the stiffness degradation of adhesive elements imposed by a damage parameter. Both models have been implemented by means of a 3D finite element model. Mode-I fracture behavior of the bonded joints was characterized using the DCB specimen. Firstly, the… More >

C.J. Zheng¹, H.B. Chen¹, T. Matsumoto², T. Takahashi²

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.1, pp. 1-30, 2011, DOI:10.3970/cmes.2011.079.001

Abstract A fast multipole boundary element approach to the shape sensitivity analysis of three dimensional acoustic wave problems is developed in this study based on the adjoint variable method. The concept of material derivative is employed in the derivation. The Burton-Miller formula which is a linear combination of the conventional and normal derivative boundary integral equations is adopted to cope with the non-uniqueness problem when solving exterior acoustic wave problems. Constant elements are used to discretize the boundary surface so that the strongly- and hyper-singular boundary integrals contained in the formulations can be evaluated explicitly and the numerical process can be… More >

T. Matsumoto¹, T. Yamada¹, T. Takahashi¹, C.J. Zheng², S. Harada¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 77-94, 2011, DOI:10.3970/cmes.2011.078.077

Abstract Shape design and topology sensitivity formulations for acoustic problems based on adjoint method and the boundary element method are presented and are applied to shape sensitivity analysis and topology optimization of acoustic field. The objective function is assumed to consist only of boundary integrals and quantities defined at certain number of discrete points. The adjoint field is defined so that the sensitivity of the objective function does not include the unknown sensitivity coefficients of the sound pressures and particle velocities on the boundary and in the domain. Since the final sensitivity expression does not have the sensitivity coefficients of the… More >