• Structural Engineering and Mechanics
• /
• v.24 no.6
• /
• pp.709-725
• /
• 2006

In this paper the buckling and post-buckling behavior of slender bars under self-weight are studied. In order to study the post-buckling behavior of the bar, a geometrically exact formulation for the non-linear analysis of uni-directional structural elements is presented, considering arbitrary load distribution and boundary conditions. From this formulation one obtains a set of first-order coupled nonlinear equations which, together with the boundary conditions at the bar ends, form a two-point boundary value problem. This problem is solved by the simultaneous use of the Runge-Kutta integration scheme and the Newton-Raphson method. By virtue of a continuation algorithm, accurate solutions can be obtained for a variety of stability problems exhibiting either limit point or bifurcational-type buckling. Using this formulation, a detailed parametric analysis is conducted in order to study the buckling and post-buckling behavior of slender bars under self-weight, including the influence of boundary conditions on the stability and large deflection behavior of the bar. In order to evaluate the quality and accuracy of the results, an experimental analysis was conducted considering a clamped-free thin-walled metal bar. As this kind of structure presents a high index of slenderness, its answers could be affected by the introduction of conventional sensors. In this paper, an experimental methodology was developed, allowing the measurement of static or dynamic displacements without making contact with the structure, using digital image processing techniques. The proposed experimental procedure can be used to a wide class of problems involving large deflections and deformations. The experimental buckling and post-buckling behavior compared favorably with the theoretical and numerical results.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.259-267
• /
• 2007

An efficient algorithm for inverting static-shifted magnetotelluric (MT) data has been proposed to produce a three-dimensional (3D) resistivity model. In the Gauss-Newton approach, computational costs associated with construction of a full sensitivity matrix usually make 3D MT inversion impractical. This computational difficulty may be overcome by using approximate sensitivities. We use four kinds of sensitivities in particular orders in the inversion process. These sensitivities are computed 1) analytically for an initial, homogeneous earth, 2) exactly for a current model, 3) approximately by the Broyden method, and 4) approximately using the previous adjoint fields. Inversion experiments with static-shifted synthetic and field MT data indicate that inversion results are highly dependent on characteristics of data and thus applying various combinations of sensitivities is helpful in obtaining a good image of the subsurface structure with reasonable computation time.

• Journal of rehabilitation welfare engineering & assistive technology
• /
• v.7 no.1
• /
• pp.21-27
• /
• 2013

Manual wheelchair propulsion can lead to pain and injuries of users due to mechanical inefficiency of wheelchair propulsion motion. The kinetic analysis of the upper limbs during manual wheelchair propulsion needs to be studied. A two dimensional inverse dynamic model of upper limbs was developed to compute the joint torque during manual wheelchair propulsion. The model was composed of three segments corresponding to upper arm, lower arm and hand. These segments connected in series by revolute joints constitute open chain mechanism in sagittal plane. The inverse dynamic method is based on Newton-Euler formalism. The model was applied to data collected in experiments. Kinematic data of upper limbs during wheelchair propulsion were obtained from three dimensional trajectories of markers collected by a motion capture system. Kinetic data as external forces applied on the hand were obtained from a dynamometer. The joint rotation angles and joint torques were computed using the inverse dynamic model. The developed model is for upper limbs biomechanics and can easily be extended to three dimensional dynamic model.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.40 no.1
• /
• pp.88-95
• /
• 2015

This paper presents the study of the efficient localization algorithm for non-linear least square estimation. Although non-linear least square(NLS) estimation algorithms are more accurate algorithms than linear least square(LLS) estimation, NLS algorithms have more computation loads because of iterations. This study proposed the efficient algorithm which reduced complexity for small accuracy loss in NLS estimation. Simulation results show the accuracy and complexity of the localization system compared to the proposed algorithm and conventional schemes.

• Journal of Astronomy and Space Sciences
• /
• v.4 no.1
• /
• pp.11-23
• /
• 1987

In dward nova we assume the primary star as a white dwarf and the secondary as the late type star which filled Roche lobe. Mass flow from the secondary star leads to the formation of thin accretion disk around the white dwarf. We use the $\alpha$ parameter as viscosity to maintain the disk form and propose that the outburst in dwarf nova cause the step increase of source term. With these assumptions we solve the basic equations of stellar structure using New-Raphson method. We show the physical parameters like temperature, density, pressure, opacity, surface density, height and flux to the radius of disk. Changing the value of $\alpha$, we compare several parameters when mass flow rate is constant with those of when luminosity of disk is brightest. At the same time, we obtain time-dependent variations of luminosity and mass of disk. We propose the suitable range of $\alpha$is 0.15-0.18 to the difference of luminosity. We compare several parameters of disk with those of the normal late type stars which have the same molecular weight of disk. These show the temperature and pressure of disk are similar to those of normal stars but the density of disk is lower. Maybe the outburst in dwarf nova is due to the variation of the $\alpha$ value instead of increment of mass flow from the secondary star.

