# # 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()
# # 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)
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)
# # 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()