• 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 Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.42-57
• /
• 1995

An efficient and accurate scheme has been constructed by taking advantages of the hi-quadratic spline scheme and the higher-order boundary element method selectively depending on computation domains. Boundary surfaces are represented by 8-node boundary elements to describe curved surfaces of a ship and its neighboring free surface more accurately. The variation of the velocity potential complies with the characteristics of the 8-node element on the body surface. But on the free surface, it is assumed to follow that of the hi-quadratic spline scheme. By which, the free surface solution is free from numerical damping and has better numerical dispersion property. As numerical examples, steady and unsteady Neumann-Kelvin problems are considered. Numerical results for a submerged spheroid, Series 60($C_B=0.6$) and a modified support the proposed method. Finally, a new upstream radiation condition is derived using a wave equation operator in order to deal with problems for subcritical reduced frequency. The relevance of this operator has been confirmed in the case of unsteady Kelvin source potential.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.5
• /
• pp.619-628
• /
• 2010

Shape Memory Alloy (SMA) has recently received attention in developing implantable surgical equipments and it is expected to lead the future medical device market by adequately imitating surgeons' flexible and delicate hand movement. However, SMA actuators have not been used widely because of their nonlinear behavior called hysteresis, which makes their control difficult. Hence, we propose a parameter, $t_1$, which is necessary for temperature control, by analyzing the open-loop step response between current and temperature and by comparing it with the values of linear differential equations. $t_1$ is a pole of the transfer function in the invariant linear model in which the input and output are current and temperature, respectively; hence, $t_1$ is found to be related to the state variable used for temperature control. When considering the parameter under heat treatment conditions, $T_{max}$ was found to assume the lowest value, and $t_1$ was irrelevant to the heat treatment.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.81-94
• /
• 2010

This paper deals with the fracture simulation of UHPFRC girder with the interface type model. Based on the existing numerical simulation of quasi-brittle fracture in normal strength concrete, constitutive modeling for UHPFRC I-girder has been improved by including a tensile hardening at the failure surface. The finite element formulation is based on a triangular unit, constructed from constant strain triangles, with nodes along its sides and neither at the vertex nor the center of the unit. Fracture is simulated through a hardening/softening fracture constitutive law in tension, a softening fracture constitutive law in shear as well as in compression at the boundary nodes, with the material within the triangular unit remaining linear elastic. LCP is used to formulate the path-dependent hardening-softening behavior in non-holonomic rate form and a mathematical programming algorithm is employed to solve the LCP. The piece-wise linear inelastic yielding-failure/failure surface is modeled with two compressive caps, two Mohr-Coulomb failure surfaces, a tensile yielding surface and a tensile failure surface. The comparison between test results and numerical results indicates this method effectively simulates the deformation and failure of specimen.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2003.05a
• /
• pp.181-186
• /
• 2003

The satellite navigation system is widely used for identifying a user's position regardless of weather or geographic conditions and also make effect on new technology of marine LBS(Location Based Service), which has the technology of geographic information such as the ENC. Generally, there are conceivable systems of marine LBS such as ECDIS, or ECS that use the ENC itself with powerful processor in installed type on ships bridge. Since the ENC is relatively heavy structure with dummy format for data transfer between different systems, we should reduce the ENC to small and compact size in order to use it in mobile platform. In this paper, we assumed that the mobile system like PDA, or Webpad can be used for small capability of mobile platform. However, the ENC should be updated periodically by update profile data produced by HO. If we would reduce the ENC without a consideration of update, we could not get newly updated data furthermore. As summary, we studied considerations for ENC reduction with update capability. It will make the ENC be useful in many mobile platforms for various applications.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.3
• /
• pp.263-272
• /
• 2005

For the Communication, Ocean and Meteorological Satellite (COMS) which will be launched in 2008, an algorithm for finding the precise location of the sun-glint point on the ocean surface is studied. The precise locations of the sun-glint are estimated by considering azimuth and elevation angles of Sun-satellite-Earth geometric position and the law of reflection. The obtained nonlinear equations are solved by using the Newton-Raphson method. As a result, when COMS is located at $116.2^{\circ}E$ or $128.2^{\circ}E$ longitude, the sun-glint covers region of ${\pm}10^{\circ}(N-S)$ latitude and $80-150^{\circ}(E-W)$ longitude. The diurnal path of the sun-glint in the southern hemisphere is curved towards the North Pole, and the path in the northern hemisphere is forwards the south pole. The algorithm presented in this paper can be applied to predict the precise location of sun-glint region in any other geostationary satellites.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.4
• /
• pp.225-243
• /
• 2014

As the economic matters are involved, applying the WEC, which is used for controlling waves as well as utilizing the wave energy on existing breakwater, is preferred rather than installing exclusive WEC. This study examines the OWC mounted on a pressurized air chamber floating breakwater regarding the functionality of both breakwater and wave-power generation. In order to verify the performance as a WEC, the velocity of air flow from pressurized air chamber to WEC has to be evaluated properly. Therefore, numerical simulation was implemented based on BEM from linear velocity potential theory as well as Boyle's law with the state equation to analyze pressurized air flow. The validity of the obtained values can be determined by comparing the previous results from numerical analysis and empirically obtained values of different shapes. In the actual numerical analysis, properties of wave deformation around OWC system mounted on fixed type and floating type breakwaters, motions of the structure with air flow velocities are investigated. Since, the wind power generating system can be hybridized on the structure, it is expected to be applied on complex power generation system which generates both wind and wave power energy.

• Journal of Radiation Protection and Research
• /
• v.16 no.2
• /
• pp.41-54
• /
• 1991

This paper aimed to analyse dose sensitivity to the controllable parameters of indoor radon $(^{222}Rn)$ and its decay products (Rn-D) by applying the input~output linear system theory. Physical behaviors of $^{222}Rn\;&\;Rn-D$ were analyzed in terms of $(^{222}Rn)$ gas -generation, -migation and -infiltration to indoor environments, and the performance output-function, i. e. mean dose equivalent to Tracho-Bronchial (TB) lung region, was assessed to the following extented ranges of the controllable paramenters; a) the ventilation rate $constant({\lambda}_v)\;:\;0{\sim}50[h^{-l}].\;b)$ the attachment rate $constant({\lambda}_a)\;:\;0{\sim}500[h^{-l}].\;c)$ the unattached-deposition rate constant (${\lambda}^u_d)\;:\;0-50[h-l]$. A linear input-output model was reconstructed from the original models in literatures, as follows, which was modified into the matrices consisting of 111 nodal equations; a) indoor $^{222}Rn\;&\;Rn-D$ Behaviour; Jacobi-Porstendoerfer-Bruno model.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.2
• /
• pp.103-108
• /
• 2019

Non-equilibrium molecular dynamics simulation on the thermal boundary resistance(TBR) of an aluminum(Al)/silicon(Si) interface was performed in the present study. The constant heat flux across the Si/Al interface was simulated by adding the kinetic energy in hot Si region and removing the same amount of the energy from the cold Al region. The TBR estimated from the sharp temperature drop at the interface was independent of heat flux and equal to $5.13{\pm}0.17K{\cdot}m^2/GW$ at 300K. The simulation result was experimentally confirmed by the time-domain thermoreflectance technique. A 90nm thick Al film was deposited on a Si(100) wafer using an e-beam evaporator and the TBR on the film/substrate interface was measured using the time-domain thermoreflectance technique based on a femtosecond laser system. A numerical solution of the transient heat conduction equation was obtained using the finite difference method to estimate the TBR value. Experimental results were compared to the prediction and discussions on the nanoscale thermal transport phenomena were made.