E.J. Sapountzakis¹, V.G. Mokos²

CMC-Computers, Materials & Continua, Vol.10, No.1, pp. 1-40, 2009, DOI:10.3970/cmc.2009.010.001

Abstract In this paper the boundary element method is employed to develop a displacement solution for the general transverse shear loading problem of composite beams of arbitrary constant cross section. The composite beam (thin or thick walled) consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli and are firmly bonded together. The analysis of the beam is accomplished with respect to a coordinate system that has its origin at the centroid of the cross section, while its axes are not necessarily the principal bending ones. The transverse… More >

T. Nishioka¹, S. Tchouikov¹, T. Fujimoto¹

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 147-154, 2006, DOI:10.3970/cmc.2006.003.147

Abstract In this study, phenomena of multiple branching of dynamically propagating crack are investigated numerically. The complicated paths of cracks propagating in a material are simulated by moving finite element method based on Delaunay automatic triangulation (MFEM BODAT), which was extended for such problems. For evaluation of fracture parameters for propagating and branching cracks switching method of the path independent dynamic J integral was used. Using these techniques the generation phase simulation of multiple dynamic crack branching was performed. Various dynamic fracture parameters, which are almost impossible to obtain by experimental technique alone, were accurately evaluated. 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 >

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 >

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 >

Aimin Jiang^1,2, Haojiang Ding²

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 1-12, 2006, DOI:10.3970/cmc.2007.003.001

Abstract This paper presents a development of the boundary contour method (BCM) for magneto-electro-elastic media. Firstly, the divergence-free of the integrand of the magneto- electro-elastic boundary element is proved. Secondly, the boundary contour method formulations are obtained by introducing linear shape functions and Green's functions (Computers & Structures, 82(2004):1599-1607) for magneto-electro-elastic media and using the rigid body motion solution to regularize the BCM and avoid computation of the corner tensor. The BCM is applied to the problem of magneto-electro-elastic media. Finally, numerical solutions for illustrative examples are compared with exact ones and those of the conventional boundary element method (BEM). The… More >

Eduardo Divo¹, Alain J. Kassab²

CMC-Computers, Materials & Continua, Vol.2, No.4, pp. 277-288, 2005, DOI:10.3970/cmc.2005.002.277

Abstract A Dual Reciprocity Boundary Element Method is formulated to solve non-linear heat conduction problems. The approach is based on using the Kirchhoff transform along with lagging of the effective non-linear thermal diffusivity. A posteriori error estimate is used to provide effective estimates of the temporal and spatial error. A numerical example is used to demonstrate the approach. More >

Carlos J. S. Alves, Pedro R. S. Antunes¹

CMC-Computers, Materials & Continua, Vol.2, No.4, pp. 251-266, 2005, DOI:10.3970/cmc.2005.002.251

Abstract In this work we show the application of the Method of Fundamental Solutions(MFS) in the determination of eigenfrequencies and eigenmodes associated to wave scattering problems. This meshless method was already applied to simple geometry domains with Dirichlet boundary conditions (cf. Karageorghis (2001)) and to multiply connected domains (cf. Chen, Chang, Chen, and Chen (2005)). Here we show that a particular choice of point-sourcescan lead to very good results for a fairly general type of domains. Simulations with Neumann boundary conditionare also considered. More >