• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.5
• /
• pp.932-938
• /
• 2002

In structural applications, brittle materials such as soda-lime glasses and ceramics are usually subjected to multiaxial stress state. Brittle materials with cracks or damage by foreign object impacts are apt to fracture abruptly from cracks, because of their properities of very high strength and low fracture toughness. But in most cases, the residual strength of structural members with damage has been tested under uniaxial stress condition such as the 4-point bend test. Depending upon the crack pattern developed, the strength under multiaxial stress state might be different from the one under uniaxial. A comparative study was carried out to investigate the influence of stress state on the residual strength evaluation. In comparable tests, the residual strength under biaxial stress state by the ball-on-ring test was greater than that under the uniaxial one by the 4-point bend test, when a small size indendation crack was introduced. In the case that crack having an angle of 90deg. to the applied stress direction, the ratio of biaxial to uniaxial flexure strength was about 1.12. The residual strength was different from crack angles to loading direction when it was evaluated by the 4-point bend test. The ratio of residual strength of 45deg. crack to 90deg. one was about 1.20. In the case of specimen cracked by a spherical impact, it was shown that an overall decrease in flexure strength with increasing impact velocity, and the critical impact velocity for formation of a radial and/or cone crack was about 30m/s. In those cases that relatively large cracks were developed as compared with the case of indented cracks, the ratio of residual strength under biaxial stress state to one uniaxial became small.

• Computers and Concrete
• /
• v.13 no.1
• /
• pp.49-70
• /
• 2014

Steel fiber-reinforced concrete (SFRC) is known as one of the efficient modern composites that can greatly enhance the material performance of cracked concrete in tension. Such improved tensile resistance mechanism at crack interfaces in SFRC members can be heavily influenced by methodologies of treatments of crack direction. While most existing studies have focused on developing the numerical analysis model with the rotating-angle theory, there are only few studies on finite element analysis models with the fixed-angle model approach. According to many existing experimental studies, the direction of principal stress rotated after the formation of initial fixed-cracks, but it was also observed that new cracks with completely different angles relative to the initial crack direction very rarely occurred. Therefore, this study introduced the direct tension force transfer model (DTFTM), in which tensile resistance of the fibers at the crack interface can be easily estimated, to the nonlinear finite element analysis algorithm with the fixed-angle theory, and the proposed model was also verified by comparing the analysis results to the SFRC shear panel test results. The secant modulus method adopted in this study for iterative calculations in nonlinear finite element analysis showed highly stable and fast convergence capability when it was applied to the fixed-angle theory. The deviation angle between the principal stress direction and the fixed-crack direction significantly increased as the tensile stresses in the steel fibers at crack interfaces increased, which implies that the deviation angle is very important in the estimation of the shear behavior of SFRC members.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.297-304
• /
• 2019

In this research, crack propagation in a silicon anode during two-phase lithiation was evaluated using a cohesive zone model. The phase transition from crystalline silicon to lithiated silicon causes compressive yielding due to the high volume expansion rate. Li-ion diffuses from the surface of the silicon to its core, and the complex deformation mechanisms during lithiation cause tensile hoop stress along the surface. The Park-Paulino-Roesler (PPR) potential-based cohesive zone model that guarantees consistent energy dissipation in mixed-mode fracture was adopted to simulate edge crack propagation. It was confirmed that the edge crack propagation characteristics during lithiation from the FEM simulation results coincided with the real experimental results. Crack turning observed from real experiments could also be predicted by evaluating the angles of maximum tensile stress directions.

• Smart Structures and Systems
• /
• v.20 no.1
• /
• pp.75-83
• /
• 2017

Effectively monitoring the concrete cracks is an urgent question to be solved in the structural safety monitoring while cracks in hydraulic concrete structures are ubiquitous. In this paper, two experiments are designed based on the measuring principle of Pulse-Pre pump Brillouin Optical Time Domain Analysis (PPP-BOTDA) utilizing Brillouin optical fiber sensor to monitor concrete cracks. More specifically, "V" shaped optical fiber sensor is proposed to determine the position of the initial crack and the experiment illustrates that the concrete crack position can be located by the mutation position of optical fiber strain. Further, Brillouin distributed optical fiber sensor and preinstall cracks are set at different angles and loads until the optical fiber is fractured. Through the monitoring data, it can be concluded that the variation law of optical fiber strain can basically reflect the propagation trend of the cracks in hydraulic concrete structures.

• Carbon letters
• /
• v.19
• /
• pp.89-98
• /
• 2016

