• Structural Engineering and Mechanics
• /
• 제33권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
• /
• 제32권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.

• 대한수학회논문집
• /
• 제15권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.

• 대한수학회논문집
• /
• 제11권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
• /
• 제14권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.

• 대한기계학회논문집
• /
• 제11권6호
• /
• pp.1001-1004
• /
• 1987

본 연구에서는 우선 선형 탄성문제의 변분해(variational solution)가 Sobol- ev 공간[ $H^{1}$(.OMEGA.)]= $H^{1}$(.OMEGA.)* $H^{1}$(.OMEGA.)* $H^{1}$(.OMEGA.)에서 유일하게 존재함을 재 검토하고 다음으로 경계적분식의 해도 변분해와 같음을 보인다. 이것은 선형 탄 성문제의 경우 고전해(classical solution)가 존재하지 않을 경우에도 BEM을 사용하여 변분해의 수치적 근사치를 구할 수 있다는 수학적 근거가 된다. 이를 위해서 Sobol- ev 공간 내에서의 Green's formula를 적용하는데 점하중해의 특이점(singularity)때문 에 Green's formula를 적용하기가 곤란해진다. 이 문제는 적분영역 .OMEGA.를 .OMEGA.-Ｂ$_{\rho }$로 치환하고 .rho.를 0으로 접근시키는 방법으로 해결한다. 이 때 Ｂ$_{\rho}$는 특이 점에 중심을 두고 매우 작은 변경 .rho.를 갖는 구이다.ho.를 갖는 구이다.

• Steel and Composite Structures
• /
• 제25권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
• /
• 제7권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
• /
• 제24권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
• /
• 제11권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.