HIRENASD Beam FEM

This contains attempts to create BEAM FEM model.

I have started a Blog to discuss this... please put your comments there and I will attempt to keep everything updated. Correspondence between A.B. and CDW is provided to assist in developing NASTRAN Beam FEM
(through July 1)(PDF)

The original data provided by Alexander is shown in model_input-extract.txt

Additional files loaded July 21, 2011:

  • Comparison of Wing Masses, and Lengths, and Mass calculated from provided densities, and lengths
  • Comparison of Model Masses and frequencies
  • Emails dealing with offsets for bending and center of mass (between Proctor and Wieseman)
  • Spreadsheet that contains element beam properties used for debugging
  • One element beam model to check beam stiffnesses
  • Email correspondence as of July 28, 2011 - includes a cursory check of the bending and center of mass locations.

    X1,X2,X3 Inertias/Kappa Mode 1 Freq, Hz Mode 2 Freq, Hz File name (ms103_beam_) 1,0,0 modified* 26.3785483.69360 100_modi_ccs_o1_k_nsm1_nm1 1,0,0 modified*/Beam properties modified 26.3785583.69344 100_modi_ccs_o1_k2_nsm1_nm1 cosa,0,-sina as provided 26.3794983.69533 ca_0_nsa_ccs_o1_k_nsm1_nm1 cosa,0,-sina as provided 26.3794983.69533 ca_0_nsa_ccs_o1_k_nsm1_nm1

    Below are frequencies with different permutations of

  • X1,X2,X3 definitions on CBEAM card
    A) - 1,0,0 or
    B) ca_0_nsa = cos(alpha),0,-sin(alpha) or
    C)nca_0_sa = -cos(alpha),0,sin(alpha) on the beam card
  • coordinate system (CCS) is converted from original
  • O1 indicates that offsets W1-W3 were included between the node and the shear center
  • nsm1 = beam density is unchanged and includes the non-structural mass
  • nsm2 = non-structural mass is separated from the beam mass density
  • nm = calculation of offset N1,N2 for beginning and end of each beam segment is performed - different ways depending on the index
  • modi = inertias provided by Boucke were modified by the alpha.
  • _k_ = kappa(K1,K2) were included in the model - values provided by Boucke (A.B.)
  • _k2_ = kappa(K1/K2) modified along with beam properties per equations provided by A.B.

    Description 1st mode2nd mode Boucke provided frequencies25.3387954778.63642299 100_modi_ccs_o1_k2_nsm1_nm126.3785583.69344 100_modi_ccs_o1_k_nsm1_nm1 26.3785483.6936 ca_0_nsa_ccs_o1_k_nsm1_nm1 26.3794983.69533 ca_0_nsa_modi_ccs_o1_k2_nsm1_nm126.379583.69762 ca_0_nsa_modi_ccs_o1_k_nsm1_nm126.3794983.69778 nca_0_sa_ccs_o1_k 26.3746683.62444 nca_0_sa_ccs_o1_k_nsm1_nm126.3395582.98219 nca_0_sa_ccs_o1_k_nsm1_nm226.377883.67184 nca_0_sa_ccs_o1_k_nsm2_nm126.345283.07531 nca_0_sa_ccs_o1_k_nsm2_nm226.3778583.6722 nca_0_sa_modi_ccs_o1_k_nsm1_nm1 26.3395382.98429

    Matlab code has been requested that reads the input deck. This code also does some debugging and calculation. The disclaimer is that it isn't easy to read. If you find a coding error please let me know. There are two codes - one for reading that also does some distance calculations and then a code that writes the BDF file.

    • Data and Resources

    Field Value
    Groups
    • AmeriGEOSS
    • National Provider
    • North America
    Tags
    • amerigeo
    • amerigeoss
    • ckan
    • geo
    • geoss
    • national
    • north-america
    • united-states
    isopen False
    license_id us-pd
    license_title us-pd
    maintainer CAROL WIESEMAN
    maintainer_email carol.d.wieseman@nasa.gov
    metadata_created 2025-12-02T06:10:17.282057
    metadata_modified 2025-12-02T06:10:17.282062
    notes This contains attempts to create BEAM FEM model. I have started a <a href="https://c3.nasa.gov/dashlink/projects/47/blog/">Blog</a> to discuss this... please put your comments there and I will attempt to keep everything updated. <a href="https://c3.nasa.gov/dashlink/static/media/dataset/Emails_AB_CDW.pdf">Correspondence between A.B. and CDW is provided to assist in developing NASTRAN Beam FEM <br> (through July 1)(PDF)</a> The original data provided by Alexander is shown in <a href="https://c3.nasa.gov/dashlink/static/media/dataset/model_input-extract.txt">model_input-extract.txt </a> Additional files loaded July 21, 2011:<br> <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/WingMassLength_Comparisons.xls">Comparison of Wing Masses, and Lengths, and Mass calculated from provided densities, and lengths</a> <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/Compare_W_eig_20110720.xlsx.pdf">Comparison of Model Masses and frequencies</a> <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/LP_CDW_emails.pdf">Emails dealing with offsets for bending and center of mass (between Proctor and Wieseman)</a> <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/BeamProperties_modrho.xlsx">Spreadsheet that contains element beam properties</a> used for debugging <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/OneElement_TipLoad.bdf">One element beam model to check beam stiffnesses</a> <LI><a href="https://c3.nasa.gov/dashlink/static/media/dataset/Emails_20110728.pdf"> Email correspondence as of July 28, 2011</a> - includes a cursory check of the bending and center of mass locations. <P> <TABLE BORDER = 1> <TR> <TD>X1,X2,X3</TD> <TD>Inertias/Kappa</TD> <TD>Mode 1 Freq, Hz</TD> <TD>Mode 2 Freq, Hz</TD> <TD>File name (ms103_beam_)</TD> </TR><TR> <TD>1,0,0</TD> <TD>modified*</TD> <TD>26.37854</TD><TD>83.69360</TD> <TD><a href="https://c3.nasa.gov/dashlink/static/media/dataset/ms103_beam_100_modi_ccs_o1_k_nsm1_nm1.bdf"> 100_modi_ccs_o1_k_nsm1_nm1</a> </TD></TR><TR> <TD>1,0,0</TD> <TD>modified*/Beam properties modified</TD> <TD>26.37855</TD><TD>83.69344</TD> <TD><a href="https://c3.nasa.gov/dashlink/static/media/dataset/ms103_beam_100_modi_ccs_o1_k2_nsm1_nm1.bdf"> 100_modi_ccs_o1_k2_nsm1_nm1</a> </TD></TR><TR> <TD>cosa,0,-sina</TD> <TD>as provided</TD> <TD>26.37949</TD><TD>83.69533</TD> <TD><a href="https://c3.nasa.gov/dashlink/static/media/dataset/ms103_beam_ca_0_nsa_ccs_o1_k_nsm1_nm1.bdf"> ca_0_nsa_ccs_o1_k_nsm1_nm1</a> </TD></TR><TR> <TD>cosa,0,-sina</TD> <TD>as provided</TD> <TD>26.37949</TD><TD>83.69533</TD> <TD><a href="https://c3.nasa.gov/dashlink/static/media/dataset/ms103_beam_ca_0_nsa_ccs_o1_k_nsm1_nm1.bdf"> ca_0_nsa_ccs_o1_k_nsm1_nm1</a> </TD></TR> </TABLE> <BR> <P> <P> Below are frequencies with different permutations of <LI> X1,X2,X3 definitions on CBEAM card <br>A) - 1,0,0 or <br>B) ca_0_nsa = cos(alpha),0,-sin(alpha) or <br> C)nca_0_sa = -cos(alpha),0,sin(alpha) on the beam card <br> <LI> coordinate system (CCS) is converted from original <LI> O1 indicates that offsets W1-W3 were included between the node and the shear center <LI> nsm1 = beam density is unchanged and includes the non-structural mass <LI> nsm2 = non-structural mass is separated from the beam mass density <LI> nm = calculation of offset N1,N2 for beginning and end of each beam segment is performed - different ways depending on the index <LI> modi = inertias provided by Boucke were modified by the alpha. <LI> _k_ = kappa(K1,K2) were included in the model - values provided by Boucke (A.B.) <LI> _k2_ = kappa(K1/K2) modified along with beam properties per equations provided by A.B. <P> <TABLE BORDER = 1> <TR>Description<TD></TD> <TD>1st mode</TD><TD>2nd mode</TD></TR> <TR><TD>Boucke provided frequencies</TD><TD>25.33879547</TD><TD>78.63642299</TD></TR> <TR><TD>100_modi_ccs_o1_k2_nsm1_nm1</TD><TD>26.37855</TD><TD>83.69344</TD></TR> <TR><TD>100_modi_ccs_o1_k_nsm1_nm1</TD><TD> 26.37854</TD><TD>83.6936</TD></TR> <TR><TD>ca_0_nsa_ccs_o1_k_nsm1_nm1</TD><TD> 26.37949</TD><TD>83.69533</TD></TR> <TR><TD>ca_0_nsa_modi_ccs_o1_k2_nsm1_nm1</TD><TD>26.3795</TD><TD>83.69762</TD></TR> <TR><TD>ca_0_nsa_modi_ccs_o1_k_nsm1_nm1</TD><TD>26.37949</TD<TD>83.69778</TD></TR> <TR><TD>nca_0_sa_ccs_o1_k </TD><TD>26.37466</TD<TD>83.62444</TD></TR><TR> <TD>nca_0_sa_ccs_o1_k_nsm1_nm1</TD><TD>26.33955</TD><TD>82.98219</TD> </TR><TR><TD>nca_0_sa_ccs_o1_k_nsm1_nm2</TD><TD>26.3778</TD<TD>83.67184</TD> </TR><TR> <TD>nca_0_sa_ccs_o1_k_nsm2_nm1</TD><TD>26.3452</TD<TD>83.07531</TD></TR> <TR> <TD>nca_0_sa_ccs_o1_k_nsm2_nm2</TD><TD>26.37785</TD<TD>83.6722</TD> </TR><TR><TD>nca_0_sa_modi_ccs_o1_k_nsm1_nm1</TD> <TD>26.33953</TD<TD>82.98429</TD></TR></TABLE> <P> Matlab code has been requested that reads the input deck. This code also does some debugging and calculation. The disclaimer is that it isn't easy to read. If you find a coding error please let me know. There are two codes - one for reading that also does some distance calculations and then a code that writes the BDF file.
    num_resources 17
    num_tags 8
    title HIRENASD Beam FEM