Taguchi’s experimental design was employed in the melt spinning of molten mesophase pitch to produce carbon fibers. The textures of the obtained carbon fibers were radial with varied crack angles, as observed by scanning electron microscopy and polarized optical imaging. The diameter, crack angle, preferred orientation, and tensile modulus of the produced samples were examined to investigate the influence of four spinning variables. The relative importance of the variables has been emphasized for each characteristic. The results show that thicker carbon fiber can be obtained with a smaller entry angle, a higher spinning temperature, a reduced winding speed, and an increased extrusion pressure. The winding speed was found to be the most significant factor in relation to the fiber diameter. While it was observed that thicker carbon fiber generally shows improved preferred orientation, the most important variable affecting the preferred orientation was found to be the entry angle. As the entry angle decreased from 120° to 60°, the shear flow was enhanced to induce more ordered radial alignment of crystallite planes so as to obtain carbon fibers with a higher degree of preferred orientation. As a consequence, the crack angle was increased, and the tensile modulus was improved.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.1
• /
• pp.61-68
• /
• 2006

The p-version nonlinear finite element model has been developed to analyze the nonlinear behavior of simply supported RC slabs strengthened with carbon fiber reinforced plastic sheets. The shape function is adopted with integral of Legendre polynomials. The compression model of concrete is based on the Kupfer's yield criterion, hardening rule, and crushing condition. The cracking behavior is modeled by a smeared crack model. In this study, the fixed crack approach is adopted as being geometrically fixed in direction once generated. Each steel layer has a uniaxial behavior resisting only the axial force in the bar direction. Identical behavior is assumed fur tension and compression of steel according to the elastic modulus. The carbon fiber reinforced plastic sheets are considered as reinforced layers of equivalent thickness with uniaxial strength and rigidity properties in the present model. It is shown that the proposed model is able to adequately predicte the displacement and ultimate load of nonlinear simply supported RC slabs by a patch with respect to reinforcement ratio, thickness and angles of CFRP sheets.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.9
• /
• pp.859-869
• /
• 2015

In this study, we analyze the problem of wedge cracks, which are geometrically unsymmetrical in transversely piezoelectric materials. A single concentrated antiplane mechanical load and inplane electrical load are applied at the point of the wedge surface, while one concentrated antiplane load is applied at the crack surface. The crack surfaces are considered as permeable thin slits, where both the normal component of electric displacement and the electric potential are assumed to be continuous across these slits. Using Mellin transform method, the problem is formulated and the Wiener-Hopf equation is derived. By solving the equation, the solution is obtained in a closed form. The intensity factors of the stress and the electric displacement are obtained for any crack length as well as inclined and wedge angles. Based on the results, the intensity factors are independent of the applied electric loads. The electric displacement intensity factor is always dependent on that of stress intensity factor, while the electric field intensity factor is zero. In addition, the energy release rate is computed. These solutions can be used as fundamental solutions which can be superposed to arbitrary electromechanical loadings.

• Steel and Composite Structures
• /
• v.37 no.6
• /
• pp.663-678
• /
• 2020

This paper proposes a novel topology optimization method generating multiple materials for external linear plane crack structures based on the combination of IsoGeometric Analysis (IGA) and eXtended Finite Element Method (X-FEM). A so-called eXtended IsoGeometric Analysis (X-IGA) is derived for a mechanical description of a strong discontinuity state's continuous boundaries through the inherited special properties of X-FEM. In X-IGA, control points and patches play the same role with nodes and sub-domains in the finite element method. While being similar to X-FEM, enrichment functions are added to finite element approximation without any mesh generation. The geometry of structures based on basic functions of Non-Uniform Rational B-Splines (NURBS) provides accurate and reliable results. Moreover, the basis function to define the geometry becomes a systematic p-refinement to control the field approximation order without altering the geometry or its parameterization. The accuracy of analytical solutions of X-IGA for the crack problem, which is superior to a conventional X-FEM, guarantees the reliability of the optimal multi-material retrofitting against external cracks through using topology optimization. Topology optimization is applied to the minimal compliance design of two-dimensional plane linear cracked structures retrofitted by multiple distinct materials to prevent the propagation of the present crack pattern. The alternating active-phase algorithm with optimality criteria-based algorithms is employed to update design variables of element densities. Numerical results under different lengths, positions, and angles of given cracks verify the proposed method's efficiency and feasibility in using X-IGA compared to a conventional X-FEM.