• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.

• Advances in materials Research
• /
• v.5 no.1
• /
• pp.11-20
• /
• 2016

In bonded composite repair of aircraft structures, the damage of the adhesive can thus reduce significantly the efficiency and the durability of the bonded composite repair. The adhesive damage models using critical zone have proven their effectiveness due to simplicity and ap-plicability of the damage criteria in these models. The scope of this study is to analyze the effects of the patch thickness and the adhesive thickness on the damage damage in bonded composite repair of aircraft structures by using modified damage zone theory. The obtained results show that, when the thickness of adhesive increases the damage zone increases and the adhesive loses its rigidity, inversely when the patch is reduced the adhesive damage be-comes more significant.

• Proceedings of the KSR Conference
• /
• 2003.10c
• /
• pp.79-84
• /
• 2003

A fictitious crack model was used to analyze fatigue crack growth under the influence of residual stress. In the fictitious crack model, crack is represented in terms of the separation of two adjacent interfaces and the constitutive equation between the separation and traction is assumed. The effect of fatigue loading was included in the constitutive equation by considering damage accumulation in the cohesive zone. To investigate the effect of the residual stress on the fatigue crack growth, we calculated the residual stress distribution due to transient heat flux to the specimen by finite element method. Fatigue crack growth was simulated by the fictitious crack model with repeated loading. The mode-I crack growth rates were compared for the cases with and without the compressive residual stress around the crack tip. It was observed that the mode-I crack growth can be suppressed by compressive residual stress.

• Interaction and multiscale mechanics
• /
• v.5 no.1
• /
• pp.57-73
• /
• 2012

The paper presents the analysis of two groups of piles subjected to lateral loads incorporating the non-linear behaviour of soil. The finite element method is adopted for carrying out the parametric study of the pile groups. The pile is idealized as a one dimensional beam element, the pile cap as two dimensional plate elements and the soil as non-linear elastic springs using the p-y curves developed by Georgiadis et al. (1992). Two groups of piles, embedded in a cohesive soil, involving two and three piles in series and parallel arrangement thereof are considered. The response of the pile groups is found to be significantly affected by the parameters such as the spacing between the piles, the number of piles in a group and the orientation of the lateral load. The non-linear response of the system is, further, compared with the one by Chore et al. (2012) obtained by the analysis of a system to the present one, except that the soil is assumed to be linear elastic. From the comparison, it is observed that the non-linearity of soil is found to increase the top displacement of the pile group in the range of 66.4%-145.6%, while decreasing the fixed moments in the range of 2% to 20% and the positive moments in the range of 54% to 57%.

• Coupled systems mechanics
• /
• v.5 no.1
• /
• pp.41-58
• /
• 2016

The study deals with physical modeling of a typical building frame resting on pile raft foundation and embedded in cohesive soil mass using finite element based software ETABS. Both- the elements of superstructure and substructure (i.e., foundation) including soil is assumed to remain in elastic state at all the time. The raft is modelled as a thin plate and the pile and soils are treated as interactive springs. Both- the resistance of the piles as well as that of raft base - are incorporated into the model. Interactions between raft-soil-pile are computed. The proposed method makes it possible to solve the problems of uniformly and large non-uniformly arranged piled rafts in a time saving way using finite element based software ETABS. The effect of the various parameters of the pile raft foundation such as thickness of raft and pile diameter is evaluated on the response of superstructure. The response included the displacement at the top of the frame and bending moment in columns. The soil-structure interaction effect is found to increase displacement and increase the absolute maximum positive and negative moments. The effect of the soil- structure interaction is observed to be significant for the type of foundation and soil considered in the present study.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.2 no.2
• /
• pp.19-28
• /
• 1982

The development of unsteady, depth-averaged two dimensional sediment transport prediction model in estuaries and harbors by the Galerkin finite element technique is presented. The model consists of two submodels, flow induced circulation model and sediment transport model. The sediment transport submodel is formulated by incorporating sediment continuity equation and sediment diffusion equation. Numerical experiments of the model, which were carried out in one dimensional channel under different conditions for circulation and sediment transport, show the adaptability of the formulation for predicting the migration of both cohesive and noncohesive sediments. The model was applied to Busan harbor to simulate circulation and sediment transport for simplified conditions. Of the results by the model the flow pattern are shown to be similar to observed data.

• International Journal of Concrete Structures and Materials
• /
• v.18 no.3E
• /
• pp.191-198
• /
• 2006

