• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.10a
• /
• pp.25-30
• /
• 1994

A vector algorithm for calculating the stiffness matrices of reinforced concrete shell elements is presented. The algorithm is based on establishing vector lengths equal to the number of elements. The computational efficiency of the proposed algorithm is assessed on a Cray Y-MP supercomputer. It is shown that the vector algorithm achieves scalar-to-vector speedup of 1.7 to 7.6 on three inelastic problems.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1997.10a
• /
• pp.659-663
• /
• 1997

In this paper, a simplified vibration analysis model for bellows was presented to avoid excessive effort required for shell model. To reduce degree of freedom, bellows was modelled using one dimensional beam element. The equivalent mass and stiffness matrices were obtained based upon Guyan reduction process. The results were compared with the confirmed results, which were in good agreement.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1518-1525
• /
• 2006

Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.2
• /
• pp.250-255
• /
• 2008

Springs are widely utilized in machine element. To find out stiffness of frustum-shaped coil spring, the space beam theory using the finite element method is adopted in this paper. In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. To find out load vector of coil spring subjected to distributed compression. principle of virtual work is adapted. And this theory was programming using MATLAB software. To compare FEM using MATLAB software was applied MSC. Nastran software. The geometry model for MSC. Patran was produced by 3-D design modeling software. Finite element model was produced by MSC. Patran. Finite element was applied tetra (CTETRA) having 10 node. The analysis results of the MATLAB and MSC. Nastran are fairly well agreed with those of various experiments. Using MATLAB program proposed in this paper and MSC. Nastran, spring constants and stresses can be predicted by input of few factors.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.363-373
• /
• 1995

In this study, a FE model tuning method using experimental modal data was suggested after examining all the published conventional methods. The idea of this method is introducing scale factors to maintain both the structural connectivity and the consistency in the corrected stiffness matrix which makes it always possible to interpret the stiffness elements with the corresponding physical configuration of the targeting structure. The scale factors are determined to minimize the objective function of eigen-properties. The proposed method was tested to determine the joint stiffness of a T shaped beam. The test results were also compared with the tuned stiffness obtained from a probed commercial package (SYSTUNE) and found that this method is very accurate and compatible.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.169-177
• /
• 1996

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear programming formulation of the topology problem is presented. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.4 no.1
• /
• pp.1-8
• /
• 2013

The finite element analysis (FEA) is a numerical technique to find solutions of field problems. A field problem is approximated by differential equations or integral expressions. In a finite element, the field quantity is allowed to have a simple spatial variation in terms of linear or polynomial functions. This paper represents a review and an accuracy-study of the finite element method comparing the FEA results with the exact solution. The exact solutions were calculated by solid mechanics and FEA using matrix stiffness method. For this study, simple bar and cantilever models were considered to evaluate four types of basic elements - constant strain triangle (CST), linear strain triangle (LST), bi-linear-rectangle(Q4),and quadratic-rectangle(Q8). The bar model was subjected to uniaxial loading whereas in case of the cantilever model moment loading was used. In the uniaxial loading case, all basic element results of the displacement and stress in x-direction agreed well with the exact solutions. In the moment loading case, the displacement in y-direction using LST and Q8 elements were acceptable compared to the exact solution, but CST and Q4 elements had to be improved by the mesh refinement.

• Structural Engineering and Mechanics
• /
• v.31 no.2
• /
• pp.137-152
• /
• 2009

Mechanic behavior of Y-shape thin-walled box girder bridge structure is complex, so one can not exactly hold the mechanical behavior of the Y-shape thin-walled box girder bridge structure through general calculation theory and analytical method. To hold the mechanical behavior better, based on elementary beam theory, by increasing the degree of freedom analytical method, taking account of restrained torsiondistortion angledistortion warp and shearing lag effect at the same time, authors obtain a thin-walled box beam analytical element of 10 degrees of freedom of every node, derive stiffness matrix of the element, and code a finite element procedure. In addition, authors combine the obtained procedure with spatial grillage analytical method, meanwhile, they build a new analytical method that is the spatial thin-walled box girder element grillage analysis method. In order to validate the precision of the obtained analysis method, authors analyze a type Y-shape thin-walled box girder bridge structure according to the elementary beam theory analytical method, the shell theory analytical method and the spatial thin-walled box girder element grillage analysis method respectively. At last, authors test a type Y-shape thin-walled box girder bridge structure. Comparisons of the results of theory analysis with the experimental text show that the spatial thin-walled box girder element grillage analysis method is simple and exact. The research results are helpful for the knowledge of the mechanics property of these Y-shape thin-walled box girder bridge structures.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.990-1004
• /
• 1995

A modified 16-node element for composite plate has been proposed and compared with the 20-node element to check the validity of it. The fields of numerical inspection include mode analysis and specific damping analysis. By symetrizing the conventional unsymmetric damping matrix in the analysis of specific damping capacity, we could compute the specific damping capacity and make a program, effectively. In addition, we could predict the errors caused by reduction of integration order in thickness direction depending upon the number of layers.