• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• pp.541-547
• /
• 2000

This paper presents the closed-form forward kinematics of the 6-6 Stewart platform of planar base and moving platform. Based on algebraic elimination method and with one extra linear sensor, it first derives an 8th-degree univariate equation and then finds tentative solution sets out of which the actual solution is to be selected. In order to provide more exact solution despite the error between measured sensor value and the theoretical one, a correction method is also used. The overall procedure requires so little computation time that it can be efficiently used for realtime applications. In addition, unlike the iterative schemes e.g. Newton-Raphson, the algorithm does not require initial estimates of solution and is free of the problems that it does not converge to actual solution within limited time. The presented method has been implemented in C language and a numerical example is given to confirm the effectiveness and accuracy of the developed algorithm.

• Journal of Mechanical Science and Technology
• /
• 제19권5호
• /
• pp.1065-1071
• /
• 2005

In this study, residual stress distribution in multi-stacked film by MEMS (Micro-Electro Mechanical System) process is predicted using Finite Element method (FEM). We evelop a finite element program for residual stress analysis (RESA) in multi-stacked film. The RESA predicts the distribution of residual stress field in multi-stacked film. Curvatures of multi­stacked film and single layers which consist of the multi-stacked film are used as the input to the RESA. To measure those curvatures is easier than to measure a distribution of residual stress. To verify the RESA, mean stresses and stress gradients of single and multi layers are measured. The mean stresses are calculated from curvatures of deposited wafer by using Stoney's equation. The stress gradients are calculated from the vertical deflection at the end of cantilever beam. To measure the mean stress of each layer in multi-stacked film, we measure the curvature of wafer with the left film after etching layer by layer in multi-stacked film.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2005년도 추계학술대회 논문집
• /
• pp.269-272
• /
• 2005

For the three decades, the mold-filling of injection molding process was modeled as Hele-Shaw model. However, this model can not consider the 3D effect. In this paper, numerical simulations of three dimensional mold-filling during the filling phase were performed. The governing equations were discretized by segregated finite element method, which used equal order interpolation for pressure and velocity fields. The iterative linear equation solver (JCG, SOR) was employed for the solution of the momentum and pressure equations. Volume of Fluid (VOF) was employed for the melt front advancement. To check the validity of the numerical results, the results were compared with the experimental ones. The agreements between the experiment and the numerical results were found to be satisfactory.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제2권1호
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• Structural Engineering and Mechanics
• /
• 제43권5호
• /
• pp.583-605
• /
• 2012

This paper presents a spatial catenary cable element for the nonlinear analysis of cable-supported structures. An incremental-iterative solution based on the Newton-Raphson method is adopted for solving the equilibrium equation. As a result, the element stiffness matrix and nodal forces are determined, wherein the effect of self-weight and pretension are taken into account. In the case of the initial cable tension is given, an algorithm for form-finding of cable-supported structures is proposed to determine precisely the unstressed length of the cables. Several classical numerical examples are solved and compared with the other available numerical methods or experiment tests showing the accuracy and efficiency of the present elements.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1999년도 가을 학술발표회 논문집
• /
• pp.176-183
• /
• 1999

This paper describes the application of discretized continuum-type optimality criteria (DCOC) for the reinforced concrete continuous beams. The cost of construction as objective function which includes the costs of concrete, reinforced steel, formwork is minimized. The design constraints include limits on the maximum deflection in a given span, on bending and shear strengths, optimality criteria is given based on the well known Kuhn-Tucker necessary conditions, followed by an iterative procedure for designs when the design variables are the depth and the steel ratio. The self-weight of the beam is included in the equilibrium equation of the real system. Two numerical examples of reinforced concrete continuous beams with rectangular cross-section are solved to show the applicability and efficiency for the DCOC-based technique

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1999년도 가을 학술발표회 논문집
• /
• pp.359-369
• /
• 1999

The paper describes the study of concurrent subspace optimization(CSSO) for coupled multidisciplinary design optimization (MDO) techniques in mechanical systems. This method is a solution to large scale coupled multidisciplinary system, wherein the original problem is decomposed into a set of smaller, more tractable subproblems. Key elements in CSSO are consisted of global sensitivity equation(GSE), subspace optimization (SSO), optimum sensitivity analysis(OSA), and coordination optimization problem(COP) so as to inquiry valanced design solutions finally, Automatic differentiation has an ability to provide a robust sensitivity solution, and have shown the numerical numerical effectiveness over finite difference schemes wherein the perturbed step size in design variable is required. The present paper will develop the automatic differentiation based concurrent subspace optimization(AD-CSSO) in MDO. An automatic differentiation tool in FORTRAN(ADIFOR) will be employed to evaluate sensitivities. The use of exact function derivatives in GSE, OSA and COP makes Possible to enhance the numerical accuracy during the iterative design process. The paper discusses how much influence on final optimal design compared with traditional all-in-one approach, finite difference based CSSO and AD-CSSO applying coupled design variables.

• 대한기계학회논문집
• /
• 제17권10호
• /
• pp.2457-2466
• /
• 1993

A simple method is presented for determining transient responses of locally nonlinear structures using substructure eigenproperties and Lagrange multiplier technique. Although the method is based upon the mode synthesis formulation procedure, the equations of the combined whole structure are not constructed compared with the conventional methods. Lagrange multi-pliers are used to enforce the conditions of geometric compatibility between the substructure interfaces and they are treated as external forces on each substructure itself. Substructure eigenvalue problem is defined with the substructure interface free of fixed. The transient analysis is based upon the recurrence discrete-time state equations and offers the simplicity of the Euler integration method without requiring small time increment and iterative solution procedure. Numerical examples reveal that the method is very accurated and efficient in calculating transient responses compared with the direct numerical integration method.

• 전자공학회논문지A
• /
• 제31A권6호
• /
• pp.34-44
• /
• 1994

In this paper, an iterative inversion method in angular spectral domain is presented for microwave imaging of a perfectly conducting cylinder. Angular spectra are calculated from measured far-field scattered fields. And then both the propagating modes and the evanescent modes are defined. The center and initial shape of an unknown conductor may be obtained by the characteristics of angular spectra and the total scattering cross section (TSCS). Finally, the orignal shape is reconstructed by the modified Newton algorithm. By using well estimated initial shape the local minima can be avoided, which might appear when the nonlinear equation is solved with Newton algorithm. It is shown to be robust to noise in scattered fields via numerical examples by keeping only the propagating modes and filtering out the evanescent modes.

• 대한조선학회지
• /
• 제17권2호
• /
• pp.1-10
• /
• 1980

The optimal plastic design of framed structures has been treated as the minimum weight design while satisfying the limit equilibrium condition that the structure may not fail in any of the all possible collapse modes before the specified design ultimate load is reached. Conventional optimum frame designs assume that a continuous spectrum of member size is available. In fact, the vailable sections merely consist of a finite range of discrete member sizes. Optimum frame design using discrete sections has been performed by adopting the plastic collapse theory and using the Complex Method of Box. This study has presented an iterative approach to the optimal plastic design of plane structures that involves the performance of a series of minimum weight design where the limit equilibrium equation pertaining to the critical collapse mode is added to the constraint set for the next design. The critical collapse mode is found by the collapse load analysis that is formulated as a linear programming problem. This area of research is currently being studied. This study would be applied and extended to design the larger and more complex framed structures.