• Computational Structural Engineering
• /
• v.8 no.1
• /
• pp.127-136
• /
• 1995

In this paper the stochastic FE analysis considering the material and geometrical property of the plate structure is performed by the weighted integral method. To consider the stochasity of the material and geometrical property, the stochastic field is assumed respectively. The mean value of the stochastic field is 0 and the value of variance is assumed as 0.1. The characteristics of the assumed stochastic field is represented by auto-correlation function. This auto-correlation function is used in evaluating the response variability of the plate structure. In this study a new auto-correlation function is derived to concern the uncertainty of the plate thickness. The newly derived auto-correlation function is a function of auto-correlation function and coefficient of variation of the assumed stochastic field. The two results, obtained by proposed Weighted Integral method and Monte Carlo Simulation method, are coincided with each other and these results are almost equal to the theoretical result that is derived in this study. In the case of considering the variability of plate thickness, the obtained result is well coincide with those of Lawrence and Monte Carlo simulation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.4
• /
• pp.245-253
• /
• 2020

In this study, natural frequencies are compared using several pile-soil interaction (PSI) models to evaluate the effects of each model on resonance safety checks for a monopile type of wind turbine structure. Base spring, distributed spring, and three-dimensional brick-shell models represented the PSIs in the finite element model. To analyze the effects of the PSI models on a natural frequency, after a stiffness matrix calculation and Winkler-based beam model for base spring and distributed spring models were presented, respectively; natural frequencies from these models were investigated for monopiles with different geometries and soil properties. These results were compared with those from the brick-shell model. The results show that differences in the first natural frequency of the monopiles from each model are small when the small diameter of monopile penetrates hard soil and rock, while the distributed spring model can over-estimate the natural frequency for large monopiles installed in weak soil. Thus, an appropriate PSI model for natural frequency analyses should be adopted by considering soil conditions and structure scale.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.54-61
• /
• 1998

This paper presents a prediction method of natural frequencies of vertical hull girder vibration based on design sensitivity analysis in case of design modification and the variation of loading condition. The resented method premises the vibration analysis by the transfer matrix method. Governing sensitivity equation is derived from the direct differentiation of state vector and transfer matrix to parameters and its transfer over all the hull girder elements. Derivatives of natural frequencies and mode shapes are determined by two trial calculation of the governing equation. Using the derivatives, the changes of natural frequencies and mode shapes can be predicted when mass and stiffness parameter's are changed. As results, it is possible to optimize ship structure as well as to avoid troublesome calculation in hull girder vibration analysis rationally and efficiently. To verify the accuracy and efficiency of the resented method numerical results obtained by both the sensitivity analysis and the ordinary reanalysis far a real ore/bulk carrier in case of the change of mass and stiffness parameters are compared.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.719-730
• /
• 2007

This study presents a thin-walled space frame element based on the classical Timoshenko beam theory. The element is derived according to the assumed strain field in order to resolve the shear-locking phenomenon. The shape function is developed in accordance with the strain field which is assumed to be constant at a 2-noded straight frame element. In this study, the geometrically nonlinear analysis applies the Corotational procedure in order to evaluate unbalanced loads. The bowing effect is also considered faithfully. Two numerical examples are given; monosymmetric curved and nonsymmetric straight cantilever. When these example structures behave lateral-torsional bucking, the critical loads are obtained by this study and ABAQUS shell elements. Also, the post-buckling behavior is examined. The results give good agreement between this study and ABAQUS shell.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.18 no.2
• /
• pp.64-73
• /
• 2014

