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
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
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
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']) ################################### # Store in global data dictionary # ################################### from pyfem.util.dataStructures import GlobalData globdat = GlobalData(nodes, elements, dofs) ################################ # Solution procedure (Box 2.3) #
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'] ) ################################### # Store in global data dictionary # ################################### from pyfem.util.dataStructures import GlobalData globdat = GlobalData( nodes, elements, dofs ) ################################ # Solution procedure (Box 2.3) #