!  cc-pVDZ  EMSL  Basis Set Exchange Library   7/19/15 12:13 AM! Elements
                             References
! --------                             ----------
! H     : T.H. Dunning, Jr. J. Chem. Phys. 90, 1007 (1989).
! He    : D.E. Woon and T.H. Dunning, Jr. J. Chem. Phys. 100, 2975 (1994).
! Li - Ne: T.H. Dunning, Jr. J. Chem. Phys. 90, 1007 (1989).
! Na - Mg: D.E. Woon and T.H. Dunning, Jr.  (to be published)
! Al - Ar: D.E. Woon and T.H. Dunning, Jr.  J. Chem. Phys. 98, 1358 (1993).
! Sc - Zn: N.B. Balabanov and K.A. Peterson, J. Chem. Phys. 123, 064107 (2005),
! N.B. Balabanov and K.A. Peterson, J. Chem. Phys. 125, 074110 (2006)
! Ca     : J. Koput and K.A. Peterson, J. Phys. Chem. A, 106, 9595 (2002).
! $DATA
HYDROGEN
S   3
  1     13.0100000              0.0196850        
  2      1.9620000              0.1379770        
  3      0.4446000              0.4781480        
S   1
  1      0.1220000              1.0000000        
P   1
  1      0.7270000              1.0000000    

$END

$DATAと$ENDを削除、HYDROGENをH cc-pVDZのように元素名＋基底関数名(間はスペースでないとだめみたい)にする。

ちなみにhttp://classic.chem.msu.su/gran/gamess/pvxz.zipからダウンロードできる。

後は 
$BASIS EXTFIL=.T. GBASIS=CC-PVDZ $END
を追加するだけ。
 

 $CONTRL RUNTYP=rsurface NZVAR=18 COORD=ZMT mplevl=2 $END
!NZVARでIZMATの配列数（行の数）を指定する。通常は3N-6、直線分子で3N-5を指定する。

 $SYSTEM TIMLIM=500 MEMORY=30000000 $END $BASIS  GBASIS=tzv ndfunc=1 $END! To freeze coord # 12 during constrained optimization $statpt opttol=1d-5 nstep=1000 IFREEZ(1)=12 $end! To scan along coord # 12 $surf ndisp1=33 disp1=3.75 vect1(12)=1 orig1=0 $end
!vect1(i)で指定される二面角（ZMATのi番目の行）についてポテンシャル障壁を求める。今回は3.75°ずつ33点を計算する。
 $DATACH3-CH3C1CC    1 r1h    1 r2    2 a1h    1 r2    2 a1    3 d1h    1 r2    2 a1    3 d2h    2 r2    1 a1    4 dxh    2 r2    1 a1    6 d1h    2 r2    1 a1    6 d2r1=1.54r2=1.09a1=109d1=120d2=-120dx=0 $END $ZMAT   IZMAT(1)=
!結合と最適化する結合角を指定する。(code, I, J, K, L, M, N)
!code = 1: I,Jが結合する。
!code = 2: 角I,J,Kの角度
!code = 3: 角I,J,Kと角J,K,Lの二面角
         1,   2,   1,
         1,   3,   1,
         2,   3,   1,   2,
         1,   4,   1,
         2,   4,   1,   2,
         3,   4,   1,   2,   3,
         1,   5,   1,
         2,   5,   1,   2,
         3,   5,   1,   2,   3,
         1,   6,   2,
         2,   6,   2,   1,
         3,   6,   2,   1,   4,
         1,   7,   2,
         2,   7,   2,   1,
         3,   7,   2,   1,   6,
         1,   8,   2,
         2,   8,   2,   1,
         3,   8,   2,   1,   6,
 $END

２．　ツール→リレーションをクリック

３．　リーレーションウィンドウに次の式を打ち込む

Z=D/M
ha=m
hf=1.16*m
Ro=(D/2)+ha
Rd=(D/2)-hf
Ri=(D/2)*cos(PHI)
g_c=sqrt((D^2)/4-Ri^2)
theta_0=(360/(4*Z))-((g_c/Ri)*(180/pi))+atan(g_c/Ri)

(piが変数名として予約されているっぽい)

D=200(歯車のピッチ円直径)

M=10(モージュール数)

PHI=20(圧力角)

B=15(ヘリカルの角)

を初期値として入れる。

４．　適当な平面に適当な大きさの円を描く

５．　ツール→リレーションをクリック

６．　リーレーションウィンドウに次の式を打ち込む

sd0=Ro*2(sd0は適当に描いた円の直径)

７．　スケッチウィンドウを閉じる

８．　円を押し出す。(押し出し距離は100)

９．　データム→カーブで任意のカーブを描く

１０．リボンの右端にデータアイコンがあり、その下にオプションが出せるので、オプションにあるカーブ→方程式利用カーブをクリック

１１．円柱座標系を選ぶ

１２．方程式をクリック、次の式を打ち込む

gamma = sqrt(Ro^2-Ri^2)*t
r = sqrt(gamma^2+Ri^2)
theta = theta_0+((gamma/Ri)*(360/(2*pi)))-atan(gamma/Ri)
z = 0

(r,z,thetaが座標軸名として予約されている？、tは媒介変数？)

１３．下限、上限はここの媒介変数の範囲を示す？

１４．カーブ：方程式使用を終わる

１５．歯車の中心軸を描く

１６．描いた軸を選択し、スイープをクリック

１７．サーフェイスとしてスイープを選択

１８．参照→軌道の下にある詳細アイコンをクリックして終点1を方程式利用カーブが交差している面上にする。

１９．スケッチを始める

２０．中心を円周を結ぶ直線を適当に引く

２１．ツール→リレーションをクリック

２２．　リーレーションウィンドウに次の式を打ち込む

sd3=90+trajpar*1*B(sd3は基準面と直線のなす角、trajparは0から1まで変化する変数らしい)

２３．スケッチ、スイープを終了する。

２４．円周上にできた線を選択し、スイープをクリック

２５．参照→断面平面の制御：一定垂直方向、方向参照：FRONT:F3(データム平面)、水平/垂直制御：サーフェスに垂直を選ぶ

２６．スケッチを始める

２７．適当な図形を描く(この操作がないと投影が選べない)

２８．適当な大きさの円を作図モードで描く

２９．ツール→リレーションをクリック

３０．　リーレーションウィンドウに次の式を打ち込む

sd17=Rd(sd17は適当に描いた円の半径)

３１．１４．で描いたカーブをクリックして投影をクリック

３２．投影されたカーブの中心側の端点と円の中心を結ぶ直線を作図モードで引く

３３．３０．で描いた円とと３２．の直線の交点を打つ

３４．３３．の交点と中心側の端点を直線で結ぶ

３５．３４．の線分と３１．のカーブを選択してミラーをクリック、その側にある基準線をクリックする

３６．３１．のカーブと３５．のカーブの中心側と円周側の点をそれぞれ円の中心とする弧で結ぶ

３７．スケッチを終了し、材料の除去を選択した後スイープを終了する

３８．３７．のスイープを選択し、パターンをクリックする

３９．軸を選択、歯の数を入力、パターンの角度範囲を選択する

４０．歯の数だけ点が表示されたらパターンを終了する

ヘリカルギアがモデリングできているはず