# Organize the laminates into an array laminates_Lay3 = [lam1, lam2, lam3, lam4] # Create the cross-section object and mesh it xsect_Lay3 = XSect(4, af3, xdim3, laminates_Lay3, matLib, typeXSect="box", meshSize=4) # Run the cross-sectional analysis. Since this is an airfoil and for this, # symmetric airfoils the AC is at the 1/c chord, we will put the reference axis # here xsect_Lay3.xSectionAnalysis() # ref_ax=[0.25*c3,0.]) K_tmp = xsect_Lay3.K F_tmp = xsect_Lay3.F xs_tmp = xsect_Lay3.xs ys_tmp = xsect_Lay3.ys # Let's see what the rigid cross-section looks like: xsect_Lay3.plotRigid() # Print the stiffness matrix xsect_Lay3.printSummary(stiffMat=True) x1 = np.array([0, 0, 0]) x2 = np.array([0, 0, 160.5337]) # Initialize a superbeam ID SBID = 1 # Next we need to initialize the number of elements the superbeam should mesh: noe = 40 # Now let's make the superbeam sbeam1 = SuperBeam(SBID, x1, x2, xsect_Lay3, noe) model = Model() # Easy right? Now let's add the superbeam to the model. model.addElements([sbeam1]) # Now that our beam is loaded into the FEM, let's visualize it! model.plotRigidModel(numXSects=8)
laminates_Lay3,
matLib,
typeXSect='box',
meshSize=4)
# Run the cross-sectional analysis. Since this is an airfoil and for this,
# symmetric airfoils the AC is at the 1/c chord, we will put the reference axis
# here
xsect_Lay3.xSectionAnalysis() #ref_ax=[0.25*c3,0.])
K_tmp = xsect_Lay3.K
F_tmp = xsect_Lay3.F
xs_tmp = xsect_Lay3.xs
ys_tmp = xsect_Lay3.ys
# Let's see what the rigid cross-section looks like:
xsect_Lay3.plotRigid()
# Print the stiffness matrix
xsect_Lay3.printSummary(stiffMat=True)
x1 = np.array([0, 0, 0])
x2 = np.array([0, 0, 160.5337])
# Initialize a superbeam ID
SBID = 1
# Next we need to initialize the number of elements the superbeam should mesh:
noe = 40
# Now let's make the superbeam
sbeam1 = SuperBeam(SBID, x1, x2, xsect_Lay3, noe)
model = Model()
# Easy right? Now let's add the superbeam to the model.
model.addElements([sbeam1])
# Now that our beam is loaded into the FEM, let's visualize it!
model.plotRigidModel(numXSects=8)
