• Journal of Mechanical Science and Technology
• /
• 제16권8호
• /
• pp.1072-1078
• /
• 2002

The constrained motion requires the determination of constraint force acting on unconstrained systems for satisfying given constraints. Most of the methods to decide the force depend on numerical approaches such that the Lagrange multiplier method, and the other methods need vector analysis or complicated intermediate process. In 1992, Udwadia and Kalaba presented the generalized inverse method to describe the constrained motion as well as to calculate the constraint force. The generalized inverse method has the advantages which do not require any linearization process for the control of nonlinear systems and can explicitly describe the motion of holonomically and/or nongolonomically constrained systems. In this paper, an explicit equation to describe the constrained motion is derived by minimizing the performance index, which is a function of constraint force vector, with respect to the constraint force. At this time, it is shown that the positive-definite weighting matrix in the performance index must be the inverse of mass matrix on the basis of the Gauss's principle and the derived differential equation coincides with the generalized inverse method. The effectiveness of this method is illustrated by means of two numerical applications.

• 대한기계학회논문집A
• /
• 제24권1호
• /
• pp.190-195
• /
• 2000

It has been established that a crack has an important effect on the dynamic behavior of a structure. This effect depends mainly on the location and depth of the crack. To identify the location and depth of a crack in a structure. Nikolakopoulos et. al. used the intersection point of the superposed contours that correspond to the eigenfrequency caused by the crack presence. However the intersecting point of the superposed contours is not only difficult to find but also incorrect to calculate. A method is presented in this paper which uses optimization technique for the location and depth of the crack. The basic idea is to find parameters which use the structural eigenfrequencies on crack depth and location and optimization algorithm. With finite element model of the structure to calculate eigenfrequencies, it is possible to formulate the inverse problem in optimization format. Method of optimization is augmented lagrange multiplier method and search direction method is BFGS variable metric method and one dimensional search method is polynomial interpolation.

• 한국공작기계학회논문집
• /
• 제12권2호
• /
• pp.71-78
• /
• 2003

This research proposes an implementation method of linearized equations of motion for multibody systems with closed loops. The null space of the constraint Jacobian is first pre-multiplied to the equations of motion to eliminate the Lagrange multiplier and the equations of motion are reduced down to a minimum set of ordinary differential equations. The resulting differential equations are functions of all relative coordinates, velocities, and accelerations. Since the variables are tightly coupled by the position, velocity, and acceleration level coordinates, direct substitution of the relationships among these variables yields very complicated equations to be implemented. As a consequence, the reduced equations of motion are perturbed with respect to the variations of all variables, which are coupled by the constraints. The position velocity and acceleration level constraints are also perturbed to obtain the relationships between the variations of all relative coordinates, velocities, and accelerations and variations of the independent ones. The Perturbed constraint equations are then simultaneously solved for variations of all variables only in terms of the variations of the independent variables. Finally, the relationships between the variations of all variables and these of the independent ones are substituted into the variational equations of motion to obtain the linearized equations of motion only in terms of the independent variables variations.

• 응용통계연구
• /
• 제19권3호
• /
• pp.461-479
• /
• 2006

본 논문에서는 이원오차성분을 갖는 패널회귀모형에서 모형식별을 위하여 LM 검정통계량을 유도하고 검정통계량의 연산을 위하여 인공회귀방법(Double-Length Artificial Regression, DLR)을 이용한다. 모의 실험 결과, 소표본의 경 우에는 Outer-Product Gradient(OPG)에 근거한 LM 검정통계량은 유위수준이 과대기각하는 경향을 보인 반면 DLR에 근거한 LM 검정통계량은 명목유의수준을 잘 유지하고 검정력도 높게 나타났다.

• Industrial Engineering and Management Systems
• /
• 제4권2호
• /
• pp.145-157
• /
• 2005

There are a lot of raw materials, work-in-processes and finished goods in manufacturing industry. Here, the less stock of materials and work-in-processes manufacturing industry has, the worse the rate of the production is. Inversely, the more manufacturing industry has, the more expensive the cost to support them is. Thus, it is important for us to balance them efficiently. In general, inventory problems are to decide appropriate times to produce goods and to determine appropriate quantities of goods. Therefore, inventory problems require as more useful information as possible. For example, there are demand, lead time, ordering point and so on. In this paper, we deal with an optimal ordering policy on both way substitutable two-commodity inventory control system. That is, there is a problem of how to allocate the produced two kinds of goods in a factory to m areas so as to minimize the total expected inventory cost. The demand of each area is probabilistic, and we adopt the exponential distribution as a probability density function of demand. Moreover, we provide numerical examples of the problem.

• 한국항공우주학회지
• /
• 제38권2호
• /
• pp.141-149
• /
• 2010

