• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.9
• /
• pp.1287-1292
• /
• 2010

In this paper, a hybrid technique is presented. First, the isochromatic fringe data of a given set of points are calculated by the finite element method and are used as input data in complex variable formulations. Then the numerical model of the specimen with a central inclined crack is transformed from the physical plane to the complex plane by conformal mapping. The stress field is analyzed and the mixed-mode stress intensity factors are calculated for this complex plane. The stress intensity factors are calculated by the finite element method as well as by a theoretical method and compared with each other. In order to conveniently compare these values with each other, both actual and regenerated photoelastic fringe patterns are multiplied by a factor of two and sharpened by digital image processing.

• Structural Engineering and Mechanics
• /
• v.27 no.4
• /
• pp.409-424
• /
• 2007

In this paper, an identification method of impact force is proposed for composite structures. In this method, the relation between force histories and strain responses is first formulated. The transfer matrix, which relates the strain responses of sensors and impact force information, is constructed from the finite element method (FEM). Based on this relation, an optimization model to minimize the difference between the measured strain responses and numerically evaluated strain responses is built up to obtain the impact force history. The identification of force history is performed by a modified least-squares method that imposes the penalty on the first-order derivative of the force history. Moreover, from the relation of strain responses and force history, an error vector indicating the force location is defined and used for the force location identification. The above theory has also been extended into the cases when using acceleration information instead of strain information. The validity of the present method has been verified through two experimental examples. The obtained results demonstrate that the present approach works very well, even when the internal damages in composites happen due to impact events. Moreover, this method can be used for the real-time health monitoring of composite structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.311-319
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate, hollow cylinder and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.157-162
• /
• 2007

A new level set based topology optimization employing inner-front creation algorithm is presented. In the conventional level set based topology optimization, the optimum topology strongly depends on the initial level set distribution due to the incapability of inner-front creation during optimization process. In the present work, in this regard, an inner-front creation algorithm is proposed. in which the sizes. shapes. positions, and number of new inner-fronts during the optimization process can be globally and consistently identified by considering both the value of a given criterion for inner-front creation and the occupied volume (area) of material domain. To facilitate the inner-front creation process, the inner-front creation map which corresponds to the discrete valued criterion of inner-front creation is applied to the level set function. In order to regularize the design domain during the optimization process, the edge smoothing is carried out by solving the edge smoothing partial differential equation (PDE). Updating the level set function during the optimization process, in the present work, the least-squares finite element method (LSFEM) is employed. As demonstrative examples for the flexibility and usefulness of the proposed method. the level set based topology optimization considering lightweight design of 3D shell structure is carried out.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1335-1340
• /
• 2007

The media transport systems, such as printers, copy machines, facsimiles, ATMs, cameras, etc. have been widely used and being developed rapidly. In the development of those sheet-handling machineries, it is important to predict the static and dynamic behavior of the sheet with a high degree of reliability because the sheets are fed and stacked at such a high speed. Flexible media are very thin, light and flexible, so they behave in geometric nonlinearity with large displacement and large rotation but small strain. In the flexible media analysis, aerodynamic effect from the surrounding air must be included because any small force can make large deformation. In this paper, surrounding air was modeled by incompressible Navier-Stokes flow and an arbitrary Lagranigan-Eulerian(ALE) finite element method with automatic mesh-updating technique was formulated for large domain changes. In the numerical simulations, the results with consideration of the air fast decayed and converged into static results while the results without considering air oscillated continuously.

• Structural Engineering and Mechanics
• /
• v.65 no.5
• /
• pp.547-556
• /
• 2018

Based on the shaking table tests of a 1:3.52 scale one-bay and one-story ancient wooden structure, a simplified structural mechanics model was established, and the structural state equation and observation equation were deduced. Under the action of seismic waves, the damage rule of initial stiffness and yield stiffness of the joint was obtained. The force hammer percussion test and finite element calculations were carried out, and the structural response was obtained. Considering the 5% noise disturbance in the laboratory environment, the stiffness parameters of the mortise-tenon joint were identified by the partial least squares of singular value decomposition (PLS-SVD) and the Extended Kalman filter (EKF) method. The results show that dynamic and static cohesion method, PLS-SVD, and EKF method can be used to identify the damage degree of structures, and the stiffness of the mortise-tenon joints under strong earthquakes is reduced step by step. Using the proposed model, the identified error of the initial stiffness is about 0.58%-1.28%, and the error of the yield stiffness is about 0.44%-1.21%. This method has high accuracy and good applicability for identifying the initial stiffness and yield stiffness of the joints. The identification method and research results can provide a reference for monitoring and evaluating actual engineering structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.249-256
• /
• 2019

