根据规范大批量计算轴压稳定系数

有时候需大批量计算轴压稳定系数，若采用查表的方式则严重影响工作效率，且附表未给出长细比过大时的稳定系数；若仅采用公式计算，则需要进行一定量的复核，防止公式输入错误。为便于数据操作及进行复核对比，特编写python脚本对弱硬化铝合金轴心受力构件的稳定系数进行计算。脚本共采用两种计算方法，一是根据附表进行插值计算，二是根据理论公式进行计算。两种之间可以方便地对比，以验证稳定系数计算的准确性，且可改写脚本以适用于其他稳定系数的计算。

以下为脚本：

# -*- coding: utf-8 -*-
#############脚本用于计算弱硬化铝合金稳定系数#############################
################可进行修改，以适用于其他材料##############################
######################本脚本为羽落之声编写################################
#注：代码糅合了多个代码，故引入多个包，可删减
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import xlrd
import xlwt

#获取换算长细比
#读取excel表格
wb = xlrd.open_workbook('data1.xlsx')
sheet = wb.sheet_by_name('Sheet1')
lamda = []
#写入数组lamda
for a in range(sheet.nrows):
 cells = sheet.row_values(a)  #每行数据赋值给cells
 a = cells[0]                 #因为表内可能存在多列数据，0代表第一列数据，1代表第二列，以此类推
 lamda.append(a)              #把每次循环读取的数据插入到list

#读取规范数据，以便放入数组
wb1 = xlrd.open_workbook('data2.xlsx')
sheet1 = wb1.sheet_by_name('Sheet1')

#写入数组fa
rows = sheet1.nrows                        # 获得行数
cols = sheet1.ncols                        # 获得列数
fa = [[] for i in range(rows-1)] # 定义了一个二维列表，行数等于sheet的行数
for i in range(1,rows):                      # 去掉表头，从第二行开始读
 fa[i-1] = sheet1.row_values(i)[1:cols]     # 给二维列表赋值

lamda1 = []
#写入数组lamda1，规范竖向数值，用于判断计算值位置
for a in range(1,sheet1.nrows):
 cells = sheet1.row_values(a)  #每行数据赋值给cells
 a = cells[0]                 #因为表内可能存在多列数据，0代表第一列数据，1代表第二列，以此类推
 lamda1.append(a)              #把每次循环读取的数据插入到list

lamda2 = []
#写入数组lamda1，规范竖向数值，用于判断计算值位置
for a in range(1,sheet1.ncols):
 cells = sheet1.col_values(a)  #每行数据赋值给cells
 a = cells[0]                 #因为表内可能存在多列数据，0代表第一列数据，1代表第二列，以此类推
 lamda2.append(a)              #把每次循环读取的数据插入到list

完整代码见收费部分。