• Structural Engineering and Mechanics
• /
• v.61 no.6
• /
• pp.807-816
• /
• 2017

In this paper, we introduce an efficient new model reduction method, named the automated static condensation method, which is developed for the local analysis of large finite element models. The algebraic multilevel substructuring procedure is modified appropriately, and then applied to the original static condensation method. The retained substructure, which is the local finite element model to be analyzed, is defined, and then the remaining part of the global model is automatically partitioned into many omitted substructures in an algebraic perspective. For an efficient condensation procedure, a substructural tree diagram and substructural sets are established. Using these, the omitted substructures are sequentially condensed into the retained substructure to construct the reduced model. Using several large practical engineering problems, the performance of the proposed method is demonstrated in terms of its solution accuracy and computational efficiency, compared to the original static condensation method and the superelement technique.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1992.04a
• /
• pp.65-70
• /
• 1992

This study deals with the finite element analysis of the monotonic behavior of reinforced concrete beams and beam-column joint subassemblages. It is assumed that the behavior of these members can be discribed by a plane stress field. Concrete and reinforcing steel are represented by separate material models which are combined together with a model of the interaction between reinforcing bar and concrete through bond-slip to discribe the behavior of the composite reinforced concrete material. To discribe the concrete behavior, a nonlinear orthotropic model is adopted and the crack is discribed by a system of orthogonal cracks, which are rotating as the principal strain directions are changed. A smeared finite element model based on the fracture mechanics principles are used to overcome the numerical defect according to the finite element mesh size. Finally, correlation studies between analytical and experimental results and several parameter studies are conducted with the objective to estabilish the validity of the proposed model and identify the significance of various effects on the local and global response of reinforced concrete members.

• Structural Engineering and Mechanics
• /
• v.23 no.6
• /
• pp.691-711
• /
• 2006

This paper uses the finite element method to simultaneously consider the coupled cable-deck vibrations and the parametric vibrations of stay cables in dynamic analysis of a cable-stayed bridge. The stay cables are represented by some cable finite elements, which can consider the parametric vibration of the cables. In addition to modeling stay cables using multiple link cable elements, a procedure for removing the self-weight term of cable element is proposed. A eigenvalue analysis process using dynamic condensation method for sorting out the natural modes of the girder-tower vibrations and the Rayleigh damping considering element damping for damping matrix are also proposed for dynamic analyses of cable-stayed bridges. The possibilities of using cable elements and of using global and local vibrations to evaluate the parametric vibrations of stay cables in a cable-stayed bridge are confirmed, respectively.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.05a
• /
• pp.571-577
• /
• 2003

This research proposes a finite element method to determine the natural vibration characteristics of the spinning disk-spindle system in a HDD including the flexibility of supporting structure. Finite element equations of each substructure are derived with the introduction of consistent variables to satisfy the geometric compatibility at the internal boundaries. The natural frequencies and modes from the global asymmetric matrix equations of motion are determined by using the implicit restarted Arnoldi iteration method. The validity of the proposed method is verified by the experimental modal testing. It also shows that the flexibility of base plate plays an important role to determine the natural frequencies of the spinning disk-spindle system in a HDD.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.3
• /
• pp.196-203
• /
• 2007

The postbuckling behavior and failure of composite cylinders subjected to external hydrostatic pressure are investigated by a finite element method and test. A nonlinear finite element program, ACOS, is used for the postbuckling progressive failure analysis of composite cylinders. A total of 5 carbon/epoxy composite cylinders were fabricated and tested to verify the finite element results. For comparison, analyses by MSC/NASTRAN and MSC/MARC are additionally conducted. Among the softwares, the finite element program, ACOS, predicts the buckling loads the best with about 11 to 26% deviation from experimental results except for one specimen. While the finite element analysis shows global buckling modes with 4 waves in hoop direction, in the experiments the local buckling appears first and results in the final failure without global buckling.

