• Structural Engineering and Mechanics
• /
• v.28 no.4
• /
• pp.461-477
• /
• 2008

An adaptive finite element method for analyzing two-dimensional and axisymmetric nonlinear elastic fracture mechanics problems with cracks is presented. The J-integral is used as a parameter to characterize the severity of stresses and deformation near crack tips. The domain integral technique, for which all relevant quantities are integrated over any arbitrary element areas around the crack tips, is utilized as the J-integral solution scheme with 9-node degenerated crack tip elements. The solution accuracy is further improved by incorporating an error estimation procedure onto a remeshing algorithm with a solution mapping scheme to resume the analysis at a particular load level after the adaptive remeshing technique has been applied. Several benchmark problems are analyzed to evaluate the efficiency of the combined domain integral technique and the adaptive finite element method.

• Journal of Power System Engineering
• /
• v.21 no.6
• /
• pp.40-45
• /
• 2017

The objective of this study is the finite element-transfer stiffness coefficient method, which is the combination of the modeling technique of finite element method and the transfer technique of transfer stiffness coefficient method, is applied in the static analyses of two dimensional curved beam structures. To confirm the effectiveness of the applied method, two computational models are selected and analyzed by using finite element method, finite element-transfer stiffness coefficient method and exact solution. The computational results of the static analyses for two computational models using finite element-transfer stiffness coefficient method are equal to those using finite element method. When the element partition number of curved beam structure is increased, the computational results of the static analyses using both methods approach the exact solution. We confirmed that the finite element-transfer stiffness coefficient method is superior to finite element method when the number of the curved beam elements is increased from the viewpoints of the computational speed and the utility of computer memory.

• Bulletin of the Korean Mathematical Society
• /
• v.57 no.2
• /
• pp.393-406
• /
• 2020

We analyze a posteriori error estimator for the conforming P2 finite element on triangular meshes which is based on the solution of local Neumann problems. This error estimator extends the one for the conforming P1 finite element proposed in [4]. We prove that it is asymptotically exact for the Poisson equation when the underlying triangulations are mildly structured and the solution is smooth enough.

• Structural Engineering and Mechanics
• /
• v.50 no.2
• /
• pp.181-199
• /
• 2014

Locking phenomenon is a mesh problem and can be staved off with mesh refinement. If the studier is not preferred going to the solution with increasing mesh size or the computer memory can stack over flow than using higher order plate finite element or using integration techniques is a solution for this problem. The purpose of this paper is to show the shear locking phenomenon can be avoided by increase low order finite element mesh size of the plates and to study shear locking-free analysis of thick plates using Mindlin's theory by using higher order displacement shape function and to determine the effects of various parameters such as the thickness/span ratio, mesh size on the linear responses of thick plates subjected to uniformly distributed loads. A computer program using finite element method is coded in C++ to analyze the plates clamped or simply supported along all four edges. In the analysis, 4-, 8- and 17-noded quadrilateral finite elements are used. It is concluded that 17-noded finite element converges to exact results much faster than 8-noded finite element, and that it is better to use 17-noded finite element for shear-locking free analysis of plates.

• Honam Mathematical Journal
• /
• v.38 no.3
• /
• pp.661-674
• /
• 2016

In [15] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities, which is useful for the problem with known stress intensity factor. They consider the Poisson equations with homogeneous Dirichlet boundary condition, compute the finite element solution using standard FEM and use the extraction formula to compute the stress intensity factor, then they pose a PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor, which converges with optimal speed. From the solution we could get accurate solution just by adding the singular part. This approach works for the case when we have the accurate stress intensity factor. In this paper we consider Poisson equations with mixed boundary conditions and show the method depends the accrucy of the stress intensity factor by considering two algorithms.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.3
• /
• pp.11-20
• /
• 1986

The purpose of this study is to devise an efficient and economic solution algorithm for the nonlinear finite element equations. First, procedures and characteristics of the solution methods of ordinary nonlinear equations are critically reviewed and discussed. Based on the discussion, some promising nonlinear finite element analysis procedures are presented as an algorithmic form. Finally, a conceptually combined algorithm for a solution of nonlinear finite element equations is proposed and analyzed, in which the computational effort is minimized and numerical difficulties can be avoided.

• Honam Mathematical Journal
• /
• v.40 no.4
• /
• pp.785-794
• /
• 2018

In [8] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities. They consider the Poisson equations with homogeneous Dirichlet boundary condition with one corner singularity at the origin, and compute the finite element solution using standard FEM and use the extraction formula to compute the stress intensity factor, then pose a PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor, which converges with optimal speed. From the solution they could get an accurate solution just by adding the singular part. This approach uses the polar coordinate and the cut-off function to control the singularity and the boundary condition. In this paper we consider Poisson equations with multiple singular points, which involves different cut-off functions which might overlaps together and shows the way of cording in FreeFEM++ to control the singular functions and cut-off functions with numerical experiments.

• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.281-295
• /
• 2005

An analytical method is presented to solve the elastodynamic problem of finitely long hollow cylinder subjected to torsional impact often occurs in engineering mechanics. The analytical solution is composed of a solution of quasi-static equation satisfied with the non-homogeneous boundary condition and a solution of dynamic equation satisfied with homogeneous boundary condition. The quasi-static solution is obtained directly by solving the quasi-static equation satisfied with the non-homogeneous boundary condition. The solution of the non-homogeneous dynamic equation is obtained by means of finite Hankel transform on the radial variable, r, Laplace transform on time variable, t, and finite Fourier transform on axial variable, z. Thus, the solution for finitely long, hollow cylinder subjected to torsion impact is obtained. In the calculating examples, the response histories and distributions of shear stress in the finitely long hollow cylinder subjected to an exponential decay torsion load are obtained, and the results have been analyzed and discussed. Finally, a dynamic finite element for the same problem is carried out by using ABAQUS finite element analysis. Comparing the analytical solution with the finite element solution, it can be found that two kinds of results obtained by means of two different methods agree well. Therefore, it is further concluded that the analytical method and computing process presented in the paper are effective and accurate.

• Coupled systems mechanics
• /
• v.9 no.1
• /
• pp.5-28
• /
• 2020

Despite a widely-held belief that the finite element method is the method for the solution of solid mechanics problems, which has for 30 years dissuaded solid mechanics scientists from paying any attention to the finite volume method, it is argued that finite volume methods can be a viable alternative. It is shown that it is simple to understand and implement, strongly conservative, memory efficient, and directly applicable to nonlinear problems. A number of examples are presented and, when available, comparison with finite element methods is made, showing that finite volume methods can be not only equal to, but outperform finite element methods for many applications.

• Structural Engineering and Mechanics
• /
• v.69 no.6
• /
• pp.615-626
• /
• 2019

The 3-noded metric Timoshenko beam element with an offset of the internal node from the element centre is used here to demonstrate the best-fit paradigm using function space formulation under locking and mesh distortion. The best-fit paradigm follows from the projection theorem describing finite element analysis which shows that the stresses computed by the displacement finite element procedure are the best approximation of the true stresses at an element level as well as global level. In this paper, closed form best-fit solutions are arrived for the 3-noded Timoshenko beam element through function space formulation by combining field consistency requirements and distortion effects for the element modelled in metric Cartesian coordinates. It is demonstrated through projection theorems how lock-free best-fit solutions are arrived even under mesh distortion by using a consistent definition for the shear strain field. It is shown how the field consistency enforced finite element solution differ from the best-fit solution by an extraneous response resulting from an additional spurious force vector. However, it can be observed that when the extraneous forces vanish fortuitously, the field consistent solution coincides with the best-fit strain solution.