• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.3 s.234
• /
• pp.369-386
• /
• 2005

Various aspects of structural optimization techniques under transient loads are extensively reviewed. The main themes of the paper are treatment of time dependent constraints, calculation of design sensitivity, and approximation. Each subject is reviewed with the corresponding papers that have been published since 1970s. The treatment of time dependent constraints in both the direct method and the transformation method is discussed. Two ways of calculating design sensitivity of a structure under transient loads are discussed - direct differentiation method and adjoint variable method. The approximation concept mainly focuses on re- sponse surface method in crashworthiness and local approximation with the intermediate variable Especially, as an approximated optimization technique, Equivalent Static Load method which takes advantage of the well-established static response optimization technique is introduced. And as an application area of dynamic response optimization technique, the structural optimization in flexible multibody dynamic systems is re- viewed in the viewpoint of the above three themes

• Computational Structural Engineering
• /
• v.1 no.1
• /
• pp.87-94
• /
• 1988

The linerization method as one of the recursive quadratic programming method is applied for the optimal design of a large-scale structure supported by Pshenichny's proof of global convergence of the algorithm and convergence rate estimates. The linearization method transforms all constants of the design problem into an equivalent linearized constraint and employs the active-set strategy. This results in substantial computational savings by reducing the number of sate and adjoint to be solved at every design iteration. The illustrative example of plates with beams supported by columns is the typical one of a large-scale structure to give successful optimum solutions with satisfactory convergence criteria. Hopefully, the method may be applicable to all classes of optimization problems.

• 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.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.11
• /
• pp.943-950
• /
• 2013

This paper discusses a harmonic response estimation method on the L$\acute{e}$vy plate with two opposite edges simply supported and the other two edges having free boundary conditions. Since the equation of motion of the plate is not self-adjoint, the modes are not orthogonal to each other on the domain. Noting that the L$\acute{e}$vy plate can be expressed using one term sinusoidal function that is orthogonal to other sinusoidal functions, this paper suggested the calculation method that is equivalent to finding a least square error minimization solution of the finite number of algebraic equations. Example problems subjected to a distributed area loading and a distributed line loading are defined and their solutions are provided. The solutions are compared to those of the commercial code, ANSYS. According to the verification results, it is expected that the suggested method will be useful to predict the forced response on the L$\acute{e}$vy plate with the distributed area or line loading conditions.

• Steel and Composite Structures
• /
• v.30 no.1
• /
• pp.57-67
• /
• 2019

This paper deals with the static analysis of axially functionally graded rectangular plates. It is assumed that the flexural rigidity of the plate varies exponentially along one of the plate's in-plane dimensions. Both an analytical approach and a numerical method are utilized to solve the problem. The analytical solution is obtained by using the Green's function method. To employ this approach, the adjoint boundary value problem is established. Then, exact solutions for deflection of the plate for different boundary conditions are found. In another way, a finite element formulation for the problem is developed. In order to demonstrate the validity of the Authors' formulation, the results obtained via both mentioned schemes are compared with each other for functionally graded plates and with results of previously published works for homogeneous plates. The effect of plate parameters on the response of the plate is also investigated. To remind the research background, a brief review on the application of Green's function method in plates' analysis and functionally graded plates is also presented.

• Computational Structural Engineering
• /
• v.7 no.1
• /
• pp.91-98
• /
• 1994

Design sensitivity analysis method for the vibration of vehicle structure is developed using adjoint variable method. A variational approach with complex response method is used to derive sensitivity expression. To evaluate sensitivity, FEM analysis of ship deck and vehicle structure are performed using MSC/NASTRAN installed in the super computer CRAY2S, and sensitivity computation is performed by PC. The accuracy of sensitivity is verified by the results of finite difference method. When compared to structural analysis time on CRAY2S, sensitivity computation is remarkably economical. The sensitivity of vehicle frame can be used to reduce the vibration responses such as displacement and acceleration of vehicle.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1715-1725
• /
• 1998

A two-dimensional transient inverse heat conduction problem involving the estimation of the unknown location, ($X^*$, $Y^*$), and timewise varying unknown strength, $G({\tau})$, of a line heat source embedded inside a rectangular bar with insulated boundaries has been solved simultaneously. The regular conjugate gradient method, RCGM and the modified conjugate gradient method, MCGM with adjoint equation, are used alternately to estimate the unknown strength $G({\tau})$ of the source term, while the parameter estimation approach is used to estimate the unknown location ($X^*$, $Y^*$) of the line heat source. The alternate use of the regular and the modified conjugate gradient methods alleviates the convergence difficulties encountered at the initial and final times (i.e ${\tau}=0$ and ${\tau}={\tau}_f$), hence stabilizes the computation and fastens the convergence of the solution. In order to examine the effectiveness of this approach under severe test conditions, the unknown strength $G({\tau})$ is chosen in the form of rectangular, triangular and sinusoidal functions.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.10
• /
• pp.915-922
• /
• 2010

This paper presents an efficient method for determining the forced response of a spinning flexible disk-spindle system supported by fluid dynamic bearings(FDBs) in a computer hard disk drive(HDD). The spinning flexible disk-spindle system is represented by the asymmetric finite element equations of motion originating from the asymmetric dynamic coefficients of the FDBs and the gyroscopic moment of a spinning disk-spindle system. The proposed method utilizes only the right eigenvectors of the eigenvalue problem to transform the large asymmetric finite element equations of motion into a small number of coupled equations, guaranteeing the accuracy of their numerical integration. The results are then back-substituted into the equations of motion to determine the forced response. The effectiveness of the proposed method was verified by comparing it with the responses from the classical methods of mode superposition with the general eigenvalue problems, and mode superposition with modal approximation. The proposed method was shown to be effective in determining the forced response represented by the asymmetric finite element equations of motion of a spinning flexible disk-spindle system supported by FDBs.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2625-2634
• /
• 2022

Discretization errors are extremely challenging conundrums of discrete ordinates calculations for radiation transport problems with void regions. In previous work, we have presented a multi-collision source method (MCS) to overcome discretization errors, but the efficiency needs to be improved. This paper proposes a goal-oriented algorithm for the MCS method to adaptively determine the partitioning of the geometry and dynamically change the angular quadrature in remaining iterations. The importance factor based on the adjoint transport calculation obtains the response function to get a problem-dependent, goal-oriented spatial decomposition. The difference in the scalar fluxes from one high-order quadrature set to a lower one provides the error estimation as a driving force behind the dynamic quadrature. The goal-oriented algorithm allows optimizing by using ray-tracing technology or high-order quadrature sets in the first few iterations and arranging the integration order of the remaining iterations from high to low. The algorithm has been implemented in the 3D transport code ARES and was tested on the Kobayashi benchmarks. The numerical results show a reduction in computation time on these problems for the same desired level of accuracy as compared to the standard ARES code, and it has clear advantages over the traditional MCS method in solving radiation transport problems with reflective boundary conditions.