• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.102-105
• /
• 2000

Recently, applications of integrated large composite structures have been attempted to many structures of vehicles. To improve the cost performance and reliability of the integrated composite structures, it is necessary to judge structural integrity of the composite structures. For the judgement, we need fracture simulation techniques for composite structures. Many researches oil the fracture simulation method using FEM have been reported by now. Most of the researches carried out simulations considering only matrix cracking and fiber breaking as fracture modes, and did not consider delamination. Several papers have reported the delamination simulation, but all these reports require three-dimensional elements or quasi three- dimensional elements for FEM analysis. Among fracture mechanisms of composite laminates, delamination is the most important factor because it causes stiffness degradation in composite structures. It is known that onset and propagation of delamination are dominated by the strain energy release rate and interfacial moment. In this study, laminated composite has been described by using 3 dimensional finite elements. Then impact behavior of the laminated composite is simulated using FEM(ABAQUS/Explicit) with progressive failure mechanism. These results are compared with experimental results.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.911-914
• /
• 2005

Soft tissue characterization and modeling based on living tissues has been investigated in order to provide a more realistic behavior in a virtual reality based surgical simulation. In this paper, we characterize the nonlinear viscoelastic properties of intra-abdominal organs using the data from in vivo animal experiments and inverse FE parameter estimation algorithm. In the assumptions of quasi-linear-viscoelastic theory, we estimated the nonlinear material parameters to provide a physically based simulation of tissue deformations. To calibrate the parameters to the experimental results, we developed a three dimensional FE model to simulate the forces at the indenter and an optimization program that updates new parameters and runs the simulation iteratively. The comparison between simulation and experimental behavior of pig intra abdominal soft tissue are presented to provide a validness of the tissue model using our approach.

• Wind and Structures
• /
• v.5 no.6
• /
• pp.479-492
• /
• 2002

A finite element aerodynamic model that can be used to analyse flutter instability of long span bridges in the time domain is presented. This approach adopts a simplified quasi-steady formulation of the wind forces neglecting the vortex shedding effects. The governing equations used are effective only for reduced velocities $V^*$ sufficiently great: this is generally acceptable for long-span suspension bridges and, then, the dependence of the wind forces expressions of the flutter derivatives can be neglected. The procedure describes the mechanical response in an accurate way, taking into account the non-linear geometry effects (large displacements and large strains) and considering also the compressed locked coil strands instability. The time-dependence of the inertia force due to fluid structure interaction is not considered. The numerical examples are performed on the three-dimensional finite element model of the Great Belt East Bridge (DK). A mode frequency analysis is carried out to validate the model and the results show good agreement with the experimental measurements of the full bridge aeroelastic model in the wind tunnel tests. Significant parameters affecting bridge response are introduced and accurately investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.7
• /
• pp.1294-1306
• /
• 1992

In order to determine the design parameters and processing conditions in injection molding, it is very important to establish the theoretical model with scientific base. In this study, a two dimensional model has been developed for the purpose and flow simulations of filling process are carried out. The moving boundary transient flow problem along the flat plane is solved efficiently by the Iterative Boundary Pressure Reflection Method which rearranges the impinged melt front along the physical boundary in scientific manner. The two dimensional modeling of filling process is applied to two examples : a three dimensional cover with two screw holes and a two-gated flat cavity with unbalanced runners. The numerical results show good agreement with experimental short shots, especially for the weldline locations and the pressure traces at various locations. They also provide the temperature, clamp force, and velocity field in the mold at different times during filling of cavity.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.437-440
• /
• 2010

This paper represents a simple and effective calculation method to predict the orthotropic engineering constants for C/SiC woven fabric composite. The method, a quasi-analytical method using the modified equivalent laminated model, idealizes the woven fabric structure as a symmetric three-ply laminate to utilize a classical laminated plate theory. The required initial parameters are in-plane modulus from experiments and crimp ratio of the woven fabric. This study shows its feasibility by demonstrating example to calculate the engineering constants to thickness direction needed for three dimensional thermo-mechanical stress calculations.

• Structural Engineering and Mechanics
• /
• v.83 no.4
• /
• pp.435-449
• /
• 2022

