chapter_2/reliability_problem_new_policy.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "# Bearing Life(Hours)\n",
    "bearing_life = [1000, 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800, 1900]\n",
    "# Bearing Life(Hours) - Probability\n",
    "bearing_probabilities = [0.10, 0.14, 0.24, 0.14, 0.12, 0.10, 0.06, 0.05, 0.03, 0.02]\n",
    "\n",
    "# calculate cumulative probabilities\n",
    "cumulative_probabilities = [bearing_probabilities[0]]\n",
    "for i in range(1, len(bearing_probabilities)):\n",
    "    cumulative_probabilities.append(cumulative_probabilities[i-1] + bearing_probabilities[i])\n",
    "\n",
    "# generate random digit assignments\n",
    "bearing_random_digits = []\n",
    "count = 0\n",
    "for i in range(len(bearing_probabilities)):\n",
    "    bearing_random_digits.append([])\n",
    "    for j in range(int(bearing_probabilities[i]*100)): \n",
    "        bearing_random_digits[i].append(count)\n",
    "        count += 1\n",
    "        \n",
    "\n",
    "# for i in range(len(cumulative_probabilities)):\n",
    "#     print(f\"{cumulative_probabilities[i]:.2f}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "# delay Life(minutes)\n",
    "delay_life = [4 , 6 , 8]\n",
    "# delay Life(minutes) - Probability\n",
    "delay_probabilities = [0.3 , 0.6 , 0.1]\n",
    "\n",
    "# calculate cumulative probabilities\n",
    "cumulative_probabilities = [delay_probabilities[0]]\n",
    "for i in range(1, len(delay_probabilities)):\n",
    "    cumulative_probabilities.append(cumulative_probabilities[i-1] + delay_probabilities[i])\n",
    "\n",
    "# generate random digit assignments\n",
    "delay_random_digits = []\n",
    "count = 0\n",
    "for i in range(len(delay_probabilities)):\n",
    "    delay_random_digits.append([])\n",
    "    for j in range(int(delay_probabilities[i]*10)): \n",
    "        delay_random_digits[i].append(count)\n",
    "        count += 1\n",
    "\n",
    "\n",
    "\n",
    "# for i in range(len(cumulative_probabilities)):\n",
    "#     print(f\"{cumulative_probabilities[i]:.2f}\")\n",
    "\n",
    "# delay_random_digits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "import random\n",
    "\n",
    "# for bearing 1 , 2 , 3\n",
    "total_hours_of_operation = 30000\n",
    "total_clock_hour = 0\n",
    "total_delay_minutes = 0\n",
    "total_change = 0    # how many bearing have to change\n",
    "\n",
    "while total_hours_of_operation > total_clock_hour:\n",
    "    total_change += 1\n",
    "    bearing_random_num_1 = random.randint(0 , 99)\n",
    "    bearing_random_num_2 = random.randint(0 , 99)\n",
    "    bearing_random_num_3 = random.randint(0 , 99)\n",
    "    delay_random_num = random.randint(0 , 9)\n",
    "    current_clock = 0\n",
    "    current_delay = 0\n",
    "\n",
    "    # find bearing 1 clock from generated random num\n",
    "    for i in range(len(bearing_random_digits)):\n",
    "        if bearing_random_num_1 in bearing_random_digits[i]:\n",
    "            current_clock = bearing_life[i]\n",
    "            break\n",
    "    \n",
    "    # find bearing 2 clock from generated random num\n",
    "    for i in range(len(bearing_random_digits)):\n",
    "        if bearing_random_num_2 in bearing_random_digits[i]:\n",
    "            current_clock = min(bearing_life[i] , current_clock)\n",
    "            break\n",
    "\n",
    "    # find bearing 3 clock from generated random num\n",
    "    for i in range(len(bearing_random_digits)):\n",
    "        if bearing_random_num_3 in bearing_random_digits[i]:\n",
    "            current_clock = min(bearing_life[i] , current_clock)\n",
    "            break\n",
    "\n",
    "\n",
    "    total_clock_hour += current_clock\n",
    "\n",
    "\n",
    "    # find delay from generated random num \n",
    "    for i in range(len(delay_random_digits)):\n",
    "        if delay_random_num in delay_random_digits[i]:\n",
    "            current_delay = delay_life[j]\n",
    "            break\n",
    "\n",
    "\n",
    "    total_delay_minutes += current_delay\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total repairman delay time = 112 minutes\n",
      "Time spent in bearing replacement = 1120 minutes\n",
      "Total downtime = 1232 minutes\n",
      "Cost of downtime = Rs. 6160\n",
      "Cost of repairman = Rs. 466.6666666666667\n",
      "\n",
      "\n",
      "Total bearing replaced = 84\n",
      "Cost of bearing = Rs. 1680\n",
      "\n",
      "\n",
      "Total estimated cost = Rs. 8306.666666666666\n"
     ]
    }
   ],
   "source": [
    "# 20 minutes to change one bearing\n",
    "# Downtime of the machining center costs Rs. 5 per minute\n",
    "# direct on job cost of the repairman is Rs. 25 per hour \n",
    "# cost of a bearing is Rs. 20.\n",
    "\n",
    "three_bearing_change_time_minutes = 40\n",
    "cost_of_downtime_per_minutes = 5\n",
    "cost_of_repairman_per_hour = 25\n",
    "cost_of_one_bearing = 20\n",
    "\n",
    "print(\"Total repairman delay time =\" , total_delay_minutes , \"minutes\")\n",
    "print(\"Time spent in bearing replacement =\", total_change*three_bearing_change_time_minutes , \"minutes\")\n",
    "print(\"Total downtime =\" , total_delay_minutes+(total_change*three_bearing_change_time_minutes) , \"minutes\")\n",
    "\n",
    "cost_of_downtime = cost_of_downtime_per_minutes*(total_delay_minutes+(total_change*three_bearing_change_time_minutes))\n",
    "print(\"Cost of downtime = Rs.\",  cost_of_downtime)\n",
    "cost_of_repairman = total_change*three_bearing_change_time_minutes*cost_of_repairman_per_hour/60\n",
    "print(\"Cost of repairman = Rs.\",  cost_of_repairman)\n",
    "print(\"\\n\")\n",
    "print(\"Total bearing replaced =\", total_change*3)\n",
    "cost_of_bearing = cost_of_one_bearing*total_change*3\n",
    "print(\"Cost of bearing = Rs.\",  cost_of_bearing)\n",
    "print(\"\\n\")\n",
    "\n",
    "\n",
    "print(\"Total estimated cost = Rs.\",  cost_of_repairman+cost_of_bearing+cost_of_downtime)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.0"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}