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) xsect_B1.printStiffMat() t4 = time.time() print('LU Factor Time: %5.5f seconds' %(t2-t1)) print('Solve Time: %5.5f seconds' %(t3-t2)) print('Sparse Time: %5.5f seconds' %(t4-t3)) #xsect_B1.printStiffMat()
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) xsect_B1.printStiffMat() t4 = time.time() print('LU Factor Time: %5.5f seconds' % (t2 - t1)) print('Solve Time: %5.5f seconds' % (t3 - t2)) print('Sparse Time: %5.5f seconds' % (t4 - t3)) #xsect_B1.printStiffMat()
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() xsc = -xsc ysc = -ysc T1 = np.array([[1.,0.,0.,0.,0.,-ysc],[0.,1.,0.,0.,0.,xsc],[0.,0.,1.,ysc,-xsc,0.],[0.,0.,0.,1.,0.,0.],[0.,0.,0.,0.,1.,0.],[0.,0.,0.,0.,0.,1.]]) T2 = np.array([[1.,0.,0.,0.,0.,0.],[0.,1.,0.,0.,0.,0.],[0.,0.,1.,0.,0.,0.],[0.,0.,-ysc,1.,0.,0.],[0.,0.,xsc,0.,1.,0.],[ysc,-xsc,0.,0.,0.,1.]])