• Structural Engineering and Mechanics
• /
• 제53권2호
• /
• pp.205-226
• /
• 2015

Finite element method (FEM) is an effective quantitative method to solve complex engineering problems. The basic idea of FEM for a complex problem is to be able to find a solution by reducing the problem made simple. If mathematical tools are inadequate to obtain precise result, even approximate result, FEM is the only method that can be used for structural analyses. In FEM, the domain is divided into a large number of simple, small and interconnected sub-regions called finite elements. FEM has been used commonly for linear and nonlinear analyses of different types of structures to give us accurate results of plane stress and plane strain problems in civil engineering area. In this paper, FEM is used to investigate stress analysis of a shear wall which is subjected to concentrated loads and fundamental principles of stress analysis of the shear wall are presented by using matrix displacement method in this paper. This study is consisting of two parts. In the first part, the shear wall is discretized with constant strain triangular finite elements and stiffness matrix and load vector which is attained from external effects are calculated for each of finite elements using matrix displacement method. As to second part of the study, finite element analysis of the shear wall is made by ANSYS software program. Results obtained in the second part are presented with tables and graphics, also results of each part is compared with each other, so the performance of the matrix displacement method is demonstrated. The solutions obtained by using the proposed method show excellent agreements with the results of ANSYS. The results show that this method is effective and preferable for the stress analysis of shell structures. Further studies should be carried out to be able to prove the efficiency of the matrix displacement method on the solution of plane stress problems using different types of structures.

• 대한건축학회논문집:구조계
• /
• 제34권9호
• /
• pp.3-10
• /
• 2018

While computer environments have been dramatically developed in recent years, as the building structures become larger, the structural analysis models are also becoming more complex. So there is still a need to model one shear wall with one finite element. From the viewpoint of the concept of FEA, if one shear wall is modeled by one finite element, the result of analysis is not likely accurate. Shear wall may be modelled with various finite elements. Among them, considering the displacement compatibility condition with the beam element connected to the shear wall, plane stress element with in-plane rotational stiffness is preferred. Therefore, in order to analyze one shear wall with one finite element accurately, it is necessary to evaluate finite elements developed for the shear wall analysis and to develop various plane stress elements with rotational stiffness continuously. According to the above mentioned need, in this study, the theory about a plane stress element using hierarchical interpolation equation is reviewed and stiffness matrix is derived. And then, a computer program using this theory is developed. Developed computer program is used for numerical experiments to evaluate the analysis results using commercial programs such as SAP2000, ETABS, PERFORM-3D and MIDAS. Finally, the deflection equation of a cantilever beam with narrow rectangular section and bent by an end load P is derived according to the elasticity theory, and it is used to for comparison with theoretical solution.

• Coupled systems mechanics
• /
• 제6권4호
• /
• pp.501-522
• /
• 2017

This paper introduces a numerical procedure to incorporate elasto-plastic local deformation effects in the dynamic analysis of beams. The appealing feature is that simple beam type finite elements can be used for the global model which needs not to be altered by the localized elasto-plastic deformations. An overlapping local sophisticated 2D membrane model replaces the internal forces of the beam elements in the predefined region where the localized deformations take place. An iterative coupling technique is used to perform this replacement. Comparisons with full membrane analysis are provided in order to illustrate the accuracy and efficiency of the method developed herein. In this study, the membrane formulation is able to capture the elasto-plastic material behaviour based on the von Misses yield criterion and the associated flow rule for plane stress. The Newmark time integration method is adopted for the step-by-step dynamic analysis.

• 대한전기학회논문지:전기기기및에너지변환시스템부문B
• /
• 제55권6호
• /
• pp.306-313
• /
• 2006

A three dimensional adaptive finite element refinement algorithm is developed for non-linear magnetostatic field problems. In the method, the edge elements are used for finite element formulation, and the local error in each element is estimated from the fact that the tangential components of magnetic field intensity and the normal components of magnetic flux density should be continuous at the interface of the two adjacent elements. Based on the estimated error, the elements which have big error are divided into several elements using bisection method. The effectiveness of the developed algorithm is proved through numerical examples.

• Journal of Advanced Marine Engineering and Technology
• /
• 제6권1호
• /
• pp.41-48
• /
• 1982

One of the important subjects in fracture mechanics study is to analyze the stress intensity factor. In this paper, the stress intensity factor in Mode I ($K^{I}$) is determined by J-integral using the finite element method. In this investigation, the values of $K^{I}$ are computed for distorted and undistorted elements of 8-noded isoparametric finite elements. The numerical results obtained are summarized as follows. (1) Through a relatively coarse mesh, the $K^{I}$ values obtained by this method are fairly good accuracy. (2) The $K^{I}$ values for the distorted elements appear to be better than those obtained using the undistorted mesh. (3) Within the limits of these analyses, the solutions obtained through the integral paths in the medium region of elements approach to the analytical solution most closely.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제50권6호
• /
• pp.296-302
• /
• 2001

