• Structural Engineering and Mechanics
• /
• 제39권5호
• /
• pp.669-682
• /
• 2011

A Chebyshev spectral method (CSM) for the dynamic analysis of non-uniform Timoshenko beams under various boundary conditions and concentrated masses at their ends is proposed. The matrix-based Chebyshev spectral approach was used to construct the spectral differentiation matrix of the governing differential operator and its boundary conditions. A matrix condensation approach is crucially presented to impose boundary conditions involving the homogeneous Cauchy conditions and boundary conditions containing eigenvalues. By taking advantage of the standard powerful algorithms for solving matrix eigenvalue and generalized eigenvalue problems that are embodied in the MATLAB commands, chebfun and eigs, the modal parameters of non-uniform Timoshenko beams under various boundary conditions can be obtained from the eigensolutions of the corresponding linear differential operators. Some numerical examples are presented to compare the results herein with those obtained elsewhere, and to illustrate the accuracy and effectiveness of this method.

• Journal of applied mathematics & informatics
• /
• 제29권5_6호
• /
• pp.1205-1219
• /
• 2011

In this work, a stabilized characteristic finite volume method for the time-dependent Navier-Stokes equations is investigated based on the lowest equal-order finite element pair. The temporal differentiation and advection term are dealt with by characteristic scheme. Stability of the numerical solution is derived under some regularity assumptions. Optimal error estimates of the velocity and pressure are obtained by using the relationship between the finite volume and finite element methods.

• Journal of Korean Society of Steel Construction
• /
• 제12권5호통권48호
• /
• pp.569-579
• /
• 2000

An improved multi-level (IML) optimization algorithm using automatic differentiation (AD) of framed structures is proposed in this paper. For the efficiency of the proposed algorithm, multi-level optimization techniques using a decomposition method that separates both system-level and element-level optimizations, that utilizes and an artificial constraint deletion technique, are incorporated in the algorithm. And also to save the numerical efforts, an efficient reanalysis technique through approximated structural responses such as moments and frequencies with respect to intermediate variables is proposed in the paper. Sensitivity analysis of dynamic structural response is executed by AD that is a powerful technique for computing complex or implicit derivatives accurately and efficiently with minimal human effort. The efficiency and robustness of the IML algorithm, compared with a plain multi-level (PML) algorithm, is successfully demonstrated in the numerical examples.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제20권3호
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• Structural Engineering and Mechanics
• /
• 제44권3호
• /
• pp.267-288
• /
• 2012

Free vibration analysis of arbitrary quadrilateral thick plates with internal columns and elastic edge supports is presented by using the powerful pb-2 Ritz method and Reddy's third order shear deformation plate theory. The computing domain of arbitrary quadrilateral planform is mapped onto a standard square form by coordinate transformation. The versatile pb-2 Ritz functions defined by the product of a two-dimensional polynomial and a basic function are taken to be the admissible functions. Substituting these displacement functions into the energy functional and minimizing the total energy by differentiation, leads to a typical eigenvalue problem, which is solved by a standard eigenvalue solver. Stiffness and mass matrices are numerically integrated over the plate by using Gaussian quadrature. The accuracy and efficiency of the proposed method are demonstrated through several numerical examples by comparison and convergency studies. A lot of numerical results for reasonable natural frequency parameters of quadrilateral plates with different combinations of elastic boundary conditions and column supports at any locations are presented, which can be used as a benchmark for future studies in this area.

• Journal of applied mathematics & informatics
• /
• 제19권1_2호
• /
• pp.215-229
• /
• 2005

In this paper, constraint aggregation is combined with the adjoint and multiple shooting strategies for optimal control of differential algebraic equations (DAE) systems. The approach retains the inherent parallelism of the conventional multiple shooting method, while also being much more efficient for large scale problems. Constraint aggregation is employed to reduce the number of nonlinear continuity constraints in each multiple shooting interval, and its derivatives are computed by the adjoint DAE solver DASPKADJOINT together with ADIFOR and TAMC, the automatic differentiation software for forward and reverse mode, respectively. Numerical experiments demonstrate the effectiveness of the approach.

• The Journal of the Acoustical Society of Korea
• /
• 제16권1E호
• /
• pp.29-34
• /
• 1997

Numerical prediction of aerodynamic noise radiated by subsonic rotors are carried out. A computer program has been developed which incorporates both the discrete frequency noise as well as the broadband noise arising from the ingestion of turbulence. Acoustic analogy is used in conjunction with Homicz's formulation of turbulence ingestion noise. Formulation 1A of Farassat is used to enhance the numerical analysis performance of Ffowcs-Williams Hawkings equation by eliminating the numericla time differentiation. Homicz's trubulence ingestion noise prediction technique is used to understand the characteristics of broadband noise radiated by isotropic trubulence in gestion. Numerical predictions are carried out for a number of rotor configurations and compared with experimental data. Monopole consideration of transonic rotor agrees well with both the experimental data and the linear theory. Noise radiation characteristics of rotor at lifting hover are investigated utilizing simple blade loading obtained by thin wing section theory. By incorporating discrete noise prediction of steady loading with broadband spectrum, much better agreement with experimental data is obtained in the low frequency region. The contributions from different noise mechanisms can also be analyzed through this method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제25권2호
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Proceedings of the KWS Conference
• /
• 대한용접접합학회 2000년도 특별강연 및 추계학술발표대회 개요집
• /
• pp.260-263
• /
• 2000

A methodology is developed and many used to evaluate the response sensitivity of the thermal systems to variations in their design parameters. Technique for computing the sensitivity of temperature distributions to changes in processing parameters needed for deciding the more effective laser input parameters for laser surface hardening treatment are considered. In this study, a state equation governing the heat flow in laser surface treatment is analyzed using a three-dimensional finite element method and sensitivity data of the processing parameter obtained using a direct differentiation method applied for sensitivity analysis. The interesting processing parameter is taken as the laser scan velocity and characteristic beam radius( $r_{b}$) of the sensitivity of the temperature T versus v and $r_{b}$ is analyzed. And these sensitivity results obtained in another parameters are fixed condition. To verifying the numerical analysis results, hardened layer dimensions (width and depth) of the numerical analysis compared with the results of an experimental data.ata.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 한국전산구조공학회 2000년도 가을 학술발표회논문집
• /
• pp.238-245
• /
• 2000

A general formulation for shape design sensitivity analysis over a plane arch structure is developed based on a variational formulation of curved beam in linear elasticity. Sensitivity formula is derived using the material derivative concept and adjoint variable method for the stress defined at a local segment. Obtained sensitivity expression, which can be computed by simple algebraic manipulation of the solution variables, is well suited for numerical implementation since it does not involve numerical differentiation. Due to the complete description for the shape and its variation of the arch, the formulation can manage more complex design problems with ease and gives better optimum design than before. Several examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. Shape optimization is also conducted with two design problems to illustrate the excellent applicability.