Abaqus – Czit notes

Debug Abaqus cosimulation with user subroutines on Windows

przez kszyh | 1 lipca 2021 - 12:46 |1 lipca 2021 Abaqus

abaqus cse listenerport=4444 configure=cosHTV6_config interactive
abaqus input=cosHTV6-CFD job=cosHTV6-CFD csedirector=localhost:4444 interactive

"c:\Program Files (x86)\IntelSWTools\parallel_studio_xe_2019\compilers_and_libraries_2019\\windows\bin\compilervars.bat" intel64 vs2013
abaqus debug -standard -input cosHTV6-STD -job cosHTV6-STD -user sub.obj -interactive -csedirector localhost:4444

abaqus cse listenerport=4444 configure=cosHTV6_config interactive

abaqus input=cosHTV6-CFD job=cosHTV6-CFD csedirector=localhost:4444 interactive

"c:\Program Files (x86)\IntelSWTools\parallel_studio_xe_2019\compilers_and_libraries_2019\\windows\bin\compilervars.bat" intel64 vs2013

abaqus debug -standard -input cosHTV6-STD -job cosHTV6-STD -user sub.obj -interactive -csedirector localhost:4444

Delete old files (Windows):

del *.sim *.com *.lck *.dat *.sta *.tmp *.SMA* *.cid *.msg *.prt *.odb *.log

1	del .sim .com .lck .dat .sta .tmp .SMA .cid .msg .prt .odb *.log

Run cosimulation with one command:

abaqus cosimulation input=cosHTV6-CFD,cosHTV6-STD job=cosHTV6-CFD,cosHTV6-STD user=empty.obj,sub.obj cpus=9 interactive cpuratio=8,1 configure=cosHTV6_config cosimjob=cosHTV6

1	abaqus cosimulation input=cosHTV6-CFD,cosHTV6-STD job=cosHTV6-CFD,cosHTV6-STD user=empty.obj,sub.obj cpus=9 interactive cpuratio=8,1 configure=cosHTV6_config cosimjob=cosHTV6

Periodic boundary conditions in Abaqus (Periodyczne warunki brzegowe)

przez kszyh | 30 sierpnia 2015 - 16:53 |1 września 2015 Abaqus, Python

Dodaj komentarz

Założenia: zdefiniowane sety z węzłów na każdej krawędzi jako LEFT, RIGHT, TOP, BOTTOM i stworzone dwa punkty referencyjne w złożeniu oraz dwa sety z nimi o nazwach REF-PT-1 i REF-PT-2

# -*- coding: mbcs -*-
# Do not delete the following import lines
from abaqus import *
from abaqusConstants import *
import collections
import __main__


import section
import regionToolset
import displayGroupMdbToolset as dgm
import part
import material
import assembly
import optimization
import step
import interaction
import load
import mesh
import job
import sketch
import visualization
import xyPlot
import displayGroupOdbToolset as dgo
import connectorBehavior

modelName = 'COMPRESSION_XY'
periodicityLeftRight = True
periodicityBottomTop = True

#modelName = getInput('Enter model name:')

root = mdb.models[modelName].rootAssembly


# Clear all sets and equations connected to previous constaints
names = [set[0] for set in root.sets.items()]

for name in names:
    if "PERIOD-" in name:
        del root.sets[name] 


names = [c[0] for c in  mdb.models[modelName].constraints.items()]

for name in names:
    if "PERIOD-" in name:
        del  mdb.models[modelName].constraints[name] 

# Get nodes from edges. Sets need to be prepared before this step
leftSet = root.sets['LEFT'].nodes
rightSet = root.sets['RIGHT'].nodes
bottomSet = root.sets['BOTTOM'].nodes
topSet = root.sets['TOP'].nodes

# Connect Left to Right node
if periodicityLeftRight:
    idsLR = collections.OrderedDict()
    for nodeL in leftSet:
        for nodeR in rightSet:
            if nodeL.coordinates[1] == nodeR.coordinates[1]:
                idsLR[nodeL.label] = nodeR.label
        
    sortedLeft = tuple(idsLR.keys())
    sortedRight = tuple(idsLR.values())

# Connect Bottom to Top node
if periodicityBottomTop:
    idsBT = collections.OrderedDict()
    for nodeB in bottomSet:
        for nodeT in topSet:
            if nodeB.coordinates[0] == nodeT.coordinates[0]:
                idsBT[nodeB.label] = nodeT.label
        
    sortedBottom = tuple(idsBT.keys())
    sortedTop = tuple(idsBT.values())

