• Structural Engineering and Mechanics
• /
• 제6권7호
• /
• pp.727-746
• /
• 1998

Finite element models and numerical results are presented for bending and natural vibration using the unified third-order plate theory developed in Part 1 of this paper. The unified third-order theory contains the classical, first-order, and other third-order plate theories as special cases. Analytical solutions are developed using the Navier and L$\acute{e}$vy solution procedures (see Part 1 of the paper). Displacement finite element models of the unified third-order theory are developed herein. The finite element models are based on $C^0$ interpolation of the inplane displacements and rotation functions and $C^1$ interpolation of the transverse deflection. Numerical results of bending and natural vibration are presented to evaluate the accuracy of various plate theories.

• Computers and Concrete
• /
• 제3권6호
• /
• pp.423-437
• /
• 2006

Stress wave propagation through concrete is simulated by finite element analysis. The concrete medium is modeled as a homogeneous material with smeared properties to investigate and establish the suitable finite element analysis method (explicit versus implicit) and analysis parameters (element size, and solution time increment) also suitable for rigorous investigation. In the next step, finite element analysis model of the medium is developed using a digital image processing technique, which distinguishes the mortar and aggregate phases of concrete. The mortar and aggregate phase topologies are, then, directly mapped to the finite element mesh to form a heterogeneous concrete model. The heterogeneous concrete model is then used to simulate wave propagation. The veracity of the model is demonstrated by evaluating the intrinsic parameters of nondestructive ultrasonic pulse velocity testing of concrete. Quantitative relationships between aggregate size and testing frequency for nondestructive testing are presented.

• 소성∙가공
• /
• 제9권4호
• /
• pp.379-388
• /
• 2000

This study aims to examine the influence of experimental and numerical factors on the results of the test and finite element simulation to evaluate the formability of sheet metals. The stretch-forming test with a hemispherical punch is carried out to obtain the limiting dome height (LDH) and forming limit diagram (FLD) for several kinds of aluminium and steel sheet. The results of the LDH and FLD tests are analysed to find any correlation with the uniaxial tensile properties. It proves that the size of the prescribed grid has great influence on the measured value of strain. The finite element analysis of the punch stretching process is also carried out and the result is compared with the experimental data. The influence of the numerical parameters such as friction coefficient, element size and anisotropy model on the simulation results tms out to be very considerable.

• Structural Engineering and Mechanics
• /
• 제17권2호
• /
• pp.203-214
• /
• 2004

The dynamic response characteristics of a rectangular fluid container are investigated by using finite element method. The fluid is assumed to be linear-elastic, inviscid and compressible. A displacement-based fluid finite element was employed to allow for the effects of the fluid. A typical rectangular fluid container, which is used in recent studies, is considered for the numerical analysis. The North-South component of El Centro Earthquake records is used as input ground acceleration. Rigid and flexible fluid containers solutions are obtained for the chosen sample tank. Hydrodynamic pressures and sloshing motions are determined using Lagrangian fluid finite element. The results obtained from this study are compared with the results obtained by boundary-finite element method (BEM-FEM) and requirements of Eurocode-8. Based on the numerical analysis, some conclusions and discussions on the design considerations for rectangular fluid containers are presented.

• 소성∙가공
• /
• 제17권1호
• /
• pp.35-45
• /
• 2008

Electromagnetic Forming (EMF) technology such as magnetic pulse forming, which is one of the high velocity forming methods, has been used for the joining and forming process in various industry fields. This method could be derived a series of deformation of sheet metal by using a strong magnetic field. In this study, numerical approach by finite element simulation of the electromagnetic forming process was presented. A transient electromagnetic finite element code was used to obtain the numerical model of the time-varying currents that are discharged through the coil in order to obtain the transient magnetic forces. Also, the body forces generated in electromagnetic field were used as the loading condition to analyze deformation of thin sheet metal workpiece using explicit dynamic finite element code. In this study, after finite element analysis for thin sheet metal forming process with free surface configuration was performed, analytical approach for a dimpled shape by using EMF was carried out. Furthermore, the simulated results of the dimpled shape by EMF were compared with that by a conventional solid tool in view of the deformed shape. From the results of finite element analysis, it is confirmed that the EMF process could be applied to thin sheet metal forming.

