• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.255-262
• /
• 2013

In this paper, a shape design sensitivity analysis(DSA) method is presented for Mindlin plates using an isogeometric approach. The isogeometric method possesses desirable advantages; the representation of exact geometry and the higher order inter-element continuity, which lead to the fast convergence of solution as well as accurate sensitivity results. Unlike the finite element methods using linear shape functions, the isogeometric method considers the exact normal vector and curvature of the CAD geometry, taking advantages of higher order NURBS basis functions. A selective reduced integration(SRI) technique is incorporated to overcome the difficulty of 'shear locking' phenomenon. This simple technique is surprisingly helpful for the accuracy of the isogeometric shape sensitivity without complicated formulation. Through the numerical examples of plate bending problems, the accuracy of the proposed isogeometric analysis method is compared with that of finite element one. Also, the isogeometric shape sensitivity turns out to be very accurate when compared with finite difference sensitivity.

• Nuclear Engineering and Technology
• /
• v.49 no.4
• /
• pp.817-824
• /
• 2017

Residual stress is a critical element in determining the integrity of parts and the lifetime of welded structures. It is necessary to estimate the residual stress of a welding zone because residual stress is a major reason for the generation of primary water stress corrosion cracking in nuclear power plants. That is, it is necessary to estimate the distribution of the residual stress in welding of dissimilar metals under manifold welding conditions. In this study, a cascaded support vector regression (CSVR) model was presented to estimate the residual stress of a welding zone. The CSVR model was serially and consecutively structured in terms of SVR modules. Using numerical data obtained from finite element analysis by a subtractive clustering method, learning data that explained the characteristic behavior of the residual stress of a welding zone were selected to optimize the proposed model. The results suggest that the CSVR model yielded a better estimation performance when compared with a classic SVR model.

• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.129-152
• /
• 2008

This paper considers the solution of the stochastic differential equations (SDEs) with random operator and/or random excitation using the spectral SFEM. The random system parameters (involved in the operator) and the random excitations are modeled as second order stochastic processes defined only by their means and covariance functions. All random fields dealt with in this paper are continuous and do not have known explicit forms dependent on the spatial dimension. This fact makes the usage of the finite element (FE) analysis be difficult. Relying on the spectral properties of the covariance function, the Karhunen-Loeve expansion is used to represent these processes to overcome this difficulty. Then, a spectral approximation for the stochastic response (solution) of the SDE is obtained based on the implementation of the concept of generalized inverse defined by the Neumann expansion. This leads to an explicit expression for the solution process as a multivariate polynomial functional of a set of uncorrelated random variables that enables us to compute the statistical moments of the solution vector. To check the validity of this method, two applications are introduced which are, randomly loaded simply supported reinforced concrete beam and reinforced concrete cantilever beam with random bending rigidity. Finally, a more general application, randomly loaded simply supported reinforced concrete beam with random bending rigidity, is presented to illustrate the method.

• Structural Engineering and Mechanics
• /
• v.15 no.3
• /
• pp.299-314
• /
• 2003

This investigation presents an efficient method for identifying modal characteristics from the measured displacement, velocity and acceleration signals of multiple channels on structural systems. A Vector Backward Auto-Regressive model (VBAR) that describes the relationship between the output information in different time steps is used to establish a backward state equation. Generally, the accuracy of the identified dynamic characteristics can be improved by increasing the order of the Auto-Regressive model (AR) in cases of measurement of data under noisy circumstances. However, a higher-order AR model also induces more numerical modes, only some of which are the system modes. The proposed VBAR model provides a clear characteristic boundary to separate the system modes from the spurious modes. A numerical example of a lumped-mass model with three DOFs was established to verify the applicability and effectiveness of the proposed method. Finally, an offshore platform model was experimentally employed as an application case to confirm the proposed VBAR method can be applied to real-world structures.

• Smart Structures and Systems
• /
• v.18 no.6
• /
• pp.1233-1250
• /
• 2016

Cracks in plate-like structures are some of the main reasons for destruction of the entire structure. In this study, a novel two-stage methodology is proposed for damage detection of flexural plates using an optimized artificial neural network. In the first stage, location of damages in plates is investigated using curvature-moment and curvature-moment derivative concepts. After detecting the damaged areas, the equations for damage severity detection are solved via Bat Algorithm (BA). In the second stage, in order to efficiently reduce the computational cost of model updating during the optimization process of damage severity detection, multiple damage location assurance criterion index based on the frequency change vector of structures are evaluated using properly trained cascade feed-forward neural network (CFNN) as a surrogate model. In order to achieve the most generalized neural network as a surrogate model, its structure is optimized using binary version of BA. To validate this proposed solution method, two examples are presented. The results indicate that after determining the damage location based on curvature-moment derivative concept, the proposed solution method for damage severity detection leads to significant reduction of computational time compared with direct finite element method. Furthermore, integrating BA with the efficient approximation mechanism of finite element model, maintains the acceptable accuracy of damage severity detection.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1243-1250
• /
• 1995

