• Journal of the Korean Society of Combustion
• /
• v.12 no.2
• /
• pp.10-19
• /
• 2007

Flame surface area is a critical parameter determining turbulent flame speed. Three-dimensional direct numerical simulations(DNS) were conducted to figure out the evolution process of flame surface area. Fully compressible Navier-Stokes equations are solved to reproduce premixed flame embedded in isotropic decaying turbulent flow. The tangential straining and curvature of propagating surface affect development of flame area. In this study, four different turbulent intensity flows and three different Le number flames are investigated to force changes in straining and curvature effects. Consistent results are obtained for the probability density functions (PDF) of strain and curvature with previous researches. It is revealed that displacement speed, which is a speed of flame surface relative to unburnt flow, controls the balance between sink and source of flame surface area.

• Water Engineering Research
• /
• v.3 no.4
• /
• pp.247-258
• /
• 2002

This paper presents a non-isotropic turbulence modeling of flows over bedforms. The Reynolds stress model is used for the turbulence closure. In the model, Launder, Reece, and Rodi's model and Hanjalic and Launder's model are employed f3r the pressure strain correlation term and the diffusion term, respectively. The mean flow and turbulence structures are simulated and compared with profiles measured in the experiments. The numerical solutions from two-equation turbulence models are also provided for comparisons. The Reynolds stress model yields the separation length of eddy similar to the other numerical results. Using the developed model, the resistance coefficients are also estimated for the flows at different Froude numbers. Karim's (1999) relationship is used to determine the bedform geometry. It is found that the values of the form drag and the skin friction are very similar to those obtained by the other turbulence models. meaning higher values of the form drag and lower values of the skin friction compared with the empirical formulas.

• Geotechnical Engineering
• /
• v.12 no.4
• /
• pp.145-156
• /
• 1996

An implicit stress integration algorithm was formulated for implementing an aiusotorpic hardening constitutive model which has been based op the generalized isotropic hardening rule in nonlinear finite element analysis technique. the rate form of stress tensor was implicitly integrated using the generalized trapezoidal rule and the tangent stress-strain modulus was evaluated consistently with the nonlinear solution technique. As a result, it has been found that the nonlinear analysis with the anisotropic hardening constitutive model might be performed accurately and efficiently.

• Transactions of Materials Processing
• /
• v.15 no.8 s.89
• /
• pp.533-538
• /
• 2006

A simplified numerical procedure to predict drawbead restraining forces(DBRF) has been developed based on the hybrid membrane/bending method which superposes bending effects onto membrane solutions. As a semi-analytical method, the new approach is especially useful to analyze the effects of various constitutive parameters. The present model can accommodate general anisotropic yield functions along with non-linear isotropic-kinematic hardening under the plane strain condition. For the preliminary results, several sensitivity analyses for the process and material effects such as friction, drawbead depth, hardening behavior including the Bauschinger effect and yield surface shapes on the DBRF are carried out.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.2
• /
• pp.149-156
• /
• 2012

In fracture mechanics, the weight function can be used for calculating stress intensity factors. In this paper, a two-dimensional electroelastic analysis is performed on a transversely isotropic piezoelectric material with an open crack. A plane strain formulation of the piezoelectric problem is solved within the Leknitskii formalism. Weight function theory is extended to piezoelectric materials. The stress intensity factors and electric displacement intensity factor are calculated by the weight function theory.

• Earthquakes and Structures
• /
• v.13 no.5
• /
• pp.465-471
• /
• 2017

This work derives a generalized time fractional differential equation governing wave propagation in a radially vibrating non-classical cylindrical medium. The cylinder is made of a transversely isotropic hyperelastic John's material which obeys frequency-dependent power law attenuation. Employing the definition of the conformable fractional derivative, the solution of the obtained generalized time fractional wave equation is expressed in terms of product of Bessel functions in spatial and temporal variables; and the resulting wave is characterized by the presence of peakons, the appearance of which fade in density as the order of fractional derivative approaches 2. It is obtained that the transversely isotropic structure of the material of the cylinder increases the wave speed and introduces an additional term in the wave equation. Further, it is observed that the law relating the non-zero components of the Cauchy stress tensor in the cylinder under consideration generalizes the hypothesis of plane strain in classical elasticity theory. This study reinforces the view that fractional derivative is suitable for modeling anomalous wave propagation in media.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.70-81
• /
• 1996

Solutions of geotechnical engineering problems require predictions of deformation and stresses during various stages of loading. Powerful numerical methods are available to make such predictions even for complicated problems. To get accurate results, realistic stress-strain relationships of soils are dependent on a number of factors such as soil type, density, stress level and stress path. Attempts are continuously being made to develope analytical models for soils incorporating all such factors. Isotropic compression-expansion test and a series of drained conventional triaxial tests with several stress path for Baekma river sand were performed to investigate stress-strain and volume change characteristics of Lade's single work hardening model dependant on the stress path. In order to predicted of stress-strain and volumetric strain behavior were determined the values of parameters for the mode by the computer program based on the regression analysis. Predicted stress-strain behavior of triaxial compression tests and optional stress path tests for increasing confining pressure with parameters obtained conventional triaxial compression tests agreed with several test results but the prediction results for decreasing confining pressure reduced triaxial compression tests make a little difference with test results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.5 no.2
• /
• pp.122-130
• /
• 1981

Equation of equilibrium is derived and solved for an infinite isotropic solid under applied hydrostatic stress which is uniform at large distance, and disturbed by a spherical precipitate which has isotropoc elastic constants dirrerent form those of the matrix. A linear strain hardening behavior of the matrix is assumed, and an elasto-plastic sloution is obtained. The difference of the total strain energy stored inthe infinite solid with and without a precipitate is computed, and compared with that for purely elastic case. Finally the effect of the ratio of the bulk modulus of the precipitate to that of the matrix and the effct of linear strain hardening rate on the plastic zone size and the energy difference are discussed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.2056-2061
• /
• 2001

Non-proportional loading tests of Nylon 66 at room temperature exhibit path dependent behavior and plasticity-relaxation interactions. The uniaxial formulation of the viscoplasticity theory based on overstress (VBO), which has been used to reproduce the nonlinear strain rate sensitivity, relaxation, significant recovery and cyclic softening behaviors of Nylon 66, is extended to three-dimensions to predict the response in strain-controlled, comer-path tests. VBO consists of a flow law that is easily written for either the stress or the strain as the independent variable. The flow law depends on the overstress, the difference between the stress and the equilibrium stress that is a state variable in VBO. The evolution law of the equilibrium stress in turn contains two additional state variables, the kinematic stress and the isotropic stress. The simulations show that the constitutive model is competent at modeling the deformation behavior of Nylon 66 and other solid polymers.

• Computers and Concrete
• /
• v.1 no.4
• /
• pp.433-455
• /
• 2004

The paper presents results of FE-calculations on shear localizations in quasi-brittle materials during both an uniaxial plane strain compression and uniaxial plane strain extension. An elasto-plastic model with a linear Drucker-Prager type criterion using isotropic hardening and softening and non-associated flow rule was used. A non-local extension was applied in a softening regime to capture realistically shear localization and to obtain a well-posed boundary value problem. A characteristic length was incorporated via a weighting function. Attention was focused on the effect of mesh size, mesh alignment, non-local parameter and imperfections on the thickness and inclination of shear localization. Different methods to calculate plastic strain rates were carefully discussed.