建筑热过程太阳辐射负荷计算

太阳辐射负荷计算

注释

'''
'''
1、赤纬角:地球中心与太阳中心的连线与赤道平面的夹角= 23.45*sin(360*(284*n)/365) n是一年中某一天对应365天中的天数 
2、真太阳时:H = H_s +(-) (L-L_s)/15 + e/60 H_s为该地区标准时间，L，L_s分别为当地的纬度和地区标准时间位置的经度，东半球为+，西半球为+
3、将真太阳时用角度表示，称为当地太阳时角 = (H_s +(-) (L-L_s)/15 + e/60 -12) * 15
4、太阳高度角h：sinh = sin(纬度)*sin(赤纬角) + cos(纬度)*cos(赤纬角)*cos(当地太阳时角)
5、太阳方位角a：太阳偏东为-，偏西为+，sina= cos(赤纬角)*sin(当地太阳时角)/cos(太阳高度角)，当计算出sina大于1，或sina的绝对值较小时，cosa = (sin(太阳高度角)*sin(纬度) - sin(赤纬角))/(cos(太阳高度角)cos(纬度))
6、太阳位置：太阳入射角i(太阳射线与壁面法向之间的夹角):cosi = cos(壁面与水平面的夹角)*sin(太阳高度角) + sin(壁面与水平面的夹角)*cos(太阳高度角)cos(壁面太阳方位角)
7、太阳常数：单位面积单位时间内来自太阳的辐射能量，用I_0表示，1353W/平方米
8、太阳光谱：I_0 等于各个光谱辐射强度之和
9、大气透明度：阳光经过大气层时，强度按指数规律衰减，为exp(-k*x)
10、太阳位于天顶时，到达地面的法向太阳直射辐射为I_l = I_0 * p
11、太阳不位于天顶是，阳光路径为l',到达地面上的法向太阳辐射强度应为：I_DN = I_0*(exp-k*l') = I_0 * exp(-k*m*l) = I_0*p的m次方，p的m次方为大气透明度，m = 1/sinh
12、太阳直射辐射强度计算：I_Dθ = I_DNcosi = I_0*p的m次方cosi,(i是太阳入射角)
13、天空散射辐射：I_dH = 0.5 * I_0 * sinh * (1-p的m次方) / (1-1.4lnp)
14、太阳总辐射I_Sθ = I_Dθ + I_dθ


'''
'''

引入模块

import datetime
from math import sin, cos, asin, acos, log, pi
from tkinter import E

某一天对应365天中的天数

list_31 = [1, 3, 5, 7, 8, 10, 12]
list_30 = [4, 6, 9, 11]
def number_day(mouth, day):
    num = 0
    for i in range(1, mouth):
        if i in list_31:
            num = num + 31
        elif i in list_30