For the numerical solution of frictional dynamic contact problems, correct contact points and displacements are determined by iteratively reducing the displacement error vector monotonically toward zero And spurious oscillations are prevented from the solution by enforcing the velocity and acceleration compatibilities of the contact points with the corresponding error vectors. Numerical simulations are conducted to demonstrate the accuracy of the solution and the necessity of the velocity and acceleration compatibilities on the contact surface.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1993.04a
• /
• pp.61-68
• /
• 1993

A computer program for free vibration analysis of plane framed structures has been developed by object oriented programming technique using C" language. The object oriented programming concepts such as object, class, method and inheritance are represented. The static and free vibration analyses for framed structures were satisfactorily performed by this program which consists of TOP, VECTOR, MATRIX, STRU, GUI and other classes. Numerical test shows the validity and capability of the present study which can be expandable to develop a general purpose object oriented finite element analysis program of structures ,res ,

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.3 s.96
• /
• pp.251-258
• /
• 2005

The free vibration of a spinning flexible disk-spindle system supported by hydro dynamic bearings (HDB) in an HDD is analyzed by FEM. The spinning flexible disk is described using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. It is discretized by annular sector element. The rotating spindle which includes the clamp, hub, permanent magnet and yoke, is modeled by Timoshenko beam including the gyroscopic effect. The flexible supporting structure with a complex shape which includes stator core, housing, base plate, sleeve and thrust pad is modeled by using a 4-node tetrahedron element with rotational degrees of freedom to satisfy the geometric compatibility. The dynamic coefficients of HDB are calculated from the HDB analysis program, which solves the perturbed Reynolds equation using FEM. Introducing the virtual nodes and the rigid link constraints defined in the center of HDB, beam elements of the shaft are connected to the solid elements of the sleeve and thrust pad through the spring and damper element. The global matrix equation obtained by assembling the finite element equations of each substructure is transformed to the state-space matrix-vector equation, and the associated eigen value problem is solved by using the restarted Arnoldi iteration method. The validity of this research is verified by comparing the numerical results of the natural frequencies with the experimental ones. Also the effect of supporting structures to the natural modes of the total HDD system is rigorously analyzed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.11a
• /
• pp.572-578
• /
• 2003

The free vibration of a spinning flexible disk-spindle system supported by hydro dynamic bearings in a HDD is analyzed by FEM. The spinning flexible disk is described using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. It is discretized by annular sector element. The rotating spindle which includes the clamp, hub, permanent magnet and yoke, is modeled by Timoshenko beam including the gyroscopic effect. The flexible supporting structure with a complex shape which includes stator core, housing, base plate, sleeve and thrust pad is modeled by using a 4-node tetrahedron element with rotational degrees of freedom to satisfy the geometric compatibility. The dynamic coefficients of HDB are calculated from the HDB analysis program, which solves the perturbed Raynolds equation using FEM. Introducing the virtual nodes and the rigid link constraints defined in the center of HDB, beam elements of the shaft are connected to the solid elements of the sleeve and thrust pad through the spring and damper element. The global matrix equation obtained by assembling the finite element equations of each substructure is transformed to the state-space matrix-vector equation, and the associated eigenvalue problem is solved by using the restarted Arnoldi iteration method. The validity of this research is verified by comparing the numerical results of the natural frequencies with the experimental ones. Also the effect of supporting structures to the natural modes of the total HDD system is rigorously analyzed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.889-898
• /
• 2022

In this paper, parallel computing method on the three-dimensional electromagnetic field is proposed. The present electromagnetic scattering analysis is conducted based on the time-harmonic vector wave equation and the finite element method. The edge-based element and 2^nd -order absorbing boundary condition are used. Parallelization of the elemental numerical integration and the matrix assemblage is accomplished by allocating the partitioned finite element subdomain for each processor. The graph partitioning library, METIS, is employed for the subdomain generation. The large sparse matrix computation is conducted by MUMPS, which is the parallel computing library based on the multi-frontal method. The accuracy of the present program is validated by the comparison against the Mie-series analytical solution and the results by ANSYS HFSS. In addition, the scalability is verified by measuring the speed-up in terms of the number of processors used. The present electromagnetic scattering analysis is performed for a perfect electric conductor sphere, isotropic/anisotropic dielectric sphere, and the missile configuration. The algorithm of the present program will be applied to the finite element and tearing method, aiming for the further extended parallel computing performance.