In this paper, the vector finite elements are adopted to calculate eigenvalues of RF and microwave components. Simulation results show that spurious are completely avoided because of the divergence free nature of the vector elements. This paper seeks to extend these low-order elements to higher orders to improve the accuracy of numerical solution. Investigation of numerical results for a rectangular waveguide was provided. A vector finite element program was implemented to allow propagation constants and electric field distributions to be directly computed in the rectangular and circular waveguides which are partially filled with the dielectric.

• Journal of Theoretical and Applied Mechanics
• /
• 제3권1호
• /
• pp.66-80
• /
• 2002

The bilinear formulation proposed earlier by Peters and Izadpanah to develop finite elements in time to solve undamped linear systems, Is extended (and found to be readily amenable) to develop time finite elements to obtain transient responses of both linear and nonlinear, and damped and undamped systems. The formulation Is used in the h-, p- and hp-versions. The resulting linear and nonlinear algebraic equations are differentiated to obtain the first- and second-order sensitivities of the transient response with respect to various system parameters. The present developments were tested on a series of linear and nonlinear examples and were found to yield, when compared with results obtained using other methods, excellent results for both the transient response and Its sensitivity to system parameters. Mostly. the results were obtained using the Legendre polynomials as basis functions, though. in some cases other orthogonal polynomials namely. the Hermite. the Chebyshev, and integrated Legendre polynomials were also employed (but to no great advantage). A key advantage of the time finite element method, and the one often overlooked in its past applications, is the ease In which the sensitivity of the transient response with respect to various system parameters can be obtained. The results of sensitivity analysis can be used for approximate schemes for efficient solution of design optimization problems. Also. the results can be applied to gradient-based parameter identification schemes.

• 수산해양기술연구
• /
• 제35권1호
• /
• pp.54-59
• /
• 1999

High-speed electronic digital computers have enabled engineers to employ various numerical discretization techniques for solutions of complex problems. The Finite Element Method is one of the such technique. The Finite Element Method is one of the numerical analysis based on the concepts of fundamental mathematical approximation. Three dimensional plate structures used often in partition of ship, box girder and frame are analyzed by Finite Element Method. In design of structures, the static deflections, stress concentrations and dynamic deflections must be considered. However, these problem belong to geometrically nonlinear mechanical structure analysis. The analysis of each element is independent, but coupling occurs in assembly process of elements. So, to overcome such a difficulty the shell theory which includes transformation matrix and a fictitious rotational stiffness is taken into account. Also, the Mindlin's theory which is considered the effect of shear deformation is used. The Mindlin's theory is based on assumption that the normal to the midsurface before deformation is "not necessarily normal to the midsurface after deformation", and is more powerful than Kirchoff's theory in thick plate analysis. To ensure that a small number of element can represent a relatively complex form of the type which is liable to occur in real, rather than in academic problem, eight-node quadratic isoparametric elements are used. are used.

• Journal of Welding and Joining
• /
• 제21권6호
• /
• pp.77-83
• /
• 2003

The finite element analysis was performed for the cracks existing in residual stress fields in order to investigate the effects of configuration of residual stress distribution to the fatigue crack opening behaviour. And the variation of stress distributions adjacent to the crack caused by uploading was examined. The finite element model with contact elements for the crack plane and plane stress elements for the base material and the analytical method based on the superposition principle to estimate crack opening behaviour and the stress distribution adjacent to the crack subjected to uploading were used. The results of the analysis showed that crack opening behaviors and variations of stress distribution caused by uploading were changed depending on the configuration of residual stress distribution. When the crack existed in the region of compressive residual stress and the configuration of compressive residual stress distribution were inclined, a partial crack opening just behind of a crack tip occurred during uploading. Based on the above results, it was clarified that the crack opening behaviour in the residual stress field could be predicted accurately by the finite element analysis using these analytical method and model.

• Interaction and multiscale mechanics
• /
• 제3권1호
• /
• pp.95-110
• /
• 2010

In recent years, study of the static response of pavements to moving vehicle and aircraft loads has received significant attention because of its relevance to the design of pavements and airport runways. The static response of beams resting on an elastic foundation and subjected to moving loads was studied by several researchers in the past. However, most of these studies were limited to steady-state analytical solutions for infinitely long beams resting on Winkler-type elastic foundations. Although the modelling of subgrade as a continuum is more accurate, such an approach can hardly be incorporated in analysis due to its complexity. In contrast, the two-parameter foundation model provides a better way for simulating the underlying soil medium and is conceptually more appealing than the one-parameter (Winkler) foundation model. The finite element method is one of the most suitable mathematical tools for analysing rigid pavements under moving loads. This paper presents an improved solution algorithm based on the finite element method for the static analysis of rigid pavements under moving vehicular or aircraft loads. The concrete pavement is discretized by finite and infinite beam elements, with the latter for modelling the infinity boundary conditions. The underlying soil medium is modelled by the Pasternak model allowing the shear interaction to exist between the spring elements. This can be accomplished by connecting the spring elements to a layer of incompressible vertical elements that can deform in transverse shear only. The deformations and forces maintaining equilibrium in the shear layer are considered by assuming the shear layer to be isotropic. A parametric study is conducted to investigate the effect of the position of moving loads on the response of pavement.