Constitutive modeling that could reasonably predict and effectively evaluate the post-peak structural behavior while eliminating the mesh-dependency in numerical simulation remains to be developed for general engineering applications. Based on the previous work, a simple one-dimensional modeling procedure is proposed to predict and evaluate the post-peak response, as characterized by the evolution of localized strain field, of a steel member to monotonically uniaxial tension. The proposed model extends the classic one-dimensional softening with localization model as introduced by (Schreyer and Chen 1986) to account for the localization length, and bifurcation and rupture points. The new findings of this research are as follows. Two types of strain-softening functions (bilinear and nonlinear) are proposed for comparison. The new failure criterion corresponding to the constitutive modeling is formulated based on the engineering strain inside the localization zone at rupture. Furthermore, a new mathematical expression is developed, based on the strain rate inside and outside the localization zone, to describe the displacement field at which bifurcation occurs. The model solutions are compared with the experimental data on four low-carbon cylindrical steel bars of different lengths. For engineering applications, the model solutions are also compared to the experimental data of a cylindrical steel bar system (three steel bars arranged in series). It is shown that the bilinear and nonlinear softening models can predict the energy dissipation in the post-peak regime with an average difference of only 4%.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.11 no.1
• /
• pp.9-16
• /
• 1974

To contribute towards more accurate estimation of the virtual inertia coefficient for the horizontal vibration of ships, three dimensional correction factor $J_H$ for the added mass of finitely long elliptic prismatic bars in horizontal vibration in a free surface of an ideal fluid are calculated. In the problem formulation Dr. T. Kumai's quasi-finite length concept[1,11,12] is employed. Now that, in Dr. Kumai's work[1] for the horizontal vibration the mathematical model was a circular cylinder, the principal aim of the authors' work is to investigate the influence of the beam-draft ratio B/T on $J_H$. The numerical results of this work are shown in Fig.3 graphically, from which we may recognize that the influence of B/T on $J_H$ is remarkable as much as that of the length-draft ratio L/T(refer to Fig.1 also). In Fig.3 the curves for B/T=2.00 are of those based on Dr. Kumai's result[1]. On the other hand, the experimental data obtained by Burril et al.[9] for the horizontal vibration of finitely long prismatic bars of various cross-section shapes are compared with the theoretical added mass coefficients defined by combination of the authors' $J_H$ from Fig.3 and two dimensional coefficients $C_H$ obtained by Lewis form approximation for the corresponding sections. They are in reasonable correspondence with each other as shown in Fig.2. Finally, considering that the longitudinal profile of full-form ship's hull is well resembled to that of an elliptic cylinder and that the influences of other factors such as the sectional area coefficient and the shape of section contour itself can be well merged in the two dimensional added mass coefficient, the authors recommend that the data given in Fig.3 may be successfully adopted for the three dimensional correction factor the added mass in the horizontal vibration of hull-form ships.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2004.05a
• /
• pp.462-465
• /
• 2004

Since quasi-brittle materials like concrete show strain localization behavior accompanied by strain softening, a numerical drawback such as mesh sensitivity is appeared in the finite element analysis. In this paper, the so-called homogenized crack model which was introduced for three dimensional finite element analysis of fracture in concrete is studied for the mesh size dependence problem in fracture analysis. A homogenized crack element having a velocity discontinuity. is averaged to remove the mesh sensitivity in finite element analysis of concrete fracture. Numerical examples show that softening behavior of concrete fracture is successfully predicted without mesh sensitivity using the homogenized crack model.

• Journal of Korea Water Resources Association
• /
• v.33 no.S1
• /
• pp.28-36
• /
• 2000

A TIN, Triagulated Irregular Network, based topographic modeling method and a distributed rainfall-runoff model using the topographic representation is presented. In the TIN based topographic representation, a watershed basin is modeled as a set of contiguous non-overlapping triagular facets : the watershed basin is subdivided according to streamlines to deal with water movement one-dimensionally ; and each partitioned catchment is approximated to a slope element having a quasi-three-dimensional shape by using cubic spline functions. On an approximated slope element, water movement is represented by combined surface-subsurface kinematic wave equations considering a change of slope gradient and slope width. By using the distributed rainfall-runoff model, the effects of spatial variability of soil properties on runoff response are examined.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2000.05a
• /
• pp.28-36
• /
• 2000

A TIN, Triangulated Irregular Network, based topographic modeling method and a distributed rainfall-runoff model using the topographic representation is presented. In the TIN based topographic representation, a watershed basin is modeled as a set of contiguous non-overlapping triangular facets: the watershed basin is subdivided according to streamlines to deal with water movement one-dimensionally; and each partitioned catchment is approximated to a slope element having a quasi-three-dimensional shape by using cubic spline functions. On an approximated slope element, water movement is represented by combined surface-subsurface kinematic wave equations considering a change of slope gradient and slope width. By using the distributed rainfall-runoff model, the effects of spatial variability of soil properties on runoff response are examined.