• Steel and Composite Structures
• /
• v.18 no.6
• /
• pp.1465-1478
• /
• 2015

The effect of cutting off fibers on transient load in a polymeric matrix composite lamina was studied in this paper. The behavior of fibers was considered to be linear elastic and the matrix behavior was considered to be linear viscoelastic. To model the viscoelastic behavior of matrix, a three parameter solid model was employed. To conduct this research, finite difference method was used. The governing equations were obtained using Shear-lag theory and were solved using boundary and initial conditions before and after the development of break. Using finite difference method, the governing integro-differential equations were developed and normal stress in the fibers is obtained. Particular attention is paid the dynamic overshoot resulting when the fibers are suddenly broken. Results show that considering viscoelastic properties of matrix causes a decrease in dynamic load concentration factor and an increase in static load concentration factor. Also with increases the number of broken fibers, trend of increasing load concentration factor decreases gradually. Furthermore, the overshoot of load in fibers adjacent to the break in a polymeric matrix with high transient time is lower than a matrix with lower transient time, but the load concentration factor in the matrix with high transient time is lower.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.8
• /
• pp.27-36
• /
• 2013

LEM (Limit Equilibrium Method) based programs are commonly used for the designs of soil nailing as a slope reinforcement. However, there is a drawback that the interaction between ground and soil nailing is not properly reflected in those programs, which needs to be solved. For economical constructions and designs, research is also required on the support pattern of soil nailing. In this study, therefore, reinforcement effects of soil nailing were compared and analyzed by performing finite difference analyses which could properly consider the interaction between ground and soil nailing. As a result, when the angle from slope to nail is $90^{\circ}$, failure slip surface becomes the largest and thus the factor of safety becomes maximum.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.1
• /
• pp.15-24
• /
• 1986

Pressure trasients must deal with safety problem of system. For identification of physical situation that can and method of limiting surges are essential consideration in sucessful design. The finite difference equation by method of characteristics are derived from the governing equation of unsteady flow in a pipe, and solved by using boundary condition derived. A computer program which can simulate general hydraulic system is developed by using finite difference equations and boundary conditions derived. The sumulated resulted by developed computer program are in fair agreement with experiment result.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.2 no.2
• /
• pp.20-25
• /
• 2011

In this paper, two dimensional concrete slabs for a railroad bridge were analyzed by the specially orthotropic laminates theory. Both the geometrical and material property of the cross section of the slab was considered symmetrically with respect to the neutral surface so that the bending extension coupling stiffness, $B_{ij}$ = 0, and $D_{16}=D_{26}=0$ Bridge deck behaves as specially orthotropic plates. In general, the analytical solution for such complex systems is very difficult to obtain. Thus, finite difference method was used for analysis of the problem. In this paper, the finite difference method and the beam theory were used for analysis.

• Solar Energy
• /
• v.5 no.2
• /
• pp.77-83
• /
• 1985

An analysis is performed to investigate the influence of the buoyancy force and the thickness variation of melting layer in the containment that is filled with phase-change Material surrounding a cylindrical heating tube during melting process. The phase-change material is assumed to be initially solid at its phase-change temperature and the remaining solid at any given time is still at the phase-change temperature and neglecting the effect of heat transfer occuring within the solid. At the start of melting process, the thickness of melting layer is assumed to be a stefan-problem and after the starting process, the change of temperature and velocity is calculated using a two dimensional finite difference method. The governing equations for velocity and temperature are solved by a finite difference method which used SIMPLE (Semi Implicit Method Pressure linked Equations) algorithm. Results are presented for a wide range of Granshof number and in accordance with the time increment and it is founded that two dimensional fluid flow occurred by natural convection decreases the velocity of melting process at the bottom of container. The larger the radius of heating tube, the higher heat transfer is occurred in the melting layer.

• Geomechanics and Engineering
• /
• v.21 no.1
• /
• pp.11-21
• /
• 2020

