• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권1호
• /
• pp.1-21
• /
• 2015

This paper presents two algorithms based on the Jamshidian equation which is from the Black-Scholes partial differential equation. The first algorithm is for American call options and the second one is for American put options. They compute numerically free boundary and then option price, iteratively, because the free boundary and the option price are coupled implicitly. By the upwind finite-difference scheme, we discretize the Jamshidian equation with respect to asset variable s and set up a linear system whose solution is an approximation to the option value. Using the property that the coefficient matrix of this linear system is an M-matrix, we prove several theorems in order to formulate a bisection method, which generates a sequence of intervals converging to the fixed interval containing the free boundary value with error bound h. These algorithms have the accuracy of O(k + h), where k and h are step sizes of variables t and s, respectively. We prove that they are unconditionally stable. We applied our algorithms for a series of numerical experiments and compared them with other algorithms. Our algorithms are efficient and applicable to options with such constraints as r > d, $r{\leq}d$, long-time or short-time maturity T.

• 한국주조공학회지
• /
• 제22권3호
• /
• pp.130-136
• /
• 2002

The growing use of aluminum for castings over the past decade has brought with it the increased scrutiny of component properties. One area that has received much attention is the effect of in inclusions - or impurities particles held in the metal - on casting properties. A new method of electromagnetic separation for removal of inclusions in aluminum alloy melts is proposed. The principle is that as the electromagnetic force induced in metal acts on inclusions due to low electric conductivity, they are moved to the direction opposite to electromagnetic force and can be separated and removed from the melt. Experiments were carried out on A356 melt mixed alumina particles and commercial Al alloys of ADC 10 and 12. In the experiment using A356, it was proved that $Al_2O_3$ particles was separated and removed continuously from matrix melt by electromagnetic force. Based on these results, the continuous separation experiment that used ADC 10, 12 was carried and the cleanliness of melt was assessed by the amount of porosity, hydrogen contents, PoDFA and mechanical properties through tensile test. As the results of analyses, the amount of porosity and hydrogen contents decreased without variation of chemical composition in the specimen that passed the electromagnetic continuous separator. In addition, tensile strength and elongation of this specimen increased by $20{\sim}30%$ because of reduction of inclusions.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.6-15
• /
• 2008

This paper deals with the aeroelastic instability of vibrating multiple blade rows under aerodynamic coupling with each other. A model composed of three blade rows, e.g., rotor-stator-rotor, where blades of the two rotor cascades are simultaneously vibrating, is considered. The displacement of a blade vibrating under aerodynamic force is expanded in a modal series with the natural mode shape functions, and the modal amplitudes are treated as the generalized coordinates. The generalized mass matrix and the generalized stiffness matrix are formulated on the basis of the finite element concept. The generalized aerodynamic force on a vibrating blade consists of the component induced by the motion of the blade itself and those induced not only by vibrations of other blades of the same cascade but also vibrations of blades in another cascade. To evaluate the aerodynamic forces, the unsteady lifting surface theory for the model of three blade rows is applied. The so-called k method is applied to determine the critical flutter conditions. A numerical study has been conducted. The flutter boundaries are compared with those for a single blade row. It is shown that the effect of the aerodynamic blade row coupling substantially modifies the critical flutter conditions.

• Journal of Advanced Marine Engineering and Technology
• /
• 제25권5호
• /
• pp.1130-1139
• /
• 2001

Springs are widely utilized in machine element. To find out stiffness of coil spring, the space beam theory using the finite element method is adopted in this paper. In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displacements are added to coordinates of nodal points. The new stiffness matrix of the system using the new coordinates of nodal points is adopted to calculated the another increments of nodal displacements, that is, the step by step method is used in this paper. The results of the finite element method are fairly well agreed with those of various experiments. Using MATLAB program developed in this paper, spring constants can be predicted by input of few factors.

• 한국마린엔지니어링학회:학술대회논문집
• /
• 한국마린엔지니어링학회 2001년도 추계학술대회 논문집(Proceeding of the KOSME 2001 Autumn Annual Meeting)
• /
• pp.21-27
• /
• 2001

Springs are widely utilized in machine element. To find out stiffness of frustum-shaped coil spring, the space beam theory using the finite element method is adopted in this paper In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. To find out load vector of coil spring subjected to distributed compression, principle of virtual work is adapted The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displacements are added to coordinates of nodal points. The new stiffness matrix of the system using the new coordinates of nodal points is adopted to calculate the another increments of nodal displacements, that is, the step by step method is used in this paper. The results of the finite element method are fairly well agreed with those of various experiments. Using MATLAB program developed in this paper, spring constants and stresses can be predicted by input of few factors.

