• Structural Engineering and Mechanics
• /
• v.49 no.4
• /
• pp.427-448
• /
• 2014

At present, the time of easy oil and gas is over. Now, the largest part of fossil fuels is concentrated in the deepest levels of tectonic structures and in the sea shelves. One of the most cumbersome operations of their extraction is the bore-hole drilling. In connection with austere tectonic and climate conditions, their drivage every so often is associated with great and diversified technological difficulties causing emergencies on frequent occasions. As a rule, they are linked with drill string accidents. A key role in prediction of these situations should play methods of theoretical modelling. For this reason, there is a growing need for development and implementation of new numerical methods for computer simulation of critical and post-critical behavior of drill strings (DSs). In this paper, the processes of non-linear deforming of a DS in cylindrical cavity of a deep bore-hole are considered. On the basis of the theory of curvilinear flexible rods, non-linear constitutive differential equations are deduced. The effects of the longitudinal non-uniform preloading, action of torque and interaction between the DS and the bore-hole surface are taken into account. Owing to the use of curvilinear coordinates in the constraining cylindrical surface and a specially chosen concomitant reference frame, it became possible to separate the desired variables and to reduce the total order of the equation system. To solve it, the method of continuation the solution by parameter and the transfer matrix technique are applied. As a result of the completed numerical analysis, the critical states of the DS loading in the cylindrical channels of inclined bore-holes are found. It is shown that the modes of the post-critical deforming of the DS are associated with its irregular spiral curving prevailing in the zone of bottom-hole-assembly. The possibility of invariant state generation during post-critical deforming is established, condition of its bifurcation is formulated. It is shown that infinite variety of loads can correspond to one geometrical configuration of the DS. They differ each from other by contact force functions.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.99-102
• /
• 2006

Primitive chain network model for block copolymers is used here to simulate molecular dynamics in the entangled state with acceptable computational cost. It was found that i) the hooking procedure rearranging the topology of the entangled network is critical for the equilibrium structure of the system, and ii) simulations accounting for the different chemistry, i.e., with a biased hooking probability based on interaction parameter ${\chi}$ for selection of the hooked partner, generates a reasonable phase diagram.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.2
• /
• pp.137-142
• /
• 2011

In this paper we formulate a predator-prey system in two patches in which the per capita migration rate of each species is influenced only by its own density, i.e. there is no response to the density of the other one. Numerical studies show that at a critical value of the bifurcation parameter the system undergoes a Turing bifurcation, i. e. the stable constant steady state loses its stability and spatially non-constant stationary solutions, a pattern emerge.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.1
• /
• pp.36-43
• /
• 2015

Since the crack generation and its propagation caused by welding defects is one of the main hull damage patterns, the simulation of crack propagation process has an important significance for ship safety. Based on interface element method, a finite element model to simulate crack propagation is studied in the paper. A Lennard-Jones type potential function is employed to define potential energy of the interface element. Tensile tests of steel flat plates with initial central crack are carried out. Surface energy density and spring critical stress that are suitable for the simulation of crack propagation are determined by comparing numerical calculation and tests results. Based on a large number of simulation results, the curve of simulation correction parameter plotted against the crack length is calculated.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.241-260
• /
• 2017

An analysis is made to study the solute transport in a Casson fluid flow through an annulus in presence of oscillatory flow field and determine how this flow influence the solute dispersion along the annular region. Axial dispersion coefficient and the mean concentration expressions are calculated using the generalized dispersion model. Dispersion coefficient in oscillatory flow is found to be a function of frequency parameter, Schmidt number, and the pressure fluctuation component besides its dependency on yield stress of the fluid, annular gap and time in the case of steady flow. Due to the oscillatory nature of the flow, the dispersion coefficient changes cyclically and the amplitude and magnitude of the dispersion increases initially with time and reaches a non - transient state after a certain critical time. This critical value varies with frequency parameter and independent of the other parameters. It is found that the presence of inner cylinder and increase in the size of the inner cylinder inhibits the dispersion process. This model may be used in understanding the dispersion phenomenon in cardiovascular flows and in particular in catheterized arteries.

• Structural Engineering and Mechanics
• /
• v.14 no.4
• /
• pp.491-504
• /
• 2002

