• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary element.

• Journal of Korean Association for Spatial Structures
• /
• v.12 no.1
• /
• pp.109-119
• /
• 2012

For membrane structure, there are three main steps in design and construction, which are form finding, statistical load analysis, and cutting patterning. Unlike the first two stages, the step of cutting pattern involves the translation of a double-curved surface in 3D space into a 2D plane with minimal error. For economic reasons, the seam lines of generated cutting patterns rely greatly on the geodesic line. Generally, as searching regions of the seam line are plane elements in the step of shape analysis, the seam line is not a smooth curve, but an irregularly divided straight line. So, it is how we make an irregularly divided straight line a smooth curve that defines the quality of the pattern. Accordingly, in this paper, we analyzed interpolation schemes using spline, and apply these methods to cutting pattern generation on the curved surface. To generate the pattern, three types of spline functions were used, i.e., cubic spline function, B-spline, and least-square spline approximation, and simple model and the catenary-shaped membrane was adopted to examine the result of generation. The result of comparing the approximation curves by the number of elements and the number of extracted nodes of simple model revealed that the seam line for less number of extracted nodes with large number of elements were more efficient, and the least-square spline approximation provided smoother seam line than other methods.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.1
• /
• pp.37-44
• /
• 2019

A bond-based peridynamic model has been reported dynamic fracture characteristic of brittle materials through a simple constitutive model. In the model, each bond is assumed to be a simple spring operating independently. As a result, this simple bond interaction modeling restricts the material behavior having a fixed Poisson's ratio of 1/4 and not being capable of expressing shear deformation. We consider a state-based peridynamics as a generalized peridynamic model. Constitutive models in the state-based peridynamics are corresponding to those in continuum theory. In state-based peridynamics, thus, the response of a material particle depends collectively on deformation of all bonds connected to other particles. So, a state-based peridynamic theory can represent the volume and shear changes of the material. In this paper, the perfect plasticity is considered to express plastic deformation of material by the state-based peridynamic constitutive model with perfect plastic flow rule. The elastic-plastic behavior of the material is verified through the stress-strain curves of the flat plate example. Furthermore, we simulate the high-speed impact on 3D granite model with a nonlocal contact modeling. It is observed that the damage patterns obtained by peridynamics are similar to experimental observations.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.4
• /
• pp.221-234
• /
• 1999

This paper deals with the application of joint element in the finite element modeling of discontinuities encountered during rock tunneling. A nodal displacement joint element was implemented in a two dimensional finite element program GEOFE2D. The applicability of the joint element for modeling of discontinuities and the numerical stability of the implemented algorithm were examined by comparing the results of reduced small scale model tests as well as commercially available FEM program. The GEOFE2D was then used to analyze a tunnel crossed by a major discontinuity for the purpose of understanding the effect of discontinuity on the tunnel behavior. In addition, a modeling technique for the junction of discontinuity and shotcrete lining was presented. The results of analysis indicated that the stress-strain field around the tunnel is significantly altered by the presence of discontinuity, and that the stresses in the shotcrete lining considerably increase at the junction of the shotcrete lining and the discontinuity. It is therefore concluded that the major discontinuities must be carefully modeled in the finite element analysis of a tunneling problem in order to obtain more reliable results close to actual tunnel behavior.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.641-652
• /
• 2007

In this paper, a half-scaled substructure test was performed to evaluate the buckling and structural safety of an existing transmission tower subjected to wind load. A loading scheme was devised to reproduce the dead and wind loads of a prototype transmission tower, which uses a triangular jig that is mounted on the reduced model to which the similarity law of a half length was applied. As a result of the preliminary numerical analysis carried out to evaluate the stability of a specimen for the design load, is was confirmed that the calculated axial forces of tower leg members were distributed to $80{\sim}90%$ of an admissible buckling load. When the substructured transmission tower was loaded by 270% of its maximum admissible buckling load, it was failed due to the local buckling that is occurred in joints with weak constraints for out-of-plane behavior of leg members. By inspection of load-displacement curves, displacements and strains of members, it is considered that this local buckling was due to additional eccentric force by unbalanced deformation because the time that is reached to yielding stress due to the bending moment is different at each point of a same section.

