• Computational Structural Engineering
• /
• v.4 no.1
• /
• pp.73-82
• /
• 1991

Recently, the finite element method has been sucessfully extended to treat the rather complex phenomena such as nonlinear buckling problems which are of considerable practical interest. In this study, a finite element program to evaluate the elasto-plastic buckling stress is developed. The Stowell's deformation theory for the plastic buckling of flat plates, which is in good agreement with experimental results, is used to evaluate bending stiffness matrix. A bifurcation analysis is performed to compute the elasto-plastic buckling stress. The subspace iteration method is employed to find the eigenvalues. The results are compared with corresponding exact solutions to the governing equations presented by Stowell and also with experimental data due to Pride. The developed program is applied to obtain elastic and elasto-plastic buckling stresses for various loading cases. The effect of different plate aspect ratio is also investigated.

• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.559-567
• /
• 2012

A level set based topological shape optimization method for nonlinear structure considering hyper-elastic problems is developed. To relieve significant convergence difficulty in topology optimization of nonlinear structure due to inaccurate tangent stiffness which comes from material penalization of whole domain, explicit boundary for exact tangent stiffness is used by taking advantage of level set function for arbitrary boundary shape. For given arbitrary boundary which is represented by level set function, a Delaunay triangulation scheme is used for current structure discretization instead of using implicit fixed grid. The required velocity field in the actual domain to update the level set equation is determined from the descent direction of Lagrangian derived from optimality conditions. The velocity field outside the actual domain is determined through a velocity extension scheme based on the method suggested by Adalsteinsson and Sethian(1999). The topological derivatives are incorporated into the level set based framework to enable to create holes whenever and wherever necessary during the optimization.

• Structural Engineering and Mechanics
• /
• v.89 no.6
• /
• pp.557-570
• /
• 2024

Due to interfacial ageing, chemical action and interfacial damage, the interface debonding may appear in the interfaces of composite laminates. Particularly, the laminates display a side-dependent effect at small scale. In this work, a three-dimensional (3D) and anisotropic thick nanoplate model is proposed to investigate the effects of imperfect interface and nonlocal parameter on the bending deformation, vibrational response and buckling stability of one-dimensional (1D) hexagonal quasicrystal (QC) layered nanoplates. By combining the linear spring model with the transferring matrix method, exact solutions of phonon and phason displacements, phonon and phason stresses of bending deformation, the natural frequencies of vibration and the critical buckling loads of 1D hexagonal QC layered nanoplates are derived with imperfect interfaces and nonlocal effects. Numerical examples are illustrated to demonstrate the effects of the imperfect interface parameter, aspect ratio, thickness, nonlocal parameter, and stacking sequence on the bending deformation, the vibrational response and the critical buckling load of 1D hexagonal QC layered nanoplate. The results indicate that both the interface debonding and nonlocal effect can reduce the stiffness and stability of layered nanoplates. Increasing thickness of QC coatings can enhance the stability of sandwich nanoplates with the perfect interfaces, while it can reduce first and then enhance the stability of sandwich nanoplates with the imperfect interfaces. The biaxial compression easily results in an instability of the QC layered nanoplates compared to uniaxial compression. QC material is suitable for surface layers in layered structures. The mechanical behavior of QC layered nanoplates can be optimized by imposing imperfect interfaces and controlling the stacking sequence artificially. The present solutions are helpful for the various numerical methods, thin nanoplate theories and the optimal design of QC nano-composites in engineering practice with interfacial debonding.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.1
• /
• pp.195-202
• /
• 2004

The shell due to the effect of initial normal pressures on the shell surface was based on the assumption that the directions of the pressures are always normal to the undeformed shell surface, and that the change in the surface area of the shell is negligible. But the fact that the pressure are always normal to the deforming surface leads "follower force". The follower effect in the analysis can significantly alter the solution for natural frequency and buckling load as compared to the case when the direction of the pressures are assumed to be normal to the uniform shell surface. The expression for the part of strain energy contribution from normal pressure due to the effect of follower force was derived and added to the element stiffness matrix of axisymmetric shell. In the case of increasing external pressure, the natural frequencies decrease until one of them reaches zero. Theoretically the smallest applied load that reduces the frequency of any mode to zero, will have same magnitude as that of the buckling load. In order to determine the bucking load of the shell a few sets of frequencies are computed and the results considering the follower effects are well with the exact solution while the case without that are quite different. But in case of hemispherical dome, there are little difference in buckling pressure between with and without the effect of follower force.