• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.209-216
• /
• 1996

The primary objective of this paper is to grasp many characteristics of buckling behavior of latticed spherical domes under various conditions. The Arc-Length Method proposed by E.Riks is used for the computation and evaluation of geometrically nonlinear fundamental equilibrium paths and bifurcation points. And the direction of the path after the bifurcation point is decided by means of Hosono's concept. Three different nonlinear stiffness matrices of the Slope-Deflection Method are derived for the system with rigid nodes and results of the numerical analysis are examined in regard to geometrical parameters such as slenderness ratio, half-open angle, boundary conditions, and various loading types. But in case of analytical model 2 (rigid node), the post-buckling path could not be surveyed because of Newton-Raphson iteration process being diversed on the critical point since many eigenvalues become zero simultaneously.

• Nuclear Engineering and Technology
• /
• v.46 no.2
• /
• pp.263-272
• /
• 2014

Scattering source calculations using conventional spherical harmonic expansion may require lots of computation time to treat full-coupled three-dimensional photon-electron transport in a highly anisotropic scattering medium where their scattering cross sections should be expanded with very high order (e.g., $P_7$ or higher) Legendre expansions. In this paper, we introduce a modified scattering kernel approach to avoid the unnecessarily repeated calculations involved with the scattering source calculation, and used it with parallel computing to effectively reduce the computation time. Its computational efficiency was tested for three-dimensional full-coupled photon-electron transport problems using our computer program which solves the multi-group discrete ordinates transport equation by using the discontinuous finite element method with unstructured tetrahedral meshes for complicated geometrical problems. The numerical tests show that we can improve speed up to 17~42 times for the elapsed time per iteration using the modified scattering kernel, not only in the single CPU calculation but also in the parallel computing with several CPUs.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.2
• /
• pp.193-201
• /
• 2004

The aim of this study is focused on getting an almost exact solution which is the simplicity and exactness of an axisymmetrically loaded cylindrical shell. This method replaces the finite element method which is a very powerful tool for analysis of any kind of structure which has an arbitrary shape, but is still a numerical analysis. Instead, this study uses the method of distribution of end actions which is a kind of iteration technique to implement the leading matrix method. The distribution and carry-over factors of a cylinder are calculated by the theory of a differential equation of a beam on an elastic foundation. The results are satisfactory when this method is applied to a cylinder that is subjected to a concentrated load and hydrostatic pressure when compared with the BEF analogy separately.

• Journal of Magnetics
• /
• v.16 no.2
• /
• pp.181-185
• /
• 2011

To improve trapped field characteristics of a high temperature superconducting (HTS) bulk, a technique to implement artificial holes has been studied. The artificial holes, filled up with epoxy or metal, may provide better cooling channel and enhance mechanical strength of the HTS bulk. Although many useful researches based on experiments have been reported, a numerical approach is still limited because of several reasons that include: 1) highly non-linear electromagnetic properties of HTS; and 2) difficulty in modeling of randomly scattered "small" artificial holes. In this paper, a 2-D finite element method with iteration is adopted to analyze trapped field characteristics of HTS bulk with artificial holes. The validity of the calculation is verified by comparison between measurement and calculation of a trapped field in a $40{\times}40\;mm$ square and 3.1 mm thick HTS bulk having 16 artificial holes with diameter of 0.7 mm. The effects of sizes and array patterns of artificial holes on distribution of trapped field within HTS bulk are numerically investigated using suggested method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.2 s.257
• /
• pp.211-220
• /
• 2007

This work presents an iterated improved reduced system (IIRS) procedure combined with sub-structuring scheme for large structures. Iterated IRS methods are usually more efficient than others because the dynamic condensation matrix is updated repeatedly until the desired convergent values are obtained. However, using these methods simply for large structures causes expensive computational cost and even makes analyses intractable because of the limited computer storage. Therefore, the application of sub-structuring scheme is necessary. Because the large structures are subdivided into several (or more) sub-domains, the construction of dynamic condensation matrix does not require much computation cost in every iteration. This makes the present method much more efficient to compute the eigenpairs both in lower and intermediate modes. In Part I, iterated IRS method combined with sub-structuring scheme for undamped structures is presented. The validation of the proposed method and the evaluation of computational efficiency are demonstrated through the numerical examples.

• Structural Engineering and Mechanics
• /
• v.8 no.3
• /
• pp.233-242
• /
• 1999

In the present study, the elasto-plastic analysis of prismatic plate structures subjected to pure bending is carried out using the finite strip method. The end cross-sections of the structure are assumed to remain plane during deformation, and the compatibility along corner lines is ensured by choosing proper displacement functions. The effects of both the initial geometrical imperfections and residual stresses due to fabrication are included in the combined geometrically and materially nonlinear simulation. The von-Mises yield criterion and the Prandtl-Reuss flow theory of plasticity are applied in modelling the elasto-plastic behavior of material. Newton-Raphson iterations are carried out as the rotation of the end cross sections of the structure is increased step by step. The parameter representing the overall axial strain of structure is adjusted constantly during the iteration process in order to eliminate the resulting overall axial force on any cross-section of the structure in correspondence with the assumption of zero axial force in pure bending. Several numerical examples are presented to validate the present method and to investigate the effects of some material and geometrical parameters.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.2
• /
• pp.109-123
• /
• 2008

The present paper presents the CFD result for a beam wave test case. An ONR tumblehome ship model with bilge keels is used. The beam wave test is for zero forward speed and roll and heave 2DOF with wave slope $a_k=0.156$ and wavelength ${\lambda}=1.12L_{PP}$, with $L_{PP}$ the ship length. The problems is solved numerically with an unsteady Reynolds averaged Navier-Stokes approach. The free surface flow is computed using a single-phase level-set method and the motions in each time step are integrated using a predictor-corrector iteration approach which uses dynamic overset grids moving with relative ship motion. The predicted CFD results for motions and forces are compared with experimental data, showing a reasonable agreement.

• Smart Structures and Systems
• /
• v.12 no.3_4
• /
• pp.235-250
• /
• 2013

Optimal sensor placement (OSP) technique plays a key role in the structural health monitoring (SHM) of large-scale structures. According to the mathematical background and implicit assumptions made in the triaxial effective independence (EfI) method, this paper presents a novel multi-dimensional OSP method for the Canton Tower focusing on application demands. In contrast to existing methods, the presented method renders the corresponding target mode shape partitions as linearly independent as possible and, at the same time, maintains the stability of the modal matrix in the iteration process. The modal assurance criterion (MAC), determinant of the Fisher Information Matrix (FIM) and condition number of the FIM have been taken as the optimal criteria, respectively, to demonstrate the feasibility and effectiveness of the proposed method. Numerical investigations suggest that the proposed method outperforms the original EfI method in all instances as expected, which is looked forward to be even more pronounced should it be used for other multi-dimensional optimization problems.