• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.12
• /
• pp.1338-1346
• /
• 2009

We present a one-dimensional annular plate element with which the in-plane vibration of full disks can be analyzed efficiently and accurately by using the FEM. Its elementary mass matrix and stiffness matrix are derived, respectively, from the virtual work by effective forces and the virtual strain energy. The static deformation modes obtained from an integration of the differential equilibrium equations of the annular plate are used as interpolation functions of the one-dimensional annular plate element. The in-plane natural vibration characteristics of a 2-step full disk and a uniform full disk are analysed. Its results are compared with the results obtained by utilizing two-dimensional 8-node quadrilateral plane elements and cyclic symmetry of the disk. And also, by comparing with the theoretical results of previous researchers, the efficiency and accuracy of the presented element are verified.

• Structural Engineering and Mechanics
• /
• v.74 no.1
• /
• pp.1-18
• /
• 2020

For the large deformation of shear walls under vertical and horizontal loads, there are difficulties in obtaining accurate simulation results using the response analysis method, even with fine mesh elements. Furthermore, concrete material nonlinearity, stiffness degradation, concrete cracking and crushing, and steel bar damage may occur during the large deformation of reinforced concrete (RC) shear walls. Matrix operations that are involved in nonlinear analysis using the traditional finite-element method (FEM) may also result in flaws, and may thus lead to serious errors. To solve these problems, a planar four-node element was developed based on vector mechanics. Owing to particle-based formulation along the path element, the method does not require repeated constructions of a global stiffness matrix for the nonlinear behavior of the structure. The nonlinear concrete constitutive model and bilinear steel material model are integrated with the developed element, to ensure that large deformation and damage behavior can be addressed. For verification, simulation analyses were performed to obtain experimental results on an RC shear wall subjected to a monotonically increasing lateral load with a constant vertical load. To appropriately evaluate the parameters, investigations were conducted on the loading speed, meshing dimension, and the damping factor, because vector mechanics is based on the equation of motion. The static problem was then verified to obtain a stable solution by employing a balanced equation of motion. Using the parameters obtained, the simulated pushover response, including the bearing capacity, deformation ability, curvature development, and energy dissipation, were found to be in accordance with the experimental observation. This study demonstrated the potential of the developed planar element for simulating the entire process of large deformation and damage behavior in RC shear walls.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.83-90
• /
• 2004

Substructuring methods are usually used in finite element structural analyses. In this study a multi-level substructuring algorithm is developed and proposed as a possible candidate for incompressible fluid solves. Finite element formulation for incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et.al.[5]. The present algorithm consists of four stages such as a gathering stage, a condensing stage, a solving stage and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At highest level, each subdomain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each subdomain has been replaced by a sequential static condensation. The global algebraic system arising feom the assembly of each subdomains is solved using Conjugate Gradient Squared(CGS) method. In this case, pre-conditioning techniques usually accompanied by iterative solvers are not needed.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.101-113
• /
• 1996

In this paper, the finite element formulation for the thermal analysis of quasi-static, uncoupled, homogeneous, isotropic and linear viscoelastic problems is presented based on the principle of virtual work. The viscoelastic material is assumed to be thermorheologically simple, which is well known material property in a large class of high polymers. The variational formulation and the finite element equation in matrix from are derived. Effective generation and storage of the hereditary stiffness matrices are given in detail especially for the case of the steady state temperature distribution T=T(x). Some numerical examples are given and compared with published results to show the versatility of the derived finite element formulations.

• Advances in Automotive Engineering
• /
• v.1 no.1
• /
• pp.41-59
• /
• 2018

In the study, investigation of fiber- reinforced composite materials that can be an alternative to conventional steel was performed by finite element analysis with the help of software. Steel and composite materials have been studied on a four axle truck chassis model. Three-dimensional finite element model was created with software, and then analyzes were performed. The analyses were performed for static and dynamic/fatigue cases. Fatigue cases are formed with the help of design spectra model and fatigue analyses were performed as static analyses with this design spectra. First, analyses were performed for steel and after that optimization analyses were made for the AS4-PEEK carbon fiber composite and Eglass-Epoxy fiber composite materials. Optimization of composite material analyzes include determining the total laminate thickness, thickness of each ply, orientation of each ply and ply stacking sequence. Analyzes were made according to macro mechanical properties of composite, micromechanics case has not been considered. Improvements in weight reduction up to %50 provided at the end of the composite optimization analyzes with satisfying stiffness performance of chassis. Fatigue strength of the composite structure depends on various factors such as, fiber orientation, ply thickness, ply stack sequence, fiber ductility, ductility of the matrix, loading angle. Therefore, the accuracy of theoretical calculations and analyzes should be correlated by testing.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.82-91
• /
• 1996

