• Korea-Australia Rheology Journal
• /
• 제14권1호
• /
• pp.33-45
• /
• 2002

Recent results obtained for the port-pom model and the constitutive equations with time-strain separability are examined. The time-strain separability in viscoelastic systems Is not a rule derived from fundamental principles but merely a hypothesis based on experimental phenomena, stress relaxation at long times. The violation of separability in the short-time response just after a step strain is also well understood (Archer, 1999). In constitutive modeling, time-strain separability has been extensively employed because of its theoretical simplicity and practical convenience. Here we present a simple analysis that verifies this hypothesis inevitably incurs mathematical inconsistency in the viewpoint of stability. Employing an asymptotic analysis, we show that both differential and integral constitutive equations based on time-strain separability are either Hadamard-type unstable or dissipative unstable. The conclusion drawn in this study is shown to be applicable to the Doi-Edwards model (with independent alignment approximation). Hence, the Hadamardtype instability of the Doi-Edwards model results from the time-strain separability in its formulation, and its remedy may lie in the transition mechanism from Rouse to reptational relaxation supposed by Doi and Edwards. Recently in order to describe the complex rheological behavior of polymer melts with long side branches like low density polyethylene, new constitutive equations called the port-pom equations have been derived in the integral/differential form and also in the simplifled differential type by McLeish and carson on the basis of the reptation dynamics with simplifled branch structure taken into account. In this study mathematical stability analysis under short and high frequency wave disturbances has been performed for these constitutive equations. It is proved that the differential model is globally Hadamard stable, and the integral model seems stable, as long as the orientation tensor remains positive definite or the smooth strain history in the flow is previously given. However cautious attention has to be paid when one employs the simplified version of the constitutive equations without arm withdrawal, since neglecting the arm withdrawal immediately yields Hadamard instability. In the flow regime of creep shear flow where the applied constant shear stress exceeds the maximum achievable value in the steady flow curves, the constitutive equations exhibit severe instability that the solution possesses strong discontinuity at the moment of change of chain dynamics mechanisms.

• Structural Engineering and Mechanics
• /
• 제75권5호
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.

• Structural Engineering and Mechanics
• /
• 제65권1호
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Earthquakes and Structures
• /
• 제7권2호
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Structural Engineering and Mechanics
• /
• 제9권4호
• /
• pp.397-416
• /
• 2000

The present study deals with the nonlinear dynamics and stability of autonomous dissipative either imperfect potential (limit point) systems or perfect (bifurcational) non-potential ones. Through a fully nonlinear dynamic analysis, performed on two simple 2-DOF models corresponding to the classes of systems mentioned above, and with the aid of basic definitions of the theory of nonlinear dynamical systems, new important phenomena are revealed. For the first class of systems a third possibility of postbuckling dynamic response is offered, associated with a point attractor on the prebuckling primary path, while for the second one the new findings are chaos-like (most likely chaotic) motions, consecutive regions of point and periodic attractors, series of global bifurcations and point attractor response of always existing complementary equilibrium configurations, regardless of the value of the nonconservativeness parameter.

• International Journal of Control, Automation, and Systems
• /
• 제2권1호
• /
• pp.55-67
• /
• 2004

In this paper, an active vibration control of a translating steel strip in a zinc galvanizing line is investigated. The control objectives in the galvanizing line are to improve the uniformity of the zinc deposit on the strip surfaces and to reduce the zinc consumption. The translating steel strip is modeled as a moving belt equation by using Hamilton’s principle for systems with moving mass. The total mechanical energy of the strip is considered to be a Lyapunov function candidate. A nonlinear boundary control law that assures the exponential stability of the closed loop system is derived. The existence of a closed-loop solution is shown by proving that the closed-loop dynamics is dissipative. Simulation results are provided.

• 한국전산유체공학회지
• /
• 제9권3호
• /
• pp.26-38
• /
• 2004

A study has been made for the investigation of the robustness and convergence of various implicit operators of the LU-SGS scheme using linear stability analysis. It is shown that the behavior of the implicit operator is not determined by its own characteristics, but is determined relatively depending on the dissipative property of the explicit operator. It is also shown that, as the dissipation level of the implicit operator increases, the robustness of the scheme increases, but the convergence rate can be deteriorated due to the excessive dissipation. The numerical results demonstrate that the dissipation level of the impliict operator needs to be higher than that of the explicit operator for computing stiff problems.

• Structural Engineering and Mechanics
• /
• 제55권4호
• /
• pp.815-837
• /
• 2015

A new family of structure-dependent integration methods is developed to enhance with desired numerical damping. This family method preserves the most important advantage of the structure-dependent integration method, which can integrate unconditional stability and explicit formulation together, and thus it is very computationally efficient. In addition, its numerical damping can be continuously controlled with a parameter. Consequently, it is best suited to solving an inertia-type problem, where the unimportant high frequency responses can be suppressed or even eliminated by the favorable numerical damping while the low frequency modes can be very accurately integrated.

• Bulletin of the Korean Chemical Society
• /
• 제20권1호
• /
• pp.35-41
• /
• 1999

Stability characteristics of hyperbolic reaction-diffusion equations with a reversible Brusselator model are investigated as an extension of the previous work. Intensive stability analysis is performed for three important parameters, Nrd, β and Dx, where Nrd is the reaction-diffusion number which is a measure of hyperbolicity, β is a measure of reversibility of autocatalytic reaction and Dx is a diffusion coefficient of intermediate X. Especially, the dependence on Nrd of stability exhibits some interesting features, such as hyperbolicity in the small Nrd region and parabolicity in the large Nrd region. The hyperbolic reaction-diffusion equations are solved numerically by a spectral method which is modified and adjusted to hyperbolic partial differential equations. The numerical method gives good accuracy and efficiency even in a stiff region in the case of small Nrd, and it can be extended to a two-dimensional system. Four types of solution, spatially homogeneous, spatially oscillatory, spatio-temporally oscillatory and chaotic can be obtained. Entropy productions for reaction are also calculated to get some crucial information related to the bifurcation of the system. At the bifurcation point, entropy production changes discontinuously and it shows that different structures of the system have different modes in the dissipative process required to maintain the structure of the system. But it appears that magnitude of entropy production in each structure give no important information related for states of system itself.