• Geomechanics and Engineering
• /
• 제13권4호
• /
• pp.585-600
• /
• 2017

This study deals with simple solutions for a spherical or circular opening excavated in elastic-brittle plastic rock mass compatible with a linear Mohr-Coulomb (M-C) or a nonlinear Hoek-Brown (H-B) yield criterion. Based on total strain approach, the closed-form solutions of stresses and displacement are derived simultaneously for circular and spherical openings using original H-B and M-C yield criteria. Two simple numerical procedures are proposed for the solution of generalized H-B and M-C yield criteria. Based on incremental approach, the similarity solution is derived for circular and spherical openings using generalized H-B and M-C yield criteria. The classical Runge-Kutta method is used to integrate the first-order ordinary differential equations. Using three data sets for M-C and H-B models, the results of the radial displacements, the spreading of the plastic radius with decreasing pressure, and the radial and circumferential stresses in the plastic region are compared. Excellent agreement among the solutions is obtained for all cases of spherical and circular openings. The importance of the use of proper initial values in the similarity solution is discussed.

• Journal of Magnetics
• /
• 제21권3호
• /
• pp.468-475
• /
• 2016

In this article, entropy generation on MHD Williamson nanofluid over a porous shrinking sheet has been analyzed. Nonlinear thermal radiation and chemical reaction effects are also taken into account with the help of energy and concentration equation. The fluid is electrically conducting by an external applied magnetic field while the induced magnetic field is assumed to be negligible due to small magnetic Reynolds number. The governing equations are first converted into the dimensionless expression with the help of similarity transformation variables. The solution of the highly nonlinear coupled ordinary differential equation has been obtained with the combination of Successive linearization method (SLM) and Chebyshev spectral collocation method. Influence of all the emerging parameters on entropy profile, temperature profile and concentration profile are plotted and discussed. Nusselt number and Sherwood number are also computed and analyzed. It is observed that entropy profile increases for all the physical parameters. Moreover, it is found that when the fluid depicts non-Newtonian (Williamson fluid) behavior then it causes reduction in the velocity of fluid, however, non-Newtonian behavior enhances the temperature and nanoparticle concentration profile.

• 대한기계학회논문집B
• /
• 제23권8호
• /
• pp.945-956
• /
• 1999

This work reports a set of approximate analytical solutions describing the initial transient process of close-contact melting between a rectangular parallelepiped solid and a flat plate on which either constant temperature or constant heat flux is imposed. Not only relative motion of the solid block tangential to the heating plate, but also the density difference between the solid and liquid phase is incorporated in the model. The thin film approximation reduces the force balance between the solid weight and liquid pressure, and the energy balance at the melting front into a simultaneous ordinary differential equation system. The normalized model equations admit compactly expressed analytical solutions which include the already approved two-dimensional solutions as a subset. In particular, the normalized liquid film thickness is independent of all pertinent parameters, thereby facilitating to define the transition period of close-contact melting. A unique behavior of the solid descending velocity due to the density difference is also resolved by the present solution. A new geometric function which alone represents the three-dimensional effect is introduced, and its properties are clarified. One of the representative results is that heat transfer is at least enhanced at the expense of the increase in friction as the cross-sectional shape deviates from the square under the same contact area.

• 한국공작기계학회논문집
• /
• 제13권5호
• /
• pp.38-48
• /
• 2004

Recently the use of high-speed equipment in machine-tool industry has greatly increased, which requires the development of prognostics and prediction methods on the design stage. Conversion of the test/experiments stage from real to virtual reality will not only significantly reduce the design and manufacturing cost, but will also increase design quality. This paper shows how it is possible to develop the automated system for the design calculations of the air-bearing spindles. First, the general calculation method is introduced. It contains several steps, namely, geometry identification, pressure calculation, stiffiness calculation, dynamics characteristics calculation. For geometry identification reducing spindle shaft to rings was proposed, which helps to automate the calculation process. For pressure calculation the Peshti method was implemented. For stiffiness calculation the analysis was made, which shown the necessity of correct calculation step selection. Then the system of ordinary differential equations containing influence coefficients was evolved, which is used for trjectories calculation. The graphical representation of the calculation results shows the dynamic behavior of the spindle unit concerning various working conditions. Finally, this automated system is illustrated by an example of the air-bearing spindle calculation.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2007년도 춘계학술대회 논문집
• /
• pp.1319-1324
• /
• 2007

Many researches have been performed to analyze the smoke movement in tunnel fires by using field model. Recently, FDS(Fire Dynamics Simulator) v.4, which is one of the field model and developed from NIST(National Institute of Standards and Technology), is widely used. In tunnel fires, FDS can show detail results in local point, but it has difficulties in boundary condition and taking long computing time as the number of grid increases. So, there is a need to use alternative method for tunnel fire simulation. A zone model is different kind of CFD method and solves ordinary differential equation based on conservation and auxiliary equations. It shows good macroscopic view in less computing time compared to field model. In this study, therefore, to confirm the applicability of CFAST in tunnel fire analysis, numerical simulations using CFAST are conducted to analyze smoke movement in longitudinal ventilation reduced-scale tunnel fires. Then the results are compared with experimental results. The differences of temperature and critical velocity between numerical results and experimental data are over $30^{\circ}C$ and 0.9m/s, respectively. These values are out of error range. It shows that CFAST 6.0 is hard to be used for tunnel fire simulation.

