• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.113-123
• /
• 1990

The paper is concerned with the elasto-plastic and geometrically nonlinear analysis of shell structures using an improved degenerated shell element. In the formulation of the element stiffness, the combined use of three different techniques was made. They are; 1) an enhanced interpolation of transverse shear strains in the natural coordinate system to overcome the shear locking problem ; 2) the reduced integration technique in in-plane strains to avoid the membrane locking behavior ; and 3) selective addition of the nonconforming displacement modes to improve the element performances. This element is free of serious shear/membrane locking problems and undesirable compatible/commutable spurious kinematic deformation modes. In the formulation for plastic deformation, the concept of a layered element model is used and the material is assumed von Mises yield criterion. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements and rotations. The resulting non-linear equilibrium equations are solved by the Netwon-Raphson method combined with load or displacement increment. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Earthquakes and Structures
• /
• v.10 no.5
• /
• pp.1143-1179
• /
• 2016

Offshore wind turbines are considered as a fundamental part to develop substantial, alternative energy sources. In this highly flexible structures, monopiles are usually used as support foundations. Since the monopiles are large diameter (3.5 to 7 m) deep foundations, they result in extremely stiff short monopiles where the slenderness (length to diameter) may range between 5 and 10. Consequently, their elastic deformation patterns under lateral loading differ from those of small diameter monopiles usually employed for supporting structures in offshore oil and gas industry. For this reason, design recommendations (API and DNV) are not appropriate for designing foundations for offshore wind turbine structures as they have been established on the basis of full-scale load tests on long, slender and flexible piles. Furthermore, as these facilities are very sensitive to rotations and dynamic changes in the soil-pile system, the accurate prediction of monopile head displacement and rotation constitutes a design criterion of paramount importance. In this paper, the Fourier Series Aided Finite Element Method (FSAFEM) is employed for the determination of static impedance functions of monopiles for OWT subjected to horizontal force and/or to an overturning moment, where a non-homogeneous soil profile has been considered. On the basis of an extensive parametric study, and in order to address the problem of head stiffness of short monopiles, approximate analytical formulae are obtained for lateral stiffness $K_L$, rotational stiffness $K_R$ and cross coupling stiffness $K_{LR}$ for both rough and smooth interfaces. Theses expressions which depend only on the values of the monopile slenderness $L/D_p$ rather than the relative soil/monopile rigidity $E_p/E_s$ usually found in the offshore platforms designing codes (DNV code for example) have been incorporated in the expressions of the OWT natural frequency of four wind farm sites. Excellent agreement has been found between the computed and the measured natural frequencies.

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

This paper aims at investigating the exact dynamic response of the system undergoing a exponential impact force from the viewpoints of conservations of momentum and energy. The midpoint method applied in the Newmark's family algorithm is found to be identical to the case of the application of the trapezoidal method which provides conservations of momentum and energy. For the linear impact force the mid point, the trapezoidal and the (n+1) point method exactly meet the conservation characteristics independent of the size of integration interval. On the other hand, constants for the dynamic motion resulting from the nonlinear impact are underestimated or overestimated by these method mentioned above. To overcome this indispensible error, the Simpson 1/3 method as one of multi step methods whose advantages is to use longer time interval with the same number of evaluation functions is adopted for the exact conservations of momentum and energy. Moreover, the suggested method is expected to expand the similar algorithm for the general dynamic motion including finite rotations.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.2 s.257
• /
• pp.277-284
• /
• 2007

The steam generator in a nuclear power plant is a large heat exchanger that uses heat from a reactor to generate steam to drive the turbine generator. Rupture of a steam generator tube can result in release of fission products to environment outside. Therefore, an accurate integrity assessment of the steam generator tubes with cracks is of great importance for maintaining the safety of a nuclear power plant. The steam generator tubes are supported at regular intervals by support plates and rotations of the tubes are restrained. Although it has been reported that the limit load for a circumferential crack is significantly affected by boundary condition of the tube, existing limit load solutions do not consider the restraining effect of support plate correctly. In addition, there are no limit load solutions for circumferential cracks in U-bend region with the effect of the support plate. This paper provides detailed limit load solutions for circumferential cracks in top of tube sheet and the U-bend regions of the steam generator tube with the actual boundary conditions to simulate the restraining effect of the support plate. Such solutions are developed based on three dimensional finite element analyses. The resulting limit load solutions are given in a polynomial form, and thus can be simply used in practical integrity assessment of the steam generator tubes.

• Structural Engineering and Mechanics
• /
• v.40 no.5
• /
• pp.637-654
• /
• 2011