This is a basic study for the static and dynamic analysis on the elasto-plastic and elasto-viscoplastic of an axi-symmetric shell. The objective of this study was to investigate the mechanical characteristics of a nuclear reinforced concrete containment structure, which was selected as a model, by a numerical analysis using a finite element method. The structure was modeled with discrete ring elements of 8-noded isoparametric element rotating against the symmetrical axis, and the interaction between the foundation and the structure was modeled by Winkler's model. Also, the meridional tendon was modeled with 2-node truss elements, and the hoop tendon was done with point elements in two degrees of freedom. The effect of the tendon was considered without the increasement in total degree of freedom as the stiffness matrix of modeled tendon elements was assembled on the stiffness matrix of ring elements linked with the tendon. The results obtained from the analysis of an example were summarized as follows : 1. The stresses in the hoop direction on the interior and exterior surfaces of the structure were shown in changes of similar trend, and high stresses appeared on the structure wall 2. The stresses in the meridional direction on the interior and exterior surfaces were shown in change of different trend. Especially, the stresses at the junctions between the dome and the wall and between the wall and the bottom plate of the structure were very high, compared with those at other parts of the structure. 3. The stress changes in the direction of thickness on the crown of the dome were much linearly distributed. However, as the amount of tendon increased, the stresses in the upper and lower parts of the wall established with the tendon were shown stress concentration. 4. The stress changes in the direction of thickness on the center of the structure wall was linearly distributed in the all cases, and special stress due to the use of the tendon was not shown.

• Structural Engineering and Mechanics
• /
• v.43 no.5
• /
• pp.583-605
• /
• 2012

This paper presents a spatial catenary cable element for the nonlinear analysis of cable-supported structures. An incremental-iterative solution based on the Newton-Raphson method is adopted for solving the equilibrium equation. As a result, the element stiffness matrix and nodal forces are determined, wherein the effect of self-weight and pretension are taken into account. In the case of the initial cable tension is given, an algorithm for form-finding of cable-supported structures is proposed to determine precisely the unstressed length of the cables. Several classical numerical examples are solved and compared with the other available numerical methods or experiment tests showing the accuracy and efficiency of the present elements.

• Steel and Composite Structures
• /
• v.47 no.3
• /
• pp.365-374
• /
• 2023

This paper presents a new scheme for constructing locking-free finite elements in thick and thin plates, called Discrete Shear Gap element (DSG), using multiphase material topology optimization for triangular elements of Reissner-Mindlin plates. Besides, common methods are also presented in this article, such as quadrilateral element (Q4) and reduced integration method. Moreover, when the plate gets too thin, the transverse shear-locking problem arises. To avoid that phenomenon, the stabilized discrete shear gap technique is utilized in the DSG3 system stiffness matrix formulation. The accuracy and efficiency of DSG are demonstrated by the numerical examples, and many superior properties are presented, such as being a strong competitor to the common kind of Q4 elements in the static topology optimization and its computed results are confirmed against those derived from the three-node triangular element, and other existing solutions.

• Ocean Systems Engineering
• /
• v.10 no.1
• /
• pp.67-86
• /
• 2020

The main goal of this study is to investigate the free vibration analysis of a large sag catenary with application to the jumper in hybrid riser system. The equation of motion is derived by using the variational method based on the virtual work principle. The finite element method is applied to evaluate the numerical solutions. The large sag catenary is utilized as an initial configuration for vibration analysis. The nonlinearity due to the large sag curvature of static configuration is taken into account in the element stiffness matrix. The natural frequencies of large sag catenary and their corresponding mode shapes are determined by solving the eigenvalue problem. The numerical examples of a large sag catenary jumpers are presented. The influences of bending rigidity and large sag shape on the free vibration behaviors of the catenary jumper are provided. The results indicate that the increase in sag reduces the jumper natural frequencies. The corresponding mode shapes of the jumper with large sag catenary shape are comprised of normal and tangential displacements. The large sag curvature including in the element stiffness matrix increases the natural frequency especially for a case of very large sag shape. Mostly, the mode shapes of jumper are dominated by the normal displacement, however, the tangential displacement significantly occurs around the lowest point of sag. The increase in degree of inclination of the catenary tends to increase the natural frequencies.

• Advanced Composite Materials
• /
• v.18 no.2
• /
• pp.105-119
• /
• 2009

In this work, the structural response of thin-walled composite beams with pretwist angle is investigated by using a mixed beam approach that combines the stiffness and flexibility methods in a unified manner. The Reissner's semi-complimentary energy functional is used to derive the stiffness matrix that approximates the beam in an Euler-Bernoulli level for extension and bending and Vlasov level for torsion. The bending and torsion-related warpings induced by the pretwist effects are derived in a closed form. The developed theory is validated with available literature and detailed finite element structural analysis results using the MSC/NASTRAN. Pretwisted composite beams with rectangular solid and thin-walled box sections are illustrated to validate the current approach. Acceptable correlation has been achieved for cases considered in this study. The effects of pretwist and fiber orientation angles on the static behavior of pretwisted composite beams are also studied.