• Coupled systems mechanics
• /
• v.13 no.1
• /
• pp.73-93
• /
• 2024

In this work, we propose combining an advanced optimal control algorithm with a geometrically exact beam model. For simplicity, the 2D Reissner beam model is chosen to represent large displacements and rotations. The difficulty pertains to the nonlinear nature of beam kinematics affecting the tangent stiffness matrix, making it non-constant, which compromises direct use of optimal control methods for linear problems. Thus, we seek to accommodate a time varying control using linear-quadratic regulator (LQR) algorithm with the proposed geometrically nonlinear beam model. We provide a detailed theoretical formulation and its numerical implementation in a variational format form. Several illustrative numerical examples are provided to confirm an excellent performance of the proposed methodology.

• Journal of Aerospace System Engineering
• /
• v.17 no.1
• /
• pp.24-32
• /
• 2023

This paper studied HART II in descending flight using rotorcraft analysis code based on geometrically exact beam (GEB) model. The present GEB model expressed by a mixed variational formulation could capture the geometrically nonlinear behavior of the blade without arbitrary assumptions. In previous results, correlation of airloads with structural moments for HART II was not as good as blade deflections. However, in present results, predictions of airloads and structural loads are fairly correlated with measured data.

• Coupled systems mechanics
• /
• v.8 no.6
• /
• pp.537-553
• /
• 2019

In this paper we present geometrically exact Kirchhoff's initially curved planar beam model. The theoretical formulation of the proposed model is based upon Reissner's geometrically exact beam formulation presented in classical works as a starting point, but with imposed Kirchhoff's constraint in the rotated strain measure. Such constraint imposes that shear deformation becomes negligible, and as a result, curvature depends on the second derivative of displacements. The constitutive law is plasticity with linear hardening, defined separately for axial and bending response. We construct discrete approximation by using Hermite's polynomials, for both position vector and displacements, and present the finite element arrays and details of numerical implementation. Several numerical examples are presented in order to illustrate an excellent performance of the proposed beam model.

• Coupled systems mechanics
• /
• v.11 no.6
• /
• pp.525-542
• /
• 2022

In this work we present an advanced approach to the design of small wind turbine support steel structures. To this end we use an improved version of previously developed geometrically exact beam models. Namely, three different geometrically exact beam models are used, the first two are the Reissner and the Kirchhoff beam models implementing bi-linear hardening response and the third is the Reissner beam capable of also representing connections response. All models were validated in our previous research for a static response, and in this work they are extended to dynamic response. With these advanced models, we can perform analysis of four practical solutions for the installation of small wind turbines in new or existing buildings including effects of elastoplastic response to vibration problems. The numerical simulations confirm the robustness of numerical models in analyzing vibration problems and the crucial effects of elastoplastic response in avoiding resonance phenomena.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.4
• /
• pp.183-190
• /
• 2018

This paper presents optimization of the main spar of Human-Powered Aircraft (HPA) wing. Mass minimization was attempted, while considering large torsional deformation of the beam. Sequential Quadratic Programming (SQP) method was adopted as a relevant tool to conduct structural optimization algorithm. An inner diameter and ply thicknesses of the main spar were selected as the design variables. The objective function includes factors such as mass minimization, constant tip bending displacement, and constant tip twist of the beam. For estimation of bending and torsional deformation, the geometrically exact beam model, which is appropriate for large deflection, was adopted. Properties of the cross sectional area which the geometrically exact beam model requires were obtained by Variational Asymptotic Beam Sectional Analysis (VABS), which is a cross sectional analysis program. As a result, maintaining tip bending displacement and tip twist within 1.45%, optimal design that accomplished 7.88% of the mass reduction was acquired. By the stress and strain recovery, structural integrity of the optimal design and validity of the present optimization procedure were authenticated.

• Structural Engineering and Mechanics
• /
• v.22 no.3
• /
• pp.263-278
• /
• 2006

Starting with the geometrically non-linear formulation and the subsequent linearization, this paper presents a consistent formulation of the exact mechanical analysis of geometrically and materially linear three-layer continuous planar beams. Each layer of the beam is described by the geometrically linear beam theory. Constitutive laws of layer materials and relationships between interlayer slips and shear stresses at the interface are assumed to be linear elastic. The formulation is first applied in the analysis of a three-layer simply supported beam. The results are compared to those of Goodman and Popov (1968) and to those obtained from the formulation of the European code for timber structures, Eurocode 5 (1993). Comparisons show that the present and the Goodman and Popov (1968) results agree completely, while the Eurocode 5 (1993) results differ to a certain degree. Next, the analytical solution is used in formulating a general procedure for the analysis of layered continuous beams. The applications show the qualitative and quantitative effects of the layer and the interlayer slip stiffnesses on internal forces, stresses and deflections of composite continuous beams.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.2
• /
• pp.98-106
• /
• 2019

Human-powered aircraft has wings with a shape of high aspect ratio which results in large bending displacement. This paper aims to improve the structural limitation by changing an incidence angle of the wings. The tendency change of bending displacement at the wing tip is observed assuming that airfoil and cross-sectional shape of the wing is fixed, and amount of the total lift generated is satisfied. Quasi-steady lift, drag and the aerodynamic moment are distributed with regard to sections of the wing. Those are analyzed using a numerical nonlinear lifting-line method and 'geometrically exact beam' (GEB) program in EDISON. 'Variational Asymptotic Beam Sectional Analysis' (VABS) program is used to check if the present wing is structurally solid. Furthermore, the predicted tip deflections are verified by comparing with DYMORE.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.4
• /
• pp.7-22
• /
• 1997

A finite element analysis is performed for lateral buckling problems on the basis of a geometrically nonlinear formulation for a beam with small elastic strain but with possibly large rotations. The total Lagrangian formulation for a general large deformation, which involves finite rotations, is chosen and the exponential map is used to treat finite rotations from the Eulerian point of view. For lateral buckling, the point of vanishing determinant of the resulting unsymmetric tangent stiffness is traced to examine its relationship to bifurcation points. It is found that the points of vanishing determinant is not corresponding to bifurcation points for large deformations in general, which suggests that the present unsymmetric tangent stiffness is not an exact first derivative of internal forces with respect to displacement. This is illustrated through several numerical examples and followed by appropriate discussion.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.167-198
• /
• 2022

In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.5
• /
• pp.287-292
• /
• 2021

We utilized the isogeometric analysis (IGA) method that uses NURBS basis functions in CAD systems, to account for the geometric exactness of a geometrically exact beam deformation, on a new type of metamaterial, twist-translation coupled structure showing a large twist angle. A two-dimensional unit cell structure was embedded in a cylindrical wall, using free-form deformation and an appropriate interpolation scheme. A parametric study on the effects of the dimensions of the cylinder and the number of cells, on the twisting angle was performed. Furthermore, the mechanism of the twist-translation coupled metamaterial was explored through numerical examples.