• Steel and Composite Structures
• /
• v.39 no.5
• /
• pp.493-509
• /
• 2021

There is not enough mixed finite element method (MFEM) model developed for static and dynamic analysis of functionally graded material (FGM) beams in the literature. The main purpose of this study is to develop a reliable and efficient computational modeling using an efficient functional in MFEM for free vibration and static analysis of FGM composite beams subject to high order shear deformation effects. The modeling of material properties was performed using mixture rule and Mori-Tanaka scheme which are more realistic determination techniques. This method based on the assumption that a two phase composite material consisting of matrix reinforced by spherical particles, randomly distributed in the beam. To explain the displacement components of the shear deformation effects, it was accepted that the shear deformation effects change sinusoidal. Partial differential field equations were obtained with the help of variational methods and then these equations were transformed into a novel functional for FGM beams with the help of Gateaux differential derivative operator. Thanks to the Gateaux differential method, the compatibility of the field equations was checked, and the field equations and boundary conditions were reflected to the function. A MFEM model was developed with a total of 10 degrees of freedom to apply the obtained functional. In the numerical applications section, free vibration and flexure problems solutions of FGM composite beams were compared with those predicted by other theories to show the effects of shear deformation, thickness changing and boundary conditions.

• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Journal of Mechanical Science and Technology
• /
• v.18 no.9
• /
• pp.1529-1536
• /
• 2004

A structure is broken down into a number of substructures by means of the finite element method and the substructures are synthesized for the complete structure. The divided substructures take two types: fixed-free and free-free elements. The flexibility and stiffness matrices of the free-free elements are the Moore-Penrose inverse of each other. Thus, it is not easy to determine the equilibrium equations of the complete structure composed of two mixed types of substructures. This study provides the general form of equilibrium equation of the entire structure through the process of assembling the equilibrium equations of substructures with end conditions of mixed types. Applications demonstrate that the proposed method is effective in the structural analysis of geometrically complicated structures.

• Korean Journal of Computational Design and Engineering
• /
• v.1 no.1
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

• Steel and Composite Structures
• /
• v.27 no.1
• /
• pp.27-33
• /
• 2018

In this study, a mixed-interpolated tensorial component 4 nodes method (MITC4) is treated as a numerical analysis model for topology optimization using multiple materials assigned within Reissner-Mindlin plates. Multi-material optimal topology and shape are produced as alternative plate retrofit designs to provide reasonable material assignments based on stress distributions. Element density distribution contours of mixing multiple material densities are linked to Solid Isotropic Material with Penalization (SIMP) as a design model. Mathematical formulation of multi-material topology optimization problem solving minimum compliance is an alternating active-phase algorithm with the Gauss-Seidel version as an optimization model of optimality criteria. Numerical examples illustrate the reliability and accuracy of the present design method for multi-material topology optimization with Reissner-Mindlin plates using MITC4 elements and steel materials.

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

In this paper we consider the nonlinear parabolic problems with mixed boundary condition. Under comparatively mild conditions of the coefficients related to the problem, we construct the discontinuous Galerkin approximation of the solution to the nonlinear parabolic problem. We discretize spatial variables and construct the finite element spaces consisting of discontinuous piecewise polynomials of which the semidiscrete approximations are composed. We present the proof of the convergence of the semidiscrete approximations in $L^{\infty}(H^1)$ and $L^{\infty}(L^2)$ normed spaces.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.1-9
• /
• 2002

There is the high possibility of crack initiation from mechanical joints, which are widely used in aircraft fuselages, due to the development of stress concentration and contact pressure. In this paper, the mixed-mode stress intensity factors at the surface and deepest points of an inclined quarter elliptical corner crack in mechanical joints are analyzed by the weight function method. The coefficients included in the weight function are obtained by finite element analyses for reference loadings. Critical angle at which mode I stress intensity factor becomes maximum is determined by analyzing the variation of stress intensity factors along incline angle of crack and the effects of the amount of clearance and crack depth on the critical angle are investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.389-398
• /
• 2002

The mixed mode problem (I and II) of a peny-shaped interface cracks in remote tension loading on a bi-material cylinder is studied using finite element method. The energy release rates for the tip of the crack in the interface were calibrated for several different moduli combinations and crack ratios using the modified crack closure integral technique and J-integral method, with numerical results obtained from a commercial finite element program. Numerical results show that non-dimensional value of$\sqrt{G_{II}E^*}/\sqrt[p]{\pi a}$ increases as the crack size or moduli ratio increases. Meanwhile, non-dimensional value of$\sqrt{G_{I}E^*}/\sqrt[p]{\pi a}$ decreases as the moduli ratio increases, but above the moduli ratio of 3 its value decreases then increases again as the crack size increases. Reliability of the numerical analysis in this study was acquired with comparison to an analytical solution for the peny-shaped interface crack in an infinite medium.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.3
• /
• pp.79-84
• /
• 1998

It is very difficult to analyze the transient temperature distribution during quenching of the steel because of coupled effects among temperature, structures and stresses. In this paper, using Inoue's equation of evolution and mixture rule, transient temperature distribution is calculated by the finite element method considering latent heat of transformation structure and temperature dependence of physical and mechanical prperties for the 0.45% carbon cylindrical steel bar with 40mm diameter and 20mm height during end-quenching.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.7
• /
• pp.989-996
• /
• 2001

High-speed and low-speed flows are simulated numerically by flowfield-dependent mixed explicit-implicit (FDMEI) method. This algorithm depends on implicitness parameters of convection, diffusion, diffusion gradients, and source terms which are calculated from the changes of local Mach, Reynolds, Peclet, and Damkohler numbers between adjacent nodes. Convection phenomena or shock waves are resolved from Mach number-dependent implicitness parameters whereas diffusion or viscous actions are simulated by Reynolds number or Peclet number-dependent implicitness parameters. Fluctuation components of all variables are properly accommodated spatially and temporally in the FDMEI procedure. To illustrate, some benchmark example problems are presented for comparisons of the FDMEI results with other available data. These results appear to be encouraging and point toward the need for further investigations of the FDMEI theory.