Esempio n. 1
0
def InputRead(fname):

    if fname[-4:] == '.pro':
        props = fileParser(fname)
    else:
        props = fileParser(fname + '.pro')

    dataFileName = props.input

    nodes = NodeSet()
    nodes.readFromFile(dataFileName)

    elems = ElementSet(nodes, props)
    elems.readFromFile(dataFileName)

    dofs = DofSpace(elems)
    dofs.readFromFile(dataFileName)

    globdat = GlobalData(nodes, elems, dofs)

    globdat.readFromFile(dataFileName)

    globdat.active = True
    globdat.cycle = 0
    globdat.prefix = os.path.splitext(fname)[0]

    return props, globdat
Esempio n. 2
0
def InputReader( argv ):

  options, remainder = getopt.getopt( argv , 'a:k:v', ['all','author='])

  proFileName  = argv[1]
  props        = fileParser( proFileName )

  dataFileName = props.input

  nodes = NodeSet()
  nodes.readFromFile( dataFileName )
  
  elems = ElementSet( nodes , props )
  elems.readFromFile( dataFileName )
  
  dofs = DofSpace( elems )
  dofs.readFromFile( dataFileName )

  globdat = GlobalData( nodes, elems, dofs ) 

  globdat.readFromFile( dataFileName )

  globdat.active = True
  globdat.cycle  = 0
  globdat.prefix = os.path.splitext(proFileName)[0]
	
  return props,globdat
Esempio n. 3
0
def InputReader(argv):

    options, remainder = getopt.getopt(argv, 'a:k:v', ['all', 'author='])

    proFileName = argv[1]
    props = fileParser(proFileName)

    dataFileName = props.input

    nodes = NodeSet()
    nodes.readFromFile(dataFileName)

    elems = ElementSet(nodes, props)
    elems.readFromFile(dataFileName)

    dofs = DofSpace(elems)
    dofs.readFromFile(dataFileName)

    globdat = GlobalData(nodes, elems, dofs)

    globdat.readFromFile(dataFileName)

    globdat.active = True
    globdat.cycle = 0
    globdat.prefix = os.path.splitext(proFileName)[0]

    return props, globdat
Esempio n. 4
0
def InputRead( fname , dname = None , parameters = None ):

  if dname is not None:
    with open(dname, 'rb') as f:
      data = pickle.load(f)
      props = data["props"]
  
  if fname is not None:
    if fname[-4:] == '.pro':
      props        = fileParser( fname )
    else:
      props        = fileParser( fname+'.pro')
    
  if parameters is not None:  
    for p in parameters:
      x = p.split("=")
      props.store(x[0],x[1])
      
  if dname is not None:
    return props,data["globdat"]

  dataFileName = props.input

  logger = setLogger( props )

  nodes = NodeSet()
  nodes.readFromFile( dataFileName )
  logger.info(nodes)
  
  elems = ElementSet( nodes , props )
  elems.readFromFile( dataFileName )
  logger.info(elems)  
  
  dofs = DofSpace( elems )
  dofs.readFromFile( dataFileName )

  globdat = GlobalData( nodes, elems, dofs ) 

  globdat.readFromFile( dataFileName )

  globdat.active = True
  globdat.prefix = os.path.splitext(fname)[0]
  
  globdat.contact = Contact( props )
  	
  return props,globdat
Esempio n. 5
0
#NodeSet
from pyfem.fem.NodeSet import NodeSet

nodes = NodeSet()

nodes.add(1, [
    0., 0.
])  #Attachement point of spring (can be positioned anywhere along the y-axis)
nodes.add(2, [-10., 0.])  #Left support
nodes.add(3, [10., 0.])  #Right support
nodes.add(4, [0., 0.5])  #Loading point

#ElementSet
from pyfem.fem.ElementSet import ElementSet

elements = ElementSet(nodes, props)

elements.add(1, 'TrussElem', [2, 4])
elements.add(2, 'TrussElem', [3, 4])
elements.add(3, 'SpringElem', [1, 4])

#DofSpace
from pyfem.fem.DofSpace import DofSpace

dofs = DofSpace(elements)

dofs.constrain(1, ['u', 'v'])
dofs.constrain(2, ['u', 'v'])
dofs.constrain(3, ['u', 'v'])

###################################
Esempio n. 6
0
from pyfem.fem.NodeSet import NodeSet

nodes = NodeSet()
nodes.add(1, [
    0., 0.
])  #Attachement point of spring (can be positioned anywhere along the y-axis)
nodes.add(2, [-b, 0.])  #Left support
nodes.add(3, [b, 0.])  #Right support
nodes.add(4, [0., h])  #Loading point

#ElementSet
from pyfem.fem.ElementSet import ElementSet
from pyfem.elements.Truss import Truss
from pyfem.elements.Spring import Spring

elements = ElementSet(nodes)
elements.add(1, Truss([2, 4]))
elements.add(2, Truss([3, 4]))
elements.add(3, Spring([1, 4]))

#Add groups
elements.addGroup('Trusses', [1, 2])
elements.addGroup('Spring', [3])

#DofSpace
from pyfem.fem.DofSpace import DofSpace

dofs = DofSpace(elements)

dofs.constrain(1, ['u', 'v'])
dofs.constrain(2, ['u', 'v'])
Esempio n. 7
0
#############################

#NodeSet
from pyfem.fem.NodeSet import NodeSet

nodes = NodeSet()

nodes.add( 1, [  0.,0.  ] ) #Attachement point of spring (can be positioned anywhere along the y-axis)
nodes.add( 2, [-10.,0.  ] ) #Left support
nodes.add( 3, [ 10.,0.  ] ) #Right support
nodes.add( 4, [ 0. ,0.5 ] ) #Loading point

#ElementSet
from pyfem.fem.ElementSet import ElementSet

elements = ElementSet( nodes , props )

elements.add( 1, 'TrussElem'  , [2,4]  )
elements.add( 2, 'TrussElem'  , [3,4]  )
elements.add( 3, 'SpringElem' , [1,4]  )

#DofSpace
from pyfem.fem.DofSpace import DofSpace

dofs = DofSpace( elements )

dofs.constrain( 1, ['u','v'] )
dofs.constrain( 2, ['u','v'] )
dofs.constrain( 3, ['u','v'] )

###################################