K. Jayabal, A. Arockiarajan, S.M. Sivakumar¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 111-124, 2008, DOI:10.3970/cmes.2008.027.111

Abstract A three dimensional micromechanically motivated model is proposed here based on firm thermodynamics principles to capture the nonlinear dissipative effects in the polycrystal ferroelectrics. The constraint imposed by the surrounding grains on a subgrain at its boundary during domain switching is modeled by a suitable modification of the switching threshold in a subgrain. The effect of this modification in the dissipation threshold is studied in the polycrystal behavior after due correlation of the subgrain behavior with the single crystal experimental results found in literature. Taking into consideration, all the domain switching possibilities, the volume fractions More >

Jyoti Ranjan Majhi, Ranjan Ganguli¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 25-36, 2008, DOI:10.3970/cmes.2008.027.025

Abstract The flapping equation for a rotating rigid helicopter blade is typically derived by considering 1) small flap angle, 2) small induced angle of attack and 3) linear aerodynamics. However, the use of nonlinear aerodynamics can make the assumptions of small angles suspect. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived in this paper with nonlinear aerodynamics . Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are More >

Y. Nishidate¹, G. P. Nikishkov^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 91-106, 2008, DOI:10.3970/cmes.2008.026.091

Abstract Atomic-scale finite element procedure for modeling of self-positioning nanostructures is developed. Our variant of the atomic-scale finite element method is based on a meshless approach and on the Tersoff interatomic potential function. The developed algorithm is used for determination of equilibrium configuration of atoms after nanostructure self-positioning. Dependency of the curvature radius of nanostructures on their thickness is investigated. It is found that for thin nanostructures the curvature radius is considerably smaller than predicted by continuum mechanics equations. Curvature radius variation with varying orientation of crystallographic axes is also modeled and results are compared to More >

Zdenko Tonković¹, Jurica Sorić^1,2, Ivica Skozrit¹

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.026.075

Abstract An efficient numerical algorithm for modeling of cyclic elastoplastic deformation of shell structures is derived. The constitutive model includes highly nonlinear multi-component forms of kinematic and isotropic hardening functions in conjunction with von Mises yield criterion. Therein, the closest point projection algorithm employing the Reissner-Mindlin type kinematic model, completely formulated in tensor notation, is applied. A consistent elastoplastic tangent modulus ensures high convergence rates in the global iteration approach. The integration algorithm has been implemented into a layered assumed strain isoparametric finite shell element, which is capable of geometrical nonlinearities including finite rotations. Numerical examples, More >

P. Le¹, N. Mai-Duy², T. Tran-Cong³, G. Baker⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 43-68, 2008, DOI:10.3970/cmes.2008.025.043

Abstract This paper describes an integrated radial basis function network (IRBFN) method for the numerical modelling of the dynamics of strain localization due to strain softening in quasi-brittle materials. The IRBFN method is a truly meshless method that is based on an unstructured point collocation procedure. We introduce a new and effective regularization method to enhance the performance of the IRBFN method and alleviate the numerical oscillations associated with weak discontinuity at the elastic wave front. The dynamic response of a one dimensional bar is investigated using both local and non-local continuum models. Numerical results, which More >

G. Fitzgerald¹, G. Goldbeck-Wood², P. Kung¹, M. Petersen¹, L. Subramanian¹, J. Wescott²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 169-184, 2008, DOI:10.3970/cmes.2008.024.169

Abstract Molecular modeling has established itself as an important component of applied research in areas such as drug discovery, catalysis, and polymers. Algorithmic improvements to these methods coupled with the increasing speed of computational hardware are making it possible to perform predictive modeling on ever larger systems. Methods are now available that are capable of modeling hundreds of thousands of atoms, and the results can have a significant impact on real-world engineering problems. The article reviews some of the modeling methods currently in use; provides illustrative examples of applications to challenges in sensors, fuel cells, and More >

Shueei-Muh Lin¹, Sen-Yung Lee², Yu-Sheng Lin³

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 175-186, 2008, DOI:10.3970/cmes.2008.023.175

Abstract The blade of a horizontal-axis wind power turbine is modeled as a rotating beam with pre-cone angles and setting angles. Based on the Bernoulli-Euler beam theory, without considering the axial extension deformation and the Coriolis forces effect, the governing differential equations for the bending vibration of the beam are derived. It is pointed out that if the geometric and the material properties of the beam are in polynomial forms, then the exact solution for the system can be obtained. Based on the frequency relations as revealed, without tedious numerical analysis, one can reach many general More >

A.P.S. Selvadurai¹, H. Ghiabi²

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 53-74, 2008, DOI:10.3970/cmes.2008.023.053

Abstract This paper deals with the problem of the consolidation of a composite consisting of alternate layers of soft clay and a granular material. A series of experiments were conducted on components to develop the constitutive models that can be implemented in a computational approach. The constitutive response of the soft clay is represented by a poro-elasto-plastic Cam clay-based model and the granular medium by an elasto-plastic model with a Drucker-Prager type failure criterion and a non-associated flow rule. The computational poro-elasto-plastic model is used to calibrate the experimental results derived from the one-dimensional tests and More >

A.Z. Lorbiecka¹, R.Vertnik², H.Gjerkeš¹, G. Manojlovič², B.Senčič², J. Cesar², B.Šarler^1,3

CMC-Computers, Materials & Continua, Vol.8, No.3, pp. 195-208, 2008, DOI:10.3970/cmc.2008.008.195

Abstract A numerical model is developed for the simulation of solidification grain structure formation (equiaxed to columnar and columnar to equiaxed transitions) during the continuous casting process of steel billets. The cellular automata microstructure model is combined with the macroscopic heat transfer model. The cellular automata method is based on the Nastac's definition of neighborhood, Gaussian nucleation rule, and KGT growth model. The heat transfer model is solved by the meshless technique by using local collocation with radial basis functions. The microscopic model parameters have been adjusted with respect to the experimental data for steel 51CrMoV4. More >

S. Baragetti¹, F. Tordini²

Structural Durability & Health Monitoring, Vol.3, No.3, pp. 165-176, 2007, DOI:10.3970/sdhm.2007.003.165

Abstract The effect of a TiN PVD (Physical Vapor Deposition) coating on the fatigue behaviour of the titanium alloy Ti-6Al-4V was investigated. Fatigue tests were performed on coated and uncoated, both smooth and 120° V-notched, specimens in order to evaluate the influence of the coating on the substrate fatigue resistance. Numerical analyses were carried out in order to determine the stress distributions below the specimen surface and on the coating. Several coating elastic moduli were used in such calculations. The residual stress gradient induced by the coating process deposition and the substrate plasticization were also taken More >