B.B. Karki¹, D. Bhattarai¹, L. Stixrude²

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 107-118, 2006, DOI:10.3970/cmc.2006.003.107

Abstract Computer modeling of liquid phase poses tremendous challenge: It requires a relatively large simulation size, long simulation time and accurate interatomic interaction and as such, it produces massive amounts of data. Recent advances in hardware and software have made it possible to accurately simulate the liquid phase. This paper reports the details of methodology used in the context of liquid simulations and subsequent analysis of the output data. For illustration purpose, we consider the results for the liquid phases of two geophysically relevant materials, namely MgO and MgSiO₃. The simulations are performed using the parallel first-principles molecular dynamics (FPMD) technique… More >

Kunigal N. Shivakumar¹, Anil Bhargava¹, Sameer Hamoush²

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 97-106, 2006, DOI:10.3970/cmc.2006.003.097

Abstract A detailed three-dimensional finite element stress analysis was conducted on straight-shank and countersunk rivet holes in a plate subjected to tension loading. The study included a wide range of plate width to radius, thickness to radius, countersunk depth to thickness ratios and countersunk angles(θ_c). The stress concentration is maximum at or near the countersunk edge. The stress concentration depends on countersunk depth, plate thickness and width and it is nearly independent of the countersunk angle for 80° ≤ θ_c ≤ 120°. Using the finite element results and limiting conditions, an equation for stress concentration factor is developed and verified. More >

Hamid Sharifi¹ and Augustin Gakwaya¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 77-96, 2006, DOI:10.3970/cmc.2006.003.077

Abstract In this paper, an object-oriented modeling of solid material constitutive behavior using the UML notation is presented. Material properties are first classified into large and small deformation kinematical models. In the small deformation package, we keep classes such as Elastic, ElastoPlastic, ViscoElastic and ViscoPlastic. In the large deformation package, we store classes such as ElastoPlastic, HyperElastic, HyperPlastic, HyperViscoElastic, HyperViscoPlastic and so on. The hierarchical structure, the association relationships as well as key attributes and methods of these classes are presented. We used a C++ implementation of the above model for developing HyperElastic, HyperElastoPlastic and Contact applications in the Diffpack environment. More >

P. B. Wang¹, Z. H. Yao^1,2, T. Lei¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 65-76, 2006, DOI:10.3970/cmc.2006.003.065

Abstract The fast multipole method (FMM) is applied to the dual boundary element method (DBEM) for the analysis of finite solids with large numbers of microcracks. The application of FMM significantly enhances the run-time and memory storage efficiency. Combining multipole expansions with local expansions, computational complexity and memory requirement are both reduced to O(N), where N is the number of DOFs (degrees of freedom). This numerical scheme is used to compute the effective in-plane bulk modulus of 2D solids with thousands of randomly distributed microcracks. The results prove that the IDD method, the differential method, and the method proposed by Feng… More >

Heng-Cheng Lin¹, Chieh Kung², Rong-Sheng Chen¹, Gin-Tiao Liang¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 55-64, 2006, DOI:10.3970/cmc.2006.003.055

Abstract Wafer level chip scale package (WLCSP) has been recognized providing clear advantages over traditional wire-bond package in relaxing the need of underfill while offering high density of I/O interconnects. Without the underfill, the solder joint reliability becomes more critical. Adding to the reliability concerns is the safety demand trend toward "green'' products on which unleaded material, e.g. lead-free solders, is required. The requirement of lead-free solders on the packages results in a higher reflow temperature profile in the package manufacturing process, in turn, complicating the reliability issue. This paper presents an optimization study, considering the fatigue reliability, for a wafer… More >

Shen Yubin¹, Xie Huicai^1,2, Den Wei¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 49-54, 2006, DOI:10.3970/cmc.2006.003.049

Abstract Bonding an overlay of new concrete onto the damaged concrete is a usual repair method. Because of the different shrinkage rate of the new and old concrete, restrained shrinkage cracks will appear in the new concrete. The cracks will reduce durability and strength of the repaired structure. A new repair method using an interface layer of carbon fiber reinforced cement mortar between new and old concrete was developed in this paper. The new method was found to be very effective in reducing shrinkage cracking of repaired beams and slabs. Comparing with normal repaired beams, the maximum observed width of the… More >

Z. Xu¹, L. G. Zhou¹, Jian Wang¹, Timothy S. Cale², Hanchen Huang^1,3

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 43-48, 2006, DOI:10.3970/cmc.2007.003.043

Abstract Recent studies show that three-dimensional Ehrlich-Schwoebel (3D ES), or facet-facet, barriers of face-centered-cubic metals are substantially higher than other surface diffusion barriers. This paper presents the numerical results of 3D ES barriers for body-centered-cubic W, using classical molecular statics calculations and the nudged elastic band method. Results show that an adatom on W{110} has a diffusion barrier of 0.49 eV on the flat surface, 0.66 eV over a monolayer step, and 0.98 eV over a ridge to a neighboring {100} facet, which is one 3D ES barrier. More >

K.Y. Volokh ¹

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 37-42, 2006, DOI:10.3970/cmc.2007.003.037

Abstract Lagrangian or referential equilibrium equations for materials undergoing large deformations are of interest in the developing fields of mechanics of soft biomaterials and nanomechanics. The main feature of these equations is the necessity to deal with the First Piola-Kirchhoff, or nominal, stress tensor which is a two-point tensor referring simultaneously to the reference and current configurations. This two-point nature of the First Piola-Kirchhoff tensor is not always appreciated by the researchers and the total covariant derivative necessary for the formulation of the equilibrium equations in curvilinear coordinates is sometimes inaccurately confused with the regular covariant derivative. Surprisingly, the traditional continuum… More >

Lingfei Gao¹, Xiaoping Zheng^1,2, Zhenhan Yao¹

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 25-36, 2006, DOI:10.3970/cmc.2007.003.025

Abstract A general numerical approach is developed to model the elastic behaviours and failure processes of heterogeneous materials. The heterogeneous material body is assumed composed of a large number of convex polygon lattices with different phases. These phases are locally isotropic and elastic-brittle with the different lattices displaying variable material parameters and a Weibull-type statistical distribution. When the effective strain exceeds a local fracture criterion, the full lattice exhibits failure uniformly, and this is modelled by assuming a very small Young modulus value. An auto-select loading method is employed to model the failure process. The proposed hybrid approach is applied to… More >

I. Smojver¹, J. Sorić²

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 13-24, 2006, DOI:10.3970/cmc.2007.003.013

Abstract The paper addresses the problem of impact on layered fibre composites. The behaviour of composite laminates under impact loading is dependent not only on the velocity but also on the mass and geometry of the impactor. Using micromechanical Mori-Tanaka approach, mechanical properties of the laminate have been calculated utilizing the material constants of the fibre and matrix. General purpose FEM software ABAQUS has been modified by means of user written subroutines for modelling of composite laminate and rigid impactor. The kinematics of the impact has been simulated using transient dynamic analysis. Employing user defined multi point constraints, delamination zones have… More >