• Journal of the Korean Society of Industry Convergence
• /
• v.2 no.2
• /
• pp.41-49
• /
• 1999

In general, FORTRAN is used for numerical analysis and OPS5 or LISP is used for expert systems, This causes problems at the interface because the various applications require different computing languages or environments. This paper describes the approach used to take AutoCAD as a user-interface for an existing finite element structural analysis package. Some principles concerning database management related to data-exchange interface of pre- and post-processing system for FORTRAN structural analysis program are discussed, and numerical examples demonstrate the power of the combination of these programs.

• Journal of the korean Society of Automotive Engineers
• /
• v.12 no.6
• /
• pp.48-59
• /
• 1990

A finite element program was developed using line elements for simulating the stretch/draw forming operation of an arbitrarily-shaped plane-strain section. An implicit, incremental, updated Lagrangian formulation is employed, introducing a minimum plastic work path assumption for each time step. Geometric and material nonlinearities are also considered within each time step. The finite element equation is based on the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The membrane approximation is adopted under the plane stress assumption. The sheet material is assumed to obey a rigid-viscoplastic constitutive law. The developed program was tested in the die-tryout of typical automotive inner panels. In order to determine a single friction coefficient and boundary length, FEM results and measurements of thinning for a stretched section of final die were compared. After finding analysis parameters, the sheet forming operations of original and final die designs were simulated. Excellent agreement between measured and computed thickness strains was obtained and the developed program was able to identify die designs which were rejected during die tryout.

• Proceeding of EDISON Challenge
• /
• 2015.03a
• /
• pp.285-291
• /
• 2015

In this paper, a composite material analysis program based on the finite element method(FEM) is used. The purpose of this study was to verify whether the composite material analysis program which developed as part of a project of development of softwares and educational contents for structural vibration and composite material analysis that can calculate how similar the macroscopic mechanical behavior of the composite materials actually. Because composite materials are generally anisotropic, analysis of composite laminate is used for the constitutive equations of orthotropic material. For convenience, the unit is ommited in all calculations. To verify the accuracy of the finite element method based program, the deflection and stress distribution of the simply supported composite material laminated plate subjected to a uniform load distribution is compared with exact solution. Size and properties of the composite material laminate used for analysis are fixed variables, and by changing the number of elements and the total thickness of the laminate is compared with the exact solution to the resulting value, respectively.

• Structural Engineering and Mechanics
• /
• v.87 no.3
• /
• pp.243-254
• /
• 2023

Engineering structures usually suffer from cracks. The crack geometry has an influence on the structural mechanical properties and subsequent crack propagations. However, as an extensively utilized method in fracture analysis, the extended finite element method provided by Abaqus fails to output the specific location and dimensions of fractures. In this study, a technique to capture the crack geometry is proposed. The technique is based on the invariant level set method (I-LSM), which can avoid updating the level set function during crack development. The solution is achieved by an open-source plug-in programmed by Python. Three examples were performed to verify the effectiveness and robustness of the program. The result shows that the developed program can accurately output the crack geometry in both the 2D and 3D models. The open-source plug-in codes are included as supplementary material.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.128-136
• /
• 1994

In this study, dynamic boundary element method using mass matrix is derived, using fundamental solutions for the semi-infinite domain. In constituting boundary integral equations for the dynamic equilibrium condition, inertia term in the form of domain integral is transformed into boundary integral form. Corresponding system equations are derived, and a boundary element program is developed. In addition, equations for free vibration is formulated, and eigenvalue analysis is performed. The results from the dynamic boundary element analysis for a tunnel problem are compared with those from the finite element analysis. According to the comparison, boundary element method using mass matrix is consistent with the results of finite element method. Consequently, in tunnel vibration problems, it results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.8
• /
• pp.29-37
• /
• 2022

A three-dimensional analysis is required to interpret the drainage behavior of an improved ground with vertical drains, and the macroelement method enables efficient interpretation considering the three-dimensional drainage effect of vertical drains under two-dimensional plane strain condition. In this study, a novel finite element analysis program was developed by applying the macroelement method to the vacuum consolidation method used in ground improvement practice. The conventional macroelement method was used to calculate the amount of drainage from the vertical drain by setting the excess porewater pressure in the drainage material to zero; however, the program developed in this study was improved to consider negative excess porewater pressure as an actual vacuum consolidation condition. To verify the performance of the program, because of a comparison with the measurement values at the site where the vacuum consolidation method was applied, results predicted by the program and field measurement data showed similar settlement behavior.

• Journal of the korean Society of Automotive Engineers
• /
• v.18 no.4
• /
• pp.50-67
• /
• 1996

• Proceedings of the Korean Geotechical Society Conference
• /
• 1984.06a
• /
• pp.1-23
• /
• 1984

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.288-295
• /
• 2003

An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the recovery technique. In case of the recovery technique, the SPR(superconvergent patch recovery) approach has been modified for p-adaptive mesh refinement. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly. To verify the proposed algorithm, the limit value approach is proposed which utilizes the exact strain energy computed from the extrapolation equation. A new pre-processor is developed for the p-version finite element program in which the vector graphic editor is used for the automatic generation of node connection and coordinate by halfedge solid data structure according to uniform or nonuniform p-distribution. The general 2-D algorithm is also developed to generate face modes and internal modes in accordance with different mesh types. The quality of the error estimator is investigated with the help of two mumerical examples. The results show that the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.