コード例 #1
0
 def get_GenBeams(match):
     d = match.groupdict()
     sec_id = str_to_int(d["geono"])
     gen_props = GeneralProperties(
         ax=roundoff(d["area"]),
         ix=roundoff(d["ix"]),
         iy=roundoff(d["iy"]),
         iz=roundoff(d["iz"]),
         iyz=roundoff(d["iyz"]),
         wxmin=roundoff(d["wxmin"]),
         wymin=roundoff(d["wymin"]),
         wzmin=roundoff(d["wzmin"]),
         shary=roundoff(d["shary"]),
         sharz=roundoff(d["sharz"]),
         scheny=roundoff(d["shceny"]),
         schenz=roundoff(d["shcenz"]),
         sy=float(d["sy"]),
         sz=float(d["sz"]),
     )
     if sec_id in fem.parent.sections.idmap.keys():
         sec = fem.parent.sections.get_by_id(sec_id)
         sec._genprops = gen_props
         gen_props.parent = sec
     else:
         sec = Section(name=f"GB{sec_id}",
                       sec_id=sec_id,
                       sec_type="GENBEAM",
                       genprops=gen_props,
                       parent=fem.parent)
         gen_props.parent = sec
         fem.parent.sections.add(sec)
コード例 #2
0
ファイル: read_sections.py プロジェクト: Krande/adapy
def add_general_sections(match, fem) -> None:
    d = match.groupdict()
    sec_id = str_to_int(d["geono"])
    gen_props = GeneralProperties(
        Ax=roundoff(d["area"], 10),
        Ix=roundoff(d["ix"], 10),
        Iy=roundoff(d["iy"], 10),
        Iz=roundoff(d["iz"], 10),
        Iyz=roundoff(d["iyz"], 10),
        Wxmin=roundoff(d["wxmin"]),
        Wymin=roundoff(d["wymin"]),
        Wzmin=roundoff(d["wzmin"]),
        Shary=roundoff(d["shary"]),
        Sharz=roundoff(d["sharz"]),
        Shceny=roundoff(d["shceny"]),
        Shcenz=roundoff(d["shcenz"]),
        Sy=float(d["sy"]),
        Sz=float(d["sz"]),
    )

    if sec_id in fem.parent.sections.id_map.keys():
        sec = fem.parent.sections.get_by_id(sec_id)
        sec._genprops = gen_props
        gen_props.parent = sec
    else:
        stype = Section.TYPES.GENERAL
        sec = Section(name=f"GB{sec_id}",
                      sec_id=sec_id,
                      sec_type=stype,
                      genprops=gen_props,
                      parent=fem.parent)
        gen_props.parent = sec
        fem.parent.sections.add(sec)
コード例 #3
0
ファイル: read_sections.py プロジェクト: Krande/adapy
def get_flatbar(match, sect_names, fem) -> Section:
    d = match.groupdict()
    sec_id = str_to_int(d["geono"])
    return Section(
        name=sect_names[sec_id],
        sec_id=sec_id,
        sec_type=Section.TYPES.FLATBAR,
        h=roundoff(d["hz"]),
        w_top=roundoff(d["bt"]),
        w_btn=roundoff(d["bb"]),
        genprops=GeneralProperties(Sfy=float(d["sfy"]), Sfz=float(d["sfz"])),
        parent=fem.parent,
    )
コード例 #4
0
ファイル: reader.py プロジェクト: Caladeux/adapy
 def get_BoxBeams(match):
     d = match.groupdict()
     sec_id = str_to_int(d["geono"])
     return Section(
         name=sect_names[sec_id],
         sec_id=sec_id,
         sec_type="BG",
         h=roundoff(d["hz"]),
         w_top=roundoff(d["by"]),
         w_btn=roundoff(d["by"]),
         t_w=roundoff(d["ty"]),
         t_ftop=roundoff(d["tt"]),
         t_fbtn=roundoff(d["tb"]),
         genprops=GeneralProperties(sfy=float(d["sfy"]), sfz=float(d["sfz"])),
         parent=fem.parent,
     )
コード例 #5
0
ファイル: read_sections.py プロジェクト: Krande/adapy
def get_isection(match, sect_names, fem) -> Section:
    d = match.groupdict()
    sec_id = str_to_int(d["geono"])
    name = sect_names[sec_id]
    return Section(
        name=name,
        sec_id=sec_id,
        sec_type=Section.TYPES.IPROFILE,
        h=roundoff(d["hz"]),
        t_w=roundoff(d["ty"]),
        w_top=roundoff(d["bt"]),
        w_btn=roundoff(d["bb"]),
        t_ftop=roundoff(d["tt"]),
        t_fbtn=roundoff(d["tb"]),
        genprops=GeneralProperties(Sfy=float(d["sfy"]), Sfz=float(d["sfz"])),
        parent=fem.parent,
    )
