• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권2호
• /
• pp.115-138
• /
• 2019

The dual singular function method [DSFM] is a solver for corner sigulaity problem. We already construct DSFM in previous reserch to solve the Stokes equations including one singulairity at each reentrant corner, but we find out a crucial incorrection in the proof of well-posedness and regularity of dual singular function. The goal of this paper is to prove accuracy and well-posdness of DSFM for Stokes equations including two singulairities at each corner. We also introduce new applicable algorithms to slove multi-singulrarity problems in a complicated domain.

• Structural Engineering and Mechanics
• /
• 제54권4호
• /
• pp.607-622
• /
• 2015

This paper presents a comparative study of analytical method and finite element method (FEM) for analysis of a continuous contact problem. The problem consists of two elastic layers loaded by means of a rigid circular punch and resting on semi-infinite plane. It is assumed that all surfaces are frictionless and only compressive normal tractions can be transmitted through the contact areas. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. Then, finite element model of the problem is constituted using ANSYS software and the two dimensional analysis of the problem is carried out. The contact stresses under rigid circular punch, the contact areas, normal stresses along the axis of symmetry are obtained for both solutions. The results show that contact stresses and the normal stresses obtained from finite element method (FEM) provide boundary conditions of the problem as well as analytical results. Also, the contact areas obtained from finite element method are very close to results obtained from analytical method; disagree by 0.03-1.61%. Finally, it can be said that there is a good agreement between two methods.

• Structural Engineering and Mechanics
• /
• 제64권4호
• /
• pp.505-513
• /
• 2017

In this paper, a receding contact problem for an elastic layer resting on a half plane is considered. The layer is pressed by two rectangular stamps placed symmetrically. It is assumed that the contact surfaces are frictionless and only compressive traction can be transmitted through the contact surfaces. In addition the effect of body forces is neglected. Firstly, the problem is solved analytically based on theory of elasticity. In this solution, the problem is reduced into a system of singular integral equations in which half contact length and contact pressures are unknowns using boundary conditions and integral transform techniques. This system is solved numerically using Gauss-Jacobi integral formulation. Secondly, two dimensional finite element analysis of the problem is carried out using ANSYS. The dimensionless quantities for the contact length and the contact pressures are calculated under various stamp size, stamp position and material properties using both solutions. The analytic results are verified by comparison with finite element results.

• 대한기계학회논문집
• /
• 제17권10호
• /
• pp.2398-2406
• /
• 1993

Abdi's numerical method(ref.13) for representing a stress singularity by shifting the mid-side nodes of isoparametric elements is reviewed. A simple technique to obtain the optimal position of the mid-side nodes in quadratic isoparametric finite element is presented. From this technique we can directly obtain the position of the side-nodes adjacent to the crack tip. It is also observed that the present technique provides good accuracy for the expression of the opening displacement and the determination of the mid-side nodes for more wide range of material properties than that obtained by Abdicant the finite element method is applied to determine stress intensity factors for pressurized crack perpendicular to and terminating at the interface of two bonded dissimilar materials. A proper definition for stress intensity factors of a crack perpendicular to bimaterial interface is provided. It is based upon a near-tip displacement solutions on the crack surface for interface crack between two dissimilar materials. Numerical testing is carried out with the eight-node and six-node elements. The results obtained are compared with the previous solutions.

• Structural Engineering and Mechanics
• /
• 제85권5호
• /
• pp.621-633
• /
• 2023

The major objective of this paper is to study the receding contact problem between two functional graded layers under a flat indenter. The gravity is assumed negligible, and the shear moduli of both layers are assumed to vary exponentially along the thickness direction. In the absence of body forces, the problem is reduced to a system of Fredholm singular integral equations with the contact pressure and contact size as unknowns via Fourier integral transform, which is transformed into an algebraic one by the Gauss-Chebyshev quadratures and polynomials of both the first and second kinds. Then, an iterative speediest descending algorithm is proposed to numerically solve the system of algebraic equations. Both semi-analytical and finite element method, FEM solutions for the presented problem validate each other. To improve the accuracy of the numerical result of FEM, a graded FEM solution is performed to simulate the FGM mechanical characteristics. The results reveal the potential links between the contact stress/size and the indenter size, the thickness, as well as some other material properties of FGM.

• Structural Engineering and Mechanics
• /
• 제54권1호
• /
• pp.69-85
• /
• 2015

