• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.33-45
• /
• 1996

A general stiffener element which includes transverse shear deformation is formulated using the p-version finite element method. Hierarchic C/sup o/-shape functions, derived from Integrals of Legendre polynomials, are used to define the assembled stiffness matrix of the stiffener with respect to the local reference frame is transformed to the plate reference system by applying the appropriate transformation matrices in order to insure compatibility of displacements at the junction of the stiffener and plate. The transformation matrices which account for the orientation and the eccentricity effects of the stiffener with respect to the plate reference axes are used to find local behavior at the junction of the stiffener and the relative contributions of the plate and stiffener to the strength of the composite system. The results obtained by the p-version finite element method are comared with the results in literatures, especially those by the h-version finite element analysis program, MICROFEAP-II.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.6
• /
• pp.794-801
• /
• 2012

Finite Element Method is a strong tool to analyse static and dynamic problem of a structure. FEM is a good method for static problem, but for dynamic problem there are some differences between real phenomena and analyzed phenomena. Therefore some modifications are needed to identify two results. In this paper authors propose a genetic algorithm method 1) to adjust dimensions of plate for identifying natural frequencies, 2) to fit amplitude of FEM Frequency Response Function(FRF) onto it of real FRF. Analysis by raw FEM data gave questions if the results were for the same object. By only adjusting Young's modulus much better accordances were obtained, but limitation existed still. Very good agreements were achieved by shape modification and damping coefficient identification.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1289-1298
• /
• 1994

The new p-version crack model is proposed to estimate the stress intensity factors of the thick cracked plate under flexure. The proposed model is based on high order theory and $C^{\circ}$-plate element including shear deformation. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three groups such as basic mode, side mode and internal mode. The computer implementation allows arbitrary variations of p-level Up to a maximum value of 10. The stress intensity factors are computed by virtual crack extention approach. The effects of ratios of thickness to crack length(h/a), crack length to width(a/W) and boundary conditions are investigated. Very good agreement with the existing solution in the literature are shown for the uncracked plate as well as the cracked plate.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.12
• /
• pp.1165-1172
• /
• 2009

The laminar boundary layer along a streamwise corner formed by two flat plates intersecting at right angle is measured by using Particle Image Velocimetry(PIV) technique. The free stream velocity ranges from 2.96m/s to 3.0m/s. The angle of incidence of the corner is set to 1.2 degree providing slightly favourable pressure gradient to ensure a laminar flow in the corner region. A round shape leading edge is used and the length of the model is about 1000mm. In the bisector plane, the measurement data show separation type velocity profiles having an inflection point which is a typical characteristic of laminar corner boundary layers. As the distance away from the bisector plane increases, velocity profiles are found to change into the Blasius profile. The change completes around half length of the boundary layer thickness in the bisector plane away from the bisector plane along the plate. In the bisector plane, the growth characteristic of the boundary layer thickness and the approximate similarity of velocity profiles are confirmed from the measurement data.

• Journal of Advanced Marine Engineering and Technology
• /
• v.37 no.4
• /
• pp.316-323
• /
• 2013

The use of artificial roughness in various forms of shapes and sizes is the most common and effective way to improve the performance of a flat-plate solar collector. In the present study, numerical analysis on heat transfer and pressure drop was performed in a rectangular channel with various rib arrays. The uniform heat flux is applied to the channel from the upper side. The forms of ribs considered in this study were rib $90^{\circ}$, groove $90^{\circ}$, groove $60^{\circ}$, baffle $90^{\circ}$, baffle $60^{\circ}$, wave $90^{\circ}$ and wave $60^{\circ}$. Air is the working fluid, and the Reynolds number ranges from 3200 to 17800. Nusselt number and friction factor were investigated to predict the performance of the system with various type of ribs. The average Nusselt number and pressure drop were increased with the increase of velocity in all types of ribs. The highest heat transfer and pressure drop occurred for the baffle $90^{\circ}$, but highest performance factor considering heat transfer and pressure drop together occurred for the groove $60^{\circ}$. Therefore, heat transfer and pressure drop should be considered together when a flat plate solar collector is designed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.8
• /
• pp.3295-3300
• /
• 2012