This study proposes a load identification for the safety monitoring of the steel structure based on measured strain data. Instead of parameterizing the stiffness of structure in the existing system identification researches, the loads on a structure and a matrix (the unit strain matrix) defined by the relationship between strain and load on structure are parameterized in this study. The error function is defined by the difference between measured strain and strain estimated by parameters. In order to minimize this error function, the genetic algorithm which is one of the optimization algorithm is applied and the parameters are found. The loads on the structure can be identified through the founded parameters and measured strain data. When the loads are changed, the unmeasured strains are estimated based on founded parameters and measured strains on changed state of structure. To verify the load identification algorithm in this paper, the static experimental test for 3 dimensional steel frame structure was implemented and the loads were exactly identified through the measured strain data. In case of loading changes, the unmeasured strains which are monitoring targets on the structure were estimated in acceptable error range (0.17~3.13%). It is expected that the identification method in this study is applied to the safety monitoring of steel structures more practically.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.249-257
• /
• 2022

This paper presents an improved computational framework for the direct-solution-based finite element tearing and interconnecting (FETI) algorithm. The FETI-local algorithm is further improved herein, and localized Lagrange multipliers are used to define the interface among its subdomains. Selective inverse entry computation, using a property of the Boolean matrix, is employed for the computation of the subdomain interface stiffness and load, in which the original FETI-local algorithm requires a full matrix inverse computation of a high computational cost. In the global interface computation step, the original serial computation is replaced by a parallel multi-frontal method. The performance of the improved FETI-local algorithm was evaluated using a numerical example with 64 million degrees of freedom (DOFs). The computational time was reduced by up to 97.8% compared to that of the original algorithm. In addition, further stable and improved scalability was obtained in terms of a speed-up indicator. Furthermore, a performance comparison was conducted to evaluate the differences between the proposed algorithm and commercial software ANSYS using a large-scale computation with 432 million DOFs. Although ANSYS is superior in terms of computational time, the proposed algorithm has an advantage in terms of the speed-up increase per processor increase.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.95-109
• /
• 1993

In this study, the optimal configuration of arch structure has been tested by a decomposition technique. The object of this study is to provide the method of optimizing the shapes of both two hinged and fixed arches. The problem of optimal configuration of arch structures includes the interaction formulas, the working stress, and the buckling stress constraints on the assumption that arch ribs can be approximated by a finite number of straight members. On the first level, buckling loads are calculated from the relation of the stiffness matrix and the geometric stiffness matrix by using Rayleigh-Ritz method, and the number of the structural analyses can be decreased by approximating member forces through sensitivity analysis using the design space approach. The objective function is formulated as the total weight of the structures, and the constraints are derived by including the working stress, the buckling stress, and the side limit. On the second level, the nodal point coordinates of the arch structures are used as design variables and the objective function has been taken as the weight function. By treating the nodal point coordinates as design variable, the problem of optimization can be reduced to unconstrained optimal design problem which is easy to solve. Numerical comparisons with results which are obtained from numerical tests for several arch structures with various shapes and constraints show that convergence rate is very fast regardless of constraint types and configuration of arch structures. And the optimal configuration or the arch structures obtained in this study is almost the identical one from other results. The total weight could be decreased by 17.7%-91.7% when an optimal configuration is accomplished.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.329-334
• /
• 2017

The structures, which are made up with the huge number of degrees-of-freedom and the assembly of substructures, have a great complexity. In order to increase the computational efficiency, the analysis models have to be simplified. Many substructuring techniques have been developed to simplify large-scale engineering problems. The techniques are very powerful for solving nonlinear problems which require many iterative calculations. In this paper, a modal derivatives-based model order reduction method, which is able to capture the stretching-bending coupling behavior in geometrically nonlinear systems, is adopted and investigated for its performance evaluation. The quadratic terms in nonlinear beam theory, such as Green-Lagrange strains, can be explained by the modal derivatives. They can be obtained by taking the modal directional derivatives of eigenmodes and form the second order terms of modal reduction basis. The method proposed is then applied to a co-rotational finite element formulation that is well-suited for geometrically nonlinear problems. Numerical results reveal that the end-shortening effect is very important, in which a conventional modal reduction method does not work unless the full model is used. It is demonstrated that the modal derivative approach yields the best compromised result and is very promising for substructuring large-scale geometrically nonlinear problems.