The objective of this study is to examine fracture characteristics of concrete at early ages, i.g. critical stress intensity factor, critical crack-tip opening displacement, fracture energy, and bilinear softening curve based on the concepts of effective-elastic crack model and cohesive crack model. A wedge splitting test for Mode I was performed on cubic wedge specimens with a notch at the edge. By experimenting with various strengths and ages, load-crack mouth opening curves were obtained, and the results were analyzed by linear elastic fracture mechanics and FEM(finite element method). The results from the test and analysis showed that critical stress intensity factor and facture energy increased while critical crack-tip opening displacement decreased with concrete aging from 1 day to 28 days. Four parameters of bilinear softening curve from 1 day to 28 days were obtained from a numerical analysis. The obtained fracture parameters and bilinear softening curves at early ages from this study are to be used as a fracture criterion and an input data for the finite element analysis of concrete at early ages.

• Steel and Composite Structures
• /
• v.26 no.3
• /
• pp.265-272
• /
• 2018

A detailed study was carried out on the tensile properties of the single-lap joint of a steel panel bolted/bonded to a composite laminate with a flanging. Finite element model (FEM) was established to predict the strength and to analyze the damage propagation of the hybrid joints by ABAQUS/Standard, which especially adopted cohesive elements to simulate the interface between the laminate and adhesive. The strength and failure mode predicted by FEM were in good agreement with the experimental results. In addition, three influence factors including adhesive thickness, bolt preload and bolt-hole clearance were studied. The results show that the three parameters have effect on the first drop load of the load-displacement curve, but the effect of bolt-hole clearance is the largest. The bolt-hole clearance should be avoided for hybrid joints.

• Computers and Concrete
• /
• v.29 no.6
• /
• pp.393-405
• /
• 2022

Lightweight concrete is a superior material due to its light weight and high strength. There however remain significant lacunae in engineering knowledge with regards to shear failure of lightweight fiber reinforced concrete beams. The main aim of the present study is to investigate the optimum usage of steel fibers in lightweight fiber reinforced concrete (LWFRC). Multi-scale finite element model calibrated with experimental results is developed to study the effect of steel fibers on the mechanical properties of LWFRC beams. To decrease the amount of steel fibers, it is preferred to reinforce only the middle section of the LWFRC beams, where the flexural stresses are higher. For numerical simulation, a multi-scale finite element model was developed. The cement matrix was modeled as homogeneous and uniform material and both steel fibers and lightweight coarse aggregates were randomly distributed within the matrix. Considering more realistic assumptions, the bonding between fibers and cement matrix was considered with the Cohesive Zone Model (CZM) and its parameters were determined using the model update method. Furthermore, conformity of Load-Crack Mouth Opening Displacement (CMOD) curves obtained from numerical modeling and experimental test results of notched beams under center-point loading tests were investigated. Validating the finite element model results with experimental tests, the effects of fibers' volume fraction, and the length of the reinforced middle section, on flexural and residual strengths of LWFRC, were studied. Results indicate that using steel fibers in a specified length of the concrete beam with high flexural stresses, and considerable savings can be achieved in using steel fibers. Reducing the length of the reinforced middle section from 50 to 30 cm in specimens containing 10 kg/m3 of steel fibers, resulting in a considerable decrease of the used steel fibers by four times, whereas only a 7% reduction in bearing capacity was observed. Therefore, determining an appropriate length of the reinforced middle section is an essential parameter in reducing fibers, usage leading to more affordable construction costs.

• Geomechanics and Engineering
• /
• v.9 no.5
• /
• pp.645-667
• /
• 2015

Single treatment and staged treatments in vertical wells are widely applied in sandstone and mudstone thin interbedded (SMTI) reservoir to stimulate the reservoir. The keys and difficulties of stimulating this category of formations are to avoid hydraulic fracture propagating through the interface between shale and sand as well as control the fracture height. In this paper, the cohesive zone method was utilized to build the 3-dimensional fracture dynamic propagation model in shale and sand interbedded formation based on the cohesive damage element. Staged treatments and single treatment were simulated by single fracture propagation model and double fractures propagation model respectively. Study on the changes of fracture vicinity stress field during propagation is to compare and analyze the parameters which influence the interfacial induced stresses between two different fracturing methods. As a result, we can prejudge how difficult it is that the fracture propagates along its height direction. The induced stress increases as the pumping rate increasing and it changes as a parabolic function of the fluid viscosity. The optimized pump rate is $4.8m^3/min$ and fluid viscosity is $0.1Pa{\cdot}s$ to avoid the over extending of hydraulic fracture in height direction. The simulation outcomes were applied in the field to optimize the treatment parameters and the staged treatments was suggested to get a better production than single treatment.