The traditional formulations for estimation of bearing capacity in rock mechanics assume a homogeneous and isotropic rock mass. However, it is common that the rock mass consists of different layers of different rock properties or of the same rock matrix with distinct geotechnical quality levels. The bearing capacity of a heterogeneous rock is estimated traditionally through the weighted average. In this paper, the solution of the weighted average is compared to the finite difference method applied to a bilayer rock mass. The influence of different parameters such as the thickness of the layers, the rock type, the uniaxial compressive strength and the overall geotechnical quality of the rock mass on the bearing capacity of a bilayer rock mass is analyzed. A parametric study by finite difference method is carried out to develop a bearing capacity factor in function of the layer thickness and the rock mass quality expressed in terms of the geological strength index, which is presented in a form of a chart. Therefore, this correlation factor allows estimating the bearing capacity of a rock mass that is formed by two layers with distinct GSI, depending on the bearing capacity of the rock mass formed only by the upper layer and considered by that way as homogenous and isotropic rock mass.

• Journal of Industrial Technology
• /
• v.18
• /
• pp.195-202
• /
• 1998

The effect of the number of nodes on the heat loss from a rectangular fin for a finite difference method is studied. There are two ways for selecting nodes for the upper half fin in this finite difference method. In the first place, all the ${\Delta}x$ are the same and all the ${\Delta}y$ are the same for the entire upper half fin. Incremental length of x (i.e. ${\Delta}x$) is divided by two near the fin tip while all the ${\Delta}y$ are the same for another way. The results show that 1) About 30 nodes are enough to obtain the satisfactory exact analysis (relative error < 5%) on the heat loss for a given range of Biot number in case of short fin (i.e. $L{\leq}2$). 2) Under usual circumstances (Bi<0.1), the relative error of heat loss between using 30 nodes and 90 nodes is within 4% for given range of non-dimensional fin length. 3) The relative error of the calculated heat loss (the number of node=90) as compared to the expected exact heat loss is less then 1.5% for Bi=0.1 and L=10 while that is over 13% for Bi=1.0 and L=10.

• Composites Research
• /
• v.22 no.5
• /
• pp.24-29
• /
• 2009

In this paper, a post-tensioned reinforced concrete slab was analyzed by the specially orthotropic laminates theory. Both the geometrical and material property of the cross section of the slab was considered symmetrically with respect to the neutral surface so that the bending extension coupling stiffness, $B_{ij}=0$, and $D_{16}=D_{26}=0$. Reinforced concrete slab behave as specially orthotropic plates. In general, the analytical solution for such complex systems is very difficult to obtain. Thus, finite difference method was used for analysis of the problem. In this paper, the finite difference method and the beam theory were used for analysis. The result of beam analysis was modified to obtain the solution of the plate analysis.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.1920-1929
• /
• 1994

The stochastic finite element method(SFEM) based structural optimal design is presented. Random system response including uncertainties for the design variable is calculated with first order perturbation method. A method for calculating the sensitivity coefficients is developed using the equilibrium equation and first-order perturbed equation. Numerical results are presented for a truss, frame and plate structures with displacement and stress constraints. The sensitivity calculation proposed here is compared with finite difference method. A nonlinear programming technique is used to solve the problem. The procedure is easily incorporated with existing deterministic structural optimization.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1465-1477
• /
• 1993

Powder Injection Molding(PM) is an advanced and complicated technology for manufacturing ceramic or metal products making use of a conventional injection molding process, which is generally used for plastic products. Among many technologies involved in the successful PIM, injection molding process is one of the key steps to form a desired shape out of powder/binder mixtures. Thus, it is of great importance to have a numerical tool to predict the powder injection molding filling process. In this regard, a finite element analysis system has been developed for numerical simulations of filling process of powder injection molding. Powder/polymer mixtures during the filling pro cess of injection molding can be rheologically characterized as Non-Newtonian fluids with a so called yield phenomena and have a peculiar feature of apparent slip phenomena on the wall boundaries surrounding mold cavity. Therefore, in the present study, a physical modeling of the filling process of powder/polymer mixtures was developed to take into account both the yield stress and slip phenomena and a finite element formulation was developed accordingly. The numerical analysis scheme for filling simulation is accomplished by combining a finite element method with control volume technique to simulate the movement of flow front and a finite difference method to calculate the temperature distribution. The present study presents the modeling, numerical scheme and some numerical analysis results showing the effect of the yield stress and slip phenomena.