• Tribology and Lubricants
• /
• v.30 no.6
• /
• pp.337-342
• /
• 2014

The reduction of drag torque is an important research issue in terms of improving transmission efficiency. Drag torque in a wet clutch occurs because of the viscous drag generated by the transmission fluid in a narrow gap (clearance) between the friction plate and a separate plate. The objective of this paper is to observe the effects of the friction plate area and the clearance on the drag torque using finite element simulation. The two-phase flow of air and oil fluid is considered and modeled for the simulation. The simulation analysis reveals that as the rotational speed increases, the drag torque generally increases to a critical point and then decreases sharply at a high speed regime. The clearance between the two plates plays an important role in controlling drag torque peak. An increase in the clearance causes a decrease in shear stress; thus, the drag torque also decreases according to Newton's law of viscosity. An observation of the effect of the area of contact between transmission fluid and friction plate shows that the drag torque increases with the contact area. The flow vectors inside the flow channel present clear evidence that the velocity of the fluid flows is faster with a larger friction plate, that is, in the case of a larger contact area. Therefore, the optimum size of the friction plate should be determined carefully, considering both the clutch performance and drag reduction. It is expected that the results from this study can be very useful as a database for clutch design and to predict the drag torque for the initial design with respect to various clutch parameters.

• Communications of Mathematical Education
• /
• v.31 no.2
• /
• pp.139-152
• /
• 2017

In this paper, we examine how secondary school mathematics textbooks on calculus introduce the history of calculus. In order to identify the problem, we consider the Babylonian integration by trapezoidal rule, which was made to calculate the location of Jupiter in 350-50 B.C., and the integration by the method of the rotating plate of ibn al-Haytham in Egypt, about 1000 years. In conclusion, our secondary school mathematics textbooks describe Newton and Leibniz as inventing calculus and place their roots in ancient Greece. The origin of the calculus is in Babylonia and the Faṭimah Dynasty (909-1171) (Egypt) and it is desirable that the calculus is developed in Europe after the development of the power series in India, and that the value of Asia Africa is introduced in the textbooks.

• Geomechanics and Engineering
• /
• v.21 no.4
• /
• pp.357-369
• /
• 2020

Different from the conventional planar fracture and simplified Newton model, for power-law slurries with a lower water-cement ratio commonly used in grouting engineering, flow model in geological rough fractures is built based on ten standard profiles from Barton (1977) in this study. The numerical algorithm is validated by experimental results. The flow mechanism, grout superiority, and water plugging of pseudo plastic slurry are revealed. The representations of hydraulic grouting properties for JRCs are obtained. The results show that effective plugging is based on the mechanical mechanisms of the fluctuant structural surface and higher viscosity at the middle of the fissure. The formulas of grouting parameters are always variable with the roughness and shear movement, which play a key role in grouting. The roughness can only be neglected after reaching a threshold. Grouting pressure increases with increasing roughness and has variable responses for different apertures within standard profiles. The whole process can be divided into three stationary zones and three transition zones, and there is a mutation region (10 < JRCs < 14) in smaller geological fractures. The fitting equations of different JRCs are obtained of power-law models satisfying the condition of -2 < coefficient < 0. The effects of small apertures and moderate to larger roughness (JRCs > 10.8) on the permeability of surfaces cannot be underestimated. The determination of grouting parameters depends on the slurry groutability in terms of its weakest link with discontinuous streamlines. For grouting water plugging, the water-cement ratio, grouting pressure and grouting additives should be determined by combining the flow conditions and the apparent widths of the main fracture and rough surface. This study provides a calculation method of grouting parameters for variable cement-based slurries. And the findings can help for better understanding of fluid flow and diffusion in geological fractures.

• Plant Journal
• /
• v.13 no.4
• /
• pp.48-50
• /
• 2017

A program for analyzing the transient behaviors of industrial gas turbines was developed. Each component (compressor, combustor, turbine and ducts)of gas turbine is modeled as a fully module to enhance the expandability of the program. We used object-oriented programing for this purpose. The mass and energy balance equations are solved numerically by Multivariable Newton Raphson method. The characteristic maps for the compressor and turbine were used for predicting the performance of a gas turbine engine. Combustion in the combustor is assumed to be complete combustion. PID control is used to maintain constant rotational speed and turbine exhaust temperature by the control of the fuel flow rate and the changing of the compressor inlet guide vane angle at the same time. It was confirmed that stable control of the gas turbine was possible, even for a rapid load change.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.