• The Journal of Korean Physical Therapy
• /
• 제22권6호
• /
• pp.77-83
• /
• 2010

Purpose: The aim of this study was to do numerical analysis of the wavelength dependence in low level laser therapy (LLLT) using a finite element method (FEM). Methods: Numerical analysis of heat transfer based on a Pennes' bioheat equation was performed to assess the wavelength dependence of effects of LLLT in a single layer and in multilayered tissue that consists of skin, fat and muscle. The three different wavelengths selected, 660 nm, 830 nm and 980 nm, were ones that are frequently used in clinic settings for the therapy of musculoskeletal disorders. Laser parameters were set to the power density of 35.7 W/$cm^2$, a spot diameter of 0.06 cm, and a laser exposure time of 50 seconds for all wavelengths. Results: Temperature changes in tissue based on a heat transfer equation using a finite element method were simulated and were dominantly dependent upon the absorption coefficient of each tissue layer. In the analysis of a single tissue layer, heat generation by fixed laser exposure at each wavelength had a similar pattern for increasing temperature in both skin and fat (980 nm > 660 nm > 830 nm), but in the muscle layer 660nm generated the most heat (660 nm ${\gg}$ 980 nm > 830 nm). The heat generation in multilayered tissue versus penetration depth was shown that the temperature of 660 nm wavelength was higher than those of 830 nm and 980 nm Conclusion: Numerical analysis of heat transfer versus penetration depth using a finite element method showed that the greatest amount of heat generation is seen in multilayered tissue at = 660 nm. Numerical analysis of heat transfer may help lend insight into thermal events occurring inside tissue layers during low level laser therapy.

• East Asian mathematical journal
• /
• 제33권5호
• /
• pp.495-509
• /
• 2017

With the accuracy of the nonlinearity guaranteed, plenty of time and large memory space are needed when we solve the finite element numerical solution of nonlinear partial differential equations. In this paper, we use the Group Element Method (GEM) to deal with the non-linearity of the BBM-Burgers Equation with Conservation form and perform a numerical analysis for two particular initial-boundary value (the Dirichlet boundary conditions and Neumann-Dirichlet boundary conditions) problems with the Finite Element Method (FEM). Some numerical experiments are performed to analyze the error between the exact solution and the FEM solution in MATLAB.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집A
• /
• pp.324-329
• /
• 2001

A novel method for non-matching interfaces on the boundaries of the finite elements in partitioned domains is presented by introducing interface elements in this paper. The interface element method (IEM) satisfies the continuity conditions exactly through interfaces without recourse to the Lagrange multiplier technique. The moving least square (MLS) approximation in the present study is implemented to construct the shape functions of the interface elements. Alignment of the boundaries of sub-domains in the MLS approximation and integration domains provides a consistent numerical integration due to one form of rational functions in an integration domain. The compatibility of displacements on the boundaries of the finite elements and the interface elements is always preserved in this method, and the completeness of the shape functions of the interface elements guarantees the convergence of numerical solutions. The numerical examples show that the interface element method is a useful tool for the analysis of a partitioned system and for a global-local analysis.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2003년도 봄 학술발표회 논문집
• /
• pp.463-468
• /
• 2003

A finite element formulation to predict the flexural behavior of composite girder is presented in which the early-age properties of concrete are specified including maturing of elastic modulus, creep and shrinkage. The time dependent constitutive relation accounting for the early-age concrete properties is derived in an incremental format by expanding the total form of stress-strain relation by the first order Taylor series with respect to the reference time. The total potential energy of the flexural composite member is minimized to derive the time dependent finite element equilibrium equation. Numerical applications are made for the 3-span double composite steel box girders which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The numerical analysis with considering the variation of concrete elastic modulus are performed to investigate the effect of it on the early-age behavior of composite structures. The one dimensional finite element analysis results are compared with the analytical method based on the sectional analysis. Close agreement is observed among the two methods.

• Structural Engineering and Mechanics
• /
• 제40권2호
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.