• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Composites Research
• /
• v.33 no.1
• /
• pp.30-37
• /
• 2020

Stress distributions dependent on fiber arrangements are studied using the two-dimensional representative volume element (RVE) model for uni-directionally continuous fiber-reinforced composites subjected to transverse tensile loading. It is easily expected that the stresses around the fibers are concentrated mainly due to the stiffness mismatch between the fiber and matrix materials. In this presentation, it is shown that the stresses are not always increased although the distance between two fibers is shortened. The 2D RVE models, originally having a regular hexagonal fiber array, is utilized to study the effect of the fiber locations on the stress distributions. As the central fiber is relocated, the stress distributions around the fiber are obtained through finite element analysis. It is found that the stresses around the fiber are strongly dependent on the fiber distance as well as the angle between the loading direction and the line connecting two fibers.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.21 no.4
• /
• pp.425-433
• /
• 2001

A correlation between fatigue damage and acousto-ultrasonic (AU) parameters has been obtained from signals acquired during fatigue loading of the single-lap joints of a carbon-fiber reinforced plastic (CFRP) laminates and A16061 plate. The correlation showed an analogy to those representing the stiffness reduction $(E/E_0)$ of polymer matrix composites by the accumulation of fatigue damage. This has been attributed to the transmission characteristics of acoustic wave energy through bonded joints with delamination-type defects and their influence on the change of spectral content of AU signals. Another correlation between fatigue cycles and the spectral magnitude of acoustic emission (AE) signals has also been found during the final stage of fatigue loading. Both AU and AE can be applied almost in real-time to monitor the evolution of damage during fatigue loading.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.19 no.6
• /
• pp.439-447
• /
• 1999

One of the problems for practical use of fiber-reinforced plastics is the performance degradation by fatigue damage in the joints. The study is to develop a nondestructive technique for real-time evaluation of adhesively bonded composite-metal joints. From the prior study we confirmed that the bonding strength can be estimated from the correlation between the qualify of bonded parts and AUP's. We obtained a curve showing the correlation between the degree of fatigue damage and AUP's calculated from signals acquired during fatigue loading of single-lap and double-lap joints of CFRP and Al6061. The curve is an analogy to the one showing stiffness reduction ($E/E_o$) of polymer matrix composites by fatigue damage. From those facts, it is plausible to predict the degree of fatigue damage in real-time. Amplitude and AUP2 appeared to be optimal parameters to provide more reliable results for single-lap joints whereas Amplitude and AUP2 did for double-lap joints. It is recommended to select optimal parameters for different geometries in the application for real structures.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.10
• /
• pp.920-927
• /
• 2015

In this paper, numerical results of sectional analysis and stress recovery were compared with the results of VABS through the blade analysis library. The results of recovery analysis for one-dimensional model including the stiffness matrix is compared with the calculated three-dimensional stress results of three-dimensionial FEM based on the principle of virtual work. We discuss the configuration of the blade analysis library and compare verifications of numerical analysis results of VABS. Blade analysis library through dimensional reduction and stress recovery is intended to be utilized in conjunction with pre- and post-processing of the analysis program of the composite blade, high-altitude uav's wing, wind blades and tilt rotor blade.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.67-78
• /
• 2012

A three-dimensional approximate computer-based method, YSPR (Yonsei Piled Raft), was developed for analysis of behavior of piled raft foundations. The raft was modeled as a flat shell element having 6 degrees of freedom at each node and the pile was modeled as a beam-column element. The behaviors of pile head and soil were controlled by using $6{\times}6$ stiffness matrix. To model the non-linear behavior, the soil-structure interaction between soil and pile was modeled by using nonlinear load-transfer curves (t-z, q-z and p-y curves). Comparison with previous model and FEM analysis showed that YSPR gave similar load-displacement behaviors. Comparison with field measurement also indicated that YSPR gave a reasonable result. It was concluded that YSPR could be effectively used in analysis and design of piled raft foundations.

• Journal of the Society of Disaster Information
• /
• v.13 no.2
• /
• pp.249-257
• /
• 2017