A state-space method is proposed to analyze the aerodynamically coupled flutter problems of long-span bridges based on the modal coordinates of structure. The theory about complex modes is applied in this paper. The general governing equation of the system is converted into a complex standard characteristic equation in a state space format, which contains only two variables. The proposed method is a single-parameter searching method about reduced velocity, and it need not choose the participating modes beforehand and has no requirement for the form of structure damping matrix. The information about variations of system characteristics with reduced velocity and wind velocity can be provided. The method is able to find automatically the lowest critical flutter velocity and give relative amplitudes, phases and energy ratios of the participating modes in the flutter motion. Moreover, the flutter analysis of Jiangyin Yangtse suspension bridge with 1385 m main span is performed. The proposed method has proved reliable in its methodology and efficient in its use.

• Journal of Adhesion and Interface
• /
• v.2 no.3
• /
• pp.9-15
• /
• 2001

To estimate the effect of the surface state of the EVA foam on an adhesion, the surface was treated with solvents having different values of a solubility parameter and a surface tension. The morphology and the contact angle were measured by SEM and contact angle tester, respectively. The value of the critical surface tension(${\gamma}_{c,0}$) of the EVA foam calculated by Zismann plot was 27.08 dyne/cm. The surface state and the ${\gamma}_c$ of the EVA foam were changed by solvent treatment. The swelling ratio of EVA foam was influenced by solubility parameter of a solvent, on the other hand the degree of change in the surface state was dependent on a surface tension of a solvent. When a surface tension of solvent. When a surface tension of solvent is lower than ${\gamma}_{c,0}$ of the EVA foam, the wettability of a solvent was good and the surface state was greatly changed. The adhesion strength of EVA foam was greatly improved to above 600% as compared to that of the untreated one. As the ${\gamma}_c$ of EVA was closed to the surface tension of the primer, the good adhesion strength was obtained.

• Structural Engineering and Mechanics
• /
• v.85 no.5
• /
• pp.679-691
• /
• 2023

Many codes allow designers to use the bending moment diagram computed by elastic analysis and modify it by a certain amount of moment redistribution (MR) to account for plastic behaviour of continuous beams. However, several researchers indicated that the MR at the ultimate limit state (𝛽_u) for some beams deviate significantly from the specified values of various codes. This paper examines the applicability of the provisions on 𝛽_u in ACI 318-19 and Eurocode 2 through numerical investigations and comprehensively explores the influencing factors. The results show that some parameters not considered in those codes influence 𝛽_u to a certain extent, where the ratio of tensile reinforcement ratio at intermediate support to tensile reinforcement ratio at midspan (𝜌_s1/𝜌_s2) and load type are crucial parameters to consider. The specific combination of these two parameters may make the codes overestimate or significantly underestimate the 𝛽_u. On the other hand, the yield state of both critical sections is found to have an important influence on the influence degree of each parameter on 𝛽_u. The yield conditions are investigated, and an empirical judgment equation is proposed. In addition, the influence laws of the critical parameters on 𝛽_u have been further proved by theoretical derivation. Finally, due to 𝜀_t is found to have a better linear correlation with 𝛽_u than x_u/d, equations as a function of 𝜀_t for predicting the 𝛽_u of continuous beams under the two loads are proposed, respectively.

• Transactions on Electrical and Electronic Materials
• /
• v.17 no.6
• /
• pp.311-316
• /
• 2016

Batteries are critical components of electric vehicles and energy storage systems. The connection of a battery to the power grid for charge and discharge greatly affects energy storage. Therefore, an accurate and easy-to-observe battery model should be established to achieve systematic design, simulation, and SOC (state of charge) estimations. In this review, several equivalent circuit models of representative significance are explained, and their respective advantages and disadvantages are compared to determine and outline their reasonable applications to Li-ion batteries. Numerous commonly used model parameter identification principles are summarized as well, and basic model verification methods are briefly introduced for the convenient use of such models.

• Proceedings of the KIPE Conference
• /
• 2016.11a
• /
• pp.75-76
• /
• 2016

The relationship of widely used Open circuit Voltage (OCV) versus State of Charge (SOC) is critical for any reliable SOC estimation technique. However, the hysteresis existing in all type of battery which has been come to the market leads this relationship to a complicated one, especially in Lithium Iron Phosphate (LiFePO4) battery. An accurate model for hysteresis phenomenon is essential for a reliable SOC identification. This paper aims to investigate and propose a method for hysteresis modeling. The SOC estimation is done by using Extended Kalman Filter (EKF), the parameter of the battery is modeled by Auto Regressive Exogenous (ARX) and estimated by using Recursive Least Square (RLS) filter to tract each element of the parameter of the model.