In this paper, a receding contact problem for an elastic layer resting on two quarter planes is considered. The layer is pressed by a stamp and distributed loads. It is assumed that the contact surfaces are frictionless and only compressive traction can be transmitted through the contact surfaces. In addition the effect of body forces are neglected. Firstly, the problem is solved analytically based on theory of elasticity. In this solution, the problem is reduced into a system of singular integral equations in which contact areas and contact stresses are unknowns using boundary conditions and integral transform techniques. This system is solved numerically using Gauss-Jacobi integral formulation. Secondly, two dimensional finite element analysis of the problem is carried out using ANSYS. The dimensionless quantities for the contact areas and the contact pressures are calculated under various distributed load conditions using both solutions. It is concluded that the position and the magnitude of the distributed load have an important role on the contact area and contact pressure distribution between layer and quarter plane contact surface. The analytic results are verified by comparison with finite element results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제22권2호
• /
• pp.101-113
• /
• 2018

The dual singular function method(DSFM) is a numerical algorithm to get optimal solution including corner singularities for Poisson and Helmholtz equations. In this paper, we apply DSFM to solve heat equation which is a time dependent problem. Since the DSFM for heat equation is based on DSFM for Helmholtz equation, it also need to use Sherman-Morrison formula. This formula requires linear solver n + 1 times for elliptic problems on a domain including n reentrant corners. However, the DSFM for heat equation needs to pay only linear solver once per each time iteration to standard numerical method and perform optimal numerical accuracy for corner singularity problems. Because the Sherman-Morrison formula is rather complicated to apply computation, we introduce a simplified formula by reanalyzing the Sherman-Morrison method.

• Structural Engineering and Mechanics
• /
• 제31권4호
• /
• pp.439-451
• /
• 2009

The paper deals with the problem related to the modelling of riveted assemblies for crashworthiness analysis of full-scale complete aircraft structures. Comparisons between experiments and standard FE computations on high-energy accidental situations onto aluminium riveted panels show that macroscopic plastic strains are not sufficiently localised in the FE shells connected to rivet elements. The main reason is related to the structural embrittlement caused by holes, which are currently not modelled. Consequently, standard displacement FE models do not succeed in initialising and propagating the rupture in sheet metal plates and along rivet rows as observed in the experiments. However, the literature survey show that it is possible to formulate super-elements featuring defects that both give accurate singular strain fields and are compatible with standard displacement finite elements. These super-elements can be related to the displacement model of the hybrid-Trefftz principle of the finite element method, which is a kind of domain decomposition method. A feature of hybrid-Trefftz finite elements is that they are mainly used for elastic computations. It is thus proposed to investigate the possibility of formulating a hybrid displacement finite element, including the effects of a hole, dedicated to crashworthiness analysis of full-scale aeronautic structures.

• 연구논문집
• /
• 통권28호
• /
• pp.123-135
• /
• 1998

A shape function for element free Galerkin method is carved from Shepard interpolant of singular weight and consistency condition. Thus present shape function is an interpolation and has no singularities. The shape function is applied to cantilever bending problems and gives good results in comparison with beam theory. Finally it is shown that the coupling with finite element method is made easily without any additional treaties.

• 한국항공우주학회지
• /
• 제34권8호
• /
• pp.41-46
• /
• 2006

충격응답함수와 조정법(regularization methods)을 이용하여 항공기 날개의 충격하중 복원 가능성을 검토하였다. 충격하중에 대한 구조의 응답을 내타낼 수 있는 충격응답함수를 날개 유한요소모델의 강성과 질량 자료로 유도하였다. 일반적으로 부적합(ill-posed) 특성을 지닌 충격응답함수의 역행렬은 반복 Tikhonov 조종법(Iterative Tikhonov Regularization Method)과 일반화 특이치 분해법(Generalized Singular Value Decomposition Method)을 사용하여 구하였다. 수치적 입증을 위하여 전투기급 주익을 사용하였다. 해당 주익의 유한요소해석을 통하여 임의의 충격하중에 대한 변위와 변형률을 계산하였으며, 이를 충격응답함수로 계산한 결과와 비교하였다. 또한, 유한요소해석에서 계산된 변형률을 사용하여 충격하중을 복원하였다. 수치적 입증 결과 항공기 구조의 충격하중 모니터링이 본 방법으로 가능할 수 있음을 보여주었다.