コード例 #6
0
ファイル: read_sections.py プロジェクト: Krande/adapy
def get_tubular_section(match, sect_names, fem) -> Section:
    d = match.groupdict()
    sec_id = str_to_int(d["geono"])
    if sec_id not in sect_names:
        sec_name = f"TUB{sec_id}"
    else:
        sec_name = sect_names[sec_id]
    t = d["t"] if d["t"] is not None else (d["dy"] - d["di"]) / 2
    return Section(
        name=sec_name,
        sec_id=sec_id,
        sec_type=Section.TYPES.TUBULAR,
        r=roundoff(float(d["dy"]) / 2),
        wt=roundoff(t),
        genprops=GeneralProperties(Sfy=float(d["sfy"]), Sfz=float(d["sfz"])),
        parent=fem.parent,
    )
コード例 #7
0
ファイル: reader.py プロジェクト: Caladeux/adapy
 def get_gpipe(match):
     d = match.groupdict()
     sec_id = str_to_int(d["geono"])
     if sec_id not in sect_names:
         sec_name = f"TUB{sec_id}"
     else:
         sec_name = sect_names[sec_id]
     t = d["t"] if d["t"] is not None else (d["dy"] - d["di"]) / 2
     return Section(
         name=sec_name,
         sec_id=sec_id,
         sec_type="TUB",
         r=roundoff(float(d["dy"]) / 2),
         wt=roundoff(t),
         genprops=GeneralProperties(sfy=float(d["sfy"]), sfz=float(d["sfz"])),
         parent=fem.parent,
     )
コード例 #8
0
    def interpret_section(profile_name, sec_type, props):
        props_clean = [
            roundoff(x)
            for x in filter(lambda x: x.strip() != "", props.split(","))
        ]
        if sec_type.upper() == "BOX":
            b, h, t1, t2, t3, t4 = props_clean
            return Section(
                profile_name,
                "BG",
                h=h,
                w_btn=b,
                w_top=b,
                t_w=t1,
                t_fbtn=t4,
                t_ftop=t2,
                parent=fem,
            )
        elif sec_type.upper() == "CIRC":
            return Section(profile_name, "CIRC", r=props_clean[0], parent=fem)
        elif sec_type.upper() == "I":
            (
                l,
                h,
                b1,
                b2,
                t1,
                t2,
                t3,
            ) = props_clean
            return Section(
                profile_name,
                "IG",
                h=h,
                w_btn=b1,
                w_top=b2,
                t_w=t3,
                t_fbtn=t1,
                t_ftop=t2,
                parent=fem,
            )
        elif sec_type.upper() == "L":
            b, h, t1, t2 = props_clean
            return Section(profile_name,
                           "HP",
                           h=h,
                           w_btn=b,
                           t_w=t2,
                           t_fbtn=t1,
                           parent=fem)
        elif sec_type.upper() == "PIPE":
            r, t = props_clean
            return Section(profile_name, "TUB", r=r, wt=t, parent=fem)
        elif sec_type.upper() == "TRAPEZOID":
            # Currently converts Trapezoid to general beam
            b, h, a, d = props_clean
            # Assuming the Abaqus trapezoid element is symmetrical
            c = (b - a) / 2

            # The properties were quickly copied from a resource online. Most likely it contains error
            # https: // www.efunda.com / math / areas / trapezoidJz.cfm
            genprops = GeneralProperties(
                ax=h * (a + b) / 2,
                ix=h * (b * h**2 + 3 * a * h**2 + a**3 + 3 * a * c**2 +
                        3 * c * a**2 + b**3 + c * b**2 + a * b**2 + b * c**2 +
                        2 * a * b * c + b * a**2),
                iy=(h**3) * (3 * a + b) / 12,
                iz=h * (a**3 + 3 * a * c**2 + 3 * c * a**2 + b**3 + c * b**2 +
                        a * b**2 + 2 * a * b * c + b * a**2) / 12,
            )
            return Section(profile_name,
                           "GENBEAM",
                           genprops=genprops,
                           parent=fem)
        else:
            raise ValueError(
                f'Currently unsupported section type "{sec_type}"')