• Geotechnical Engineering
• /
• v.9 no.1
• /
• pp.21-30
• /
• 1993

To know the importance of soil paramters which are used to estimate the deformation and porepressure of soil, the sensitivity for soil parameters in elastic analysis is analyzed. Using the consolidation teat results of several cohesive soils, soil parameters are estimated by back analysis method, and from the parameters the deformations and porepressures of the soil are estimated by elastic analysis, In elastic analysis for soil-deformation and porepressure, the sensitivity for the Young's modulus is large, and the esimation of Young's modulus is more important in pro- portion to the size of stress. Using the measured results during initial short period in small stress, the soil parameters can be correctly estimated by back analysis method. To decrease the iteration number in back analysis and to get the better paramters, the initial measurements in more nodes are required and the more accurate initial measurements are required.

• Proceeding of KASS Symposium
• /
• 2008.05a
• /
• pp.107-112
• /
• 2008

This paper investigates the variation of post-buckling behaviours of 2-dimensional shallow arch type truss after sizing optimization. The mathematical programming technique is used to produce the optimum member size of 2D arch truss against a central point load. Total weight of structure is considered as the objective function to be minimized and the displacement occurred at loading point and member stresses of truss are used as the constraint functions. The finite difference method is used to calculate the design sensitivity of objective function with respect to design variables. The postbuckling analysis carried out by using the geometrically nonlinear finite element analysis code ISADO-GN. It is found to be that there is a huge change of post-buckling behaviour between the initial structure and optimum structure. Numerical results can be used as useful information for future research of large spatial structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.3
• /
• pp.141-154
• /
• 1987

The algorithm proposed utilizes the tow-levels technique. In the first level which consists of teeatment only the applied load and design stress as the random variables whose parent distribution has the normal distribution, the cross-sectional areas of the truss members such that the their probabilities of failure have the preseribed failure probabilites are optimized by transforming the nonlinear problem into SUMT, and solving it utilizing modified Newton-Raphson method. In the second level, the geometric shape of truss structure is optimized by utilizing the unidirectional search technique of Powell method which makes it possible to minimize only the objective function. The algorithm proposed is numerically tested for the several truss structures with various shapes and loading conditions. The numerical analysis shows that the rate of decreasing the weight of truss structures is dependent on the prescribed failure probability of the each member of truss structure and the covariance of the applied load and design stress.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.5
• /
• pp.279-286
• /
• 2020

In this study, we investigate the accuracy of higher order derivatives in the moving least square (MLS) difference method. An interpolation function is constructed by employing a Taylor series expansion via MLS approximation. The function is then applied to the mixed variational theorem in which the displacement and stress resultants are treated as independent variables. The higher order derivatives are evaluated by solving simply supported beams and cantilevers. The results are compared with the analytical solutions in terms of the order of polynomials, support size of the weighting function, and number of nodes. The accuracy of the higher order derivatives improves with the employment of the mean value theorem, especially for very high-order derivatives (e.g., above fourth-order derivatives), which are important in a classical asymptotic analysis.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2008.04a
• /
• pp.249-252
• /
• 2008

A study for the nonlinear analysis method of RC(Reinforced Concrete) beams with unbonded prestressing CFRP plate is presented. The cross-section of RC beam element is idealized as an assemblage of concrete and reinforcing steel fibers in order to account for varied material properties within the cross-section of the element. The unbonded CFRP plate is modeled as a series of the CFRP plate segments each of which is linked to the RC beam element, but slips without any resistance to simulate the unbonded behavior of the CFRP plate. The stress of each CFRP plate segment is redistributed iteratively using the force equilibrium relationship at each common node until it reaches at the same stress level. To evaluate the validity of the proposed analysis method, the results of ultimate analysis of the reinforced concrete beams with unbonded prestressing CFRP plates are compared with the experimental results obtained from other investigators. The proposed analysis method is found to predict ultimate behaviors of these beams fairly well.