\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={The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears},
 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{2020-01-15 (revised: 15 January 2020)}
\def\TheID{\makeatother }
\def\TheDate{2020-01-15}
\title{The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears}
\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]{Jian-Hua  Xue}

             \author[2]{Zhen-Hua  Zhang}

             \author[3]{Huan-rui  Wang}

\renewcommand\Authands{ and }

\date{\small \em Received: 9 December 2019 Accepted: 4 January 2020 Published: 15 January 2020}

\maketitle


\begin{abstract}
        


A non-uniform model of the load per unit of length distribution of helical gear with modification and misalignment was proposed based on the meshing stiffness, transmission error, and load-balanced equation. The distribution of unit-line load, transmission error (TE), and contact press of any point on the contact plane were calculated by the numerical method. The feature coordinate system was put forward to implement the helical preliminary design and strength rating. The thermal elastohydrodynamic lubrication (EHL) model of helical gear was established, and the pressure, film, and temperature fields were obtained from the thermal EHL model. The maximum contact temperature and minimum film thickness solved by thermal EHL were applied to check the scuffing load capacity. The highest flash temperature and thinnest film occur in the dedendum of the pinion. The thermal EHL method to evaluate the scuffing load capacity is effective.

\end{abstract}


\keywords{helical gear; meshing stiffness; load distribution; scuffing load capacity.}

