• Structural Engineering and Mechanics
• /
• v.33 no.5
• /
• pp.621-633
• /
• 2009

In this paper we compare the solution of the one-dimensional beam model and the numerical solution of the two-dimensional linearized elasticity problem for rectangular domain of the beam-like form. We first derive the beam model starting from the two-dimensional linearized elasticity, the same way it is derived from the three-dimensional linearized elasticity. Then we present the numerical solution of the two-dimensional problem by finite element method. As expected the difference of two approximations becomes smaller as the thickness of the beam tends to zero. We then analyze the applicability of the one-dimensional model and verify the main properties of the beam modeling for thin beams.

• 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.

• Communications of the Korean Mathematical Society
• /
• v.15 no.1
• /
• pp.143-153
• /
• 2000

In [8], Franca and Stenberg developed several Galerkin least squares methods for the solution of the problem of linear elasticity. That work concerned itself only with the error estimates of the method. It did not address the related problem of finding effective methods for the solution of the associated linear systems. In this work, we propose the block diagonal preconditioners. The preconditioned conjugate residual method is robust in that the convergence is uniform as the parameter, v, goes to $\sfrac{1}{2}$. Computational experiments are included.

• Communications of the Korean Mathematical Society
• /
• v.11 no.3
• /
• pp.867-880
• /
• 1996

Numerical approximations to the linear elasticity are traditionally based on the finite element method. In this paper we propose a new formulation based on the cubic spline collocation method for linear elastic problem on the unit square. We present several numerical results for the eigenvalues of the matrix represented by cubic collocation method and preconditioner matrix which is preconditioned by FEM and FDM. Finally we present the numerical solution for some example equation.

• East Asian mathematical journal
• /
• v.14 no.2
• /
• pp.399-410
• /
• 1998

In [3], Franca and Stenberg developed several Galerkin least squares methods for the solution of the problem of linear elasticity. That work concerned itself only with the error estimates of the method. It did not address the related problem of finding effective methods for the solution of the associated-linear systems. In this work, we present computational experiments of W-cycle multigrid method. Computational experiments show that the convergence is uniform as the parameter, $\nu$, goes to 1/2.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.6
• /
• pp.1001-1004
• /
• 1987

In this study mathematical properties of variational solution and solution of the boundary integral equation of the linear elasticity problem are studied. It is first reviewed that a variational solution for the three-dimensional linear elasticity problem exists in the Sobolev space [ $H^{1}$(.OMEGA.)]$^{3}$ and, then, it is shown that a unique solution of the boundary integral equation is identical to the variational solution in [ $H^{1}$(.OMEGA.)]$^{3}$. To represent the boundary integral equation, the Green's formula in the Sobolev space is utilized on the solution domain excluding a ball, with small radius .rho., centered at the point where the point load is applied. By letting .rho. tend to zero, it is shown that, for the linear elasticity problem, boundary integral equation is valid for the variational solution. From this fact, one can obtain a numerical approximatiion of the variational solution by the boundary element method even when the classical solution does not exist.exist.

• Steel and Composite Structures
• /
• v.25 no.2
• /
• pp.127-139
• /
• 2017

In this paper, an exact analytical solution for simply supported sandwich plate which considers the permeation effect of adhesives is presented. The permeation layer is described as functionally graded material (FGM), the elastic modulus of which is assumed to be graded along the thickness following the exponential law. Based on the exact three-dimensional (3-D) elasticity theory, the solution of stresses and displacements for each layer is derived. By means of the recursive matrix method, the solution can be efficiently obtained for plates with many layers. The present solution obtained can be used as a benchmark to access other simplified solutions. The comparison study indicates that the finite element (FE) solution is close to the present one when the FGM layer in the FE model is divided into a series of homogeneous layers. However, the present method is more efficient than the FE method, with which the mesh division and computation are time-consuming. Moreover, the solution based on Kirchhoff-Love plate theory is greatly different from the present solution for thick plates. The influence of the thickness of the permeation layer on the stress and displacement fields of the sandwich plate is discussed in detail. It is indicated that the permeation layer can effectively relieve the discontinuity stress at the interface.

• Structural Engineering and Mechanics
• /
• v.7 no.2
• /
• pp.181-202
• /
• 1999

Dynamic response of axisymmetric arbitrary laminated composite cylindrical shell of finite length, using three-dimensional elasticity equations are studied. The shell is simply supported at both ends. The highly coupled partial differential equations are reduced to ordinary differential equations (ODE) with variable coefficients by means of trigonometric function expansion in axial direction. For cylindrical shell under dynamic load, the resulting differential equations are solved by Galerkin finite element method, In this solution, the continuity conditions between any two layer is satisfied. It is found that the difference between elasticity solution (ES) and higher order shear deformation theory (HSD) become higher for a symmetric laminations than their unsymmetric counterpart. That is due to the effect of bending-streching coupling. It is also found that due to the discontinuity of inplane stresses at the interface of the laminate, the slope of transverse normal and shear stresses aren't continuous across the interface. For free vibration analysis, through dividing each layer into thin laminas, the variable coefficients in ODE become constants and the resulting equations can be solved exactly. It is shown that the natural frequency of symmetric angle-ply are generally higher than their antisymmetric counterpart. Also the results are in good agreement with similar results found in literatures.

• 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.

• Advances in materials Research
• /
• v.11 no.1
• /
• pp.41-57
• /
• 2022

This paper is focused on delamination analysis of a multilayered inhomogeneous viscoelastic beam subjected to linear creep under constant applied stress. The viscoelastic model that is used to treat the creep consists of consecutively connected units. Each unit consists of one spring and two dashpots. The number of units in the model is arbitrary. The modulus of elasticity of the spring in each unit changes with time. Besides, the modulii of elasticity and the coefficients of viscosity change continuously along the thickness, width and length in each layer since the material is continuously inhomogeneous in each layer of the beam. A time-dependent solution to the strain energy release rate for the delamination is derived. A time-dependent solution to the J-integral is derived too. A parametric analysis of the strain energy release rate is carried-out by applying the solution derived. The influence of various factors such as creep, material inhomogeneity, the change of the modulii of elasticity with time and the number of units in the viscoelastic model on the strain energy release rate are clarified.