[20.6e6, 1.42e6, .34, .42, .87e6, 0.], .005)

# Box Configuration 2
c2 = 0.53
xdim2 = [-0.8990566037735849, 0.8990566037735849]
af2 = Airfoil(c2, name='box')

# B1 Box beam (0.5 x 0.923 in^2 box with laminate schedule [15]_6)
n_i_B1 = [6]
m_i_B1 = [2]
th_B1 = [-15]
lam1_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam2_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam3_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam4_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam1_B1.printPlies()
laminates_B1 = [lam1_B1, lam2_B1, lam3_B1, lam4_B1]
xsect_B1 = XSect(af2, xdim2, laminates_B1, matLib, typeXsect='box', meshSize=2)

import time
t1 = time.time()
# With lu factorization
xsect_B1.xSectionAnalysis()
xsect_B1.printStiffMat()
t2 = time.time()
# With linalg.solve
xsect_B1.xSectionAnalysis(linalgSolve=True)
xsect_B1.printStiffMat()
t3 = time.time()
# With cholesky
xsect_B1.xSectionAnalysis(sparse=True)
m_i_1 = [1,1,1,1,1,1]
th_1 = [-15,-15,-15,-15,-15,-15]
lam1 = Laminate(n_i_1, m_i_1, matLib, th=th_1)
n_i_2 = [1,1,1,1,1,1]
m_i_2 = [1,1,1,1,1,1]
th_2 = [-15,15,-15,15,-15,15]
lam2 = Laminate(n_i_2, m_i_2, matLib, th=th_2)
n_i_3 = [1,1,1,1,1,1]
m_i_3 = [1,1,1,1,1,1]
th_3 = [15,15,15,15,15,15]
lam3 = Laminate(n_i_3, m_i_3, matLib, th=th_3)
n_i_4 = [1,1,1,1,1,1]
m_i_4 = [1,1,1,1,1,1]
th_4 = [-15,15,-15,15,-15,15]
lam4 = Laminate(n_i_4, m_i_4, matLib, th=th_4)
lam1.printPlies()
lam2.printPlies()
lam3.printPlies()
lam4.printPlies()
laminates_Lay3 = [lam1,lam2,lam3,lam4]

xsectDict = {}

for i in range(len(cs)):
    xsectDict['xsect'+str(i)] = XSect(Airfoil(cs[i],name='box'),xdim,laminates_Lay3,matLib,typeXsect='box',meshSize=4)
    xsectDict['xsect'+str(i)].xSectionAnalysis()
    #xsectDict['xsect'+str(i)].strn2dspl(strn,figName='xsect'+str(i),contour_Total_T=True)
'''
xsect1 = XSect(af1,xdim,laminates_Lay3,matLib,typeXsect='box',meshSize=4)
xsect2 = XSect(af2,xdim,laminates_Lay3,matLib,typeXsect='box',meshSize=4)
xsect3 = XSect(af3,xdim,laminates_Lay3,matLib,typeXsect='box',meshSize=4)
# Box Configuration 2
c2 = 0.53
xdim2 = [-0.8990566037735849,0.8990566037735849]
af2 = Airfoil(c2,name='box')



# B1 Box beam (0.5 x 0.923 in^2 box with laminate schedule [15]_6)
n_i_B1 = [6]
m_i_B1 = [2]
th_B1 = [-15]
lam1_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam2_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam3_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam4_B1 = Laminate(n_i_B1, m_i_B1, matLib, th=th_B1)
lam1_B1.printPlies()
laminates_B1 = [lam1_B1,lam2_B1,lam3_B1,lam4_B1]
xsect_B1 = XSect(af2,xdim2,laminates_B1,matLib,typeXsect='box',meshSize=2)

import time
t1 = time.time()
# With lu factorization
xsect_B1.xSectionAnalysis()
xsect_B1.printStiffMat()
t2 = time.time()
# With linalg.solve
xsect_B1.xSectionAnalysis(linalgSolve=True)
xsect_B1.printStiffMat()
t3 = time.time()
# With cholesky
xsect_B1.xSectionAnalysis(sparse=True)
m_i_1 = [1, 1, 1, 1, 1, 1]
th_1 = [-15, -15, -15, -15, -15, -15]
lam1 = Laminate(n_i_1, m_i_1, matLib, th=th_1)
n_i_2 = [1, 1, 1, 1, 1, 1]
m_i_2 = [1, 1, 1, 1, 1, 1]
th_2 = [-15, 15, -15, 15, -15, 15]
lam2 = Laminate(n_i_2, m_i_2, matLib, th=th_2)
n_i_3 = [1, 1, 1, 1, 1, 1]
m_i_3 = [1, 1, 1, 1, 1, 1]
th_3 = [15, 15, 15, 15, 15, 15]
lam3 = Laminate(n_i_3, m_i_3, matLib, th=th_3)
n_i_4 = [1, 1, 1, 1, 1, 1]
m_i_4 = [1, 1, 1, 1, 1, 1]
th_4 = [-15, 15, -15, 15, -15, 15]
lam4 = Laminate(n_i_4, m_i_4, matLib, th=th_4)
lam1.printPlies()
lam2.printPlies()
lam3.printPlies()
lam4.printPlies()
laminates_Lay3 = [lam1, lam2, lam3, lam4]

xsectDict = {}

for i in range(len(cs)):
    xsectDict['xsect' + str(i)] = XSect(Airfoil(cs[i], name='box'),
                                        xdim,
                                        laminates_Lay3,
                                        matLib,
                                        typeXsect='box',
                                        meshSize=4)
    xsectDict['xsect' + str(i)].xSectionAnalysis()
Esempio n. 5
0
xsect_B1 = XSect(af2,xdim2,laminates_B1,matLib,typeXsect='box',meshSize=2)
xsect_B1.xSectionAnalysis()
xsect_B1.printStiffMat()
strn = np.array([0.,0.,0.,0.,0.,1.0])
xsect_B1.strn2dspl(strn,figName='Validation Case B1',contour_Total_T=True)
'''

# Layup 1 Box beam (0.5 x 0.923 in^2 box with laminate schedule [0]_6)
n_i_Lay1 = [6]
m_i_Lay1 = [2]
th_Lay1 = [0]
lam1_Lay1 = Laminate(n_i_Lay1, m_i_Lay1, matLib, th=th_Lay1)
lam2_Lay1 = Laminate(n_i_Lay1, m_i_Lay1, matLib, th=th_Lay1)
lam3_Lay1 = Laminate(n_i_Lay1, m_i_Lay1, matLib, th=th_Lay1)
lam4_Lay1 = Laminate(n_i_Lay1, m_i_Lay1, matLib, th=th_Lay1)
lam1_Lay1.printPlies()
laminates_Lay1 = [lam1_Lay1,lam2_Lay1,lam3_Lay1,lam4_Lay1]
xsect_Lay1 = XSect(af2,xdim2,laminates_Lay1,matLib,typeXsect='box',meshSize=2)
xsect_Lay1.xSectionAnalysis()
xsect_Lay1.printStiffMat()
#xsect_Lay1.strn2dspl(strn,figName='Validation Case Layup 1',contour_Total_T=True)

xsect_Lay1_off = XSect(af2,xdim2_off,laminates_Lay1,matLib,typeXsect='box',meshSize=2)
xsect_Lay1_off.xSectionAnalysis()

Ktmp = xsect_Lay1_off.K

print('\n\nThe cross-section stiffness matrix is:')
print(tabulate(np.around(Ktmp,decimals=2),tablefmt="fancy_grid"))

xsc,ysc = xsect_Lay1_off.genXYSC()