\documentclass[11pt,twoside]{article}\makeatletter \IfFileExists{xcolor.sty}% {\RequirePackage{xcolor}}% {\RequirePackage{color}} \usepackage{colortbl} \usepackage{wrapfig} \usepackage{ifxetex} \ifxetex \usepackage{fontspec} \usepackage{xunicode} \catcode`⃥=\active \def⃥{\textbackslash} \catcode`❴=\active \def❴{\{} \catcode`❵=\active \def❵{\}} \def\textJapanese{\fontspec{Noto Sans CJK JP}} \def\textChinese{\fontspec{Noto Sans CJK SC}} \def\textKorean{\fontspec{Noto Sans CJK KR}} \setmonofont{DejaVu Sans Mono} \else \IfFileExists{utf8x.def}% {\usepackage[utf8x]{inputenc} \PrerenderUnicode{–} }% {\usepackage[utf8]{inputenc}} \usepackage[english]{babel} \usepackage[T1]{fontenc} \usepackage{float} \usepackage[]{ucs} \uc@dclc{8421}{default}{\textbackslash } \uc@dclc{10100}{default}{\{} \uc@dclc{10101}{default}{\}} \uc@dclc{8491}{default}{\AA{}} \uc@dclc{8239}{default}{\,} \uc@dclc{20154}{default}{ } \uc@dclc{10148}{default}{>} \def\textschwa{\rotatebox{-90}{e}} \def\textJapanese{} \def\textChinese{} \IfFileExists{tipa.sty}{\usepackage{tipa}}{} \fi \def\exampleFont{\ttfamily\small} \DeclareTextSymbol{\textpi}{OML}{25} \usepackage{relsize} \RequirePackage{array} \def\@testpach{\@chclass \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else \ifnum \@lastchclass=7 5 \else \ifnum \@lastchclass=8 \tw@ \else \ifnum \@lastchclass=9 \thr@@ \else \z@ \ifnum \@lastchclass = 10 \else \edef\@nextchar{\expandafter\string\@nextchar}% \@chnum \if \@nextchar c\z@ \else \if \@nextchar l\@ne \else \if \@nextchar r\tw@ \else \z@ \@chclass \if\@nextchar |\@ne \else \if \@nextchar !6 \else \if \@nextchar @7 \else \if \@nextchar (8 \else \if \@nextchar )9 \else 10 \@chnum \if \@nextchar m\thr@@\else \if \@nextchar p4 \else \if \@nextchar b5 \else \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi} \gdef\arraybackslash{\let\\=\@arraycr} \def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}} \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \def\abbr{} \def\corr{} \def\expan{} \def\gap{} \def\orig{} \def\reg{} \def\ref{} \def\sic{} \def\persName{}\def\name{} \def\placeName{} \def\orgName{} \def\textcal#1{{\fontspec{Lucida Calligraphy}#1}} \def\textgothic#1{{\fontspec{Lucida Blackletter}#1}} \def\textlarge#1{{\large #1}} \def\textoverbar#1{\ensuremath{\overline{#1}}} \def\textquoted#1{‘#1’} \def\textsmall#1{{\small #1}} \def\textsubscript#1{\@textsubscript{\selectfont#1}} \def\textxi{\ensuremath{\xi}} \def\titlem{\itshape} \newenvironment{biblfree}{}{\ifvmode\par\fi } \newenvironment{bibl}{}{} \newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par } \newenvironment{citbibl}{}{\ifvmode\par\fi } \newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi } \newenvironment{docDate}{}{\ifvmode\par\fi } \newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi } \newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par } \newenvironment{msHead}{\vskip 6pt}{\par} \newenvironment{msItem}{\vskip 6pt}{\par} \newenvironment{rubric}{}{} \newenvironment{titlePart}{}{\par } \newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}} \newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}} \newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}} \newcolumntype{P}[1]{){\arraybackslash}p{#1}} \newcolumntype{B}[1]{){\arraybackslash}b{#1}} \newcolumntype{M}[1]{){\arraybackslash}m{#1}} \definecolor{label}{gray}{0.75} \def\unusedattribute#1{\sout{\textcolor{label}{#1}}} \DeclareRobustCommand*{\xref}{\hyper@normalise\xref@} \def\xref@#1#2{\hyper@linkurl{#2}{#1}} \begingroup \catcode`\_=\active \gdef_#1{\ensuremath{\sb{\mathrm{#1}}}} \endgroup \mathcode`\_=\string"8000 \catcode`\_=12\relax \usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry} \usepackage{framed} \definecolor{shadecolor}{gray}{0.95} \usepackage{longtable} \usepackage[normalem]{ulem} \usepackage{fancyvrb} \usepackage{fancyhdr} \usepackage{graphicx} \usepackage{marginnote} \renewcommand{\@cite}[1]{#1} \renewcommand*{\marginfont}{\itshape\footnotesize} \def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif} \pagestyle{fancy} \usepackage[pdftitle={Analysis of Effect of MOV on Chaotic Ferroresonant Oscillations in unloaded Transformers by Chaos Theory}, pdfauthor={}]{hyperref} \hyperbaseurl{} \paperwidth210mm \paperheight297mm \def\@pnumwidth{1.55em} \def\@tocrmarg {2.55em} \def\@dotsep{4.5} \setcounter{tocdepth}{3} \clubpenalty=8000 \emergencystretch 3em \hbadness=4000 \hyphenpenalty=400 \pretolerance=750 \tolerance=2000 \vbadness=4000 \widowpenalty=10000 \renewcommand\section{\@startsection {section}{1}{\z@}% {-1.75ex \@plus -0.5ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large\bfseries}} \renewcommand\subsection{\@startsection{subsection}{2}{\z@}% {-1.75ex\@plus -0.5ex \@minus- .2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large}} \renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}% {-1.5ex\@plus -0.35ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\large}} \renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}% {-1ex \@plus-0.35ex \@minus -0.2ex}% {0.5ex \@plus .2ex}% {\reset@font\normalsize}} \renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}% {1.5ex \@plus1ex \@minus .2ex}% {-1em}% {\reset@font\normalsize\bfseries}} \def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt} \@tempdima 1.5em \begingroup \parindent \z@ \rightskip \@pnumwidth \parfillskip -\@pnumwidth \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par \endgroup} \def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}} \def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}} \def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}} \def\l@subparagraph{\@dottedtocline{5}{10em}{5em}} \@ifundefined{c@section}{\newcounter{section}}{} \@ifundefined{c@chapter}{\newcounter{chapter}}{} \newif\if@mainmatter \@mainmattertrue \def\chaptername{Chapter} \def\frontmatter{% \pagenumbering{roman} \def\thechapter{\@roman\c@chapter} \def\theHchapter{\roman{chapter}} \def\thesection{\@roman\c@section} \def\theHsection{\roman{section}} \def\@chapapp{}% } \def\mainmatter{% \cleardoublepage \def\thechapter{\@arabic\c@chapter} \setcounter{chapter}{0} \setcounter{section}{0} \pagenumbering{arabic} \setcounter{secnumdepth}{6} \def\@chapapp{\chaptername}% \def\theHchapter{\arabic{chapter}} \def\thesection{\@arabic\c@section} \def\theHsection{\arabic{section}} } \def\backmatter{% \cleardoublepage \setcounter{chapter}{0} \setcounter{section}{0} \setcounter{secnumdepth}{2} \def\@chapapp{\appendixname}% \def\thechapter{\@Alph\c@chapter} \def\theHchapter{\Alph{chapter}} \appendix } \newenvironment{bibitemlist}[1]{% \list{\@biblabel{\@arabic\c@enumiv}}% {\settowidth\labelwidth{\@biblabel{#1}}% \leftmargin\labelwidth \advance\leftmargin\labelsep \@openbib@code \usecounter{enumiv}% \let\p@enumiv\@empty \renewcommand\theenumiv{\@arabic\c@enumiv}% }% \sloppy \clubpenalty4000 \@clubpenalty \clubpenalty \widowpenalty4000% \sfcode`\.\@m}% {\def\@noitemerr {\@latex@warning{Empty `bibitemlist' environment}}% \endlist} \def\tableofcontents{\section*{\contentsname}\@starttoc{toc}} \parskip0pt \parindent1em \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \newenvironment{reflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemsep}{0pt}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\itshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{sansreflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\itemsep}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\upshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{specHead}[2]% {\vspace{20pt}\hrule\vspace{10pt}% \phantomsection\label{#1}\markright{#2}% \pdfbookmark[2]{#2}{#1}% \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}% }{} \def\TheFullDate{2011-12-00 2011-12-24 (revised: 0 December 2011 24 December 2011)} \def\TheID{\makeatother } \def\TheDate{2011-12-00 2011-12-24} \title{Analysis of Effect of MOV on Chaotic Ferroresonant Oscillations in unloaded Transformers by Chaos Theory} \author{}\makeatletter \makeatletter \newcommand*{\cleartoleftpage}{% \clearpage \if@twoside \ifodd\c@page \hbox{}\newpage \if@twocolumn \hbox{}\newpage \fi \fi \fi } \makeatother \makeatletter \thispagestyle{empty} \markright{\@title}\markboth{\@title}{\@author} \renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@ \belowdisplayskip \abovedisplayskip \abovedisplayshortskip \z@ plus2\p@ \belowdisplayshortskip 4\p@ plus2\p@ minus2\p@ \def\@listi{\leftmargin\leftmargini \topsep 2\p@ plus1\p@ minus1\p@ \parsep 2\p@ plus\p@ minus\p@ \itemsep 1pt} } \makeatother \fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm} \fancyhf{} \setlength{\headheight}{14pt} \fancyhead[LE]{\bfseries\leftmark} \fancyhead[RO]{\bfseries\rightmark} \fancyfoot[RO]{} \fancyfoot[CO]{\thepage} \fancyfoot[LO]{\TheID} \fancyfoot[LE]{} \fancyfoot[CE]{\thepage} \fancyfoot[RE]{\TheID} \hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true} \fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}} \date{} \usepackage{authblk} \providecommand{\keywords}[1] { \footnotesize \textbf{\textit{Index terms---}} #1 } \usepackage{graphicx,xcolor} \definecolor{GJBlue}{HTML}{273B81} \definecolor{GJLightBlue}{HTML}{0A9DD9} \definecolor{GJMediumGrey}{HTML}{6D6E70} \definecolor{GJLightGrey}{HTML}{929497} \renewenvironment{abstract}{% \setlength{\parindent}{0pt}\raggedright \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt \textcolor{GJBlue}{\large\bfseries\abstractname\space} }{% \vskip8pt \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt } \usepackage[absolute,overlay]{textpos} \makeatother \usepackage{lineno} \linenumbers \begin{document} \author[1]{hamid reza abbasi} \affil[1]{ iran university of science and technology} \renewcommand\Authands{ and } \date{\small \em Received: 22 November 2011 Accepted: 12 December 2011 Published: 24 December 2011} \maketitle \begin{abstract} In this paper the effect of a parallel metal oxide surge arrester on the ferroresonance oscillations of transformers is studied. It is proved that ferroresonance phenomenon in transforemer can be classified in chaotic dynamics systems. In this contribution chaos occurs in system from a sequence of period doubling bifurcation (PDB). Analysis of dynamics of ferroresonant circuit is carried out using bifurcation theory. It is expected that the arresters generally cause ferroresonance drop out. Simulation has been done on a three phase power transformer with one open phase. Effect of varying input voltage is studied. The simulation results reveal that connecting the arrester to transformers poles, exhibits a great mitigating effect on ferroresonant over voltages. Phase plane along with bifurcation diagrams are also presented. Significant effect on the onset of chaos, the range of parameter values that may lead to chaos and magnitude of ferroresonant voltages is obtained, shown and tabulated. \end{abstract} \keywords{component; Power Transformer, Phase Plane Diagram, Bifurcation Diagram, Chaotic Ferroresonance, Chaos Theory, Nonlinear Core loss Model, Metal Oxide A} \begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) \textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark} \end{textblock*} \begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) \textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}} \end{textblock*} \let\tabcellsep& \section[{Analysis of Effect of MOV on Chaotic Ferroresonant Oscillations in unloaded}]{Analysis of Effect of MOV on Chaotic Ferroresonant Oscillations in unloaded}\par Transformers by Chaos Theory erroresonance is a complex nonlinear electrical phenomenon that can cause dielectric \& thermal problems to components power system. Electrical systems exhibiting ferroresonant behaviour are categorized as nonlinear dynamical systems. Therefore conventional linear solutions cannot be applied to study ferroresonance. The prediction of ferroresonance is achieved by detailed modeling using a digital computer transient analysis program \hyperref[b0]{[1]}. Ferroresonance should not be confused with linear resonance that occurs when inductive and capacitive reactance of circuit is equal. In linear resonance the current and voltage are linearly related and are frequency dependent. In the case of ferroresonance it is characterized by a sudden jump of voltage or current from one stable operating state to another one. The relationship between voltage and current is depends not only on frequency but also on other factors such as system voltage magnitude, initial magnetic flux condition of transformer iron core, total loss in the ferroresonant circuit and moment of switching \hyperref[b1]{[2]}.\par Ferroresonance may be initiated by contingency switching operation, routine switching, or load shedding involving a high voltage transmission line. It can result in Unpredictable over voltages and high currents. The prerequisite for ferroresonance is a circuit containing iron core inductance and a capacitance. Such a circuit is characterized by simultaneous existence of several steady-state solutions for a given set of circuit parameters. The abrupt transition or jump from one steady state to another is triggered by a disturbance, switching action or a gradual change in values of a parameter. Typical cases of ferroresonance are reported in \hyperref[b0]{[1]}, \hyperref[b1]{[2]}, \hyperref[b2]{[3]} and \hyperref[b3]{[4]}. Although analyzing methods such as harmonics balance method can be use for analyzing nonlinear differential equations, but solving these equations lead to a set of complex algebraic equations \hyperref[b2]{[3]}. Thus, scientists should use other methods to solve nonlinear dynamic equations. One of these methods is bifurcation theory which some articles use from this method \hyperref[b4]{[5,}\hyperref[b5]{6,}\hyperref[b6]{7]}. Bifurcation theory enables us to describe and analyze qualitative properties of solutions (fixed points) when system parameters change. Studying ferroresonance by bifurcation theory has been carried out \hyperref[b7]{[8,} {\ref 9,}\hyperref[b8]{10,}\hyperref[b9]{11]}. But there are some problems in these articles. For example method used in \hyperref[b13]{[15]} is valid only in limited cases while creating a bifurcation diagram by a continuation method can be more systematic and save computational effort \hyperref[b2]{[3]}. The samples of ferroresonace in power system have been described in \hyperref[b10]{[12,}\hyperref[b3]{4,}\hyperref[b11]{13]}. Analyzing chaotic ferroresonant behavior in power transformer and dependence of this behavior on system parameters such as amplitude of voltage source, capacitance and resistance of system, core loss, initial conditions and e and effect of neutral resistance in damping ferroresosnant oscillations and change in system behavior from chaotic to multi frequency in \hyperref[b2]{[3,}\hyperref[b6]{7,}\hyperref[b12]{14,}\hyperref[b13]{15,}\hyperref[b14]{16,}\hyperref[b15]{17]} have been studied. Evaluation of route to chaos in transformer with modeling and solving equations in conditions that defined model for core loss is considered linear and effect of complexity of circuits breaker models in transmission and distribution lines with considering effect of damping in system and elimination of caused harmonics in \hyperref[b8]{[10,}\hyperref[b16]{18]} have been studied. Theory of nonlinear dynamics has been found to provide deeper insight into the phenomenon. \hyperref[b17]{[19]}, \hyperref[b9]{[11]}, \hyperref[b18]{[20]} and \hyperref[b19]{[21]} are among the early investigations in applying theory of bifurcation and chaosto ferroresonance. The susceptibility of a ferroresonant circuit to a quasi-periodic and frequency locked oscillations are presented in \hyperref[b20]{[22]}. The effect of initial conditions is also investigated. The effect of transformer modeling on the predicted ferroresonance oscillations has been studied in \hyperref[b21]{[23]}. Using a linear model, authors of \hyperref[b22]{[24]} have indicated the effect of core loss in damping ferroresonance oscillations. The importance of treating core loss as a nonlinear function of voltage is highlighted in \hyperref[b20]{[22]}. An algorithm for calculating core loss from no-load characteristics is given in \hyperref[b23]{[25]}. Evaluation of chaos in transformer, effect of resistance of key on the chaotic behavior transformer and subharmonics that produced with ferroresonance in this type transformer and quantification of the chaotic behavior of ferroresonant transformer circuits are studied in \hyperref[b18]{[20]}, \hyperref[b23]{[25]} and \hyperref[b24]{[26]}.\par Transformer is assumed to be connected to the Power System while one of the three switches are open and only two phases of it are energized, which produces induced voltage in the open phase. This voltage, back feeds the distribution line. Ferroresonance will occur if the distribution line is highly capacitive. System involves the nonlinear magnetizing reactance of the transformer's open phase and resulted shunt and series capacitance of the distribution line.\par Base system model is adopted from \hyperref[b2]{[3]} with the MOV arrester connected across the transformer winding which is showed in Fig. {\ref 1} Linear approximation of the peak current of the magnetization reactance can be presented by Eq. ( {\ref 1}):\par However, for very high currents, the iron core might be saturated where the flux-current characteristic becomes highly nonlinear. The li characteristic of the transformer can be demonstrated by the polynomial in Eq. ( {\ref 2})\par Arrester can be expressed by the Eq. ( {\ref 3}):\par (3) V represents resistive voltage drop, I represents arrester current and K is constant and is nonlinearity constant. The differential equation for the circuit in Fig. {\ref 1} can be derived as follows:\par Presenting in the form of state space equations, and p will be state variables as follows:\par Multiple Scales Method By using the multiple scales method one obtains a first order approximation for the solution of Eq. ( {\ref 7}) as:\par The parameters ?, a and k are of ?. Further, the frequency independent of system is such that Where ? is named external detuning. By using the multiple scales method, we seek of first order uniform expansion of Eq. ( {\ref 7}) in the form: Where T 0 = t and T 1 = ?T 0 . In term of T 1 the time derivative becomes: Substituting Eq. ( {\ref 10}) and Eq. ( {\ref 11}) into Eq. ( {\ref 36}) and equating coefficient of like power of ?, we obtain:\par The solution of Eq. ( {\ref 12}) can be expressed as: Substituting Eq. ( {\ref 14}) in Eq. ( {\ref 13}):\par Where cc is complex conjugate of preceding terms and the prime indicates the derivation with respect to T 1 . Using Eq. ( {\ref 7}) in eliminating the lead to secular terms in from Eq. ( {\ref 12}), we obtain:\par If A is defined in the polar form , where ?, ? are functions of T with separating real and imaginary part in Eq. ( {\ref 13}):\par 1 From Eq. ( {\ref 17}), we obtain Eq. ( {\ref 18}) and Eq. ( {\ref 19}):\par With multiplyingsin? in Eq. ( {\ref 18}) and cos? in Eq. ( {\ref 19}) we have: With multiplying cos? in Eq. ( {\ref 18}) and sin? in Eq. ( {\ref 19}) we have:\par Setting ??= 0 and ??= 0 in Eq. ( {\ref 20}) and ( {\ref 21}) we find that their fixed points are given by: Squaring and adding Eq. ( {\ref 22}) and ( {\ref 23} Bifurcation theory describes and studies behavior of system with change in one or more parameters of system and discusses in case of stability and instability of fixed points in the values of system parameters. Suppose system is defined as Eq. (32):\par Where x is a state vector. In fact flux and voltage in terminal of transformer are state variables. ? is a parameter of system that can be value of series capacitance or amplitude of input voltage. for ? = ? c at which the vector field f losses its structural stability is called a bifurcation point and ?c the value of bifurcation. For analyzing and studying in bifurcation diagram we use of jacobian matrix, J = Df as the linearization of f at (x0, ?0 ) which points x0 are fixed points.\par Saddle node bifurcation When J is nonhyperbolic, i.e. J has a zero eigenvalue and no other eigen value with zero real part, saddle node bifurcation (SNB) occurs. SNB is caused with changes in the number of fixed points. Indeed, one stable fixed point and unstable fixed point cause SNB. Necessary and enough conditions for SNB are: Necessary conditions: \section[{Enough conditions:}]{Enough conditions:}\par If eigen values of jacobian matrix are considered as when real 0 jacobian matrix is hyperbolic and other wise nonhyperbolic.\par The stability of the fixed points depends on the eigenvalues of the jacobian matrix ( {\ref 22}), \hyperref[b21]{(23)}; that is, the eigevalue of:\par Determinant of [?I ? A yields eigenvalues:\par Where ? is eigenvalue. Substituting the polar form of A into Eq. ( {\ref 11}) and substituting result into Eq. ( {\ref 12}), we find that, to first approximation x1is given by: For nontrivial solutions, ? ? 0 and it follows from Eq. ( {\ref 28}) that: Substituting Eq. (29) into Eq. ( {\ref 27}), we find that to the first approximation, the free oscillations of Eq. ( {\ref 7}) are given by: Where ? is given by Eq. ( {\ref 28}), which has the normal form of a supercritical pitchfork bifurcation. Equation of eigenvalues introduce as the following equation:\par We obtain first order approximation of Eq. ( {\ref 8}) by multiple scale method and by using the chaos theory we discuss in case of stability. \section[{Global ( F ) 2011 November}]{Global ( F ) 2011 November}\par Pitch fork or transcritical bifurcation points appoint necessary conditions, too. For more detail see {\ref [30]}. Hopf bifurcation If J has a pair of complex conjugate on the imaginary axis and other eigenvalues lying off the imaginary axis, hopf bifurcation (HB) occurs. if periodic solutions are unstable, bifurcation is said to be subcritical and supercritical if stable. Thus, connects fixed points to periodic solutions. SNB and HB are stationary point. Periodic solutions that are caused by a HB can increase bifurcations and complexity of system behavior, themselves. Limit cycles which are caused by a HB can involve system into chaotic region and global bifurcation occurs. Suppose x be a small perturbation to the periodic solution to Eq. (30) and 1. We obtain: Thus, periodic in T, A (t) is periodic, too. Thus: be equal with A (t) Because is Stability of periodic solutions is determined by its characteristic.\par Now, we define the mondoromy matrix M to be ? (t). Eigenvalues of M are multipliers, denoted by Mi , i = 1, n. If all eigenvalues of M lie in the unit circle, we find out, For a periodic solution one of multipliers is equal to +1, with corresponding eigenvector tangential to the periodic orbit at x.\par Stability of a limit cycle is determined by its multipliers and depending on the way in which multipliers venter or leave the unit circle. \section[{C Cyclic fold bifurcation}]{C Cyclic fold bifurcation}\par If one of multipliers enters or leaves the unit circle along the positive real axis cyclic fold bifurcation (CFB) occurs. In Fig. {\ref 2 (a}) is example of this bifurcation. \section[{Fig. 2. Multiplier crossings of the unit circle}]{Fig. 2. Multiplier crossings of the unit circle}\par If one of the multipliers leaves unit circle along the negative real axis bifurcation is said to be a period doubling bifurcation (PDB) (b in Fig. {\ref 2}). This bifurcation causes new solutions with period 2T. If this behavior continues, causes solutions with infinitive period. These solutions are aperiodic which are called chaotic solutions. If one of the lyapunov exponents be positive for systems of ODEs, is representing of chaotic behavior in system. \section[{Torus bifurcation}]{Torus bifurcation} \section[{Period doubling}]{Period doubling}\par If a complex conjugate pair of multipliers with Re\{mi \} ? 0 leaves the unit circle, causes quasi-periodic solution.\par C indicates this behavior. Quasi periodic solutions have period that is equal to in commensurate of main period T. In phase plane diagram these solutions create figures in form of torus. \section[{Routes to chaos}]{Routes to chaos}\par Chaotic solutions are aperiodic and unstable solutions. These solutions depend on initial conditions. In this section we imply 4 routes to chaos: PDB Crises Intermittency Torus bifurcation Intermittency is a route to chaos. In this route oscillation in regular mode occasionally interrupted by turbulent burst of aperiodic oscillations at irregular intervals and chaos emerges in system. In case of torus bifurcation if stable periodic solution undergoes to a supercritical secondary hopf bifurcation with changes in parameter of system. This causes two quasi periodic solutions with two in commensurate frequencies.\par When parameter increases torus is destroyed and system becomes chaotic. Sudden changes in parameter of system cause crises and system becomes chaotic. When crises occur chaotic attractor enters unstable periodic solutions or saddle points.\par Crises have different types. Some of these types are:\par Fundamental matrix for Eq. ( {\ref 75}) is (T) , such as: If lyapunov exponent be positive, routes will repel other routes and other wise will attract other route. In case of stability of fixed points, when all lyapunov exponents are negative, these points are stable and in limit cycle lyapunov exponent is zero. Necessary and enough condition for chaotic behavior system are one or more positive lyapunov exponents. For more details, see \hyperref[b25]{[28]}.If\par exponent:\par Typical values for various system parameters considered for simulation are as given below \hyperref[b4]{[5]}: Time domain simulations were performed using the MATLAB programs which are similar to EMTP simulation \hyperref[b2]{[3]}. For cases including arrester, it can be seen that ferroresonant drop out will be occurred.Fig. {\ref 3} show the phase plane plot of system states without arrester for E=1 p.u.\par Figure {\ref .} 3 Phase plane diagram for E=1, q=11 without MOV When Vin increases system is entered into saturation section of magnetization curve and ferroresonance occurs. In Figs. \hyperref[b3]{4,}\hyperref[b4]{5,}\hyperref[b5]{6,}\hyperref[b6]{7} this phenomenon is shown. Behavior of system is single frequency but PDB has occurred. Magnetization curve in Fig. {\ref 4} and phase plane diagram in Fig. {\ref 5} and voltage and flux waveforms are shown in Figs. {\ref 6 and 7}. These figures are gained whenVin = 3.5, q = 7. Phase plane diagram shows this reality that behavior of system is a single frequency behavior. But voltage and flux waveforms show that behavior of system has an undesirable effect on system insulation and maybe damage it. With consideration to Fig. {\ref 13}, Fig. {\ref 14} and Fig. {\ref 15} MOV makes a mitigation in ferroresonance chaotic behavior in transformer that in down value of q the chaotic region are removed and the behavior will be periodic, for greater value of q for example for q=11 independent chaotic regions which can be created under MOV nominal voltage have survived so chaotic behavior has been eliminated. Figs. \hyperref[b14]{16,}\hyperref[b15]{17} show that chaotic region mitigates by applying MOA surge arrester. The system shows a greater tendency for chaos for saturation characteristics with lower knee Considering to Fig. {\ref 17} MOA makes a mitigation in ferroresonance chaotic behavior in the transformer that in down value of q the chaotic region are removed and the behavior will be periodic, for greater value of q such as q=11 independent chaotic regions which can be created under MOA nominal voltage have survived so chaotic behavior has been eliminated. Tendency to chaos exhibited by the system increases while q increases too.\par Chaotic ferroresonant oscillations of unloaded transformer nonlinear core loss model have been described. The presence of the arrester results in clamping the Ferroresonant over voltages in the studied system. The arrester successfully suppresses or eliminates the chaotic behaviour of proposed model. Consequently, the system shows less sensitivity to initial conditions in the presence of the arrester. It is seen from the bifurcation diagram that chaotic ferroresonant behavior depends on parameter q. MOV makes a mitigation in ferroresonance chaotic behavior in transformer that in down value of q the chaotic region are removed and the behavior will be periodic. System stability increased with decreasing q and chaotic regions are eliminiated. It is found when q=11 at vin=4 p.u. beahavior of system is chaotic while for q=7 in the same value of vin system is in subharmonic mode and its stability is more than case that q=11. It was shown that chaos occurs in transformer from a sequence of PDB. It was found that nonlinear magnetization curve has a great influence on bifurcation diagrams and domains of ferroresonance occurrence. nonlinear core loss model has been used in dynamics equations. It was found that the nonlinear core loss model causes the mitigation and delay in chaotic ferroresonant oscillations. Also presence of nonlinear term in core loss function causes PDBs become more regular. \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-2.png} \caption{\label{fig_0}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-3.png} \caption{\label{fig_1}F}\end{figure} \begin{figure}[htbp] \noindent\textbf{121}\includegraphics[]{image-4.png} \caption{\label{fig_2}Figure. 1 2 ) ( 1 )}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-5.png} \caption{\label{fig_4}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-6.png} \caption{\label{fig_5}}\end{figure} \begin{figure}[htbp] \noindent\textbf{44576778}\includegraphics[]{image-7.png} \caption{\label{fig_6}Fig. 4 : 4 Fig. 5 : 7 Fig. 6 : 7 Fig. 7 :Fig. 8 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{8}\includegraphics[]{image-8.png} \caption{\label{fig_7}Figure 8 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{111213}\includegraphics[]{image-9.png} \caption{\label{fig_8}Figure10:Figure 11 Figure 12 :Fig. 13 ,}\end{figure} \begin{figure}[htbp] \noindent\textbf{131415}\includegraphics[]{image-10.png} \caption{\label{fig_9}Figure 13 :Figure 14 :Figure 15 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{1617}\includegraphics[]{image-11.png} \caption{\label{fig_10}Fig. 16 :Fig. 17 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-12.png} \caption{\label{fig_11}Global}\end{figure} \begin{figure}[htbp] \noindent\textbf{(} \par \begin{longtable}{P{0.23640625\textwidth}P{0.61359375\textwidth}} \tabcellsep (B) behaviour of system with mov for e= 4, 5, 6\\ 2011\tabcellsep \\ November\tabcellsep Initial conditions:\\ 26\tabcellsep \\ Volume XI Issue VIII Version I\tabcellsep Table 1. (A) Behaviour of System Without Mov For E= 1, 2, 3\\ F )\tabcellsep \\ (\tabcellsep \\ Researches in Engineering\tabcellsep (B) Behaviour of System without Mov For E= 4, 5, 6\\ Journal of\tabcellsep \\ Global\tabcellsep Table 2. (A) Behaviour of System with Mov for E= 1, 2, 3\end{longtable} \par \caption{\label{tab_0}Table ( 1}\end{figure} \footnote{© 2011 Global Journals Inc. (US)} \footnote{© 2011 Global Journals Inc. (US) ©2011 Global Journals Inc. (US)} \backmatter \begin{bibitemlist}{1} \bibitem[Zare et al. (2007)]{b17}\label{b17} ‘Analysis of Ferroresonance Modes in Power Transformers using Preisach-Type Hysteretic Magnetizing Inductance’. Rezaei Zare , A Sanaye Pasand , M Mohseni , H Farhangi , S Iravani , R . \textit{IEEE Transactions on Power Delivery} April 2007. 22 (2) . \bibitem[Anbari et al. (2001)]{b10}\label{b10} ‘Analysis of nonlinear phenomena in MOV connected transformer’. K Anbari , R Ramanjam , T Keerthiga , K Kuppusamy . \textit{Gen., Trans and Dist} Nov. 2001. 148 p. . (IEE proc.) \bibitem[Kavasseri ()]{b23}\label{b23} \textit{Analysis of subharmonic oscillations in a ferroresonant circuit}, G Kavasseri . 2005. Elsevier journal. p. . Electrical Power and Energy Systems \bibitem[Rajesh and Kavasseri (2006)]{b6}\label{b6} \textit{Analysis of subharmonics oscillations in a ferroresonant circuit}, G Rajesh , Kavasseri . March. 2006. Elsevier journal. 28 p. . Electrical Power and Energy Systems \bibitem[Sakarung and Chatratana (2005)]{b5}\label{b5} ‘Application of PSCAD/EMTDC and Chaos Theory to Power System Ferroresonance Analysis’. P Sakarung , S Chatratana . \textit{International Conference on Power Systems Transients (IPST)}, June. 2005. p. . \bibitem[Ben-Tal et al. (2001)]{b15}\label{b15} ‘Banded Chaos in Power Systems’. A Ben-Tal , V Kirk , G Wake . \textit{IEEE Transactions on Power Delivery} Jan. 2001. 16 (1) . \bibitem[Escudero et al. (2007)]{b14}\label{b14} \textit{Characterization of ferroresonant modes in HV substation with CB grading capacitors}, M V Escudero , I Dudurych , M A Redfern . Jan. 2007. Elsevier Electric Power Systems Research. p. . \bibitem[Mozaffari et al. ()]{b2}\label{b2} ‘Effect of initial conditions on chaotic ferroresonance in power transformers’. S Mozaffari , M Sameti , A C Soudack . \textit{IEE Proceedings/Generation, Transmission and Distribution}, 1997. p. . \bibitem[Al-Anbarri et al. ()]{b3}\label{b3} \textit{Effect of iron core loss nonlinearity on chaotic ferroresonance in power transformers}, K Al-Anbarri , R Ramanujam , R Saravanaselvan , K Kuppusamy . 2003. p. . Electric Power Systems Research Elsevier journal \bibitem[Al-Anbarri et al. ()]{b13}\label{b13} \textit{Effect of iron core loss nonlinearity on chaotic ferroresonance in power transformers}, K Al-Anbarri , R Ramanujam , R Saravanaselvan , K Kuppusamy . 2003. p. . Electric Power Systems Research Elsevier journal \bibitem[Abbasi et al. (2010)]{b12}\label{b12} ‘Effect of Metal Oxide Arrester o Chaotic Behavior of Power Transformers’. A Abbasi , M Rostami , S H Fathi , H R Abbasi , H Abdollahi . \textit{Energy and Power Engineering(EPE journal)} Nov. 2010. 2 p. . \bibitem[Abbasi et al. ()]{b8}\label{b8} \textit{Elimination of Chaotic Ferroresonance in power transformers including Nonlinear Core Losses applying of Neutral Resistance}, A Abbasi , H Radmanesh , M Rostami , H Abbasi ; Eeeic09 . 2009. Poland. \bibitem[Abbasi et al. ()]{b9}\label{b9} \textit{Evaluation of Chaotic Ferroresonance in power transformers including Nonlinear Core Losses}, A Abbasi , M Rostami , H Radmanesh , H R Abbasi , ; Eeeic09 . 2009. Poland. \bibitem[Pattanapakdee and Banmongkol (2007)]{b11}\label{b11} ‘Failure of Riser Pole Arrester due to Station Service Transforme Ferroresonance’. K Pattanapakdee , C Banmongkol . \textit{International Conference on Power Systems Transients (IPST)}, (France) June. 2007. \bibitem[Pattanapakdee and Banmongkol (2007)]{b16}\label{b16} ‘Failure of Riser Pole Arrester due to Station Service Transformer Ferroresonance’. K Pattanapakdee , C Banmongkol . \textit{International Conference on Power Systems Transients (IPST)}, June. 2007. \bibitem[Emin and Tong ()]{b0}\label{b0} ‘Ferroresonance experience in UK: simulations and measurements’. Z Emin , K Y Tong . \textit{International Conference on Power Systems Transients (IPST)}, 2001. \bibitem[Araujo et al. (1993)]{b18}\label{b18} ‘Ferroresonance in power systems: chaotic behaviour’. A E A Araujo , A C Soudack , J R Marti . \textit{IEE Proc.-C} May. 1993. 140 (3) p. . \bibitem[Mork ()]{b20}\label{b20} ‘Five-legged wound/core transformer model: derivation, parameters, implementation, and evaluation’. B A Mork . \textit{IEEE Transactions on Power Delivery} 1999. 14 p. . \bibitem[Chkravarthy and Nayar ()]{b19}\label{b19} ‘Frequency locked and quasi periodic (QP) oscillations in power systems’. S K Chkravarthy , C V Nayar . \textit{IEEE Transactions on Power Delivery} 1997. 13 p. . \bibitem[Mukerjee et al. ()]{b1}\label{b1} ‘Indices for ferroresonance performance assessment in power distribution network’. R N Mukerjee , B Tanggawelu , E A Ariffin , M Balakrishnan . \textit{International Conference on Power System Transients (IPST)}, 2003. \bibitem[Abbasi et al. ()]{b7}\label{b7} ‘Investigation and Control of Unstable Chaotic Behavior Using of Chaos Theory in Electrical Power Systems’. H R Abbasi , A Gholami , M Rostami , A Abbasi . \textit{Series ferroresonance in power systems}, S K Chakravarthy, C V Nayar (ed.) March 2011. 1995. 7 p. . (Electr. Power Energy Sys) \bibitem[Neves and Dommel ()]{b22}\label{b22} ‘On modeling iron core nonlinearities’. W L A Neves , H Dommel . \textit{IEEE Transactions on Power Systems} 1993. 8 p. . \bibitem[Zare et al. ()]{b21}\label{b21} \textit{Performance of various magnetic core models in comparison with the laboratory test results of a ferroresonance test on a 33 kV voltage transformer}, Rezaei Zare , A Mohseni , H Sanaye Pasand , M Farhangi , S Iravani , R . 2006. IEEE. \bibitem[Emin et al. (2001)]{b24}\label{b24} ‘Quantification of the chaotic behavior of ferroresonant voltage transformer circuits’. Z Emin , B A T Zahawi , Y K Tong , M Ugur . \textit{A. H. Nayfeh, APPLIED NONLINEAR DYNAMICS Analytical, Computational and Experimental Methods} June. 2001. 27. 2004. WILEY-VCH Verlag GmbH \& Co. KGaA. 48 (6) . (IEEE Transactions on Circuits and Systems, Fundamental Theory and Applications) \bibitem[Ben-Tal et al. ()]{b4}\label{b4} ‘Studying ferroresonance in actual power systems by bifurcation diagram’. A Ben-Tal , D Shein , S Zissu . \textit{Electric Power Systems Research}, 1999. 49 p. . \bibitem[Thompson and Stewart ()]{b25}\label{b25} J M T Thompson , H B Stewart . \textit{Nonlinear dynamic and chaos}, (Ltd, England) 2002. John Wiley \& Sons. \end{bibitemlist} \end{document}