Пример #1
0
#
# You should have received a copy of the GNU General Public License
# along with code_aster.  If not, see <http://www.gnu.org/licenses/>.
# --------------------------------------------------------------------

import code_aster
from code_aster.Commands import *

code_aster.init("--test")

test = code_aster.TestCase()

mesh = code_aster.Mesh()
mesh.readMedFile("test001f.mmed")

model = code_aster.Model(mesh)
model.addModelingOnMesh(code_aster.Physics.Mechanics,
                        code_aster.Modelings.Tridimensional)
model.build()

young = 200000.0
poisson = 0.3

acier = DEFI_MATERIAU(ELAS=_F(
    E=young,
    NU=poisson,
), )

affe_mat = code_aster.MaterialField(mesh)
affe_mat.addMaterialsOnMesh(acier)
affe_mat.buildWithoutExternalVariable()
Пример #2
0
#
# You should have received a copy of the GNU General Public License
# along with code_aster.  If not, see <http://www.gnu.org/licenses/>.
# --------------------------------------------------------------------

import code_aster
from code_aster.Commands import *

code_aster.init("--test")

test = code_aster.TestCase()

monMaillage = code_aster.Mesh()
monMaillage.readMedFile("test001f.mmed")

monModel = code_aster.Model(monMaillage)
monModel.addModelingOnMesh(code_aster.Physics.Mechanics,
                           code_aster.Modelings.Tridimensional)
monModel.build()

YOUNG = 200000.0
POISSON = 0.3

acier = DEFI_MATERIAU(ELAS=_F(
    E=YOUNG,
    NU=POISSON,
), )
acier.debugPrint(6)

affectMat = code_aster.MaterialField(monMaillage)
affectMat.addMaterialsOnMesh(acier)
Пример #3
0
test = code_aster.TestCase()

rank = code_aster.getMPIRank()
print("Nb procs", code_aster.getMPINumberOfProcs())
print("Rank", code_aster.getMPIRank())

pMesh2 = code_aster.ParallelMesh()
pMesh2.readMedFile("xxParallelMesh001a/%d.med"%rank, True)
pMesh2=DEFI_GROUP(reuse =pMesh2,MAILLAGE=pMesh2,CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI'))
del pMesh2

pMesh = code_aster.ParallelMesh()
pMesh.readMedFile("xxParallelMesh001a/%d.med"%rank, True)
pMesh.debugPrint(rank+30)

model = code_aster.Model(pMesh)
test.assertEqual(model.getType(), "MODELE_SDASTER")
model.addModelingOnMesh(code_aster.Physics.Mechanics,
                           code_aster.Modelings.Tridimensional)
model.build()

testMesh = model.getMesh()
test.assertEqual(testMesh.getType(), "MAILLAGE_P")

model.debugPrint(rank+30)

acier = DEFI_MATERIAU(ELAS = _F(E = 2.e11,
                                NU = 0.3,),)
acier.debugPrint(8)

affectMat = code_aster.MaterialField(pMesh)
Пример #4
0
#
# You should have received a copy of the GNU General Public License
# along with code_aster.  If not, see <http://www.gnu.org/licenses/>.
# --------------------------------------------------------------------

import code_aster
from code_aster.Commands import *

code_aster.init("--test")

test = code_aster.TestCase()

MA = code_aster.Mesh()
MA.readMedFile("xxContact001a.mmed")

MO2 = code_aster.Model(MA)
MO2.addModelingOnMesh(code_aster.Physics.Acoustic,
                      code_aster.Modelings.Tridimensional)
MO2.build()

load = code_aster.AcousticLoad(MO2)
load.addImposedPressureOnGroupOfCells(["FONDATION"], 1. + 2.j)
load.addImposedNormalSpeedOnGroupOfCells(["FONDATION"], 3. + 4.j)
load.addImpedanceOnMesh(5. + 6.j)
load.build()
load.debugPrint(8)

# Test trivial
test.assertTrue(True)

FIN()