Recently, a cable disconnection accident occurred due to a lightning strike at the Seohae Bridge located in Dangjin-Pyeongtaek City. This is a natural occurrence, but it is a recall that it is very important to review the safety issues due to the disconnection of cable bridges. In other words, the role of cables in cable bridges has a profound effect on the safety of the structure, and it has become necessary to grasp the effect on the entire structural system. The cable bridge is an economic bridge that builds the main tower and supports the bottom plate by cable. The influence of the cable is the main member, which is a big influence on the safety of the whole bridge system. In the cable-stayed bridge, the cables exhibit nonlinear behavior because of the change in sag, due to the dead weight of the cable, which occurs with changing tension in the cable resulting from the movement of the end points of the cable as the bridge is loaded. Modal analysis is conducted using the deformed dead-load tangent stiffness matrix. A new concept was presented by using divided a cable into several elements in order to study the effect of the cable vibration (both in-plane and swinging) on the overall bridge dynamics. The result of this study demonstrates the importance of cable vibration on the overall bridge dynamics.

• Structural Engineering and Mechanics
• /
• v.13 no.5
• /
• pp.557-568
• /
• 2002

For the evaluation of the capability of a tubular member of an offshore structure to absorb the collision energy, a simple method can be employed for the collision analysis without performing the detailed analysis. The most common simple method is the rigid-plastic method. However, in this method any characteristics for horizontal movement and rotation at the ends of the corresponding tubular member are not included. In a real structural system of an offshore structure, tubular members sustain a certain degree of elastic support from the adjacent structure. End fixity has influences in the behaviors of a tubular member. Three-dimensional FEM analysis can include the effect of end fixity fully, however in viewpoints of the inherent computational complexities of the 3-D approach, this is not the recommendable analysis at the initial design stage. In this paper, influence of end fixity on the behaviors of a tubular member is investigated, through a new approach and other approaches. A new analysis approach that includes the flexibility of the boundary points of the member is developed here. The flexibility at the ends of a tubular element is extracted using the rational reduction of the modeling characteristics. The property reduction is based on the static condensation of the related global stiffness matrix of a model to end nodal points of the tubular element. The load-displacement relation at the collision point of the tubular member with and without the end flexibility is obtained and compared. The new method lies between the rigid-plastic method and the 3-demensional analysis. It is self-evident that the rigid-plastic method gives high strengthening membrane effect of the member during global deformation, resulting in a steeper slope than the present method. On the while, full 3-D analysis gives less strengthening membrane effect on the member, resulting in a slow going load-displacement curve. Comparison of the load-displacement curves by the new approach with those by conventional methods gives the figures of the influence of end fixity on post-yielding behaviors of the relevant tubular member. One of the main contributions of this investigation is the development of an analytical rational procedure to figure out the post-yielding behaviors of a tubular member in offshore structures.

• Structural Engineering and Mechanics
• /
• v.4 no.5
• /
• pp.553-568
• /
• 1996

Finite element stiffness matrix methods are presented for finding natural frequencies (or buckling loads) and modes of repetitive structures. The usual approximate finite element formulations are included, but more relevantly they also permit the use of 'exact finite elements', which account for distributed mass exactly by solving appropriate differential equations. A transcendental eigenvalue problem results, for which all the natural frequencies are found with certainty. The calculations are performed for a single repeating portion of a rotationally or linearly (in one, two or three directions) repetitive structure. The emphasis is on rotational periodicity, for which principal advantages include: any repeating portions can be connected together, not just adjacent ones; nodes can lie on, and members along, the axis of rotational periodicity; complex arithmetic is used for brevity of presentation and speed of computation; two types of rotationally periodic substructures can be used in a multi-level manner; multi-level non-periodic substructuring is permitted within the repeating portions of parent rotationally periodic structures or substructures and; all the substructuring is exact, i.e., the same answers are obtained whether or not substructuring is used. Numerical results are given for a rotationally periodic structure by using exact finite elements and two levels of rotationally periodic substructures. The solution time is about 500 times faster than if none of the rotational periodicity had been used. The solution time would have been about ten times faster still if the software used had included all the substructuring features presented.