• Structural Engineering and Mechanics
• /
• v.37 no.3
• /
• pp.285-296
• /
• 2011

The plane crack-contact problem for an infinite elastic layer with two symmetric rectangular rigid stamps on its upper and lower surfaces is considered. The elastic layer having an internal crack parallel to its surfaces is subjected to two concentrated loads p on its upper and lower surfaces trough the rigid rectangular stamps and a pair of uniform compressive stress $p_0$ along the crack surface. It is assumed that the contact between the elastic layer and the rigid stamps is frictionless and the effect of the gravity force is neglected. The problem is reduced to a system of singular integral equations in which the derivative of the crack surface displacement and the contact pressures are unknown functions. The system of singular integral equations is solved numerically by making use of an appropriate Gauss-Chebyshev integration formula. Numerical results for stress-intensity factor, critical load factor, $\mathcal{Q}_c$, causing initial closure of the crack tip, the crack surface displacements and the contact stress distribution are presented and shown graphically for various dimensionless quantities.

• Structural Engineering and Mechanics
• /
• v.83 no.1
• /
• pp.67-77
• /
• 2022

In this study, frictionless continuous and discontinuous contact problems of a magneto-electro-elastic layer in the presence of the body force were discussed. The layer was indented by a rigid cylindrical insulating punch and supported by a rigid substrate without bond. Applying the Fourier integral transform technique, the general expressions of the problem were derived in the presence of body force. Thanks to the boundary conditions, the singular integral equations were obtained for both the continuous and the discontinuous contact cases. Gauss-Chebyshev integration formulas were used to transform the singular integral equations into a set of nonlinear equations. Contact width under the punch, initial separation distance, critical load, separation regions and contact stress under the punch and between the layer, and substrate were given as a result.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.207-216
• /
• 2023

Three time-discontinuous Galerkin quadrature element methods (TDGQEMs) are developed for structural dynamic problems. The weak-form time-discontinuous Galerkin (TDG) statements, which are capable of capturing possible displacement and/or velocity discontinuities, are employed to formulate the three types of quadrature elements, i.e., single-field, single-field/least-squares and two-field. Gauss-Lobatto quadrature rule and the differential quadrature analog are used to turn the weak-form TDG statements into a system of algebraic equations. The stability, accuracy and numerical dissipation and dispersion properties of the formulated elements are examined. It is found that all the elements are unconditionally stable, the order of accuracy is equal to two times the element order minus one or two times the element order, and the high-order elements possess desired high numerical dissipation in the high-frequency domain and low numerical dissipation and dispersion in the low-frequency domain. Three fundamental numerical examples are investigated to demonstrate the effectiveness and high accuracy of the elements, as compared with the commonly used time integration schemes.

• Communications of the Korean Mathematical Society
• /
• v.38 no.3
• /
• pp.755-772
• /
• 2023

Our aim is to establish certain image formulas of the (p, q)-extended modified Bessel function of the second kind M_ν,p,q(z) by employing the Marichev-Saigo-Maeda fractional calculus (integral and differential) operators including their composition formulas and using certain integral transforms involving (p, q)-extended modified Bessel function of the second kind M_ν,p,q(z). Corresponding assertions for the Saigo's, Riemann-Liouville (R-L) and Erdélyi-Kober (E-K) fractional integral and differential operators are deduced. All the results are represented in terms of the Hadamard product of the (p, q)-extended modified Bessel function of the second kind M_ν,p,q(z) and Fox-Wright function _rΨ_s(z).

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.1
• /
• pp.1-8
• /
• 1992

A method analizing contact pressure between plate and elastic half space is presented by using F.E.M. With the method, the pressure intensities at surface nodes of half space cae be directly calculated by using flexibility matrix of half space. The method is originally presented by Y.K. Cheung et al.(3) Insted of Y.K. Cheung's method, which use a conception of equi-contact pressure area around each surface nodes of half space in the noded rectanqular element area. We use the equi-contact pressure area around the Gaussian integration points of half space surface in the noded isoparametric element area. Numarical examples are presented and compared with other's studies.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.211-219
• /
• 2008

A three-dimensional inviscid flow solver has been developed based on unstructured meshes for the investigation of the effect of the external stores on the tail surfaces of a generic fighter aircraft. The numerical method is based on a vertex-centered finite-volume discretization and an implicit point Gauss-Seidel time integration. The calculations were made for a steady flow and the computed results were compared with experimental data to validate the flow solver. An unsteady time-accurate computation of the generic fighter aircraft with external stores at transonic flight conditions showed that the external stores cause unsteady loading on the horizontal tail surface due to the mutual interference between their wake and the horizontal tail surface. It was shown that downward deflection of the trailing edge flap significantly reduces the undesirable interference effect.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.3
• /
• pp.279-293
• /
• 2004

The analysis of linear static and free vibration problems of isotropic and laminated composite plates and shells is performed by the improved 9-node shell element with the new strain displacement relationship. In that relationship, the effect of new additional terms between the bending strain and displacement has been investigated in the warping problem. Natural co ordinate based strains, stresses and constitutive equations are used. The assumed natural strain method is used to alleviate both membrane and shear locking behavior from the element. The Lanczos method is employed in the calculation of the eigenvalues of laminated composite structures and the Gauss integration rule is adopted to evaluate the mass matrix. The numerical examples are compared with the analytical solutions to validate the current formulation and the results presented could be useful for the understanding of the behaviour of laminates under free vibration conditions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.463-474
• /
• 2009

Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.