Ejemplo n.º 1
0
class MEBODKT(Element):
    """Boundary element for DKT model : SEG2"""
    meshType = MT.SEG2
    elrefe =(
            ElrefeLoc(MT.SE2, gauss = ('RIGI=FPG2',),),
        )
    calculs = (

        OP.CARA_SECT_POUT3(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT3.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), ),
            para_out=((SP.PCASECT, ECASECT), ),
        ),

        OP.CARA_SECT_POUT4(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT4.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), (SP.PORIGIN, LC.CGEOM3D),
                     ),
            para_out=((SP.PVECTU1, MVECTUR), (SP.PVECTU2, MVECTUR),
                     ),
        ),

        OP.CARA_SECT_POUT5(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT5.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), (SP.PNUMMOD, LC.CNUMMOD),
                     (SP.PORIGFI, LC.CGEOM3D), (SP.PORIGIN, LC.CGEOM3D),
                     ),
            para_out=((SP.PVECTU1, MVECTUR), (SP.PVECTU2, MVECTUR),
                     (SP.PVECTU3, MVECTUR), ),
        ),

        OP.CHAR_MECA_FF1D3D(te=34,
            para_in=((SP.PFF1D3D, MFORCEF), (SP.PGEOMER, MGEOMER),
                     (SP.PTEMPSR, LC.MTEMPSR), ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),

        OP.CHAR_MECA_FR1D3D(te=34,
            para_in=((SP.PFR1D3D, MFORCER), (SP.PGEOMER, MGEOMER),
                     ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),

        OP.TOU_INI_ELEM(te=99,
            para_out=((OP.TOU_INI_ELEM.PGEOM_R, LC.CGEOM3D), ),
        ),

        OP.TOU_INI_ELGA(te=99,
            para_out=((OP.TOU_INI_ELGA.PGEOM_R, EGEOMER), ),
        ),

        OP.TOU_INI_ELNO(te=99,
            para_out=((OP.TOU_INI_ELNO.PGEOM_R, MGEOMER), ),
        ),

    )
Ejemplo n.º 2
0
class MEBOCQ3(Element):
    """Please document this element"""
    meshType = MT.SEG3
    elrefe =(
            ElrefeLoc(MT.SE3, gauss = ('RIGI=FPG4',),),
        )
    calculs = (

        OP.CARA_SECT_POUT3(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT3.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), ),
            para_out=((SP.PCASECT, ECASECT), ),
        ),

        OP.CARA_SECT_POUT4(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT4.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), (SP.PORIGIN, LC.CGEOM3D),
                     ),
            para_out=((SP.PVECTU1, MVECTUR), (SP.PVECTU2, MVECTUR),
                     ),
        ),

        OP.CARA_SECT_POUT5(te=570,
            para_in=((SP.PCACOQU, CCACOQU), (OP.CARA_SECT_POUT5.PCAORIE, LC.CGEOM3D),
                     (SP.PGEOMER, MGEOMER), (SP.PNUMMOD, LC.CNUMMOD),
                     (SP.PORIGFI, LC.CGEOM3D), (SP.PORIGIN, LC.CGEOM3D),
                     ),
            para_out=((SP.PVECTU1, MVECTUR), (SP.PVECTU2, MVECTUR),
                     (SP.PVECTU3, MVECTUR), ),
        ),

        OP.CHAR_MECA_FF1D3D(te=418,
            para_in=((SP.PFF1D3D, MFORCEF), (SP.PGEOMER, MGEOMER),
                     (SP.PTEMPSR, LC.MTEMPSR), ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),

        OP.CHAR_MECA_FR1D3D(te=418,
            para_in=((SP.PFR1D3D, MFORCER), (SP.PGEOMER, MGEOMER),
                     ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),

        OP.TOU_INI_ELEM(te=99,
            para_out=((OP.TOU_INI_ELEM.PGEOM_R, LC.CGEOM3D), ),
        ),

        OP.TOU_INI_ELGA(te=99,
            para_out=((OP.TOU_INI_ELGA.PGEOM_R, EGEOMER), ),
        ),

        OP.TOU_INI_ELNO(te=99,
            para_out=((OP.TOU_INI_ELNO.PGEOM_R, MGEOMER), ),
        ),

    )
Ejemplo n.º 3
0
class MECA_ARETE2(Element):
    """Please document this element"""
    meshType = MT.SEG2
    elrefe = (ElrefeLoc(
        MT.SE2,
        gauss=('RIGI=FPG2', ),
    ), )
    calculs = (
        OP.CHAR_MECA_FF1D3D(
            te=469,
            para_in=(
                (SP.PFF1D3D, MFORCEF),
                (SP.PGEOMER, NGEOMER),
                (SP.PTEMPSR, LC.MTEMPSR),
            ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),
        OP.CHAR_MECA_FR1D3D(
            te=469,
            para_in=(
                (SP.PFR1D3D, MFORCE),
                (SP.PGEOMER, NGEOMER),
            ),
            para_out=((SP.PVECTUR, MVECTUR), ),
        ),
        OP.COOR_ELGA(
            te=478,
            para_in=((SP.PGEOMER, NGEOMER), ),
            para_out=((OP.COOR_ELGA.PCOORPG, EGGEOP_R), ),
        ),
        OP.TOU_INI_ELEM(
            te=99,
            para_out=((OP.TOU_INI_ELEM.PGEOM_R, LC.CGEOM3D), ),
        ),
        OP.TOU_INI_ELGA(
            te=99,
            para_out=(
                (OP.TOU_INI_ELGA.PDEPL_R, EDEPLPG),
                (OP.TOU_INI_ELGA.PGEOM_R, EGGEOM_R),
                (OP.TOU_INI_ELGA.PNEUT_F, EGNEUT_F),
                (OP.TOU_INI_ELGA.PNEUT_R, EGNEUT_R),
            ),
        ),
        OP.TOU_INI_ELNO(
            te=99,
            para_out=(
                (OP.TOU_INI_ELNO.PGEOM_R, NGEOMER),
                (OP.TOU_INI_ELNO.PNEUT_F, LC.ENNEUT_F),
                (OP.TOU_INI_ELNO.PNEUT_R, LC.ENNEUT_R),
            ),
        ),
    )