/*                                                       */
/*  電気学会精度検証用モデル 表面磁石構造電動機 '01.6.5  */
/*                                                       */

/*             min < pal?(defolt)< max -> パラメータ項目         */
pal1 = 3.50 /* 1.00 < pal1 (3.50) < 5.00 -> (1)ステータティースの厚み */
pal2 = 1.50 /* 0.20 < pal2 (1.50) < 5.00 -> (2)ティース先端空隙部の幅 */
pal3 = 0.50 /* 0.20 < pal3 (0.50) < 1.00 -> (3)ｘｙ軸上の空隙部の厚み */
pal4 = 45.0 /* 0.00 < pal4 (45.0) < 75.0 -> (4)磁石の切欠き部分の角度 */
pal5 = 0.50 /* 0.00 < pal5 (0.50) < 1.50 -> (5)永久磁石の角部分の厚み */
pal6 = 3.20 /* 1.00 < pal6 (3.20) < 4.40 -> (6)ステータボルト穴の直径 */

/* 変数の定義 */
r0 = 26.7
r1 = 19.8
r2 = 13.15
r3 = 12.65
r4 = 12.15
r5 = 11.35
r6 = 7.85
r7 = (r1*4+r2)/5
r9 = (r4+r5*5)/6

/* 解析領域設定 */
bound = circle{0,0,r0,15,0,90}

/* ステータコア */
fig[1] = circle{0,0,r0,15,0,90}

/* ==================== ループ開始 ==================== */
/* コイルとステータ側空隙部をfor文を用いて作成する */
for (c0=0,60,30) {
	c1 = c0+asin(pal1/2,r1)
	c2 = c0+30-asin(pal1/2,r1)
	c3 = c0+asin(pal1/2,r2)
	c4 = c0+15-asin(pal2/2,r3)
	c5 = c0+15+asin(pal2/2,r3)
	c6 = c0+30-asin(pal1/2,r2)
	c7 = c0+15-asin(pal2/2,r4)
	c8 = c0+15+asin(pal2/2,r4)
	
	x1 = r1*cos(c1)		 y1 = r1*sin(c1)
	x2 = r1*cos(c2)		 y2 = r1*sin(c2)
	x3 = r2*cos(c3)		 y3 = r2*sin(c3)
	x4 = r3*cos(c4)		 y4 = r3*sin(c4)
	x5 = r3*cos(c5)		 y5 = r3*sin(c5)
	x6 = r2*cos(c6)		 y6 = r2*sin(c6)
	x7 = (r4-.1)*cos(c7) y7 = (r4-.1)*sin(c7)
	x8 = (r4-.1)*cos(c8) y8 = (r4-.1)*sin(c8)

	x10 = (r2+((r1-r2)/4)*1)*cos((c3-((c3-c1)/4)*1))
	y10 = (r2+((r1-r2)/4)*1)*sin((c3-((c3-c1)/4)*1))
	x11 = (r2+((r1-r2)/4)*2)*cos((c3-((c3-c1)/4)*2))
	y11 = (r2+((r1-r2)/4)*2)*sin((c3-((c3-c1)/4)*2))
	x12 = (r2+((r1-r2)/4)*3)*cos((c3-((c3-c1)/4)*3))
	y12 = (r2+((r1-r2)/4)*3)*sin((c3-((c3-c1)/4)*3))

	x13 = (r2+((r1-r2)/4)*1)*cos((c6-((c6-c2)/4)*1))
	y13 = (r2+((r1-r2)/4)*1)*sin((c6-((c6-c2)/4)*1))
	x14 = (r2+((r1-r2)/4)*2)*cos((c6-((c6-c2)/4)*2))
	y14 = (r2+((r1-r2)/4)*2)*sin((c6-((c6-c2)/4)*2))
	x15 = (r2+((r1-r2)/4)*3)*cos((c6-((c6-c2)/4)*3))
	y15 = (r2+((r1-r2)/4)*3)*sin((c6-((c6-c2)/4)*3))

/* コイル */
	fig[5+c0/30] = arcpoly{
	 x1,  y1, r1,  4,
	 x2,  y2,  0,  1,
	x15, y15,  0,  1,
	x14, y14,  0,  1,
	x13, y13,  0,  1,
	 x6,  y6,  0,  2, 
	 x5,  y5,-r3,  4,
	 x4,  y4,  0,  2,
	 x3,  y3,  0,  1,
	x10, y10,  0,  1,
	x11, y11,  0,  1,
	x12, y12,  0,  1
	}

/* ステータ側空隙部 */
	fig[8] = arcpoly{
	 x7, y7,  0,  1,
	 x4, y4, r3,  4,
	 x5, y5,  0,  1,
	 x8, y8,  0,  1
	}
}
/* ==================== ループ終了 ==================== */

/* ステータ側空隙部 */
fig[8] = circle{0,0,r4,90,0,90}

/* ロータ側空隙部 */
fig[9] = circle{0,0,r9,90,0,90}

/* 変数定義 */
pala1 = pal3
pala2 = asin(pala1,r6)
pala3 = (r5-r6)/2+r6+pal5
pala4 = asin(pala1,pala3)
pala5 = pal4
alpa  = pala3*cos(pala4)
bata  = pala3*sin(pala4)

/* 円と直線との交点(ptx,pty)を定式化 */
ptx = sqrt(((tan(pala5)*tan(pala5)*(bata-tan(pala5)*alpa)
        *(bata-tan(pala5)*alpa))/(1+tan(pala5)*tan(pala5))
        +r5*r5-((bata-tan(pala5)*alpa)*(bata-tan(pala5)*alpa)))
        /(1+tan(pala5)*tan(pala5)))-(tan(pala5)*(bata-tan(pala5)
        *alpa))/(1+tan(pala5)*tan(pala5))
pty = tan(pala5)*ptx+(bata-tan(pala5)*alpa)

/* 永久磁石 */
fig[4] = arcpoly{
   r6*cos(pala2),    r6*sin(pala2),  0, 1,
pala3*cos(pala4), pala3*sin(pala4),  0, 1,
             ptx,              pty, r5,35,
             pty,              ptx,  0, 1,
pala3*sin(pala4), pala3*cos(pala4),  0, 1,
   r6*sin(pala2),    r6*cos(pala2),-r6, 8
}

/* ロータコア */
fig[2] = arcpoly{
            0,             0,  0,  1,
           r6,             0,  0,  1, 
r6*cos(pala2), r6*sin(pala2), r6,  8, 
r6*sin(pala2), r6*cos(pala2),  0,  1,
            0,            r6,  0,  1
}

/* ボルト穴 */
fig[10] = circle{23.4*cos(60),23.4*sin(60),pal6/2,8}

/* 積分路指定 */
gg = (r4-r5)/6
path[1] = circle{0,0,(r4+r5)/2,90,0,90}
path[2] = circle{0,0,r4-gg*1,90,0,90}
path[3] = circle{0,0,r4-gg*2,90,0,90}
path[4] = circle{0,0,r4-gg*4,90,0,90}

/* 材質指定 */
 mat[1] = {3000,3000,7e6}
 mat[2] = {3000,3000,7e6}
 mat[4] = {0.6364,0.6364,-1}
 mat[5] = {1,1,0}
 mat[6] = {1,1,0}
 mat[7] = {1,1,0}
 mat[8] = {1,1,0}
 mat[9] = {1,1,0}
mat[10] = {1,1,0}

coil = {5,0}
coil = {6,0}
coil = {7,0}