• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2004.11a
• /
• pp.186-194
• /
• 2004

Experiments were conducted to determine the structural static and dynamic characteristics of air foil bearings. The housing of the bearing on the journal was driven by an impact hammer which was used to simulate dynamic forces acting on the bump loll with various leading condition. Two different bump foils (Cu-coated bump and viscoelastic bump) were tested and the static and dynamic coefficients of two bump foils compared, based on the experimental measurements for a wide range of operating conditions. The static and dynamic characteristics of air foil bearings were extracted 0rpm the frequency response function by least square method and IV(Instrumental Variable) method. The experiment was tested at 0rpm and $10,000\~16,000rpm$, and loaded on $50\~150N$. From the test results, the possibility of the application of high load and high speed condition is suggested.

• Structural Engineering and Mechanics
• /
• v.26 no.4
• /
• pp.441-457
• /
• 2007

A new numerical model based on the spline finite strip method is presented here for the analysis of buckling of built-up columns with and without end stay plates. The channels are modelled with spline finite strips while the connecting elements are represented by a 3D beam finite element, for which the stiffness matrix is modified in order to ensure complete compatibility with the strips. This numerical model has the advantage to give all possible failure modes of built-up columns for different boundary conditions. The end stay plates are also taken into account in this method. To validate the model a comparative study was carried out. First, a general procedure was chosen and adopted. For each numerical analysis, the lowest buckling loads and modes were calculated. The basic or "pure" buckling modes were identified and their critical loads were compared with solutions obtained using analytical methods and/or other numerical methods. The results showed that the proposed numerical model can be used in practice to study the elastic buckling of built-up columns. This model is considered accurate and efficient for the local buckling of short columns and global buckling for slender columns.

• Structural Engineering and Mechanics
• /
• v.58 no.6
• /
• pp.949-966
• /
• 2016

The second order analysis taking place due to non-linear behavior of the structures under the mechanical and geometric factors through implementing exact and approximate methods is an indispensible issue in the analysis of such structures. Among the exact methods is the slope-deflection method that due to its simplicity and efficiency of its relationships has always been in consideration. By solving the differential equations of the modified slope-deflection method in which the effect of axial compressive force is considered, the stiffness matrix including trigonometric entries would be obtained. The complexity of computations with trigonometric functions causes replacement with their Maclaurin expansion. In most cases only the first two terms of this expansion are used but to obtain more accurate results, more elements are needed. In this paper, the effect of utilizing higher order terms of Maclaurin expansion on reducing the number of required elements and attaining more rapid convergence with less error is investigated for the Bernoulli beam with various boundary conditions. The results indicate that when using only one element along the beam length, utilizing higher order terms in Maclaurin expansion would reduce the relative error in determining the critical buckling load and kinematic parameters in the second order analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.169-179
• /
• 2005

An improved method for evaluating effective buckling lengths of beam-column members in plane frames is newly proposed based on system inelastic buckling analysis. To this end, the tangent stiffness matrix of be am-column elements is first calculated using stability functions and then the inelastic buckling analysis method is presented. The scheme for determining effective length of individual members is also addressed. Design examples and numerical results ?uc presented to show the validity of the proposed method.

• Structural Engineering and Mechanics
• /
• v.40 no.1
• /
• pp.29-48
• /
• 2011

Under large storm loads sections of a long pipeline on the seabed can be uplifted. Numerically this loss of contact is extremely difficult to simulate, but accounting for uplift and any subsequent recontact behaviour is a critical component in pipeline on-bottom stability analysis. A simple method numerically accounting for this uplift and reattachment, while utilising efficient force-resultant models, is provided in this paper. While force-resultant models use a plasticity framework to directly relate the resultant forces on a segment of pipe to the corresponding displacement, their historical development has concentrated on precisely modelling increasing capacity with penetration. In this paper, the emphasis is placed on the description of loss of penetration during uplifting, modelled by 'strain-softening' of the force-resultant yield surface. The proposed method employs uplift and reattachment criteria to determine the pipe uplift and recontact. The pipe node is allowed to become free, and therefore, the resistance to the applied hydrodynamic loads to be redistributed along the pipeline. Without these criteria, a localised failure will be produced and the numerical program will terminate due to singular stiffness matrix. The proposed approach is verified with geotechnical centrifuge results. To further demonstrate the practicability of the proposed method, a computational example of a 1245 m long pipeline subjected to a large storm in conditions typical of offshore North-West Australia is discussed.

