• Journal of Welding and Joining
• /
• v.16 no.3
• /
• pp.141-147
• /
• 1998

Design sensitivity analysis of a vehicle system is an essential tool for design optimization and trade-off studies. Most optimization algorithms require the derivatives of cost and constraint function with respect to design in order to calculate the next improved design. This paper presents an efficient algorithm application for the design sensitivity analysis, using the direct differentiation method. A mounting area of suspension that welded on chassis frame is analyzed to show the validity and the efficiency of the proposed method. A mounting area of suspension that welded on chassis frame is analyzed to show the validity and the efficiency of the proposed method.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.4
• /
• pp.13-23
• /
• 1984

The purpose of this paper is to develop a numerical model which can be used to compute the suspended load concentration of which the two-dimensional unsteady diffusion equation is able to be solved by the finite difference method using the implicit scheme. The pick-up rate formula from the bottom used in the open channel as a sink source term and the Coleman's empirical formula for the diffusion coefficient were taken, and especially the hindered settling velocity and the vertical velocity of flow due to the periodical tidal motion were taken into account, while the effects of the variables, such as the horizontal and vertical velocities of flow, tidal range, the settling velocity and hindered settling, on the suspened load concentrations have been discussed, comparing the results obtained from the different cases in the simulation conditions.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.10 no.3
• /
• pp.237-244
• /
• 2002

Analyzing acoustic-structural systems such as automobiles and aircraft, the FRF-based substructuring (FBS) method is one of the most powerful tools. In this paper, a general procedure for the parametric sensitivity analysis of vibro-acoustic problems has been presented using the multi-domain FRF-based substructuring formulation. For an acoustic-structural system sub-structured by multiple domains, the substructuring formulation gives the reaction farces on the interface boundaries. The design sensitivity formula is obtained from the direct differentiation of the reaction force expression with respect to the design vector. As a practical application, the proposed design sensitivity formula is applied to an engine mount system of passenger car. An objective of the problem is to identify the most effective engine mounts and bushes in minimizing the interior noise over the concerned rpm range. The comparison of the sensitivity results with those of the finite difference method shows excellent agreement. In addition, stiffness modifications of the mounts and bushes identified through the design sensitivity analysis lead to a successful decrease of the interior noise. This results show usefulness of the present method very well.

• Nuclear Engineering and Technology
• /
• v.43 no.6
• /
• pp.531-546
• /
• 2011

The backward differentiation formula (BDF) method is applied to a three-dimensional reactor kinetics calculation for efficient yet accurate transient analysis with adaptive time step control. The coarse mesh finite difference (CMFD) formulation is used for an efficient implementation of the BDF method that does not require excessive memory to store old information from previous time steps. An iterative scheme to update the nodal coupling coefficients through higher order local nodal solutions is established in order to make it possible to store only node average fluxes of the previous five time points. An adaptive time step control method is derived using two order solutions, the fifth and the fourth order BDF solutions, which provide an estimate of the solution error at the current time point. The performance of the BDF- and CMFD-based spatial kinetics calculation and the adaptive time step control scheme is examined with the NEACRP control rod ejection and rod withdrawal benchmark problems. The accuracy is first assessed by comparing the BDF-based results with those of the Crank-Nicholson method with an exponential transform. The effectiveness of the adaptive time step control is then assessed in terms of the possible computing time reduction in producing sufficiently accurate solutions that meet the desired solution fidelity.

• 한국해양학회지
• /
• v.30 no.1
• /
• pp.13-26
• /
• 1995

In this paper three turbulence models including two-equation model by Blumberg and Mellor (1987), one-equation model with mixing length formula of Blackadar's (1962), and zero-equation model of Prandtl's (1925) were compared in homogeneous, unstratified channel flows. Steady flows which a steep-sided trapezoidal trench with uniform discharge, tidal flow and steady wind-driven flow in finite channels are considered in detail. Steady flows in a trench and tidal flows in a finite channel were reproduced fairly accurately and there was virtually no difference among results of three turbulence models. However, In case of steady wind-driven flow only two-equation model reproduced the important features of experimental data. the other two models underestimated the surface velocity. In tidal and wind-driven flows with negligibly small adjective and diffusive effects, the two-equation model gives rise to parabolic profile of eddy viscosity with maximum at the mid0depth, and the one and zero equation model based on Blackadar formula linear profile with maximum at the surface.

• Structural Engineering and Mechanics
• /
• v.15 no.2
• /
• pp.159-180
• /
• 2003

In this paper spacing and width of flexural cracks in reinforced concrete beams are determined using two-dimensional finite element analysis. At early loading stages on the beam the primary crack spacing is based on the slip length, which is the development length required to resist the steel stress increment that occurs at a cracked section on the formation of the first flexural crack. A semi-empirical formula is presented in this paper for the determination of the slip length for a given beam. At higher load levels, the crack spacing is based on critical crack spacing, which is defined as the particular crack spacing that would produce a concrete tensile stress equal to the flexural strength of concrete. The resulting crack width is calculated as the relative difference in extensions of steel reinforcement and adjacent concrete evaluated at the cracked section. Finally a comparative study is undertaken, which indicates that the spacing and width of cracks calculated by this method agree well with values measured by other investigators.

• Proceedings of the Korea Institute of Applied Superconductivity and Cryogenics Conference
• /
• 2002.02a
• /
• pp.274-276
• /
• 2002

To enhance the critical current of superconducting coil, the magnetic field experienced by superconductors in a coil should be minimized. This is true for both low $T_{c}$ and high $T_{c}$ superconductors, and the difference between the two lays in their isotropic/an-isotropic characteristics. In this paper, we propose a shape optimization algorithm to reduce radial magnetic field components in HTS solenoid to enhance the critical current of a solenoid. In the algorithm, finite element method and continuum shape design sensitivity formula were employed. The objective function is to minimize the maximum radial magnetic fields in a solenoid with a constraint of constant solenoid volume condition. In this paper, the details on algorithm are introduced and the calculated optimized shapes are presented.

• Proceedings of the KIEE Conference
• /
• 2008.07a
• /
• pp.651-652
• /
• 2008

This research presents a multi-objective optimal design employing topological approach to maximize magnetic energy while minimizing structural deformation which is caused by magnetic reluctance force and magnetostriction. A design sensitivity formula is derived by employing the adjoint variable method (AVM) to avoid numerous sensitivity evaluations for a coupled magneto-mechanical analysis. The sensitivity analysis is verified using the finite difference method (FDM) in a C-shape actuator. A linear actuator used in a home appliance is examined for optimal design and demonstrates the strength of the proposed topology optimization approach.