Pythonでボールねじとモーターとカップリングの選定をしてみる 改良編

Share on:

以前作ったボールねじ、モーター、カップリング選定プログラムを改良しました。

改良点

  • 計算式ミスを修正
  • 安全率を指定して計算結果の出力を絞り込む

座屈荷重や許容荷重などの軸方向荷重に関する計算式がミスっていたので修正しました。前回はそれぞれの計算を関数化して計算していましたが、わざわざ関数に分ける必要もないのでcalculate()という一つの関数で計算するようにしました。その他には、ストロークを指定してストロークを満たしたボールねじのみ抽出して計算して、計算結果から安全率以下を省き、処理は20000通り毎に分割するように変更したことで、処理を限定的に行うようにしました。

ストロークの限定

1bs = bs[( (bs['l0min'] < stroke) & (bs['l0max'] > stroke) )]
2bs['l1'] = bs['l7'] + bs['l8'] + stroke
3bs['l0'] = bs['l1'] + bs['l2'] + bs['l3'] + bs['l4'] + bs['l5']

計算結果を安全率で絞込み

計算結果と計算条件をすべてdfに入れて、計算した安全率sf_sl~sf_mt_d0と比較すると一致した行がTrue、一致しなければFalseになる配列が得られるので、df[]でそのブーリアンの配列を指定すると、安全率で絞り込まれた計算結果が得られます。

 1flags  = (df['sf_sl'] >= sf)
 2flags &= (df['sf_bl'] >= sf)
 3flags &= (df['sf_yl'] >= sf)
 4flags &= (df['sf_cs'] >= sf)
 5flags &= (df['sf_ir'] >= sf)
 6flags &= (df['sf_t'] >= sf)
 7flags &= (df['sf_mt_n'] >= sf)
 8flags &= (df['sf_bs_d3'] >= sf)
 9flags &= (df['sf_mt_d0'] >= sf)
10df = df[flags]

計算結果

ボールねじのパラメータ、カップリングのパラメータ、モーターのパラメータ、安全率というcsvです。

所感

計算関連の部分はこのくらいで良さげです。次は簡単に操作できるようにGUIを作ってみたいと思いますが、その前にPyQt5のQTableViewでPandasのDataFrameを表示するプログラムを作成したいと思います。