• Structural Engineering and Mechanics
• /
• v.39 no.1
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

• Structural Engineering and Mechanics
• /
• v.56 no.2
• /
• pp.169-188
• /
• 2015

There are two types of nonlinear analysis methods for building frameworks depending on the method of modeling the plastification of members including lumped plasticity and distributed plasticity. The lumped plasticity method assumes that plasticity is concentrated at a zero-length plastic hinge section at the ends of the elements. The distributed plasticity method discretizes the structural members into many line segments, and further subdivides the cross-section of each segment into a number of finite elements. When a reinforced concrete member experiences inelastic deformations, cracks tend to spread form the joint interface resulting in a curvature distribution. The program IDARC includes a spread plasticity formulation to capture the variation of the section flexibility, and combine them to determine the element stiffness matrix. In this formulation, the flexibility distribution in the structural elements is assumed to be the linear. The main objective of this study is to evaluate the accuracy of linear flexibility distribution assumed in the spread inelasticity model. For this purpose, nonlinear analysis of two reinforced concrete frames is carried out and the linear flexibility models used in the elements are compared with the real ones. It is shown that the linear flexibility distribution is incorrect assumption in cases of significant gravity load effects and can be lead to incorrect nonlinear responses in some situations.

• Tribology and Lubricants
• /
• v.36 no.6
• /
• pp.324-331
• /
• 2020

Wedging in a frictional elastic system is defined if the state of stick exists after the external loading on the system is removed. This paper presents a method to determine the critical coefficient of wedging for an elastic frictional system by considering the separation state. Wedging is always possible if the coefficient of friction exceeds a critical value known as the critical wedging coefficient. This method requires two concepts: a necessary and sufficient condition for wedging, which can be interpreted as positive spanning sets of constraint vectors existing in the wedged system, and the minimal positive basis that enables a minimum wedging coefficient. The algorithm based on the positive spanning concept is repeatedly executed after eliminating nodes from the contact stiffness matrix, for which the separation states are impending. The simulation results show that once a node enters the separation state, it never returns to the contact state again and the critical wedging coefficient reduces during repeated algorithm execution. The benefit of this method is that the computation time permits handling models with large numbers of contact nodes. The algorithm can also numerically find the critical wedging coefficient, thereby contributing to fastening and assembly performance improvements in mechanical systems.

• Advances in nano research
• /
• v.15 no.6
• /
• pp.551-565
• /
• 2023

Coupled porous curved beams, due to their low weight and high flexibility, have many applications in engineering. This study investigates the vibration behavior of coupled porous curved beams in different boundary conditions. The system consists of two curved beams connected by a mid-layer of elastic springs. These beams are made of various materials, such as homogenous steel foam, and composite materials with PMMA (polymethyl methacrylate) and SWCNT (single-walled carbon nanotube) used as the matrix and nanofillers, respectively. To obtain equivalent material properties, the role of mixture (RoM) was employed, followed by the implementation of the porosity function. The system's governing equations were obtained by employing FSDT and Hamilton's law. To investigate thermal vibration, temperature was implemented as a load in the governing equations. The GDQ method was used to solve these equations. To demonstrate the applicability of the GDQ method in calculating the frequencies of the system and the correctness of the developed program, a validation study was conducted. After validation, numerous examples were presented to investigate the behavior of single and coupled curved beams in various material properties and boundary conditions. The results indicate that the frequencies of the curved beams and the system depend highly on the amount of porosity (n) and the distribution pattern. The system frequencies decreased with an increase in the porosity coefficient. The stiffness of the springs had no effect on the first mode frequency but increased frequencies of other modes in a specific range. The frequencies of the system decreased with an increase in environmental temperature.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.3
• /
• pp.439-452
• /
• 1994

A method to optimize the cost of R/C frames and an algorithm of the optimal limit state design for R/C frames subjected to dynamic loads are presented. The modal superposition method was used to find the dynamic responses of the frames. Each member of R/C frame is made up of more than two elements and the stiffness matrix and consistent mass matrix of three d.o.f in the node of each element was used to include axial, shear and flexural effects. The objective function to be minimized formulated the cost of materials, steel and concrete, and optimised to satisfy the behaviors of R/C frame and each constraint imposed by the limit state requirements. Both objective function and each constraint are derived in terms of design variables which include the effective depth, beam width, compression and tension steel area, and column shear steel area. A few applications are presented which demonstrate the feasibility, the validity and efficiency of the algorithm for automated optimum design of R/C frames where dynamic behavior is to be considered.