• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.6
• /
• pp.367-374
• /
• 2022

This papter presents the use of the automatic differential method based on the backpropagation method to obtain the design sensitivity and its application to topology optimization considering the stress constraints. Solving topology optimization problems with stress constraints is difficult owing to singularities, the local nature of stress constraints, and nonlinearity with respect to design variables. To solve the singularity problem, the stress relaxation technique is used, and p-norm for stress constraints is applied instead of local stresses for global stress measures. To overcome the nonlinearity of the design variables in stress constraint problems, it is important to analytically obtain the exact design sensitivity. In conventional topology optimization, design sensitivity is obtained efficiently and accurately using the adjoint variable method; however, obtaining the design sensitivity analytically and additionally solving the adjoint equation is difficult. To address this problem, the design sensitivity is obtained using a backpropagation technique that is used to determine optimal weights and biases in the artificial neural network, and it is applied to the topology optimization with the stress constraints. The backpropagation technique is used in automatic differentiation and can simplify the calculation of the design sensitivity for the objectives or constraint functions without complicated analytical derivations. In addition, the backpropagation process is more computationally efficient than solving adjoint equations in sensitivity calculations.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.7
• /
• pp.2858-2864
• /
• 2012

The differential quadrature method (DQM) is applied to computation of the eigenvalues of in-plane buckling of the curved beams. Critical moments and loads are calculated for the beam subjected to equal and opposite bending moments and uniformly distributed radial loads with various end conditions and opening angles. Results are compared with existing exact solutions where available. The DQM gives good accuracy even when only a limited number of grid points is used. More results are given for two sets of boundary conditions not considered by previous investigators for in-plane buckling: clamped-clamped and simply supported-clamped ends.

• Journal of the Korean Society of Safety
• /
• v.17 no.4
• /
• pp.189-195
• /
• 2002

The differential quadrature method (DQM) is applied to computation of the eigenvalues of out-of-plane bucking of curved beams. Critical moments including the effect of radial stresses are calculated for a single-span wide-flange beam subjected to equal and opposite in-plane bending moments with various end conditions, and opening angles. Results are compared with existing exact solutions where available. The differential quadrature method gives good accuracy even when only a limited number of grid points is used. New results are given for two sets of boundary conditions not previously considered for this problem: clamped-clamped and clamped-simply supported ends.

• Computational Structural Engineering
• /
• v.3 no.1
• /
• pp.83-95
• /
• 1990

Many structural optimization problems are solved by numerical algorithms since these are complicated and nonlinear. To provide a wider base and popular it to structual design optimization, reliable, accurate and superlinearly convergent nonlinear programming algorithm with active-set strategy have been developed. One of these is RQP(recursive quadratic programming method). This algorithm has several parameters and its performance is influenced by variations of these key parameters. Therefore, an RQP algorithm is selected to enhance its numerical performances by choosing proper parameters. The paper persents these influences on its numerical performance. For comparison of performances, a structural design software for minimum weight of truss subjected to displacement, stress, and lower and upper bounds on design variables is also implemented.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.3
• /
• pp.19-28
• /
• 1990

The main purpose of this paper is to present both fundamental and some higher natural frequencies of catenary arcs. The differential equations governing planar free vibrations for these arcs are derived, in which the rotatory inertia is included, as non-dimensional forms and solved numerically to obtain frequencies and mode shapes. The hinged-hinged and clamped-clamped end constraints are applied in numerical examples. The lowest four natural frequencies are reported as the functions of non -dimensional system parameters; the slenderness ratio and the rise to span length ratio. The effects of rotatory inertia on natural frequencies are reported and some typical mode shapes are also presented.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.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 Korean Society of Computer Information Conference
• /
• 2013.01a
• /
• pp.3-6
• /
• 2013

반도체 공정에서 소자의 제조 비용 감소를 위해 제조 공정 검증을 위한 시뮬레이션을 수행하게 된다. 이 시뮬레이션은 반도체 소자 내부의 물리량 계산을 통해 반도체 소자 내부의 불순물의 거동을 해석하게 된다. 이를 위해 사용되는 알고리즘으로 3차원적 형상을 표현하는 물리적 미분 미분방정식을 계산하게 되는데, 정확한 계산을 위해 유한 차분 시간 영역법(이하 FDTD)과 같은 수치해석 기법을 이용한다. 실제적으로 반도체 공정의 시뮬레이션에서 FDTD연산의 실행 시간은 90% 이상을 소요하게 된다. 이러한 연산에서 더욱 빠른 성능을 확보하기 위해 본 논문에서는 기존의 CUDA(Compute Unified Device Architecture)로 구현된 FDTD알고리즘을 OpenMP를 통한 다중 GPU제어를 이용하여 연산 수행시간을 감소하고, 그 결과물을 통하여 성능 향상도를 측정한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.391-399
• /
• 2008

In this paper, we further generalize the work of Lin and Abel to the case of the first and the second order derivatives of energy release rates for two-dimensional, multiply cracked systems. The direct integral expressions are presented for the energy release rates and their first and second order derivatives. The salient feature of this numerical method is that the energy release rates and their first and second order derivatives can be computed in a single analysis. It is demonstrated through a set of examples that the proposed method gives expectedly decreasing, but acceptably accurate results for the energy release rates and their first and second order derivatives. The computed errors were approximately 0.5% for the energy release rates, $3\sim5%$ for their first order derivatives and $10\sim20%$ for their second order derivatives for the mesh densities used in the examples. Potential applications of the present method include a universal size effect model and a probabilistic fracture analysis of cracked structures.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.27 no.6
• /
• pp.782-786
• /
• 1994

Integrated finite difference method (IFDM) is used to solve one dimensional free oscillation problem in the coastal area. To evaluate the solution accuracy of IFDM in free oscillation analysis, two finite difference equations based on area discretization method and point discretization method are derived from the governing equations of free oscillation, respectively. The difference equations are transformed into a generalized eigenvalue problem, respectively. A numerical example is presented, for which the analytical solution is available, for comparing IFDM to conventional finite difference equation (CFDM), qualitatively. The eigenvalue matrices are solved by sub-space iteration method. The numerical results of the two methods are in good agreement with analytical ones, however, IFDM yields better solution than CFDM in lower modes because IFDM only includes first order differential operator in finite difference equation by Green's theorem. From these results, it is concluded that IFDM is useful for the free oscillation analysis in the coastal area.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.1
• /
• pp.75-81
• /
• 2002

In the re-entry comtrol system, errors apt to induce because the time deriviative of drag acceleration is analytically estimated. Still more, the difficulty of estimation of th exact drag coefficient in hypersonic velocity and the non-reality of the scale height cause a steady-state drag errer. In the Space-Shuttle, a steady-state drag error is reduced by the addition of the integral term of drag acceleation error into the control system. This method, however, induces a difficulties in respect to the modern controller composition due to the multi-poles in a closed-loop system. Thus, this paper proposes the additional method of the disturbance observer. This reduces the steady-state drag error according to the following by the analytic calculation, and then creates the new drag acceleration time derivative using the estimated error. The performance of the re-entry control system is verified about 32 refernce trajectories.