• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• 제36권1호
• /
• pp.22-30
• /
• 2008

This paper deals with the nonlinear optimal control approach to helicopter maneuver problems using the indirect method. We apply a penalty function to the deviation from a prescribed trajectory to convert the system optimality to an unconstrained optimal control problem. The resultant two-point boundary value problem has been solved by using the multiple-shooting method. This paper focuses on the effect of the number of shooting nodes and initialization methods on the numerical solution in order to define the minimum number of shooting nodes required for numerical convergence and to provide a method increasing convergence radius of the indirect method. The results of this study can provide an approach to improve numerical stability and convergence of the indirect method when we solve the optimal control problems of an inherently unstable helicopter system.

• Journal of the Korean Mathematical Society
• /
• 제57권3호
• /
• pp.747-775
• /
• 2020

In this work we investigate the nonlocal elliptic equation with critical Hardy-Sobolev exponents as follows, $$(P)\;\{(-{\Delta}_p)^su={\lambda}{\mid}u{\mid}^{q-2}u+{\frac{{\mid}u{\mid}^{p{^*_s}(t)-2}u}{{\mid}x{\mid}^t}}{\hspace{10}}in\;{\Omega},\\u=0{\hspace{217}}in\;{\mathbb{R}}^N{\backslash}{\Omega},$$ where Ω ⊂ ℝ^N is an open bounded domain with Lipschitz boundary, 0 < s < 1, λ > 0 is a parameter, 0 < t < sp < N, 1 < q < p < p^∗_s where $p^*_s={\frac{N_p}{N-sp}}$, $p^*_s(t)={\frac{p(N-t)}{N-sp}}$, are the fractional critical Sobolev and Hardy-Sobolev exponents respectively. The fractional p-laplacian (-∆_p)^su with s ∈ (0, 1) is the nonlinear nonlocal operator defined on smooth functions by $\displaystyle(-{\Delta}_p)^su(x)=2{\lim_{{\epsilon}{\searrow}0}}\int{_{{\mathbb{R}}^N{\backslash}{B_{\epsilon}}}}\;\frac{{\mid}u(x)-u(y){\mid}^{p-2}(u(x)-u(y))}{{\mid}x-y{\mid}^{N+ps}}dy$, x ∈ ℝ^N. The main goal of this work is to show how the usual variational methods and some analysis techniques can be extended to deal with nonlocal problems involving Sobolev and Hardy nonlinearities. We also prove that for some α ∈ (0, 1), the weak solution to the problem (P) is in C^1,α(${\bar{\Omega}}$).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제34권3호
• /
• pp.151-157
• /
• 2021

This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• 제30권6호
• /
• pp.270-285
• /
• 2018

The behavior of wave-induced pore water pressure inside the rubble mound and seabed, and the resultant structure failure are investigated, which are used in design of the composite breakwater representing the coastal and harbor structures. Numerical simulation techniques have been widely used to assess these behaviors through linear and nonlinear methods in many researches. While the combination of strongly nonlinear analytical method and turbulence model have not been applied yet, which can simulate these characteristics more accurately. In this study, olaFlow model considering the wave-breaking and turbulent phenomena is applied through VOF and LES methods, which gives more exact solution by using the multiphase flow analytical method. The verification of olaFlow model is demonstrated by comparing the experimental and numerical results for the interactions of regular waves-seabed and regular waves-composite breakwater-seabed. The characteristics of the spatial distributions of horizontal wave pressure, excess-pore-water pressure, mean flow velocity and mean vorticity on the upright caisson, and inside the rubble mound and seabed are discussed, as well as the relation between the mean distribution of vorticity size and mean turbulent kinetic energy. And the stability of composite breakwater are also discussed.

• Journal of the Earthquake Engineering Society of Korea
• /
• 제12권2호
• /
• pp.69-80
• /
• 2008

In this study, improved capacity spectrum method (CSM) is proposed. The method can account for higher mode contribution to the seismic response of MDOF systems. The CSM has been conveniently used for determining maximum roof displacement using both demand spectrum and capacity curve of equivalent SDOF system. Unlike the conventional CSM, the maximum roof displacement is determined without iteration using inelastic displacement ratio and R factor calculated from demand spectrum and capacity curve. Three moment resisting steel frames of 3-, 9- and 20-stories are considered to test the accuracy of the proposed method. Nonlinear response history analysis (NL-RHA) for three frames is also conducted, which is considered as an exact solution. SAC LA 10/50 and 2/50 sets of ground motions are used. Moreover, this study estimates maximum story drift ratios (IDR) using ATC-40 CSM and N2-method and compared with those from the proposed method and NL-RHA. It shows that the proposed CSM estimates the maximum IDR accurately better than the previous methods.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• 제48권9호
• /
• pp.691-701
• /
• 2020

An actuator mixer design using convex optimization technique situation where the propulsion system of an electric VTOL UAV during vertical take-off and landing maneuvers is proposed. The attainable control set to analyze the impact from failure of each motor and propeller can be calculated and illustrated using the properties of the convex function. The control allocation can be defined as a convex function optimization problem to obtain an optimal solution in real time. The mixer is implemented using a convex optimization solver, and the performance of the control allocation methods is compared to the attainable control set. Finally, the proposed mixer is compared with other techniques with nonlinear sux degree-of-freedom simulation.

• Structural Engineering and Mechanics
• /
• 제40권1호
• /
• pp.85-103
• /
• 2011

Integration of finite element analysis (FEA) software into various software platforms is commonly used in coupling systems such as systems involving structural control, fluid-structure, wind-structure, soil-structure interactions and substructure method in which FEA is used for simulating the structural responses. Integrating an FEA program into various other software platforms in an efficient and simple way is crucial for the development and performance of the entire coupling system. The lack of simplicity of the existing integration methods makes this integration difficult and therefore entails the motivation of this study. In this paper, a novel practical technique, namely CS technique, is presented for integrating a general FEA software framework OpenSees into other software platforms, e.g., Matlab-$Simulink^{(R)}$ and a soil-structure interaction (SSI) system. The advantage of this integration technique is that it is efficient and relatively easy to implement. Instead of OpenSees, a cheap client handling TCL is integrated into the other software. The integration is achieved by extending the concept of internet based client-server concept, taking advantage of the parameterization framework of OpenSees, and using a command-driven scripting language called tool command language (TCL) on which the OpenSees' interface is based. There is no need for any programming inside OpenSees. The presented CS technique proves as an excellent solution for the coupling problems mentioned above (for both linear and nonlinear problems). Application examples are provided to validate the integration method and illustrate the various uses of the method in the civil engineering.

• Computational Structural Engineering
• /
• 제7권4호
• /
• pp.57-67
• /
• 1994

It is usual that underground structures are constructed within a multi-layered medium. In this paper, an efficient numerical modelling of multi-layered structural systems is studied using coupled analysis of finite elements and boundary elements. The finite elements are applied to the area in which the material nonlinearity dominates, and the boundary elements are applied to the far field where the nonlinearity is relatively weak. In the boundary element modelling of the multi-layered medium, fundamental solutions are not readily available. Thus, methods which can utilize existing Kelvin solutions are sought for the interior multi-layered domain problem. The interior domain problem which has piecewise homogeneous layers is analyzed using boundary elements with Kelvin solution, by discretizing each homogeneous subdomain and enforcing compatibility and equilibrium conditions between interfaces. Developed methodology is verified by comparing its results with those from the finite element analysis and it is concluded that coupled analysis using boundary elements and finite elements can be reasonable and efficient.

• Advances in environmental research
• /
• 제3권3호
• /
• pp.231-251
• /
• 2014

The present study examines the phosphate adsorption potential and behavior of mixture of Ground Burnt Patties (GBP), a solid waste generated from cooking fuel used in earthen stoves and Red Soil (RS), a natural substance in fixed bed column mode operation. The characterization of adsorbent was done by Proton Induced X-ray Emission (PIXE), and Proton Induced ${\gamma}$-ray Emission (PIGE) methods. The FTIR spectroscopy of spent adsorbent reveals the presence of absorbance peak at $1127cm^{-1}$ which appears due to P = O stretching, thus confirming phosphate adsorption. The effects of bed height (10, 15 and 20 cm), flow rate (2.5, 5 and 7.5 mL/min) and initial phosphate concentration (5 and 15 mg/L) on breakthrough curves were explored. Both the breakthrough and exhaustion time increased with increase in bed depth, decrease in flow rate and influent concentration. Thomas model, Yoon-Nelson model and Modified Dose Response model were used to fit the column adsorption data using nonlinear regression analysis while Bed Depth Service Time model followed linear regression analysis under different experimental condition to evaluate model parameters that are useful in scale up of the process. The values of correlation coefficient ($R^2$) and the Sum of Square Error (SSE) revealed the Modified Dose Response model as the best fitted model to the experimental data. The adsorbent mixture responded effectively to the desorption and reusability experiment. The results of this finding advocated that mixture of GBP and RS can be used as a low cost, highly efficient adsorbent for phosphate removal from aqueous solution.