• 한국강구조학회 논문집
• /
• 제9권3호통권32호
• /
• pp.333-340
• /
• 1997

본 연구에서는 복합재료로 구성된 복합적층판 및 쉘에 대하여 3차 전단 변형이론을 이용한 변위를 가정하여 단순지지 경계조건을 만족하는 변위형상함수를 퓨리예급수로 전개하고 동적 평형 방정식을 유도하여 뉴마크의 수치적분법을 사용하여 단면특성계수, 재료의 특성, 층의 배열에 따른 복합적층판 및 쉘의 비감쇄 동적응답특성을 연구하였다.

• Structural Engineering and Mechanics
• /
• 제60권5호
• /
• pp.891-902
• /
• 2016

The beam structure with breathing cracks subjected to harmonic excitations was modeled by FEM based on Euler-Bernoulli theory, and a piecewise dynamical system was deduced. The precise integration method (PIM) was employed to propose an algorithm for analyzing the dynamic responses of the deduced system. This system was first divided into linear sub-systems, between which there are switching points resulted from the breathing cracks. The inhomogeneous terms due to the external excitations were tackled by introducing auxiliary variables to express the harmonic functions, hence the sub-systems are homogeneous. The PIM was then applied to solve the homogeneous sub-systems one by one. During the procedures, a predictor-corrector algorithm was presented to determine the switching points accurately. The presented method can provide solutions with an accuracy to a magnitude of $10^{-12}$ compared with exact solutions obtained by the theories of ordinary differential equations. The PIM results are much more accurate than Newmark ones with the same time step. Moreover, it is found that the PIM can maintain a high level of accuracy even when the time step increases within a relatively wide range.

• 대한수학회지
• /
• 제47권5호
• /
• pp.985-1000
• /
• 2010

Sufficient conditions for the existence of positive solutions for a coupled system of nonlinear nonlocal boundary value problems of the type -x"(t) = f(t, y(t)), t $\in$ (0, 1), -y"(t) = g(t, x(t)), t $\in$ (0, 1), x(0) = y(0) = 0, x(1) = ${\alpha}x(\eta)$, y(1) = ${\alpha}y(\eta)$, are obtained. The nonlinearities f, g : (0,1) $\times$ (0, $\infty$ ) $\rightarrow$ (0, $\infty$) are continuous and may be singular at t = 0, t = 1, x = 0, or y = 0. The parameters $\eta$, $\alpha$, satisfy ${\eta}\;{\in}\;$ (0,1), 0 < $\alpha$ < $1/{\eta}$. An example is provided to illustrate the results.

• 대한기계학회논문집
• /
• 제16권9호
• /
• pp.1711-1721
• /
• 1992

본 연구에서는 Nayfeh등과 이원경과 소강영은 조화가진하의 핀과 꺽쇠로 고정 된 보(hinged-clamped beam)의 강제 진동 해석을 통하여 점근적으로 안정한 정상상태 응답이 둘 이상 존재할 수 있음을 알게 되었다. 본 연구에서는 이 경우에 접근적으 로 안정한두 정상상태 응답을 구하고 이들 안정한 해로 각각 흡인되는 초기조건들의 집합인 흡인영역을 보간사상법(interpolated mapping method)과 직접 수치적분에 의해 구한 후 서로 비교하였다.

• Journal of Biosystems Engineering
• /
• 제37권2호
• /
• pp.82-89
• /
• 2012

Purpose: The simulation model of a pneumatic conveying drying (PCD) for sawdust was developed and verified with the experiments. Method: The thermal behavior and mass transfer of a PCD were modeled and investigated by comparing the experimental results given by a reference (Kamei et al. 1952) to validate the model. Momentum, energy and mass balance, one dimensional first order ordinary differential equations, were coded and solved into Matlab V. 7.1.0 (2009). Results: The simulation results showed that the moisture content reduced from 194% to 40% (dry basis), air temperature decreased from $512^{\circ}C$ to $128^{\circ}C$ with the particle residence time of 0.7 seconds. The statistical indicators, root mean square error and R-squared, were calculated to be 0.079, and 0.998, respectively, between the measured and predicted values of moisture content. The relative error between the measured and predicted values of the final pressured drop, air temperature, and air velocity were only 8.96%, 0.39% and 1.05% respectively. Conclusions: The predicted moisture content, final temperature, and pressure drop values were in good agreement with the experimental results. The developed model can be used for design and estimation of PCD system for drying of wood dust particles.