• Structural Engineering and Mechanics
• /
• v.38 no.4
• /
• pp.529-547
• /
• 2011

This paper presents the global-local finite cover method (GL-FCM) that is capable of analyzing structures involving local heterogeneities and propagating cracks. The suggested method is composed of two techniques. One of them is the FCM, which is one of the PU-based generalized finite element methods, for the analysis of local cohesive crack growth. The mechanical behavior evaluated in local heterogeneous structures by the FCM is transferred to the overall (global) structure by the so-called mortar method. The other is a method of mesh superposition for hierarchical modeling, which enables us to evaluate the average stiffness by the analysis of local heterogeneous structures not subjected to crack propagation. Several numerical experiments are conducted to validate the accuracy of the proposed method. The capability and applicability of the proposed method is demonstrated in an illustrative numerical example, in which we predict the mechanical deterioration of a reinforced concrete (RC) structure, whose local regions are subjected to propagating cracks induced by reinforcement corrosion.

• Structural Engineering and Mechanics
• /
• v.12 no.5
• /
• pp.507-526
• /
• 2001

In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

• International Journal of CAD/CAM
• /
• v.8 no.1
• /
• pp.29-36
• /
• 2009

Domain mapping is a bijective transformation of one domain to another, usually from a complicated general domain to a chosen convex domain. This is directly useful in many application problems like shape modeling, morphing, texture mapping, shape matching, remeshing, path planning etc. A new approach considering the domain as made up of structural elements, like membranes or trusses, is developed and implemented using the nonlinear finite element formulation. The mapping is performed in two stages, boundary mapping and inside mapping. The boundary of the 3-D domain is mapped to the surface of a convex domain (in this case, a sphere) in the first stage and then the displacement/distortion of this boundary is used as boundary conditions for mapping the interior of the domain in the second stage. This is a general method and it develops a bijective mapping in all cases with judicious choice of material properties and finite element analysis. The consistent global parameterization produced by this method for an arbitrary genus zero closed surface is useful in shape modeling. Results are convincing to accept this finite element structural approach for domain mapping as a good method for many purposes.

• Structural Engineering and Mechanics
• /
• v.69 no.6
• /
• pp.615-626
• /
• 2019

The 3-noded metric Timoshenko beam element with an offset of the internal node from the element centre is used here to demonstrate the best-fit paradigm using function space formulation under locking and mesh distortion. The best-fit paradigm follows from the projection theorem describing finite element analysis which shows that the stresses computed by the displacement finite element procedure are the best approximation of the true stresses at an element level as well as global level. In this paper, closed form best-fit solutions are arrived for the 3-noded Timoshenko beam element through function space formulation by combining field consistency requirements and distortion effects for the element modelled in metric Cartesian coordinates. It is demonstrated through projection theorems how lock-free best-fit solutions are arrived even under mesh distortion by using a consistent definition for the shear strain field. It is shown how the field consistency enforced finite element solution differ from the best-fit solution by an extraneous response resulting from an additional spurious force vector. However, it can be observed that when the extraneous forces vanish fortuitously, the field consistent solution coincides with the best-fit strain solution.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.172-176
• /
• 2008

In this paper, three dimensional linear conforming variable-finite elements are presented with the aid of a smoothed integration (a class of stabilized conforming nodal integration), for mnltiscale mechanics problems. These elements meet the desirable properties of an interpolation such as the Kronecker delta condition, the partition of unity condition and the positiveness of interpolation function. The necessary condition of linear exactness is fully relaxed by employing the smoothed integration, which renders us to meet the linear exactness in a straightforward manner. This novel element description extend the category of the conventional finite elements space to ration type function space and give the flexibility on the number of nodes of element which are fixed in the conventional finite elements. Several examples are provided to show the convergence and the accuracy of the proposed elements, and to demonstrate their potential with emphasis on the multiscale mechanics problems such as global/local analysis, nonmatching contact problems, and modeling of composite material with defects.