• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.4
• /
• pp.367-380
• /
• 2019

A 1kW-class horizontal axis wind turbine (HAWT) rotor blade is taken into account to investigate elastic characteristics in 2-D. The elastic blade field is composed of symmetric cross-ply laminated composite material. Blade element momentum theory is applied to obtain the boundary conditions pressuring the blade, and the plane stress elasticity problem is formulated in terms of two displacement parameters with mixed boundary conditions. For the elastic characteristics a fair of differential equations are derived based on the elastic theory. The domain is divided by triangular and rectangular elements due to the complexity of the blade configuration, and a finite element method is developed for the governing equations to search approximate solutions. The results describe that the elastic behavior is deeply influenced by the layered angle of the middle laminate and the stability of the blade can be improved by controlling the layered angle of laminates, which can be evaluated by the mathematical approach.

• Proceedings of the KSR Conference
• /
• 2000.05a
• /
• pp.439-446
• /
• 2000

In this study, 7 dof (Including warping) beam element was developed to estimate the effects of wagner effects and load height effects on the lateral buckling strength of mono-symmetric I beam. Finite element buckling analysis of mono-symmetric I-shaped girders subjected to transverse loading applied at different heights on the cross-section were conducted. Linear moment gradient were considered, too. In these cases, girders are subjected to both single-curvature and Reverse-curvature bending. An applicability of current LRFD C$\sub$b/ on the mono-symmetric I beam was studied from the finite element results. The problems of current LRFD C$\sub$b/ occurring from load height effects and reverse curvature bending in unbraced length when applied on the mono-symmetric I beam were studied. Solutions to these problems are also presented.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.199-211
• /
• 2024

The loading capacity of engineering structures/components reduces after the initiation and propagation of crack eventually leads to the final failure. Hence, it becomes essential to deal with the crack and its effects at the design and simulation stages itself, by detecting the prone area of the fracture. The phase-field (PF) method has been accepted widely in simulating fracture problems in complex geometries. However, most of the PF methods are formulated with second order continuity theoryinvolving C⁰ continuity. In the present study, PF method based on fourth-order (i.e., higher order) theory, maintaining C¹ continuity has been proposed for ductile fracture simulation. The formulation includes fourth-order derivative terms of phase field variable, varying between 0 and 1. Applications of fourth-order PF theory to ductile fracture simulation resulted in novelty in this area. The proposed formulation is numerically solved using a two-dimensional finite element (FE) framework in 3-layered manner system. The solutions thus obtained from the proposed fourth order theory for different benchmark problems portray the improvement in the accuracy of the numerical results and are well matched with experimental results available in the literature. These results are also compared with second-order PF theory and a comparison study demonstrated the robustness of the proposed model in capturing ductile behaviour close to experimental observations.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.179-188
• /
• 1997

A 4-node element for efficient finite element analysis of plate bending is presented in this paper. This element is formulated based on Mindlin plate theory to take account of shear deformation. To overcome the overestimation of shear stiffness in thin Mindlin plate element, especially in the lower order element, five nonconforming displacement modes are added to the original displacement fields. The proposed nonconforming element does not possess spurious zero-energy mode and does not show shear locking phenomena in very thin plate even for distorted mesh shapes. It was recognized from benchmark numerical tests that the displacement converges to the analytical solutions rapidly and the stress distributions are very smooth. The element also provides good results for the case of high aspect ratio.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.

• Journal of Welding and Joining
• /
• v.18 no.4
• /
• pp.104-110
• /
• 2000

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, an integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in plane finite bodies. The method was developed for in-plane (modes I and II) loadings only. In this paper, the method is formulated and applied to mode III problems involving smooth or kinked cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.implicity.

• Computers and Concrete
• /
• v.26 no.5
• /
• pp.421-427
• /
• 2020

In this study, the two-dimensional generalized finite integral transform(FIT) approach was extended for new accurate thermal buckling analysis of fully clamped orthotropic thin plates. Clamped-clamped beam functions, which can automatically satisfy boundary conditions of the plate and orthogonality as an integral kernel to construct generalized integral transform pairs, are adopted. Through performing the transformation, the governing thermal buckling equation can be directly changed into solving linear algebraic equations, which reduces the complexity of the encountered mathematical problems and provides a more efficient solution. The obtained analytical thermal buckling solutions, including critical temperatures and mode shapes, match well with the finite element method (FEM) results, which verifies the precision and validity of the employed approach.

• Geomechanics and Engineering
• /
• v.4 no.1
• /
• pp.19-37
• /
• 2012

This paper presents the results of analytical and numerical analyses of the effects of performing a pressuremeter test or driving a pile in clay. The geometry of the problem has been simplified by the assumptions of plane strain and axial symmetry. Pressuremeter testing or installation of driven piles has been modelled as an undrained expansion of a cylindrical cavity. Stresses, pore water pressures, and deformations are found by assuming that the clay behaves like normally consolidated modified Cam clay. Closed-form solutions are obtained which allow the determination of the principal effective stresses and the strains around the cavity. The analysis which indicates that the intermediate principal stress at critical state is not equal to the mean of the other two principal stresses, except when the clay is initially isotropically consolidated, also permits finding the limit expansion and excess pore water pressures by means of the Almansi finite strain approach. Results are compared with published data which were determined using finite element and finite difference methods.

• Steel and Composite Structures
• /
• v.18 no.4
• /
• pp.947-970
• /
• 2015

An equivalent friction element is proposed to simulate the friction in cable-strut joints. Equivalent stiffness matrixes and load vectors of the friction element are derived and are unified into patterns for FEM by defining a virtual node specially to store internal forces. Three approaches are described to verify the rationality of the new equivalent friction element: applying the new element in a cable-roller model, and numerical solutions match well with experimental results; applying the element in a continuous sliding cable model, and theoretical values, numerical and experimental results are compared; and the last is applying it in truss string structures, whose results indicate that there would be a great error if the cable of cable supported structures is simulated with discontinuous cable model which is usually adopted in traditional finite element analysis, and that the prestress loss resulted from the friction in cable-strut joints would have adverse effect on the mechanical performance of cable supported structures.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2390-2398
• /
• 2002

For integrity analysis of nuclear reactor pressure vessel, including the Pressurized thermal shock analysis, the fast and accurate calculation of the stress intensity factor at the crack tip is needed. For this, a simple approximation scheme is developed and the resulting stress intensity factors for axial semi-elliptical cracks in cylindrical vessel under various loading conditions are compared with those of the finite element method and other approximation methods, such as Raju-Newman's equation and ASME Sec. Xl approach. For these, three-dimensional finite-element analyses are performed to obtain the stress intensity factors for various surface cracks with t/R = 0.1. The approximation methods, incorporated in VINTIN (Vessel INTegrity analysis-INner flaws), utilizes the influence coefficients to calculate the stress intensity factor at the crack tip. This method has been compared with other solution methods including 3-D finite clement analysis for internal pressure, cooldown, and pressurized thermal shock loading conditions. The approximation solutions are within $\pm$2.5% of the those of FEA using symmetric model of one-forth of a vessel under pressure loading, and 1-3% higher under pressurized thermal shock condition. The analysis results confirm that the VINTIN method provides sufficiently accurate stress intensity factor values for axial semi-elliptical flaws on the surface of the reactor pressure vessel.