The present paper studies the variation of the natural frequencies and mode shapes of rectangular plates carrying a three degree-of-freedom spring-mass system (subsystem), when the subsystem changes (stiffness, mass, moment of inertia, location). An analytical approach based on Lagrange multipliers as well as a finite element formulation are employed and compared. Numerically reliable results are presented for the first time, illustrating the convenience of using the present analytical method which requires only the solution of a linear eigenvalue problem. Results obtained through the variation of the mass, stiffness and moment of inertia of the 3-DOF system can be understood under the effective mass concept or Rayleigh's statement. The analysis of frequency values of the whole system, when the 3-DOF system approaches or moves away from the center, shows that the variations depend on each particular mode of vibration. When the 3-DOF system is placed in the center of the plate, "new" modes are found to be a combination of the subsystem's modes (two rotations, traslation) and the bare plate's modes that possess the same symmetry. This situation no longer exists as the 3-DOF system moves away from the center of the plate, since different bare plate's modes enable distinct motions of the 3-DOF system contributing differently to the "new' modes as its location is modified. Also the natural frequencies of the compound system are nearly uncoupled have been calculated by means of a first order eigenvalue perturbation analysis.

• Journal of KIISE:Software and Applications
• /
• v.29 no.1_2
• /
• pp.20-29
• /
• 2002

The composite stock cutting problem is to allocate rectangular and/or irregular patterns onto a large composite stock sheet of finite dimensions in such a way that the resulting scrap will be minimized. In this paper, the distributed simulated annealing with the new cost error tolerant spatial decomposition is applied to the composite stock cutting problem in MPI environments. The cost error tolerant scheme relaxes synchronization and chooses small perturbations on states asynchronously in a dynamically changed stream length to keep the convergence property of the sequential annealing. This paper proposes the efficient data structures for representation of patterns and their affinity relations and also shows how to determine move generations, annealing parameters, and a cost function. The spatial decomposition method is addressed in detail. This paper identifies that the final quality is not degraded with almost linear speedup. Composite stock shapes are not constrained to convex polygons or even regular shapes, but the rotations are only allowed to 2 or 4 due to its composite nature.

• Structural Engineering and Mechanics
• /
• v.77 no.1
• /
• pp.137-149
• /
• 2021

In engineering design, the axial equivalent elastic modulus of laminated FRP pipe was mostly calculated by the average elastic modulus method or the classical laminated plate theory method, which are based on relatively simplified assumptions, and may be not accurate enough sometimes. A new analytical calculation method for the axial equivalent elastic modulus of laminated FRP pipe was established based on three-dimensional stress state. By comparing the results calculated by this method with those by the above two traditional analytical methods and the finite element method, it is found that this method for the axial equivalent elastic modulus fits well not only for thin-walled pipes with orthotropic layers, but also for thick-walled pipes with arbitrary layers. Besides, the influence of the layer stacking on the axial equivalent elastic modulus was studied with this method. It is found that a proper content of circumferential layer is beneficial for improving the axial equivalent elastic modulus of the laminated FRP pipe with oblique layers, and then can reduce its material quantity under the premise that its axial stiffness remains unchanged. Finally, the meso-mechanical mechanism of this effect was analyzed. The improving effect of circumferential layer on the axial equivalent elastic modulus of the laminated FRP pipe with oblique layers is mainly because that, the circumferential fibers can restrain the rigid body rotations of the oblique fibers, which tend to cause the significant deformations of the pipe wall units and the relatively low axial equivalent elastic modulus of the pipe.

• Steel and Composite Structures
• /
• v.44 no.3
• /
• pp.437-450
• /
• 2022

When the vertical load-bearing members in high-rise structures fail locally, the beam-column joints play an important role in the redistribution of the internal forces. In this paper, a static laboratory test of three full-scale flush flange beam-reinforced connections with side and cover plates (CP-FBSP connection) with double half-span steel beams and single L-shaped columns composed of concrete-filled steel tubes (L-CFST columns) was conducted. The influence of the side plate width and cover plate thickness on the progressive collapse resistance of the substructure was thoroughly analyzed. The failure mode, vertical force-displacement curves, strain variation, reaction force of the pin support and development of internal force in the section with the assumed plastic hinge were discussed. Then, through the verified finite element model, the corresponding analyses of the thickness and length of the side plates, the connecting length between the steel beam flange and cover plate, and the vertical-force eccentricity were carried out. The results show that the failure of all the specimens occurred through the cracking of the beam flange or the cover plate, and the beam chord rotations measured by the test were all greater than 0.085 rad. Increasing the length, thickness and width of the side plates slightly reduced the progressive collapse resistance of the substructures. The vertical-force eccentricity along the beam length reduced the progressive collapse resistance of the substructure. An increase in the connecting length between the beam flange and cover plate can significantly improve the progressive collapse resistance of substructures.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the