This paper presents an optimization framework of electrical impedance tomography for characterizing electrical conductivity profiles of concrete structures in two dimensions. The framework utilizes a partial-differential-equation(PDE)-constrained optimization approach that can obtain the spatial distribution of electrical conductivity using measured electrical potentials from several electrodes located on the boundary of the concrete domain. The forward problem is formulated based on a complete electrode model(CEM) for the electrical potential of a medium due to current input. The CEM consists of a Laplace equation for electrical potential and boundary conditions to represent the current inputs to the electrodes on the surface. To validate the forward solution, electrical potential calculated by the finite element method is compared with that obtained using TCAD software. The PDE-constrained optimization approach seeks the optimal values of electrical conductivity on the domain of investigation while minimizing the Lagrangian function. The Lagrangian consists of least-squares objective functional and regularization terms augmented by the weak imposition of the governing equation and boundary conditions via Lagrange multipliers. Enforcing the stationarity of the Lagrangian leads to the Karush-Kuhn-Tucker condition to obtain an optimal solution for electrical conductivity within the target medium. Numerical inversion results are reported showing the reconstruction of the electrical conductivity profile of a concrete specimen in two dimensions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.41-47
• /
• 2024

This paper utilizes computational asymptotic analysis to compute the boundary layer solution for composite beams and validates the findings through a comparison with ANSYS results. The boundary layer solution, presented as a sum of the interior solution and pure boundary layer effects, necessitates a mathematically rigorous formalization for both interior and boundary layer aspects. Computational asymptotic analysis emerges as a robust technique for addressing such problems. However, the challenge lies in connecting the boundary layer and interior solutions. In this study, we systematically separate the principles of virtual work and the principles of Saint-Venant to tackle internal and boundary layer issues. The boundary layer solution is articulated by calculating the Papkovich-Fadle eigenfunctions, representing them as linear combinations of real and imaginary vectors. To address warping functions in the interior solutions, we employed a least squares method. The computed solutions exhibit excellent agreement with 2D finite element analysis results, both quantitatively and qualitatively. This validates the effectiveness and accuracy of the proposed approach in capturing the behavior of composite beams.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.3
• /
• pp.210-218
• /
• 2018

In order to generate meshes automatically for finite element analysis of complex structures such as aircraft, a large number of triangular elements are typically created. However, triangular elements are less accurate than rectangular elements, so it is difficult to obtain a reliable solution. This problem can be improved through the meshfree method using the back cell integration. However, this method also causes some problems such as over-use of the integration points and inefficiency of the integral domain. In order to improve these problems, a method of performing integration by setting the integral area based on a node basis has been proposed, but in the case of incompressible material problems, the numerical accuracy deteriorates due to the vibration phenomenon of the solution. Therefore, in this paper, the modified meshfree method is proposed which sets the integral domain as an element domain instead of the nodal domain, and the proposed method improves the numerical instability caused by the conventional meshfree method without decreasing the accuracy regardles of the shape of integral domain. The effectiveness of the modified meshfree method is verified by using 2-D examples.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.23 no.3
• /
• pp.313-319
• /
• 2017

Cylindrical shells are often used in ship structures at deck plating with a camber, side shell plating at fore and aft parts, and bilge structure part. It has been believed that such curved shells can be modelled fundamentally by a part of a cylinder under axial compression. From the estimations with the usage of cylinder models, it is known that, in general, curvature increases the buckling strength of a curved shell subjected to axial compression, and that curvature is also expected to increase the ultimate strength. We conduct series of elasto-plastic large deflection analyses in order to clarify the fundamentals in buckling and plastic collapse behaviour of cylindrical shells under axial compression. From the numerical results, we derive design formula for predicting the ultimate strength of cylindrical shell, based on a series of the nonlinear finite element calculations for all edges, simply supporting plating, varying the slenderness ratio, curvature and aspect ratio, as well as the following design formulae for predicting the ultimate strength of cylindrical shell. From a number of analysis results, fitting curve can be developed to use parameter of slenderness ratio with implementation of the method of least squares. The accuracy of design formulae for evaluating ultimate strength has been confirmed by comparing the calculated results with the FE-analysis results and it has a good agreement to predict their ultimate strength.