• Structural Engineering and Mechanics
• /
• 제56권5호
• /
• pp.871-897
• /
• 2015

This paper presents a novel method based on sensitivity of structural response for identifying both the system parameters and input excitation force of a bridge. This method, referred to as "Adjoint Variable Method", is a sensitivity-based finite element model updating method. The computational cost of sensitivity analyses is the main concern associated with damage detection by these methods. The main advantage of proposed method is inclusion of an analytical method to augment the accuracy and speed of the solution. The reliable performance of the method to precisely indentify the location and intensity of all types of predetermined single, multiple and random damages over the whole domain of moving vehicle speed is shown. A comparison study is also carried out to demonstrate the relative effectiveness and upgraded performance of the proposed method in comparison to the similar ordinary sensitivity analysis methods. Moreover, various sources of error including the effects of noise and primary errors on the numerical stability of the proposed method are discussed.

• 한국전산구조공학회논문집
• /
• 제27권4호
• /
• pp.297-303
• /
• 2014

페리다이나믹스 이론과 이진분해 기법의 병렬연산을 이용하여 동적 균열진전 문제에 대한 애조인 형상 설계민감도 해석법을 개발하였다. 페리다이나믹스에서는 균열의 연속적인 분기를 다룰 수 있으며, Explicit 시간적분법을 채택한다. 설계민감도 해석은 애조인 변수법은 경로의존성 문제에는 적합하지 않으나 여기서는 응답해석의 경로를 이미 알고 있으므로 채택하여 사용할 수 있었다. 얻어진 해석적 설계민감도는 유한차분과 비교하여 그 정확성을 검증하였다. 유한차분법은 설계섭동량에 민감하여 비선형성이 강한 페리다이나믹스 문제에서 부정확한 설계민감도를 제시할 수 있다. 정확한 설계민감도 해석을 위해서는 이산화과정에서 $C^1$ 연속성을 가지는 체적율이 필요함을 알 수 있었다.

• 한국전산구조공학회논문집
• /
• 제35권6호
• /
• pp.367-374
• /
• 2022

본 논문에서는 역전파 방법 기반 자동미분법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용하였다. 응력제한조건이 있는 위상최적화문제는 특이점(singularity)과 응력의 국부성(local nature of stress constraint)문제, 그리고 설계 변수에 대한 비선형성의 문제를 포함하고 최적해를 얻기가 매우 힘들다. 특이점 문제를 해결하기 위해서 응력 완화(stress relaxation) 기법을 사용하였고, 응력의 국부성을 해결하기 위해 p-norm을 이용한 전역 응력치를 제한조건에 사용하였다. 설계 변수에 대한 비선 형성을 극복하기 위해 해석적인 방법으로 정확한 설계민감도를 구하는 것이 중요하다. 위상최적설계에서 기존에는 보조변수방법 (adjoint variable method)을 사용하여 빠르고 정확한 설계민감도를 구했지만, 설계민감도를 해석적으로 구해야 하고, 보조평형방정식을 추가로 풀어야 하는 어려움이 있다. 이를 해결하기 위해서 인공신경망에서 최적 가중치(weights)와 편차(biases)를 구할 때 쓰이는 역전파 기법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용하였다. 역전파 기법은 자동미분에 쓰이는 기법으로 목적함수나 제한조건에 대한 설계민감도를 별도의 수식유도 없이 간단하게 구할 수 있는 장점이 있다. 또한, 미분값을 구하는 역전파의 과정이 보조평형방정식을 푸는 것보다 계산시간이 빠르고 해석적 방법으로 구한 설계민감도와 같은 정확도를 보여준다.

• 대한원격탐사학회:학술대회논문집
• /
• 대한원격탐사학회 1999년도 Proceedings of International Symposium on Remote Sensing
• /
• pp.221-226
• /
• 1999

A data assimilation system for a 1-dimensional mixed layer model has been constructed using the adjoint method. The classical adjoint method does not work well for the mixed layer variabilities due to the occurrence of spikes in the gradient of the cost function. To solve this problem, the two techniques of scaling the cost function and optimization in the frequency space are used. As a result, the heat flux can be reliably estimated with an accuracy of 8Wm$^{-2}$ rms error in the identical twin experiments. We then applied this system to the tropical Pacific TOGA-TAO buoy data. The air-sea heat flux as well as the mixed layer variability were estimated in close approximation to the buoy data, particularly on time scales longer than the seasonal one.

• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1172-1180
• /
• 2017

The application of Monte Carlo (MC) to large-scale fixed-source problems has recently become possible with new hybrid methods that automate generation of parameters for variance reduction techniques. Two common variance reduction techniques, weight windows and source biasing, have been automated and popularized by the consistent adjoint-driven importance sampling (CADIS) method. This method uses the adjoint solution from an inexpensive deterministic calculation to define a consistent set of weight windows and source particles for a subsequent MC calculation. One of the motivations for source consistency is to avoid the splitting or rouletting of particles at birth, which requires computational resources. However, it is not always possible or desirable to implement such consistency, which results in inconsistent source biasing. This paper develops an original framework that mathematically expresses the coupling of the weight window and source biasing techniques, allowing the authors to explore the impact of inconsistent source sampling on the variance of MC results. A numerical experiment supports this new framework and suggests that certain classes of problems may be relatively insensitive to inconsistent source sampling schemes with moderate levels of splitting and rouletting.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 추계학술대회논문집
• /
• pp.1022-1027
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered metal bearing and flexible structures by using the finite element method and the mode superposition method considering the asymmetry of the gyroscopic effect and sintered metal bearing. The eigenvalues and eigenvectors are calculated by solving the eigenvalue problem and the adjoint eigenvalue problem by using the restarted Arnoldi iteration method. The decoupled equations of motion can be obtained from global finite element motion equations by using the orthogonal relation between the right eigenvectors and left eigenvectors. The decoupled equations of motion are used to analyze the unbalance response of a high speed polygon mirror scanner motor. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results.

• 대한조선학회논문집
• /
• 제47권1호
• /
• pp.58-66
• /
• 2010

Topology design optimization is a method to determine the optimal distribution of material that yields the minimal compliance of structures, satisfying the constraint of allowable material volume. The method is easy to implement and widely used so that it becomes a powerful design tool in various disciplines. In this paper, a large-scale topology design optimization method is developed using the efficient adjoint sensitivity and optimality criteria methods. Parallel computing technique is required for the efficient topology optimization as well as the precise analysis of large-scale problems. Parallelized finite element analysis consists of the domain decomposition and the boundary communication. The preconditioned conjugate gradient method is employed for the analysis of decomposed sub-domains. The developed parallel computing method in topology optimization is utilized to determine the optimal structural layout of human powered vessel.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1999년도 춘계학술발표회요약집
• /
• pp.15-15
• /
• 1999

• International Journal of Aeronautical and Space Sciences
• /
• 제8권1호
• /
• pp.95-104
• /
• 2007

A new design approach of complex geometries such as wing/body configuration is arranged by using overset mesh techniques under large scale computing environment. For an in-depth study of the flow physics and highly accurate design, several special overlapped structured blocks such as collar grid, tip-cap grid, and etc. which are commonly used in refined drag prediction are adopted to consider the applicability of the present design tools to practical problems. Various pre- and post-processing techniques for overset flow analysis and sensitivity analysis are devised or implemented to resolve overset mesh techniques into the design optimization problem based on Gradient Based Optimization Method (GBOM). In the pre-processing, the convergence characteristics of the flow solver and sensitivity analysis are improved by overlap optimization method. Moreover, a new post-processing method, Spline-Boundary Intersecting Grid (S-BIG) scheme, is proposed by considering the ratio of cell area for more refined prediction of aerodynamic coefficients and efficient evaluation of their sensitivities under parallel computing environment. With respect to the sensitivity analysis, discrete adjoint formulations for overset boundary conditions are derived by a full hand-differentiation. A smooth geometric modification on the overlapped surface boundaries and evaluation of grid sensitivities can be performed by mapping from planform coordinate to the surface meshes with Hicks-Henne function. Careful design works for the drag minimization problems of a transonic wing and a wing/body configuration are performed by using the newly-developed and -applied overset mesh techniques. The results from design applications demonstrate the capability of the present design approach successfully.

• Journal of Mechanical Science and Technology
• /
• 제14권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.