• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.221-238
• /
• 2008

Numerical solution to buckling analysis of beams and columns are obtained by the method of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for various support conditions considering the variation of flexural rigidity. The solution technique is applied to find the buckling load of fully or partially embedded columns such as piles. A simple semi- inverse method of DQ or HDQ is proposed for determining the flexural rigidities at various sections of non-prismatic column ( pile) partially and fully embedded given the buckling load, buckled shape and sub-grade reaction of the soil. The obtained results are compared with the existing solutions available from other numerical methods and analytical results. In addition, this paper also uses a recently developed technique, known as the differential transformation (DT) to determine the critical buckling load of fully or partially supported heavy prismatic piles as well as fully supported non-prismatic piles. In solving the problem, governing differential equation is converted to algebraic equations using differential transformation methods (DT) which must be solved together with applied boundary conditions. The symbolic programming package, Mathematica is ideally suitable to solve such recursive equations by considering fairly large number of terms.

• Proceedings of the KIEE Conference
• /
• 2000.07e
• /
• pp.99-103
• /
• 2000

The finite element method (FEM) is suitable for the analysis of a complicated region that includes nonlinear materials, whereas the boundary element method (BEM) is naturally effective for analyzing a very large region with linear characteristics. Therefore, considering the advantages in both methods, a novel algorithm for the alternate application of the FEM and BEM to magnetic field problems with the open boundary is presented. This approach avoids the disadvantages of the typical numerical methods with the open boundary problem such as a great number of unknown values for the FEM and non-symmetric matrix for the Hybrid FE-BE method. The solution of the overall problems is obtained by iterative calculations accompanied with the new acceleration method.

• Earthquakes and Structures
• /
• v.13 no.3
• /
• pp.279-288
• /
• 2017

Having the ability to keep on yielding stable solutions in problems involving high potential of instability, composite time integration methods have become very popular among scientists. These methods try to split a time step into multiple sub-steps so that each sub-step can be solved using different time integration methods with different behaviors. This paper proposes a new composite time integration in which a time step is divided into two sub-steps; the first sub-step is solved using the well-known Newmark method and the second sub-step is solved using Simpson's Rule of integration. An unconditional stability region is determined for the constant parameters to be chosen from. Also accuracy analysis is perform on the proposed method and proved that minor period elongation as well as a reasonable amount of numerical dissipation is produced in the responses obtained by the proposed method. Finally, in order to provide a practical assessment of the method, several benchmark problems are solved using the proposed method.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• Communications for Statistical Applications and Methods
• /
• v.12 no.2
• /
• pp.443-451
• /
• 2005

The studied in this paper is a new algorithm for searching the maximum likelihood estimate(MLE) in which probability density function is not explicitly expressed. Newton-Raphson's root-finding routine and a nonlinear numerical optimization algorithm with constraint (so-called feasible sequential quadratic programming) are used. This algorithm is applied to the Wakeby distribution which is importantly used in hydrology and water resource research for analysis of extreme rainfall. The performance comparison between maximum likelihood estimates and method of L-moment estimates (L-ME) is studied by Monte-carlo simulation. The recommended methods are L-ME for up to 300 observations and MLE for over the sample size, respectively. Methods for speeding up the algorithm and for computing variances of estimates are discussed.

• Journal of Mechanical Science and Technology
• /
• v.20 no.6
• /
• pp.782-794
• /
• 2006

Classical finite element programs are not well suited to the design of composite structures, because they are primarily analysis tools and need much time for the data input and as well as for the interpretation of the results. The aim of this paper is to develop a program which allows very fast analyses and reanalyses for design process, thanks to a fast reanalysis method with changes of data and conditions. Speed in the analysis Is obtained by simplification of the analysed structure and limitations in its geometrical generality and improvements in numerical methods. The use of the program is made easy with interactive user-friendly facilities.

• Communications for Statistical Applications and Methods
• /
• v.26 no.3
• /
• pp.315-323
• /
• 2019

Panel data sets have recently been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Often a dichotomous dependent variable occur in survival analysis, biomedical and epidemiological studies that is analyzed by a generalized linear mixed effects model (GLMM). The most common estimation method for the binary panel data may be the maximum likelihood (ML). Many statistical packages provide ML estimates; however, the estimates are computed from numerically approximated likelihood function. For instance, R packages, pglm (Croissant, 2017) approximate the likelihood function by the Gauss-Hermite quadratures, while Rchoice (Sarrias, Journal of Statistical Software, 74, 1-31, 2016) use a Monte Carlo integration method for the approximation. As a result, it can be observed that different packages give different results because of different numerical computation methods. In this note, we discuss the pros and cons of numerical methods compared with the exact computation method.

• Journal of the Architectural Institute of Korea Structure & Construction
• /
• v.34 no.3
• /
• pp.3-10
• /
• 2018

It's not practical to collect all information at the entire degrees of freedom of finite element model. The incomplete measurements should be expanded for subsequent analysis and damage detection. This work presents the analytical methods to expand the incomplete static or dynamic response data. Using the expanded data, introducing the concept of residual force, and minimizing the performance index expressed as the stiffness matrix and its difference before and after damage, the variation in stiffness matrix is derived. Based on the difference in the stiffness matrix, the damage detection method of structures is also provided. The validity of the proposed methods is illustrated in a numerical application, the numerical results are analyzed for applications, and the applicability of both methods is investigated.

• Nonlinear Functional Analysis and Applications
• /
• v.26 no.2
• /
• pp.347-371
• /
• 2021

A plethora of applications from mathematical programmings, such as minimax, mathematical programming, penalization and fixed point problems can be framed as variational inequality problems. Most of the methods that used to solve such problems involve iterative methods, that is why, in this paper, we introduce a new extragradient-like method to solve pseudomonotone variational inequalities in a real Hilbert space. The proposed method has the advantage of a variable step size rule that is updated for each iteration based on previous iterations. The main advantage of this method is that it operates without the previous knowledge of the Lipschitz constants of an operator. A strong convergence theorem for the proposed method is proved by letting the mild conditions on an operator 𝒢. Numerical experiments have been studied in order to validate the numerical performance of the proposed method and to compare it with existing methods.

• Journal of Computing Science and Engineering
• /
• v.3 no.2
• /
• pp.73-87
• /
• 2009

A hybrid system is a dynamical system in which states can be changed continuously and discretely. Simulation based on numerical methods is the widely used technique for analyzing complicated hybrid systems. Numerical simulation of hybrid systems, however, is subject to two types of numerical errors: truncation error and round-off error. The effect of such errors can make an impossible transition step to become possible during simulation, and thus, to generate a simulation behavior that is not allowed by the model. The possibility of an incorrect simulation behavior reduces con.dence in simulation-based analysis since it is impossible to know whether a particular simulation trace is allowed by the model or not. To address this problem, we define the notion of Instrumented Hybrid Automata (IHA), which considers the effect of accumulated numerical errors on discrete transition steps. We then show how to convert Hybrid Automata (HA) to IRA and prove that every simulation behavior of IHA preserves the discrete transition steps of some behavior in HA; that is, simulation of IHA is sound with respect to HA.