i = 1
for n1, n2 in zip(sortedLeft, sortedRight):
    leftSetName = "PERIOD-LEFT"+str(i)
    rightSetName = "PERIOD-RIGHT"+str(i)
    
    root.SetFromNodeLabels(name=leftSetName, nodeLabels=(('PART-1-1', (n1,)), ))
    root.SetFromNodeLabels(name=rightSetName, nodeLabels=(('PART-1-1', (n2,)), ))
    eqName = "PERIOD-1-"+str(i)
    mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 1), (-1.0, rightSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1)))
    # eqName = "PERIOD-2-"+str(i)
    # mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 2), (-1.0, rightSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2)))
    i += 1
    
for n1, n2 in zip(sortedBottom, sortedTop):
    bottomSetName = "PERIOD-BOTTOM"+str(i)
    topSetName = "PERIOD-TOP"+str(i)
    root.SetFromNodeLabels(name=bottomSetName, nodeLabels=(('PART-1-1', (n1,)), ))
    root.SetFromNodeLabels(name=topSetName, nodeLabels=(('PART-1-1', (n2,)), ))
    eqName = "PERIOD-2-"+str(i)
    mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 2), (-1.0, topSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2)))
    # eqName = "PERIOD-1-"+str(i)
    # mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 1), (-1.0, topSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1)))
    i += 1

100

# -*- coding: mbcs -*-

# Do not delete the following import lines

from abaqus import *

from abaqusConstants import *

import collections

import __main__

import section

import regionToolset

import displayGroupMdbToolset as dgm

import part

import material

import assembly

import optimization

import step

import interaction

import load

import mesh

import job

import sketch

import visualization

import xyPlot

import displayGroupOdbToolset as dgo

import connectorBehavior

modelName = 'COMPRESSION_XY'

periodicityLeftRight = True

periodicityBottomTop = True

#modelName = getInput('Enter model name:')

root = mdb.models[modelName].rootAssembly

# Clear all sets and equations connected to previous constaints

names = [set[0] for set in root.sets.items()]

for name in names:

if "PERIOD-" in name:

del root.sets[name]

names = [c[0] for c in mdb.models[modelName].constraints.items()]

for name in names:

if "PERIOD-" in name:

del mdb.models[modelName].constraints[name]

# Get nodes from edges. Sets need to be prepared before this step

leftSet = root.sets['LEFT'].nodes

rightSet = root.sets['RIGHT'].nodes

bottomSet = root.sets['BOTTOM'].nodes

topSet = root.sets['TOP'].nodes

# Connect Left to Right node

if periodicityLeftRight:

idsLR = collections.OrderedDict()

for nodeL in leftSet:

for nodeR in rightSet:

if nodeL.coordinates[1] == nodeR.coordinates[1]:

idsLR[nodeL.label] = nodeR.label

sortedLeft = tuple(idsLR.keys())

sortedRight = tuple(idsLR.values())

# Connect Bottom to Top node

if periodicityBottomTop:

idsBT = collections.OrderedDict()

for nodeB in bottomSet:

for nodeT in topSet:

if nodeB.coordinates[0] == nodeT.coordinates[0]:

idsBT[nodeB.label] = nodeT.label

sortedBottom = tuple(idsBT.keys())

sortedTop = tuple(idsBT.values())

i = 1

for n1, n2 in zip(sortedLeft, sortedRight):

leftSetName = "PERIOD-LEFT"+str(i)

rightSetName = "PERIOD-RIGHT"+str(i)

root.SetFromNodeLabels(name=leftSetName, nodeLabels=(('PART-1-1', (n1,)), ))

root.SetFromNodeLabels(name=rightSetName, nodeLabels=(('PART-1-1', (n2,)), ))

eqName = "PERIOD-1-"+str(i)

mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 1), (-1.0, rightSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1)))

# eqName = "PERIOD-2-"+str(i)

# mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 2), (-1.0, rightSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2)))

i += 1

for n1, n2 in zip(sortedBottom, sortedTop):

bottomSetName = "PERIOD-BOTTOM"+str(i)

topSetName = "PERIOD-TOP"+str(i)

root.SetFromNodeLabels(name=bottomSetName, nodeLabels=(('PART-1-1', (n1,)), ))

root.SetFromNodeLabels(name=topSetName, nodeLabels=(('PART-1-1', (n2,)), ))

eqName = "PERIOD-2-"+str(i)

mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 2), (-1.0, topSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2)))

# eqName = "PERIOD-1-"+str(i)

# mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 1), (-1.0, topSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1)))

i += 1

Archiwum kategorii: Abaqus

Debug Abaqus cosimulation with user subroutines on Windows

Periodic boundary conditions in Abaqus (Periodyczne warunki brzegowe)

Ostatnie wpisy

Najnowsze komentarze

Kategorie

Meta