• Structural Engineering and Mechanics
• /
• 제61권6호
• /
• pp.765-773
• /
• 2017

The quadratic B-spline finite element method yields mass and stiffness matrices which are half the size of matrices obtained by the conventional finite element method. We solve the free vibration problem of a rotating Rayleigh beam using the quadratic B-spline finite element method. Rayleigh beam theory includes the rotary inertia effects in addition to the Euler-Bernoulli theory assumptions and presents a good mathematical model for rotating beams. Galerkin's approach is used to obtain the weak form which yields a system of symmetric matrices. Results obtained for the natural frequencies at different rotating speeds show an accurate match with the published results. A comparison with Euler-Bernoulli beam is done to decipher the variations in higher modes of the Rayleigh beam due to the slenderness ratio. The results are obtained for different values of non-uniform parameter ($\bar{n}$).

• East Asian mathematical journal
• /
• 제29권5호
• /
• pp.503-510
• /
• 2013

A quadratic B-spline finite element method for the spatial variable combined with a Newton method for the time variable is proposed to approximate a solution of Benjamin-Bona-Mahony-Burgers (BBMB) equation. Two examples were considered to show the efficiency of the proposed scheme. The numerical solutions obtained for various viscosity were compared with the exact solutions. The numerical results show that the scheme is efficient and feasible.

• 대한수학회보
• /
• 제57권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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제27권4호
• /
• pp.245-271
• /
• 2023

In this paper we discuss some recovery of H(div)-conforming flux approximations from the equilibrated fluxes of Ainsworth and Oden for quadratic finite element methods of second-order elliptic problems. Combined with the hypercircle method of Prager and Synge, these flux approximations lead to a posteriori error estimators which provide guaranteed upper bounds on the numerical error. Furthermore, we prove some superconvergence results for the flux approximations and asymptotic exactness for the error estimator under proper conditions on the triangulation and the exact solution. The results extend those of the previous paper for linear finite element methods.

• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제18권3호
• /
• pp.269-277
• /
• 2014

A numerical scheme is proposed to control the BBMB (Benjamin-Bona-Mahony-Burgers) equation, and the scheme consists of three steps. Firstly, BBMB equation is converted to a finite set of nonlinear ordinary differential equations by the quadratic B-spline finite element method in spatial. Secondly, the controller is designed based on the linear quadratic regulator (LQR) theory; Finally, the system of the closed loop compensator obtained on the basis of the previous two steps is solved by the backward Euler method. The controlled numerical solutions are obtained for various values of parameters and different initial conditions. Numerical simulations show that the scheme is efficient and feasible.

• 한국정밀공학회지
• /
• 제13권7호
• /
• pp.66-77
• /
• 1996

Optimal topology design is to search the optimal configuration of a structure which can be used as a shape at the conceptual design stage. Our objective is to maximize the stiffness of the structures and ribs under a material usage constraintl. The density of each finite element is the design variable and its relationship with Young's modulus is expressed by quadratic form. The configuration is represented by the entire density distribution, the structural analysis is performed by finite element method and the optimiza- tion is performed by Feasible Direction Method. Feasible Direction Method can handle various problems simultaneously, that is, mult-objectives and multi-constraints. Total computation time can be reduced by the quadratic relationship between the density and the material property and fewer design variables than Homogenization Method. Toplogy optimization technique developed in this research is applied to design the shapes of the ribs.

• Journal of applied mathematics & informatics
• /
• 제27권1_2호
• /
• pp.301-313
• /
• 2009

This paper presents a Finite Element Method (FEM) application to thermal study of natural three layers of human dermal parts of varying properties. This paper carries out investigation of temperature distributions in these layers namely epidermis, dermis and under lying tissue layer. It is assumed that the outer skin is exposed to atmosphere and the loss of heat due to convection, radiation and evaporation of water have also been taken into account. The computations are carried out for one dimensional unsteady state case and the shape functions in dermal parts have been considered to be quadratic. A Finite Element scheme that uses the Crank-Nicolson method is used to solve the problem and the results computed have been exhibited graphically.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1988년도 전기.전자공학 학술대회 논문집
• /
• pp.87-90
• /
• 1988

This paper reports a simple posteriori error estimate method for adaptive finite element mesh generation using quadratic shape function especially for the magnetic field problems. The elements of quadratic shape function have more precise solution than those of linear shape function. Therefore, the difference of two solutions gives error quantity. The method uses the magnetic flux density error as a basis for refinement. This estimator is tested on two dimensional problem which has singular points. The estimated error is always under estimated but in same order as exact error, and this method is much simpler and more convenient than other methods. The result shows that the adaptive mesh gives even better rate of convergence in global error than the uniform mesh.

• Nuclear Engineering and Technology
• /
• 제49권1호
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.