• 한국해양공학회지
• /
• 제33권5호
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• 한국주조공학회지
• /
• 제15권4호
• /
• pp.388-396
• /
• 1995

This study has been conducted with the purpose of examining the fatigue crack growth characteristics of $Al_2O_3$ short fiber reinforced aluminum matrix composites made by squeeze casting process with different applied pressure and binder amount. Fatigue crack growth experiments have been performed under constant load amplitude method with a fixed load ratio. The rate of crack propagation was decreased with binder amount as well as applied pressure. Also fatigue crack growth path in matrix was changed from flat to rough mode with an increase of applied pressure. In the composites, fatigue crack was propagated to interface between matrix and reinforcement at 10MPa, but it was propagated to reinforcement at 20MPa. The major reason of thee result was considered that interfacial bonding force and microstructure of matrix were improved due to an increase of applied pressure. Localized ductile striation in the composites was observed at low growth rate region and such a phenominon was remarkable with an increase of applied pressure. At high growth rate region, the propensity of fracture appearance was changed from interfacial debonding to reinforcement fracture with an increase of applied pressure.

• 대한기계학회논문집A
• /
• 제24권8호
• /
• pp.1931-1939
• /
• 2000

Operating excitation forces of the linear vibratory system are normally determined by direct measurement techniques using load cells, strain gauges, etc. But, hydraulic forces of the rotating turbomachinery such as centrifugal pumps are exerted on an impeller due to asymmety of the flow by the interaction between pump impeller and volute. So, investigations of wide range of hydraulic designs and geometric deviations are difficult by direct method. This paper presents a hybrid approach for fourier transformed operational excitation forces, which uses pseudo-inverse matrix of the transfer matrix for the system and the measured vibrational data with standard installed pump. The determination of the transfer function matrix is based on a linear rotor/stationary system and steady state harmonic response in finite element analysis. And, vibrational data is collected in both vertical and horizontal directions at inboard and outboard bearing housings. The results of the process may be enhanced by making acceleration measurements at many more locations than there are forces to be determined.

• Steel and Composite Structures
• /
• 제39권5호
• /
• pp.543-564
• /
• 2021

In this study free vibration analysis of a cracked Goland composite wing is investigated. The wing is modelled as a cantilevered beam based on Euler- Bernoulli equations. Also, composite material is modelled based on lamina fiber-reinforced. Edge crack is modelled by additional boundary conditions and local flexibility matrix in crack location, Castigliano's theorem and energy release rate formulation. Governing differential equations are extracted by Hamilton's principle. Using the separation of variables method, general solution in the normalized form for bending and torsion deflection is achieved then expressions for the cross-sectional rotation, the bending moment, the shear force and the torsional moment for the cantilevered beam are obtained. The cracked beam is modelled by separation of beam into two interconnected intact beams. Free vibration analysis of the beam is performed by applying boundary conditions at the fixed end, the free end, continuity conditions in the crack location of the beam and dynamic stiffness matrix determinant. Also, the effects of various parameters such as length and location of crack and fiber angle on natural frequencies and mode shapes are studied. Modal analysis results illustrate that natural frequencies and mode shapes are affected by depth and location of edge crack and coupling parameter.

• 드라이브 ㆍ 컨트롤
• /
• 제14권4호
• /
• pp.45-50
• /
• 2017

The gas spring is a hydropneumatic adjusting element, consisting of a pressure tube, a piston rod, a piston and a connection fitting. The gas spring is filled with compressed nitrogen within the cylinder. The filling pressure acts on both sides of the piston and because of area difference it produces an extension force. Therefore, a gas spring is similar in function compare to mechanical coil spring. Conversely, optimization is a process of finding the best set of parameters to reach a goal while not violating certain constraints. The AMESim software provides NLPQL (Nonlinear Programming by Quadratic Lagrangian) and GA (genetic algorithm) for optimization. The NLPQL method builds a quadratic approximation to the Lagrange function and linear approximations to all output constraints at each iteration, starting with the identity matrix for the Hessian of the Lagrangian, and gradually updating it using the BFGS method. On each iteration, a quadratic programming problem is solved to find an improved design until the final convergence to the optimum design. In this study, we conducted optimization design of the gas spring reaction force with NLPQL.