\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[{Introduction}]{Introduction}\par
helical gear is an common transmission device and has been widely used in all the fields, especially the machine under high speeds and heavy loads. The load distribution is the foundation of the gear preliminary design and strength rating process. It is known that the load distribution depends on the meshing stiffness of the tooth pair, and the load per unit of length is different at any point in the contact plane. The simple equations in standard ISO \hyperref[b0]{[1,}\hyperref[b1]{2]} to describe the load distribution, which is not in good agreement with experimental results. The contact lines of a helical gear are not parallel with the axial line, and the length of contact lines is dynamic changing in the meshing process.\par
Some studies on the meshing stiffness and load distribution of involute gears can be found in technical literature. Z. Chen et al. \hyperref[b2]{[3]}\hyperref[b3]{[4]}\hyperref[b4]{[5]}\hyperref[b5]{[6]} studied the tooth mesh stiffness and transmission error via the finite element method (FEM). However, this method is feasible but has the problem no generality of the obtained results. Afterward, J.I.Pedreroet al. \hyperref[b6]{[7]}\hyperref[b7]{[8]}\hyperref[b8]{[9]} proposed a method to calculated non-uniform load distribution along the line of contact from the minimum elastic criterion potential, which depends on the transverse contact ratio. Through this method, the author analyzed the bending strength and pitting load capacity of helical gears. But the balanced load equation was not considered in this method. Thus it can't provide the load distribution of any point on the contact plane, and it is hard to locate the maximum value of the load.\par
The heavy load and high-speed gear generate a lot of heat and temperature rise. High contact temperature of lubricant and tooth surfaces at the instantaneous contact position may lead to the breakdown of the lubricant film at the contact interface. The scuffing failure is unpredictable and fatal for the gears system. Therefor the scuffing load capacity is of great importance in the process of helical gear system preliminary design and strength rating, especially for the heavy load and high-speed gear system. ISO \hyperref[b9]{[10,}\hyperref[b10]{11]} provides two methods, namely the flash temperature method and integral temperature method, to evaluate the scuffing load capacity of the gear system, nevertheless the load distribution use the simplified form and the flash temperature calculated based on the Blok flash temperature equation \hyperref[b11]{[12]}, which can't get the accurate flash temperature and need to attach a large safety factor to amend it.\par
The thermal elastohydrodynamic lubrication (EHL) is also a hot research topic. So far, most studies focused on the thermal EHL of spur gears. Wang and Cheng obtained a comprehensive research on a numerical simulation of the contact conditions of straight spur gear pairs \hyperref[b12]{[13,}\hyperref[b13]{14]}. L.M. Li proposed an inverse approach to establish the pressure, temperature rise, and apparent viscosity distribution in an EHL line contact \hyperref[b14]{[15]}. Some methods and beneficial work to study on the EHL of spur gear, and a mass of cases were obtained, but these research mostly stay in the theoretical and ignore the practical application \hyperref[b15]{[16]}\hyperref[b16]{[17]}\hyperref[b17]{[18]}\hyperref[b18]{[19]}. Besides the spur gears, the EHL research of helical gears is rare, no matter under the isothermal or thermal conditions. Recently, P.Yang and P.R.Yang used the multilevel multi-integration method to study the thermal elastohydrodynamic lubrication of tapered rollers in the opposite orientation; this model can be applied to the helical gear system \hyperref[b19]{[20]}. These technical literature mentioned above mostly focus on the calculation method and directly offer the load; they can express the thermal EHL characteristic in theory but can't apply to A Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version I the actual conditions. The scuffing load capacity can be evaluated by the maximum contact temperature and minimum film thickness, which can be solved by the thermal EHL method. The literature of thermal EHL mostly focus on the temperature and film thickness of some single points \hyperref[b15]{[16]}\hyperref[b16]{[17]}\hyperref[b17]{[18]}\hyperref[b18]{[19]}\hyperref[b19]{[20]}. The literature which makes the thermal EHL theory to design and check gear is absent.\par
This paper proposes a method to study the load per unit of length distribution of all the points on contact plane of helical gears accurately based on the balanced Load equation, transmission error, and meshing stiffness. The feature coordinate to simplify the preliminary design and strength check process of helical gears is established. Based on the load distribution, the thermal EHL model, which corresponds more to actual conditions, is put forward. The hydrodynamic pressure, film thickness, and contact temperature, as well as the flash temperature, to check the scuffing load capacity are calculated via the numerical method. 
\section[{II.}]{II.} 
\section[{The Load Distribution Model of Helical Gear}]{The Load Distribution Model of Helical Gear}\par
a) The contact model of the helical gear under heavy load For operating helical gear pair under heavy load, even though the driving pinion rotates at a constant speed, the gear as well as fall behind the angle ? ? than the theoretical location because of the deformation of driving and driven gear along the action line. As Fig. \hyperref[fig_2]{1} shows, the meshing condition viewing from helical gear transverse direction, the profile of the solid line is the theoretical location, and the dashed line is the actual location under deformation. N 1 N 2 is the theoretical action line. Thus for any point K at the action line, the deformation along the action line is   The analysis model of the helical gear is shown as Fig. \hyperref[fig_3]{2}, the contact plane N 1 N 2 N 3 N 4 is the tangent plane of two gear base circle. K 1 K 2 is one of the contact line. The actual action line is A 1 E 1 .The transverse contact ratio calculated by Equ.2.  is the single contact tooth region. In double contact tooth regions, the contact lines always occur double in the same location. The coordinate system is established to describe the actual contact plane A 1 A 2 E 1 E2. The axial direction and action line direction are described by B and Î?" [10]. Î?" is the dimensionless parameter defined as follow: For helical gears, the contact lines at the same time always more than one (depend on the total contact ratio), so setting L denotes the sum of the length of all the contact lines, the sum load of the helical gear is the integral of w k , at any moment, it should be balanced to the extern load, the balanced load equation as follow:C N C N K N 1 1 1 / ) ( ? = Î?"sum b b L L K W r n P dl r l k dl w = = ? = ? ? 1 1 2 0 0 9549 ) ( ? (1)\par
Where C is the pitch point, K the contact point, thus ] , [1 1 E A Î?" Î?" ? Î?" , ] , 0 [ b B ? . 
\section[{Fig. 3: Analysis model of helical real action plane}]{Fig. 3: Analysis model of helical real action plane}\par
Now assuming A 1 is the gear pair approach point (begin meshing), the contact line will move along with the line A 1 E 2 to the recess action point E 2 . Axial contact ratio is expressed by the ? ? .  The length distribution of contact lines with the parameters listed in Table \hyperref[tab_0]{1} shows as Fig. \hyperref[fig_6]{5}. For different gear pairs, the distribution of load per unit of length is different. Compare the gear pair one and two or three and four, the face width doubled, and the length of the contact line almost double as well. The wider the tooth width, the longer the total length. But the value that the maximum minus the minimum has the upper limit value. We are letting the length ratio max min / L L = ? , the length ratio under different helix angle is show in Fig.  {\ref 6}. It is obvious that when 2 < ? ? , the length ratio is only 0.5, and the fluctuation of length is greatly, the maximum is as double as the minimum length. When2 > ? ?\par
, the length ratio is close to 1. 
\section[{Fig. 6: The length ratio varying with the total contact ratio c) The tooth meshing stiffness}]{Fig. 6: The length ratio varying with the total contact ratio c) The tooth meshing stiffness}\par
The gear profile is a complex graphics and can be simply as the combination of a rectangle and a trapezium, as the Fig. \hyperref[fig_7]{7} shows. The gear deformation includes the bending deformation, shear deformation, and contact deformation. The total is the sum of the rectangle deformation and trapezium deformation. The total of the contact deformation point along the action line can be expressed as the sum of the bending deformation B ? , shear deformation, S ? and lean deformation G ? . ? ? ? ? ? ? + ? = 3 ) ( cos 12 3 3 2 r r x r x f i x Br h h h h h S E w ? ? ? ? ? ? ? ? ? ? ? + + = x i r i r i r f i x S h h h h h h h S E w ln ) ( cos ) 1 ( 2 2 1 ? ? ? 2 2 cos 24 f i x x G S E wh ? ? ? =\par
Where h is the height of the tooth, h x the height of contact point, h r the height of the rectangle, x ? the pressure angle of the contact point,a f a r f i S S S h hS h ? ? = ,v i\par
notes the Poisson's ratio of pinion, E i the elastic modulus of pinion or gear, w the load per unit of length.\par
When a pair of tooth meshing, the sum deformation along the action line can be described as the Equ.6?PV ? ? ? ? + + = ? 2 1 (6) Where E v v w PV ? ? )] 1 ( ) 1 [( 2 2 2 2 1 ? + ? =\par
is the contact deformation of the meshing point, so the stiffness of this point can be expressed as Equ.7.? = ? / w k (7)\par
Fig. \hyperref[fig_8]{8} shows the distribution of stiffness along the action line; the maximum occurs close to the pitch point. The regulation of stiffness is similar to the inverse unitary potential ) (? v in Ref. [7-9], which all reflect the capacity that the tooth bears the deformation. ? ? ? = sum L sum dl TE l k W 0 ) (\textbf{(8)}\par
For the convenience of numerical calculation, the line A 1 E 1 is divided into N nodes, and) 1 /( 1 1 ? = ? N E A l\par
, so the total load can be expressed as Equ.9, the number of the contact line is m.b x ij m j n i ij sum l TE k W ? sin / ) 0 ), max(( 1 1 ? = = ? ? ? = ?? (9)\par
The load distribution calculation process of helical gear is shown as Fig.  {\ref 9}. According to the equations above, we can calculate the load distribution of every point on the contact plane. Firstly, letting L=L 0 , k=k c , and TE 0 will be calculated; secondly, from the Equ.9, the sum load was obtained, then modified the TE according to Equ.10, until the sum load satisfied the Equ.11, then repeat this process until the results of all points were obtained.  r i r i x i r i x i r i h h h h S E w ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ?\par
The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears© 2020 Global Journals\par
Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version The contact model of the helical gear can be regarded as two tapered rollers \hyperref[b19]{[20]}. In this model, the actual geometric model is two oval in oxz coordinate systems, and the geometric parameters as follows: According to the Hertz contact theory, the maximum contact press of the helical gear is expressed as Equ.14. 
\section[{R Ew}]{R Ew}H ? ? ? 2 4 = (14)\par
Where  f) The load distribution calculation of helical gear without modification In the calculation process, the modification, setting ij ? to 0, then the transmission and the unit-line load of contact plane can be obtained. The threedimension unit-line load distribution is shown in Fig. \hyperref[fig_13]{11}. The load distribution on the contact plane is not only depends on the transmission error, but also depends on the stiffness distribution. The value of the load is small in the dedendum and addendum region of helical gear, and large close to pitch point. For the different helix angles and different face width, the load distribution is different. As we know, the sliding speed is maximum in the dedendum or addendum. From the 4 cases, the biggest occurs in the begin meshing and engaging-out point of, the helical contact plane.2 1 2 1 R R R R R + = the curvature sum, E is equivalent to Young's modulus, ? ? ? ? ? ? ? ? ? + = 2 2 2 1 2 1 1 - 1 2 1 1 E E E ? ? .\par
The transmission errors of helical gear pairs shown in Fig. \hyperref[fig_14]{12}, its distribution law corresponding to the distribution of length. The fluctuation of transmission error is the main indicator of the vibration and dynamic load. It is significant to choose the advisable helix angle and face width to make the fluctuation of transmission error minimum.\par
The contact press distribution of the helical gears is shown in Fig. \hyperref[fig_2]{13}. It is similar to the unit-line load.\par
But the contact stress of foot is greater than that of top of pinion. It is because the sum radius of the curvature of the root is smaller than that of the tooth top. For a helical cylindrical gear without modification, the maximum contact press is located at the engagement of the pinion root. As shown in Fig. \hyperref[fig_2]{14}, it is a fatigue pinion of an electric axle after the loading bench test, fatigue pitting at the meshing point of the root of the pinion. The contact stress with the modification and misalignment In practical application, the tooth surface of helical gear needs to be modified. For one reason, the maximum contact press located at the engagement of the pinion root. For the second reason, the shafts, bearings, and the housing will be deformed under the heavy load.\par
The modification includes the profile modification and helix modification. Profile modification includes profile crowning, pressure angle modification, tip relief, and root relief. The helix modification includes lead crowning, helix angel modification, and end relief. According to the calculation method in Fig.  {\ref 9}, the transmission error(TE), unit-line load distribution, and contact press distribution of helical gear could be obtained. For the four gear pairs in table 1, the modification parameters are shown as table \hyperref[tab_1]{2}.   In a transmission system, due to the deformation of the transmission shaft, the machining error of the gearbox, and the changes of the bearing stiffness, the gears pair will operate with misalignment. The contact state with misalignment can be simulated by the helix angle modification. Take gear pair 4 for example, the contact press distribution with different helix modification is shown in Fig. \hyperref[fig_18]{17}. With the increase of tooth inclination deviation, the contact stress inclines to one end, and the maximum stress is also increasing. The feature coordinate system From the above analysis, all the results in the contact plane were calculated by the numerical method. But in the preliminary design and strength rating process of helical gear, our focus is the maximum stress or the highest contact temperature of the contact plane. For the helical gear, the sliding speed is large in the addendum and dedendum region, and most scuffing failure occurs in these regions. The sliding speed approximately zero in the region close to the pitch point, so it is safer than the addendum and dedendum region. The 3-dimensional load distribution covers all the information about the contact plane, but it is timewasting and not intuitionistic. So the feature coordinate system should be established. 
\section[{Fig. 18: Unit-linear load and transmission error distribution}]{Fig. 18: Unit-linear load and transmission error distribution}\par
In Fig. \hyperref[fig_0]{3}, the line A 1 E 2 is recommended as the feature coordinate system, marked by ? , defined the same as Î?" . The A 1 is the approach point (begin meshing), and E 2 is the recess point. The parameters in feature coordinate cover both the tooth profile and axial information. The unit-line load distribution along with the feature coordinate is shown in Fig. \hyperref[fig_8]{18}.Compared the feature coordinate and 3-dimension coordinate, the maximum is the same in the addendum and dedendum region, the value in the feature coordinate system is a little lower than the three-dimension coordinate close to the pitch point, because the dangerous region of helical gear is the addendum and dedendum, so that the feature coordinate can satisfy the need of design and strength ratting. Furthermore, the feature coordinate is more succinct than the three-dimension. From the load distribution of gear pair 1 and 3, or gear pair 2 and 4, the helix angle has influenced a lot on the load distribution, in the case, the helix angle from the 9.8 to 20.2, and the maximum load per unit of length from the In the Ref.10, the load distribution of narrow helical gear is given as the Fig. \hyperref[fig_19]{19}, the buttressing effect near the end points A and E of the line of action, compared the results in this paper, the load distribution is similar with the gear pair 1,but is different from the gear pair 2, so the ISO can reflect the characteristic to some extent. Still it can't adapt to all the conditions. In this way, the method proposed by this paper is effective, and the point AU and EU is given by Equ. 15.b EU E A AU ? sin 2 . 0 = Î?" ? Î?" = Î?" ? Î?" (15)\par
In the ISO standard theory, the helical gear with the total contact ratio When gears pair subjected to heavy loads, the lubricating film may not separate the surfaces adequately, leading to localized damage of the tooth surface. This type of failure is known as scuf fing and it is checked by maximum contact temperature and film thickness. Scuffing failure may occur at any stage during the lifetime of a set of gears, and it may lead to failure in a few of hours if the contact temperature is too high. So it is significant to research the scuffing load capacity of the gears.\par
The high contact temperature is the main reason for scuffing failure, and the interfacial contact temperature conceived as the sum of two components, namely the bulk and flash temperature \hyperref[b9]{[10]}. The bulk temperature fields are constant with time. Nevertheless, the flash temperature varies with time, and only appears in the interface of pinion and gear. In the thermal steady state, the bulk temperature of the tooth is higher than the ambient oil temperature.f M B t t t + = (16)\par
Where t B is the contact temperature, t M the bulk temperature, and t f the flash temperature. The bulk temperature can be calculated by the finite element method (FEM). For the pinion of the pair 1, the bulk temperature is 360K (87?), the ambient oil temperature is 60?.\par
The internal energy of oil film will increase due to the compression and viscous damping, and the temperature of oil film and tooth will rise because of the heat convection and conduction. After some time, the whole system will be on a thermal steady state. The bulk temperature of gear was higher than the ambient oil temperature. Therefore, in the inlet zone, the tooth profile is the heat source to heat the lubrication film, and the lubrication film temperature will rise instantly as high as the bulk temperature because that the lubricating film thickness is so thin. Then the temperature of the lubrication film will get higher due to the compression and viscous damping. In this region, the lubricating film is the heat source to heat the gear by convection and conduction. The ambient temperature was chosen as the inlet temperature in the technical literature \hyperref[b15]{[16]}\hyperref[b16]{[17]}\hyperref[b17]{[18]}\hyperref[b18]{[19]}\hyperref[b19]{[20]}. By this means, the interface temperatures of the inlet contact region solved by the numerical method will lower than the interface bulk temperature; It is not in good agreement with the fact. In this paper, the interface bulk temperature was taken as the inlet temperature, which is more reasonable. 
\section[{b) The thermal elastohydrodynamic lubrication equations i. Reynolds equation}]{b) The thermal elastohydrodynamic lubrication equations i. Reynolds equation}\par
For the thermal steady state, the Reynolds equation of helical gear can be expressed in the following form as:? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? 0 1 1 2 2 2 ) ( ) ( F F x u u x h u x p F x ? ?\textbf{(17)}\par
The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears © 2020 Global Journals 
\section[{Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version}]{Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version}\par
Wheredz F h ? = 0 0 / 1 ? , dz z F h ? = 0 1 /? , dz F F z z F h ? ? ? = 0 0 1 2 ) / ( /? ? .\par
The sliding speed of contact point:30 / 2 , 1 i i i i R n u ? = = ?\textbf{(18)}\par
In this equation, the mass density and the viscosity of lubricant can be described by Dowson and Higginson temperature-pressure-density relationship \hyperref[b20]{[21]} and Roelands temperature-pressure-viscosity relationship \hyperref[b21]{[22]} can be expressed as: ? ? ? ? ? ? ? ? ? ? × + × + = ? ? ) (? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + ? = ? 0 138 138 ) 1 ( 1 exp ) , ( 0 2 1 0 s z t t p A A t p ? ?\textbf{(20)}\par
Where 67 . 9 ln 0 1+ = ? A , 9 - 2 10 1 . 5 × = A , ) /( 2 1 A A z ? = , )) 138 ( /( 0 1 0 ? = t A s ? .\par
ii. The film thickness is given asds s x x p E R x h h b a x x ? ? ? + = ln ) ( 4 2 2 0 ? (21)\par
Whereh 0 is film thickness constant depend on the applied load. 
\section[{iii. The balanced load equation}]{iii. The balanced load equation}\par
The hydrodynamic pressure of the contact region must be in equilibrium with the input load, which calculated by the numerical method in section 2.Thus the pressure equilibrium is expressed by the following equation:? = b a x x dx x p w ) (\textbf{(22)}\par
The three-dimensional temperature distribution of oil film was calculated by the energy equation as follows:? ? ?? ? ?? ?? ? ? ? ? ?? ?? ? ion FluidFrict Expation n Compressio Conduction H p z u x p u t t z t x t u c 2 / eat 2 2 Convection ) ( ) ( ? ? ? ? ? ? ? ? = ? ? ? ? + ? ? ? ? ? ? ? ? ? ? ? (23)\par
The boundary condition for Eq.23 is the thermal interface equation: The equation governing the shear stress in the lubricant film is a momentum equation which for the Newtonian fluid is as follow:? ? ? ? ? ? ? + ? ? ? = + ? ? ? = ? ? ? ? = ? ? = x b h z x a z t s x ds z t u c h x t t s x ds z t u c x t? ? ? ? ? ? ? ? ? ? = ? ? z u z x p ? (24) iv. Numerical method\par
All equations and boundary conditions mentioned above about thermal were written into dimensionless forms to facilitate the numerical calculation. There are pressure loop and temperature loop in the entire calculation process. In the pressure loop, assuming the temperature of all the nodes is the ambient temperature (360K). In the contact coordinate in the rolling direction ( x direction), the whole contact region (from inlet position to outlet position) was divided into 81 non-equidistant nodes; the interval is large in the inlet region and small in the secondary pressure peak region. The distribution of pressure and the film thickness can be calculated by the Newton-Raphson method depending on the Equ.17-24.Divided the film gap (obtained in the pressure loop, z direction) into 21 equidistant nodes. The temperature of all the nodes can be calculated depend on the equ.20-21 and equ.24-25.Then repeat these two processes until the pressure error less than 5 -10 .The lubricant properties are shown in Table \hyperref[tab_1]{2}.  N/m shown in Fig. \hyperref[fig_2]{14}. For the gear transimission system under heavy load, the film center pressure is higher; in this case, the maximum pressure is as high as 1.51GPa, the minimum film thickness is 0.42 m µ . Because of the heavy load, the curve of oil pressure and oil thick has't the typical characteristic of thermal EHL. The secondary pressure peak close to the outlet and inconspicuous. The distribution of the pressure is similar to Hertzian pressure distribution. The maximum occurs the location x=0 and the heavier of the load, the higher of the maximum. The heavier, the wider the contact zone. Furthermore, it has little influence on the minimum oil thickness. The three-dimension and boundaries temperature distribution of the contact region are shown in Fig. \hyperref[fig_23]{21}. The film center temperature is much higher than two boundaries; this is because the heat conducts from the film to all the teeth to maintain the bulk temperature. The distribution of film center temperature is similar to the film center pressure and peaked in the Hertzian maximum, but two boundaries are different from the film center; the temperature increase until close to the outlet and slightly decrease. The reason is in the rear part of the contact zone, the compression effect is no longer generate the heat, so the heat from the film center to the two boundaries via conduction effect and to the inner tooth through convection. 
\section[{b) Comparison with the Blok's flash temperature}]{b) Comparison with the Blok's flash temperature}\par
The contact and flash temperature is the main reason for the gear scuffing failure, Blok derived the flash temperature equation in 1937 \hyperref[b11]{[12]} as follow:0 2 2 2 2 1 1 1 1 2 1 ) ( 7858 . 0 b u c u c u u fw t f ? ? ? ? + ? = (25)\par
Where f=0.06 is coefficient of friction, w the load per unit of length,  In the thermal EHL results, the maximum temperature rise of two boundaries is its flash temperature. The flash temperature along the feature coordinate was calculated by the thermal EHL method and Blok equation as show in Fig. \hyperref[fig_27]{22}. In the dedendum region of the pinion (close to approach point), the thermal EHL flash temperature is higher than the Blok result, and in the addendum region of the pinion(close to recess point), the Blok temperature is lower than thermal EHL flash temperature. In theory, the temperature-pressure-density and temperaturepressure-viscosity effect and the compression of oil film were considered in the thermal EHL theory. The Blok flash temperature equation is very concise and corresponding to the experiment to some extent. However, it loses sight of the influence of the bulk temperature, so it is not a good agreement of the experiment under the heavy load. Thus the thermal EHL theory is even close to the actual condition. The safety factor defined as follow \hyperref[b9]{[10]}: Where t s is the scuffing temperature, t oil the ambient temperature, t cmax the maximum contact temperature. When gear teeth are separated completely by a full fluid film of lubricant, there is no contact between the asperities of tooth surfaces, and usually, there is no scuffing and wear. For the heavy load helical gear, the thickness of oil film is very small, and incidental asperity contact takes place. As the minimum film thickness decreases, the number of contacts increases, abrasive wear, adhesive wear, and scuffing became possible. So the minimum thickness of lubricant film is a property of scuffing load capacity, especially for the gear pair under heavy load. The minimum thickness equation was given by Dowson and Higginson \hyperref[b21]{[22]} as the Equ.27. The safe factor can be measured by the thickness ratio. When the thickness ratio ) 4 , 1 ) ? ? ,it is considered as mixed friction. When the thickness ratio is less than 1, the scuffing failure probably takes place to a great extent. \begin{figure}[htbp]
\noindent\textbf{3}\includegraphics[]{image-2.png}
\caption{\label{fig_0}The 3 mw}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-3.png}
\caption{\label{fig_1}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1}\includegraphics[]{image-4.png}
\caption{\label{fig_2}Fig. 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2}\includegraphics[]{image-5.png}
\caption{\label{fig_3}Fig. 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-6.png}
\caption{\label{fig_4}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-7.png}
\caption{\label{fig_5}The.}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5}\includegraphics[]{image-8.png}
\caption{\label{fig_6}Fig. 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{7}\includegraphics[]{image-9.png}
\caption{\label{fig_7}Fig. 7 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{8}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Fig. 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{910}\includegraphics[]{image-11.png}
\caption{\label{fig_11}Fig. 9 :Fig. 10 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-12.png}
\caption{\label{fig_12}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{11}\includegraphics[]{image-13.png}
\caption{\label{fig_13}Fig. 11 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{12}\includegraphics[]{image-14.png}
\caption{\label{fig_14}Fig. 12 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1314}\includegraphics[]{image-15.png}
\caption{\label{fig_15}Fig. 13 :Fig. 14 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{15}\includegraphics[]{image-16.png}
\caption{\label{fig_16}Fig. 15 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{16}\includegraphics[]{image-17.png}
\caption{\label{fig_17}Fig. 16 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{17}\includegraphics[]{image-18.png}
\caption{\label{fig_18}Fig. 17 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{19}\includegraphics[]{image-19.png}
\caption{\label{fig_19}Fig. 19 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-20.png}
\caption{\label{fig_20}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-21.png}
\caption{\label{fig_21}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{220}\includegraphics[]{image-22.png}
\caption{\label{fig_22}? and 2 ?Fig. 20 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{21}\includegraphics[]{image-23.png}
\caption{\label{fig_23}Fig. 21 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-24.png}
\caption{\label{fig_24}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-25.png}
\caption{\label{fig_25}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{22}\includegraphics[]{image-26.png}
\caption{\label{fig_27}Fig. 22 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{23}\includegraphics[]{image-27.png}
\caption{\label{fig_28}Fig. 23 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.47458333333333336\textwidth}P{0.09208333333333334\textwidth}P{0.09208333333333334\textwidth}P{0.095625\textwidth}P{0.095625\textwidth}}
\tabcellsep Gear Pair 1\tabcellsep Gear Pair 2\tabcellsep Gear Pair 3\tabcellsep Gear Pair 4\\
Number of teeth (pinion/gear)\tabcellsep 23/30\tabcellsep 23/30\tabcellsep 23/30\tabcellsep 23/30\\
Normal pressure angle\tabcellsep 20\tabcellsep 20\tabcellsep 20\tabcellsep 20\\
Normal module, mm\tabcellsep 3\tabcellsep 3\tabcellsep 3\tabcellsep 3\\
Face width, mm\tabcellsep 20\tabcellsep 40\tabcellsep 20\tabcellsep 40\\
Input power(kw)\tabcellsep 50\tabcellsep 50\tabcellsep 50\tabcellsep 50\\
Modulus of elasticity ?? Helix angle,°\tabcellsep 9.8\tabcellsep 9.8\tabcellsep 20.2\tabcellsep 20.2\end{longtable} \par
  {\small\itshape [Note: 1 /?? 2 , GPa 206/206 206/206 206/206 206/206 Poisson ratio, ?? 1 /?? 2 0.3/0.3 0.3/0.3 0.3/0.3 0.3/0.3]} 
\caption{\label{tab_0}Table 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2} \par 
\begin{longtable}{P{0.735672514619883\textwidth}P{0.11432748538011694\textwidth}}
Modification form\tabcellsep Modification/um\\
profile crowning(barreling) ?? ??\tabcellsep 10\\
Lead crowning ?? ??\tabcellsep 10\\
tip relief\tabcellsep 5\\
End relief\tabcellsep 0\\
Helix angle modification fH?\tabcellsep 0\\
Pressure angle modification?fH?\tabcellsep 0\end{longtable} \par
 
\caption{\label{tab_1}Table 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2} \par 
\begin{longtable}{P{0.775438596491228\textwidth}P{0.07456140350877193\textwidth}}
Symbol,unit\tabcellsep Value\\
Ambient viscosity of lubricant , 0 ? ,Ns/ 2 m\tabcellsep 0.08\\
Specific heat of lubricant, c, J/kgK\tabcellsep 2000\\
Specific heats of solids,J/kgK\tabcellsep 470\\
Thermal conductivity of lubricant, ? , W/mK\tabcellsep 0.14\\
Thermal conductivities of solids a and b, 1\tabcellsep \end{longtable} \par
 
\caption{\label{tab_3}Table 2 :}\end{figure}
 			\footnote{b r base radius of gear, m ba} 			\footnote{© 2020 Global Journals} 			\footnote{x i f x Bt h h h h h h h h h h} 			\footnote{© 2020 Global Journals g)} 			\footnote{© 2020 Global Journals h)} 		 		\backmatter  			 \par
The minimum film thickness solved by thermal EHL and Dowson equation is shows as Fig.  {\ref 23}. The minimum thickness curves are similar along the feature coordinate of gear pair one and two. The minimum film thick is lowest in the dedendum of the pinion, so the dedendum of the pinion is dangerous. The results are consistent with the maximum flash temperature method. The value from thermal EHL is lower than that from the Dowson equation. The Dowson equation didn't consider the influence of temperature on film thickness.\par
In contrast to the fatigue damage, a single momentary overload may initiate scuffing failure, so the scuffing capacity is crucial for the heavy load and highspeed helical gear system. The minimum film thickness and maximum contact temperature, as well as the flash temperature obtained by the thermal EHL theory, are applied to check the scuffing strength. The scuffing capacity can be checked by the flash temperature or minimum thickness method.\par
V. 
\subsection[{Conclusions}]{Conclusions}\par
In this paper, a model of non-uniform along the contact line of helical gear teeth, obtained from the meshing stiffness, transmission error, and balanced load equation, has been proposed. The feature coordinate system was established to simply preliminary design and strength rating process. The thermal effect of helical gear lubrication will affect the scuffing load capacity. The thermal EHL method is applied to check the scuffing load capacity.\par
The following conclusions can be drawn:\par
1. The length of contact lines depends on the face width and basic helix angle. And the length is oline is not paralleled with the axial line. The fluctuation is evident when the total contact ratio less than 2.The distribution of transmission error is similar to the sum length. The distribution of load per unit of length is dependent on the meshing stiffness and transmission, it is similar to the meshing stiffness along the contact line and similar to the transmission error along the Baxis. The modification of gears can greatly improve the distribution of contact stress. The maximum contact stress is concentrated in the center of the tooth surface to avoid the contact between the tooth surface. The contact stress inclines to one end, and the maximum stress is also increasing under with misalignment. 2. The feature coordinate system established in this paper can cover bother the tooth profile and axial information. It can satisfy the need for preliminary design and strength check of helical gears. The weakest strength region is the dedendum of the pinion. A load of the dedendum region in the feature coordinate system is as same as the threedimension coordinates system.\par
3. The thermal EHL theory can provide more information to check the scuffing load capacity. The highest contact temperature and minimum film thickness method were applied to evaluate the scuffing load capacity. The pressure distribution under heavy load is similar to the Hertzian pressure distribution, and the secondary presser peak close the outlet and inconspicuous. The temperature of the film center is much higher than two boundaries. The film center temperature distribution is similar to the pressure distribution. The two boundaries temperature increase until close to the outlet and slightly decrease. The thermal EHL flash temperature is corresponding to the Blok flash temperature but higher in the dedendum region and lower in the addendum region of the pinion. The minimum thickness film is smaller than the Dowson equation. So the scuffing load capacity that solved by thermal EHL is more practical than the traditional method.			 			  				\begin{bibitemlist}{1}
\bibitem[Fernandez Del Rincon et al. ()]{b3}\label{b3} 	 		\textit{A model for the study of meshing stiffness in spur gear transmissions. Mechanism and Machine Theory},  		 			A Fernandez Del Rincon 		,  		 			F Viadero 		,  		 			M Iglesias 		.  		2013. p. .  	 
\bibitem[Wang and Cheng ()]{b13}\label{b13} 	 		‘A Numerical solution to the dynamic load, film thickness and surface temperatures in spur gears, part II results’.  		 			K L Wang 		,  		 			H S Cheng 		.  	 	 		\textit{Journal of Mechanical Design}  		1981. 103 p. .  	 
\bibitem[Zhang and Fang ()]{b4}\label{b4} 	 		\textit{Analysis of tooth contact and load distribution of helical gears with crossed axes. Mechanism and Machine Theory},  		 			Y Zhang 		,  		 			Z Fang 		.  		1999. 34 p. .  	 
\bibitem[Yang and Yang ()]{b19}\label{b19} 	 		‘Analysis on the thermal elastohydrodynamic lubrication of tapered rollers in opposite orientation’.  		 			Ping Yang 		,  		 			Peiran Yang 		.  	 	 		\textit{Tribology International}  		2007. 40 p. .  	 
\bibitem[Wang and Cheng ()]{b12}\label{b12} 	 		‘Anumerical solution to the dynamic load, film thickness and surface temperatures in spur gears, part I analysis’.  		 			K L Wang 		,  		 			H S Cheng 		.  	 	 		\textit{Journal of Mechanical Design}  		1981. 103 p. .  	 
\bibitem[Calculation of Load Capacity of Spur and Helical Gears -Part 2: Calculation of Surface Durability (Pitting), International Organization for Standardization ()]{b0}\label{b0} 	 		\textit{Calculation of Load Capacity of Spur and Helical Gears -Part 2: Calculation of Surface Durability (Pitting), International Organization for Standardization},  		 ISO Standard 6336-2:1996.  		1996. Geneva, Switzerland.  	 
\bibitem[Calculation of Load Capacity of Spur and Helical Gears -Part 3: Calculation of Tooth Bending Strength International Organization for Standardization ()]{b1}\label{b1} 	 		‘Calculation of Load Capacity of Spur and Helical Gears -Part 3: Calculation of Tooth Bending Strength’.  		 ISO Standard 6336-3:1996.  	 	 		\textit{International Organization for Standardization}  		1996.  	 
\bibitem[Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Flash temperature method ISO Standard ()]{b9}\label{b9} 	 		‘Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Flash temperature method’.  		 13989-1:2000.  	 	 		\textit{ISO Standard},  				 (London)  		2000. The United Kingdom.  	 
\bibitem[Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Integral temperature method ()]{b10}\label{b10} 	 		\textit{Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Integral temperature method},  		 ISO Standard 13989-2:2000.  		2000. London: The United Kingdom.  	 
\bibitem[Roelands ()]{b20}\label{b20} 	 		\textit{Correlation aspects of Viscositytemperature-pressure Relationship of Lubrication oils},  		 			Cja Roelands 		.  		1966. Netherlands.  		 			Delft Univeristy of Technology 		 	 
\bibitem[Pedrero et al. ()]{b7}\label{b7} 	 		‘Critical stress and load conditions for bending calculations of spur and helical gear’.  		 			J I Pedrero 		,  		 			M Pleguezuelos 		,  		 			M Muñoz 		.  	 	 		\textit{International Journal of Fatigue}  		2013. 48 p. .  	 
\bibitem[Pedrero et al. ()]{b8}\label{b8} 	 		‘Critical stress and load conditions for pitting calculations of spur and helical gear’.  		 			J I Pedrero 		,  		 			M Pleguezuelos 		,  		 			M Muñoz 		.  	 	 		\textit{Mechanism and Machine Theory}  		2011. 46 p. .  	 
\bibitem[Anuradha and Kumar ()]{b16}\label{b16} 	 		‘Effect of lubricant selection on EHL performance of involute spur gears’.  		 			P Anuradha 		,  		 			Punit Kumar 		.  	 	 		\textit{Tribology International}  		2012. 50 p. .  	 
\bibitem[Dowson and Higginson ()]{b21}\label{b21} 	 		\textit{Elastohydrodynamic Lubrication},  		 			D Dowson 		,  		 			G R Higginson 		.  		1977. Oxford: Pergamon Press.  	 
\bibitem[Yesilyurt et al. ()]{b5}\label{b5} 	 		‘Gear tooth stiffness reduction measurement using modal analysis and its use in wear fault severity assessment of spur gears’.  		 			Yesilyurt 		,  		 			Andrew D Gu 		,  		 			Ball 		.  	 	 		\textit{NDT\&E International}  		2003. 36 p. .  	 
\bibitem[Chu and Hsu ()]{b14}\label{b14} 	 		‘Jaw-Ren Lin, dtal. Inverse approach for calculating temperature in EHL of line contacts’.  		 			Li-Ming Chu 		,  		 			Hsiang-Chen Hsu 		.  	 	 		\textit{Tribology International}  		2009. 42 p. .  	 
\bibitem[Pedrero et al. ()]{b6}\label{b6} 	 		\textit{Load distribution model along the line of contact for involute external gears, Mechanism and Machine Theory},  		 			J I Pedrero 		,  		 			M Pleguezuelos 		,  		 			M Artés 		,  		 			J A Antona 		.  		2010. 45 p. .  	 
\bibitem[Chen and Shao ()]{b2}\label{b2} 	 		\textit{Mesh stiffness calculation of a spur gear pair with tooth pro file modification and tooth root crack. Mechanism and Machine Theory},  		 			Zaigang Chen 		,  		 			Yimin Shao 		.  		2013. 62 p. .  	 
\bibitem[Almqvist and Larsson ()]{b15}\label{b15} 	 		‘The Navier-Stokes approach for thermal EHL line contact solutions’.  		 			T Almqvist 		,  		 			R Larsson 		.  	 	 		\textit{Tribology International}  		2002. 35 p. .  	 
\bibitem[Blok ()]{b11}\label{b11} 	 		‘Theoretical study of temperature rise at surfaces of actual contact under oiliness lubricating conditions’.  		 			H Blok 		.  	 	 		\textit{Proceedings of the General Discussion on Lubrication and Lubricants[J]. Imech E},  				 (the General Discussion on Lubrication and Lubricants[J]. Imech E)  		1937. 2 p. .  	 
\bibitem[Bobach et al. ()]{b17}\label{b17} 	 		‘Thermal elastohydrodynamic simulation of involute spur gears incorporating mixed friction’.  		 			Lars Bobach 		,  		 			Ronny Beilicke 		,  		 			Dirk Bartel 		.  	 	 		\textit{Tribology International}  		2012. p. .  	 
\bibitem[Wang et al. ()]{b18}\label{b18} 	 		‘Transient thermo elastohydrodynamic lubrication analysis’.  		 			P W Wang 		,  		 			H Q Li 		,  		 			J W Tong 		.  	 	 		\textit{Tribology International}  		2004. 37 p. .  	 
\end{bibitemlist}
 			 		 	 
\end{document}
