• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.1 no.1
• /
• pp.105-125
• /
• 1997

In this paper we consider multigrid algorithms for solving elasticity problems by using Wilson's nonconforming finite element method. We consider two types of intergrid transfer operators which is needed to define the multigrid algorithm and prove convergence of $\mathcal{W}$-cycle mutigrid algorithm and uniform condition number estimates for the variable $\mathcal{V}$-cycle multigrid preconditioner.

• Journal of Ship and Ocean Technology
• /
• v.8 no.3
• /
• pp.50-60
• /
• 2004

In this paper, a continuum-based design sensitivity analysis (DSA) method is developed for the weakly coupled thermo-elasticity problems. The temperature and displacement fields are described in a common domain. Boundary value problems such as an equilibrium equation and a heat conduction equation in steady state are considered. The direct differentiation method of continuum-based DSA is employed to enhance the efficiency and accuracy of sensitivity computation. We derive design sensitivity expressions with respect to thermal conductivity in heat conduction problem and Young's modulus in equilibrium equation. The sensitivities are evaluated using the finite element method. The obtained analytical sensitivities are compared with the finite differencing to yield very accurate results. Extensive developments of this method are useful and applicable for the optimal design problems incorporating welding and thermal deformation problems.

• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.761-781
• /
• 2014

An exponentially accurate nonconforming spectral element method for elasticity systems with discontinuities in the coefficients and the flux across the interface is proposed in this paper. The method is least-squares spectral element method. The jump in the flux across the interface is incorporated (in appropriate Sobolev norm) in the functional to be minimized. The interface is resolved exactly using blending elements. The solution is obtained by the preconditioned conjugate gradient method. The numerical solution for different examples with discontinuous coefficients and non-homogeneous jump in the flux across the interface are presented to show the efficiency of the proposed method.

• Structural Engineering and Mechanics
• /
• v.26 no.5
• /
• pp.591-615
• /
• 2007

An improvement is introduced to solve the plane problems of linear elasticity by reciprocal theorem for orthotropic materials. This method gives an integral equation with complex kernels which will be solved numerically. An artificial boundary is defined to eliminate the singularities and also an algorithm is introduced to calculate multi-valued complex functions which belonged to the kernels of the integral equation. The chosen sample problem is a plate, having a circular or elliptical hole, stretched by the forces parallel to one of the principal directions of the material. Results are compatible with the solutions given by Lekhnitskii for an infinite plane. Five different orthotropic materials are considered. Stress distributions have been calculated inside and on the boundary. There is no boundary layer effect. For comparison, some sample problems are also solved by finite element method and to check the accuracy of the presented method, two sample problems are also solved for infinite plate.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.347-354
• /
• 2004

Using an efficient adjoint variable method, we develop a unified design sensitivity analysis (DSA) method considering both steady state nonlinear heat conduction and geometrical nonlinear elasticity problems. Design sensitivity expressions with respect to thermal conductivity and Young's modulus are derived. Beside the temperature and displacement adjoint equations, another coupled one is defined regarding the obtained adjoint displacement field as the adjoint load in temperature field. The developed DSA method is shown to be very efficient and further extended to a topology design optimization method for the nonlinear weakly coupled thermo-elasticity problems using a density approach.

• Interaction and multiscale mechanics
• /
• v.1 no.3
• /
• pp.329-337
• /
• 2008

This paper presents the derivation of the generalized governing equations in theory of elasticity, their weak forms and the some applications in the numerical analysis of structural mechanics. Unlike the differential equations in classical elasticity theory, the generalized equations of the equilibrium and compatibility equations presented here take the form of integral equations, and the generalized equilibrium equations contain the classical differential equations and the boundary conditions in a single equation. By using appropriate test functions, the weak forms of these generalized governing equations can be established. It can be shown that various variational principles in structural analysis are merely the special cases of these weak forms of generalized governing equations in elasticity. The present weak forms of elasticity equations extend greatly the choices of the trial functions for approximate solutions in the numerical analysis of various engineering problems. Therefore, the weak forms of generalized governing equations in elasticity provide a powerful modeling tool in the computational structural mechanics.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.3
• /
• pp.187-195
• /
• 2024

This paper presents two-dimensional boundary value problems of the stress-based gradient elasticity within the smoothed finite element method (S-FEM) framework. Gradient elasticity is introduced to address the limitations of classical elasticity, particularly its struggle to capture size-dependent mechanical behavior at the micro/nano scale. The Ru-Aifantis theorem is employed to overcome the challenges of high-order differential equations in gradient elasticity. This theorem effectively splits the original equation into two solvable second-order differential equations, enabling its incorporation into the S-FEM framework. The present method utilizes a staggered scheme to solve the boundary value problems. This approach efficiently separates the calculation of the local displacement field (obtained over each smoothing domain) from the non-local stress field (computed element-wise). A series of numerical tests are conducted to investigate the influence of the internal length scale, a key parameter in gradient elasticity. The results demonstrate the effectiveness of the proposed approach in smoothing stress concentrations typically observed at crack tips and dislocation lines.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.16-23
• /
• 1997

The Reproducing Kernel Particle Method (RKPM), one of the popular meshless methods, is developed and applied to stress concentration problems. Since the meshless methods require only a set of particles (or nodes) and the description of boundaries in their formulation, the adaptivity can be implemented with much more ease than finite element method. In addition, due to its intrinsic property of multiresolution, the shape function of RKPM provides us a new criterion for adaptivity. Recently, this multiple scale Reproducing Kernel Particle Method and its adaptive procedure have been formulated for large deformation problems by the authors. They are also under development for damage materials and localization problems. In this paper the multiple scale RKPM for linear elasticity is presented and the adaptive procedure is applied to stress concentration problems. Therefore, this work may be regarded as the edition of linear elasticity in the complete framework of multiple scale RKPM and the associated adaptivity.

• Korean Journal of Mathematics
• /
• v.7 no.2
• /
• pp.245-269
• /
• 1999

Multigrid method and acceleration by conjugate gradient method for first-order system least squares (FOSLS) using bilinear finite elements are developed for various boundary value problems of planar linear elasticity. They are two-stage algorithms that first solve for the displacement flux variable, then for the displacement itself. This paper focuses on solving for the displacement flux variable only. Numerical results show that the convergence is uniform even as the material becomes nearly incompressible. Computations for convergence factors and discretization errors are included. Heuristic arguments to improve the convergences are discussed as well.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.1
• /
• pp.93-102
• /
• 2020

The concept of temperature distribution in inhomogeneous semi-infinite solids is examined by making use of direct integration method. The analysis is done on the solution of the in-plane steady state heat conduction problem under certain boundary conditions. The method of direct integration has been employed, which is then reduced to Volterra integral equation of second kind, produces the explicit form analytical solution. Using resolvent- kernel algorithm, the governing equation is solved to get present solution. The temperature distribution obtained and calculated numerically and the relation with distribution of heat flux generated by internal heat source is shown graphically.