// torquemodel.cpp: implementation of the torquemodel class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "2level.h"
#include "torquemodel.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

torquemodel::torquemodel()
{
	
	
	
	
}

torquemodel::~torquemodel()
{

}
double torquemodel::torque(int i,global x)
{
	de_drag=x.de_drag[i];
	f=x.f[i];
	Tb=x.Tb[i];
	Tf=x.Tf[i];
	m=20;
	C8=sqrt((1-EE*EE)/YY);
	H=x.H[i];
	h=x.h[i];
	R=x.R[i];
	C15=0.5;
	D2=0.58;
	D3=0.42;
	D1=0.98;
	B=x.B[0];
	double P=x.P[i];
	double v=x.v[i];
	double torque1;  //張應力產生的力矩
	double torque2;  //軋制力產生的力矩
    double torque3;  //損耗功率
	double DG;
	
	double dh=H-h;
	double RR=R*(1+AA*P/B/dh);
	double lc=sqrt(RR*dh);
	double torque=0;
	double r=(H-h)/H;
	DG=1.05+(0.07+1.32*r)*f*sqrt(RR/H)-0.85*r;
	torque1=B*de_drag*R*(H-h)*DG;
	torque2=B*R*h*(Tb-Tf);
	torque3=m*P*v/6000;
	torque=torque1+torque2+torque3;   //unit:KN*mm,N*m;
	
	return torque;
}