• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.217-223
• /
• 1997

This paper deals with the basic concepts of sensitivity analysis in a time domain for the transient response of a lumped mass system. Sensitivity analysis methods in thme domain for determining the effects of parameter changes on the response of a dynamic system by external excitation are presented. The parametric sensitivity of a lumped mass system in time domain can be investigated using different types of sensitivity functions, including first order standard and percentage sensitivity functions. These sensitivity functions are determined as a function of partial derivatives of system variables taken with respect to system parameters. In addition, we compared the results of the analytical method by direct method and those of numerical methods.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.12 s.117
• /
• pp.1192-1200
• /
• 2006

Viscoelastic damping materials are widely used to reduce noise and vibration because of its low cost and easy implementation, for examples, on the body structure of passenger cars, air planes, electric appliances and ships. To design the damped structures, the material property such as elastic modulus and loss factor is essential information. The four-parameter fractional derivative model well describes the dynamic characteristics of the viscoelastic damping materials with respect to both frequency and temperature. However, the identification procedure of the four-parameter is very time-consuming one. In this study a new identification procedure of the four-parameters is proposed by using an FE model and a gradient-based numerical search algorithm. The identification procedure goes two sequential steps to make measured frequency response functions(FRF) coincident with simulated FRFs: the first one is a peak alignment step and the second one is an amplitude adjustment step. A numerical example shows that the proposed method is useful in identifying the viscoelastic material parameters of fractional derivative model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.299-306
• /
• 2006

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.7
• /
• pp.634-641
• /
• 2008

Automatic Differentiation(AD) is a tool for generating sensitivities, such as gradient or Jacobian, automatically. AD tools provide mathematically exact sensitivities for the given source code. In this paper applications of automatic differentiation are studied. Derivative codes are generated with AD tools for structural analysis code and flow analysis code. How to apply AD tools is explained and the accuracy of sensitivities is compared with the finite difference. Sensitivities of generated derivative code accord well with finite difference, but the calculation time of derivative code increases. It was found that the calculation time can be decreased by additional modification of derivative code.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.176-182
• /
• 2000

This paper presents methods for dynamic stall control utilizing an optimization approach. Unsteady aerodynamic sensitivity code is developed using a direct diffentiation method from a discrete two-dimensional unsteady compressible Navier-Stokes code including a two-equation turbulence model. Dynamic stall control is conducted by minimizing an objective function defined at an instant instead of integrating a response for a period of time. Unsteady sensitivity derivative of the objective function is calculated by the sensitivity code, and optimization is conducted using a linear line search method at every physical time level. Examples of dynamic stall control utilizing airfoils nose radius or maximum thickness variation show very satisfactory results.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.4
• /
• pp.435-447
• /
• 1999

It is main objective of this approach to present a method to analyse stochastic design sensitivity for problems of structural dynamics with randomness in design parameters. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element approach. An alternative form of the constant functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The terminal problem of the adjoint system is solved using equivalent homogeneous equations excited by initial velocities. The numerical procedures are shown to be much more efficient when based on the fold superposition method: the generalized co-ordinates are normalized and the correlated random variables are transformed to uncorrelated variables, whereas the secularities are eliminated by the fast Fourier transform of complex valued sequences. Numerical algorithms have been worked out and proved to be accurate and efficient : they can be readily adapted to fit into the existing finite element codes whose element derivative matrices can be explicitly generated. The numerical results of two cases -2 dimensional portal frame for the comparison with reference and 3-dimensional frame structure - for the deterministic sensitivity analysis are presented.

• Proceedings of the KIEE Conference
• /
• 1992.07b
• /
• pp.567-569
• /
• 1992

Design sensitivity analysis is proposed for the optimal shape design of three dimensional magnetostatic problems. The direct differentiation method is introduced for design sensitivity analysis and the boundary element method with reduced magnetic scalar potential as the state variable is used to analyze the magnetic characteristics. In the direct differentiation method, the design sensitivity, defined as the total derivative of the objective function with respect to the design variables, is calculated based on the variation of the state variable with respect to the design variable. And the variation of He state variable is calculated by differentiating the both sides of the system matrix equation obtained by applying boundary element method. Through the numerical example with simple electromagnet, the usefullness is proved.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.41 no.8
• /
• pp.850-857
• /
• 1992

Design sensitivity analysis is proposed for the optimal shape design of three-dimensional magnetostatic problems. The direct differentiation method is introduced for design sensitivity analysis and the boundary element method with reduced magnetic scalar potential as the state variable is used to analyze the magnetic characteristics. In the direct differentiation method, the design sensitivity, defined as the total derivative of the objective function with respect to the design variables, is calculated based on the variation of the state variable with respect to the design variable. And the variation of the state variable is calculated by differetiating the both sides of the system matrix equation obtained by applying boundary element method. Through the numerical example with simple electromagnet, the usefulness is proved.

• Korean Journal of Optics and Photonics
• /
• v.25 no.5
• /
• pp.273-278
• /
• 2014

We present an algorithm for surface reconstruction from the second-derivative data for free-form aspherics, which uses a subaperture scanning system that measures the local surface profile and determines the three second-derivative values at those local sampling points across the free-form surface. The three second-derivative data were integrated to get a map of x- and y-slopes, which went through a second Southwell integration step to reconstruct the surface profile. A synthetic free-form surface 200 mm in diameter was simulated. The simulation results show that the reconstruction error is 19 nm RMS residual difference. Finally, the sensitivity to noise is diagnosed for second-derivative Gaussian random noise with a signal to noise ratio (SNR) of 16, the simulation results proving that the suggested method is robust to noise.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.