Three Dimensional finite element method (FEM) was used to obtain the stress intensity factor for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model. Nodes are generated by the bubble packing, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. To examine accuracy and efficiency of the present system, the stress intensity factor for a semi-elliptical surface crack in a plate subjected to uniform tension is calculated, and compared with Raju-Newman's solutions. Then the system is applied to analyze interaction effects of two dissimilar semi-elliptical cracks in a plate subjected to uniform tension.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.189-196
• /
• 2009

This paper presents the subparametric finite element model formulated by partial-linear layerwise theory for the analysis of laminate composites. The proposed model is based on refined approximations of two dimensional plane for orthotropic thick laminate plate as well as thin case. Three dimensional problem can be reduced to two dimensional case by assuming piecewise linear variation of in-plane displacement and a constant value of out-of-plane displacement across the thickness. The integrals of Legendre polynomials are chosen to define displacement fields and Gauss-Lobatto numerical integration is implemented in order to directly obtain maximum values occurred at the nodal points of each layer without other extrapolation techniques. The validity and characteristics of the proposed model have been tested by using orthotropic multilayered plate problem as compared to the values available in the published references. In this study, the convergence test has been carried out to determine the optimal layer model in terms of central deflection and stresses. Also, the distribution of displacements and stresses across the thickness has been investigated as the number of layer is increased.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.1
• /
• pp.1-11
• /
• 2002

The suitable turbulence model is found to be required in the course of establishing a proper analysis methodology for thermal stripping phenomena which are shown in strong temperature variation area such as reactors and propulsion devices. Three different turbulence models of $\kappa$-$\varepsilon$ model, modified $\kappa$-$\varepsilon$ model, and full Reynolds stress(FRS) model, are applied to analyze unsteady turbulent flows with temperature variation. Three test cases are selected for verification. These are vertical jet flows with water and sodium, and parallel jet flow with sodium. Analysis yields the conclusion that 3-D computation with FRS betters others. However, modified modeling is required to improve its heat transfer characteristic analysis. Further analysis is performed to find momentum variation effects on temperature distribution. It is found that the momentum increase results increase of fluid mixing and magnitude of temperature variation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.3
• /
• pp.291-299
• /
• 2010

Aspect-ratio-based numerical analysis is carried out to investigate the air-side convective heat transfer characteristics in compact heat exchangers with circular tubes and flat tubes-plate fins. The RNG $k-{\varepsilon}$ model is adopted for turbulence analysis. The numerical analysis is carried out for aspect ratios ranging from 3.06 to 5.44 and for Reynolds numbers ranging from 1,000 to 10,000. The calculated results indicate a correlation between the friction factor and Colburn j factor in the compact heat exchanger system for the range of aspect ratios under consideration. The results obtained for circular tubes and flat tubes-plate fins in this study can be utilized to realize the optimal design of an air conditioning system.

• Composites Research
• /
• v.18 no.6
• /
• pp.9-18
• /
• 2005

With the wide applications of fiber-reinforced composite material in aero-structures and mechanical parts, the design of composite joints have become a very important research area because the joints are often the weakest areas in composite structures. This paper presents an analytical study of the stress distributions in mechanically single-fastened and multi-fastened composite laminates. The finite element models which treat the pin and hole contact problem using a contact stress analysis are described. A dimensionless stress concentration factor is used to compare the stress distributions in composite laminates quantitatively In the case of single-pin loaded composite laminate, the effects of stacking sequence, the ratio of a hole diameter and the width of a laminate (W/D ratio), the ratio of hole diameter and distance from edge to hole (E/D ratio), friction coefficient and clamping force are considered. In the case of multi-pin loaded composite laminate, the influence of the number of pins, pitch distance, number of rows, row spacing and hole pattern are considered. The results show that P/D ratio and E/D ratio affect more on stress distributions near the hole boundary than the other factors. In the case of multi-pin loaded composite laminate, the stress concentration in the double column case is better than the other cases of multi-pin loaded composite laminate.