위성 편대 비행 시 주위성과 부위성 간의 연료 소비 균형을 고려한 임펄스 기동에 관한 연구를 수행하였다. 위성 간 사용가능한 연료량을 비교하여 가중치(weight)를 두고 가격함수(cost function)를 설계하여 라그랑지 승수법을 통해 필요한 임펄스를 획득하였다. 상대궤도 발산 방지를 위해 에너지 매칭 기법을 사용하였고, 임펄스 기동 후 상대 거리 구속이 이루어짐을 시뮬레이션을 통해 확인하였다. 시뮬레이션은 지구 중력 외의 외란이 없는 경우와 대기 항력이 외란으로 존재하는 상황으로 시나리오를 나누어 수행하였다. 본 논문 결과는 이후 실제로 위성을 편대로 사용한 위성 군집 비행 시, 상대 궤도 구속 요건을 만족하고 각 위성의 연료량을 비교한 임펄스 기동이 요구될 때 사용가능할 것으로 기대된다.

• 전자공학회논문지CI
• /
• 제38권2호
• /
• pp.16-26
• /
• 2001

본 논문에서는 주어진 필터와 구현 복잡도에 대하여 최대 코딩이득을 내는 부 주파수 분활을 가진 서브밴드 부호화기를 구축하는 고속 알고리듬을 제안한다. 이를 위하여 본 논문에서는 직교 기저 및 비 직교 기저와 임의의 부 주파수 분할에 대하여 적용할 수 있는 통합적인 코딩이득의 식을 유도한 다음, 부 주파수 대역수에 대하여 코딩이득이 단순 증가 함수임을 증명한다, 이를 바탕으로 복잡도에 대하여 최대 코딩이득을 내는 최적화 된 부 주파수 분할을 찾아내기 위하여 그 단순 증가 함수를 부 주파수 대역 수에 따른 왜곡 함수로 다룬다. 이 왜곡 함수을 목적함수로 두고 Lagrange 방법에 근거하여 최적화 된 해를 고속으로 제공하는 알고리듬을 개발한다.

• 대한기계학회논문집
• /
• 제19권3호
• /
• pp.639-646
• /
• 1995

The objective of this study is to minimize the weight of a damped anisotropic roto-bearing system considering whirl natural frequency and stability. The system is modeled as an assemblage of rigid disks, flexible shafts and discrete bearings. The system design variables are the crosssectional areas of shaft elements and the properties of bearings. To analyze the system, the polynomial method which is derived by rearranging the calculations performed by a transfer matrix method is adopted. For the optimization, the optimization software IDOL (Integrated Design Optimization Library) which is based on the Augmented Lagrange Multiplier (ALM) method is employed. Also, an analytical design sensitivity analysis of the system is used for high accuracy and efficiency. To demonstrate the usefulness of the proposed optimal design program incorporating analysis, design sensitivity analysis, and optimization modules, a damped anisotropic rotor-bearing system is optimized to obtain 34$ weight reduction.

• Wind and Structures
• /
• 제22권3호
• /
• pp.349-368
• /
• 2016

A preconditioning technique is presented for a simultaneous solution to wind-membrane interaction. In the simultaneous equations, a linear elastic model was employed to deal with the fluid-structure data transfer at the interface. A Lagrange multiplier was introduced to impose the specified boundary conditions at the interface and strongly coupled simultaneous equations are derived after space and time discretization. An initial linear elastic model preconditioner and modified one were derived by treating the linearized elastic model equation as a saddle point problem, respectively. Accordingly, initial and modified fluid-structure interaction (FSI) preconditioner for the simultaneous equations were derived based on the initial and modified linear elastic model preconditioners, respectively. Wind-membrane interaction analysis by the proposed preconditioners, for two and three dimensional membranous structures respectively, was performed. Comparison was made between the performance of initial and modified preconditioners by comparing parameters such as iteration numbers, relative residuals and convergence in FSI computation. The results show that the proposed preconditioning technique greatly improves calculation accuracy and efficiency. The priority of the modified FSI preconditioner is verified. The proposed preconditioning technique provides an efficient solution procedure and paves the way for practical application of simultaneous solution for wind-structure interaction computation.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.395-402
• /
• 2006

In analyzing the nano-scale phenomena or behaviors of nano devices or materials, it is often desirable to deal with more atoms than can be treated only with a full atomistic simulation. However, even now, it is advisable to apply the atomistic simulation to the narrow region where the deformation field changes rapidly but to apply the conventional continuum model to the region far from that region. This equivalent continuum model can be formulated by applying the Cauchy-Born rule to the exact atomistic potential as in the quasicontinuum method. To couple the atomistic model with the equivalent continuum model, continuum displacements are conformed to the molecular displacements at the discrete positions of the atoms within the bridging domain. To satisfy the coupling constraints, we apply the Lagrange multiplier method. The continuum model in the bridging model should be applied on the region where the deformation field changes gradually. Then we can make the nodal spacing in the continuum model be much larger than the atomic spacing. In the first step, we generate the atomic-resolution mesh with the nodal spacing equal to the atomic spacing, and then we eliminate the nodal degrees of freedom adaptively using the node deactivation techniques. We eliminate more DOFs as the regions are more far from the atomistic region. Computing time and computational resources can be greatly reduced by the present node deactivation technique in multi scale analysis.