ソースコード

  1# -*- coding: utf-8 -*-
  2import itertools
  3import math
  4import numpy as np
  5import pandas as pd
  6def calculate(_bs, _cp, _mt, sf, stroke, velocity, acceleration, buckling_load_factor, static_allowable_load_factor, 
  7        load_factor, critical_speed_factor, work_mass, guide_friction_factor, guide_friction):
  8    '''
  9    bs : ボールねじリスト(pandasデータフレーム)
 10    cp : カップリングリスト(pandasデータフレーム)
 11    mt : モーターリスト(pandasデータフレーム)
 12    sf : 安全率
 13    velocity : 速度[m/s]
 14    acceleration : 加速度[m/ss]
 15    buckling_load_factor : 座屈荷重係数 19.9, 10.0, 1.2, 5.0
 16    static_allowable_load_factor : 静的許容荷重係数 1.0~2.0, 1.5~3.0
 17    load_factor : 荷重係数 1.0~1.2, 1.2~1.5, 1.5~3.0
 18    critical_speed_factor : 危険速度係数 21.9, 15.1, 9.7, 3.4
 19    work_mass : ワーク質量[kg]
 20    guide_friction_factor : 案内面摩擦係数
 21    guide_friction : 案内面抵抗[N]
 22    '''
 23    # データフレームをコピー
 24    bs = _bs.copy()
 25    bs.insert(0, 'model', bs.index)
 26    cp = _cp.copy()
 27    cp.insert(0, 'model', cp.index)
 28    mt = _mt.copy()
 29    mt.insert(0, 'model', mt.index)
 30    # ストロークを満たす行を抽出
 31    bs = bs[( (bs['l0min'] < stroke) & (bs['l0max'] > stroke) )]
 32    bs['l1'] = bs['l7'] + bs['l8'] + stroke
 33    bs['l0'] = bs['l1'] + bs['l2'] + bs['l3'] + bs['l4'] + bs['l5']
 34    # 全通りの組み合わせリストを作成
 35    df = pd.DataFrame(itertools.product(bs.values, cp.values, mt.values))
 36    df = pd.concat([
 37        pd.DataFrame(df[0].tolist(), columns= 'bs_' + bs.columns), 
 38        pd.DataFrame(df[1].tolist(), columns= 'cp_' + cp.columns), 
 39        pd.DataFrame(df[2].tolist(), columns= 'mt_' + mt.columns)
 40    ], axis=1)
 41    
 42    df['v']  = velocity  # 速度[m/s]
 43    df['a']  = acceleration  # 加速度[m/ss]
 44    df['m']  = buckling_load_factor  # 座屈荷重係数 19.9, 10.0, 1.2, 5.0
 45    df['fs'] = static_allowable_load_factor # 静許容荷重係数 1.0~2.0, 1.5~3.0
 46    df['fw'] = load_factor # 荷重係数 1.0~1.2, 1.2~1.5, 1.5~3.0
 47    df['f']  = critical_speed_factor  # 危険速度係数 21.9, 15.1, 9.7, 3.4
 48    df['mw'] = work_mass # ワーク質量[kg]
 49    # ボールねじ計算
 50    df['n']  = df['v'] * 60000.0 / df['bs_l6'] # 回転数
 51    df['at'] = df['v'] / df['a'] # 加速時間 [sec]
 52    df['Fa'] = -1.0 * guide_friction_factor * df['mw'] * 9.80665 + guide_friction + df['mw'] * df['a'] # 軸方向最大荷重 [N] -myu * mg + ff + ma
 53    df['P1'] = df['m'] * np.power(df['bs_d1'], 4.0) * 10000.0 / np.power(( df['bs_l3'] / 2.0 + df['bs_l2'] + df['bs_l1'] - df['bs_l7'] / 2.0 ),2.0) # 座屈荷重
 54    df['P2'] = 115.45353 * np.power(df['bs_d1'], 2.0) # 降伏荷重
 55    df['P0'] = df['bs_C0a'] / df['fs'] # 静最大許容荷重
 56    df['nc'] = df['f'] * df['bs_d1'] * 10000000.0 / np.power(( df['bs_l5'] + df['bs_l1'] - df['bs_l8'] - df['bs_l7']/2.0 ),2.0) # 危険速度
 57    df['L']  = np.power( ( df['bs_Ca']/(df['Fa']*df['fw']) ) , 3.0) * 1000000.0 # 寿命回転数
 58    df['Lt'] = df['L'] / 60.0 / df['n'] # 寿命時間
 59    df['Ls'] = df['L'] * df['bs_l6'] * 0.000001 # 寿命距離
 60    # モーター計算
 61    df['ar'] = 2.0 * math.pi * df['n'] / 60.0 / df['at'] # 角加速度 [rad/ss]
 62    df['Iw'] = df['mw'] * ( df['bs_l6']*0.001 / math.pi )**2.0 /4.0 # ワークのイナーシャ m*(l/pi)^2/4
 63    df['Is'] = math.pi * 0.24375 * df['bs_l0'] * df['bs_d0']**4.0 * 0.000000000001 # ボールねじ軸のイナーシャ [kg*m**2]
 64    df['Ir'] = ( df['Is'] + df['Iw'] + df['cp_i'] ) / df['mt_i'] # イナーシャ比 = モーター以外のイナーシャ合計 / モーターのイナーシャ
 65    df['t']  = ( df['Is'] + df['Iw'] + df['cp_i'] + df['mt_i'] ) * df['ar'] # トルク = イナーシャ合計 * 角加速度
 66    
 67    # 安全率計算
 68    df['sf_sl'] = df['P0'] / df['Fa'] # 静最大許容荷重
 69    df['sf_bl'] = df['P1'] / df['Fa'] # 座屈荷重
 70    df['sf_yl'] = df['P2'] / df['Fa'] # 降伏荷重
 71    df['sf_cs'] = df['nc'] / df['n'] # 危険速度
 72    df['sf_ir'] = df['mt_ir'] / df['Ir'] # イナーシャ比
 73    df['sf_t']  = df['mt_t0'] / df['t'] # トルク
 74    df['sf_mt_n']  = df['mt_nmax'] / df['n'] # モーター回転数
 75    df['sf_bs_d3'] = 999.0 # ボールねじシャンク径
 76    df['sf_bs_d3'].where( (df['cp_d1'] <= df['bs_d3']) & (df['bs_d3'] <= df['cp_d2']), 0.0, inplace=True)
 77    df['sf_mt_d0'] = 999.0 # モーター軸径
 78    df['sf_mt_d0'].where( (df['cp_d1'] <= df['mt_d0']) & (df['mt_d0'] <= df['cp_d2']), 0.0, inplace=True)
 79    # 安全率で絞り込み
 80    flags  = (df['sf_sl'] >= sf)
 81    flags &= (df['sf_bl'] >= sf)
 82    flags &= (df['sf_yl'] >= sf)
 83    flags &= (df['sf_cs'] >= sf)
 84    flags &= (df['sf_ir'] >= sf)
 85    flags &= (df['sf_t'] >= sf)
 86    flags &= (df['sf_mt_n'] >= sf)
 87    flags &= (df['sf_bs_d3'] >= sf)
 88    flags &= (df['sf_mt_d0'] >= sf)
 89    return df[flags]
 90def main():
 91    # csv読み込み
 92    bs = pd.read_csv('ballscrew.csv',  encoding='utf-8', index_col=0).astype(float)
 93    cp = pd.read_csv('coupling.csv',  encoding='utf-8', index_col=0).astype(float)
 94    mt = pd.read_csv('motor.csv',  encoding='utf-8', index_col=0).astype(float)
 95    df = calculate(
 96            bs, # ボールねじリスト
 97            cp, # カップリングリスト
 98            mt, # モーターリスト
 99           1.0, # 安全率
100         600.0, # ストローク[mm]
101           1.0, # 速度[m/s]
102          10.0, # 加速度[m/ss]
103          19.9, # 座屈荷重係数 19.9, 10.0, 1.2, 5.0
104           2.0, # 静許容荷重係数 1.0~2.0, 1.5~3.0
105           1.5, # 荷重係数 1.0~1.2, 1.2~1.5, 1.5~3.0
106          15.1, # 危険速度係数 21.9, 15.1, 9.7, 3.4
107          20.0, # ワーク質量[kg]
108         0.003, # 案内面摩擦係数
109          15.0  # 案内面の抵抗[N]
110    )
111    df.to_csv("output.csv")
112if __name__ == '__main__':
113    main()

関連記事