• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.1
• /
• pp.67-77
• /
• 2000

In this paper we consider finite element mathods for nonlinear parabolic problems defined in ${\Omega}{\subset}R^d$ ($d{\leq}4$). A new initial approximation is taken. Optimal order error estimates in $L_p$ for $2{\leq}p{\leq}{\infty}$ are established for arbitrary order finite element. One order superconvergence in $W^{1,p}$ for $2{\leq}q{\leq}{\infty}$ are demonstrated as well.

• Steel and Composite Structures
• /
• v.18 no.6
• /
• pp.1353-1368
• /
• 2015

A finite element model is presented for the analysis of composite steel-concrete beams based on a refined high-order theory. The employed theory satisfies all the kinematic and stress continuity conditions at the layer interfaces and considers effects of the transverse normal stress and transverse flexibility. The global displacement components, described by polynomial or combinations of polynomial and exponential expressions, are superposed on local ones chosen based on the layerwise or discrete-layer concepts. The present finite model does not need the incorporating any shear correction factor. Moreover, in the present $C^1$-continuous finite element model, the number of unknowns is independent of the number of layers. The proposed finite element model is validated by comparing the present results with those obtained from the three-dimensional (3D) finite element analysis. In addition to correctly predicting the distribution of all stress components of the composite steel-concrete beams, the proposed finite element model is computationally economic.

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

• Structural Engineering and Mechanics
• /
• v.65 no.3
• /
• pp.213-222
• /
• 2018

The purpose of this paper is to study free vibration analysis of thick plates resting on Winkler foundation using Mindlin's theory with shear locking free fourth order finite element, to determine the effects of the thickness/span ratio, the aspect ratio, subgrade reaction modulus and the boundary conditions on the frequency paramerets of thick plates subjected to free vibration. In the analysis, finite element method is used for spatial integration. Finite element formulation of the equations of the thick plate theory is derived by using higher order displacement shape functions. A computer program using finite element method is coded in C++ to analyze the plates free, clamped or simply supported along all four edges. In the analysis, 17-noded finite element is used. Graphs are presented that should help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that 17-noded finite element can be effectively used in the free vibration analysis of thick plates. It is also concluded that, in general, the changes in the thickness/span ratio are more effective on the maximum responses considered in this study than the changes in the aspect ratio.

• East Asian mathematical journal
• /
• v.33 no.5
• /
• pp.511-525
• /
• 2017

We introduce an extrapolated higher order characteristic finite element method to construct approximate solutions of a Sobolev equation with a convection term. The higher order of convergence in both the temporal direction and the spatial direction in $L^2$ normed space is established and some computational results to support our theoretical results are presented.

• Journal of applied mathematics & informatics
• /
• v.5 no.1
• /
• pp.23-40
• /
• 1998

Mixed finite element method is developed to approxi-mate the solution of the initial-boundary value problem for a strongly nonlinear second-order hyperbolic equation in divergence form. Exis-tence and uniqueness of the approximation are proved and optimal-order $L\infty$-in-time $L^2$-in-space a priori error estimates are derived for both the scalar and vector functions approximated by the method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Journal of the Korean earth science society
• /
• v.38 no.2
• /
• pp.105-114
• /
• 2017

Two dimensional finite element method with quadrilateral basis functions was applied to the spherical high order filter on the spherical surface limited area domain. The basis function consists of four shape functions which are defined on separate four grid boxes sharing the same gridpoint. With the basis functions, the first order derivative was expressed as an algebraic equation associated with nine point stencil. As the theory depicts, the convergence rate of the error for the spherical Laplacian operator was found to be fourth order, while it was the second order for the spherical Laplacian operator. The accuracy of the new high order filter was shown to be almost the same as those of Fourier finite element high order filter. The two-dimension finite element high order filter was incorporated in the weather research and forecasting (WRF) model as a hyper viscosity. The effect of the high order filter was compared with the built-in viscosity scheme of the WRF model. It was revealed that the high order filter performed better than the built in viscosity scheme did in providing a sharper cutoff of small scale disturbances without affecting the large scale field. Simulation of the tropical cyclone track and intensity with the high order filter showed a forecast performance comparable to the built in viscosity scheme. However, the predicted amount and spatial distribution of the rainfall for the simulation with the high order filter was closer to the observed values than the case of built in viscosity scheme.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.7 no.2
• /
• pp.75-81
• /
• 2003

Compact schemes are shown to be effective for a class of problems including convection-diffusion equations when combined with multigrid algorithms [7, 8] and V-cycle convergence is proved[5]. We apply the multigrid algorithm for an semianalytic finite difference scheme, which is desinged to preserve high order accuracy despite of singularities.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.587-614
• /
• 2014

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.