def getNeutralFormat(self, counter, typ, mapLoadCases): retval= dict() for key in self: comb= self[key] retval[comb.name]= nld.LoadComb(counter,comb.name,self.description,typ,nld.getComponentsFromStr(comb.expr,mapLoadCases)) counter+=1 return retval
def __init__(self): self.pointLoadCounter= 1 self.surfaceLoadCounter= 1 #Loads self.loads= nld.LoadData() lg= self.loads.loadGroups lg0= nld.LoadGroup(0,'Permanent loads') lg0.permanent= True lg["LG0"]= lg0 lg["LG1"]= nld.LoadGroup(1,'Variable loads') self.mapLoadCases= dict()
def dumpLoadPattern(self, counter, lpName, loadPatterns,permanentLoadCaseNames): loadCaseId= counter loadCaseName= lpName descr= ''#lp.description newLoadCase= None if(loadCaseName in permanentLoadCaseNames): newLoadCase= nld.LoadCase(loadCaseId,loadCaseName, descr,0,1) else: newLoadCase= nld.LoadCase(loadCaseId,loadCaseName, descr,1,1) newLoadCase.actionType= 'Variable' self.loads.loadCases[loadCaseId]= newLoadCase self.mapLoadCases[loadCaseName]= newLoadCase return newLoadCase
def dumpSurfaceLoads(self, lp, destLoadCase): '''Dump loads over surfaces as free surface loads.''' domain= lp.getDomain preprocessor= lp.getDomain.getPreprocessor eLoadIter= lp.loads.getElementalLoadIter eLoad= eLoadIter.next() loadSets= [] while eLoad: setName= 'surfaceLoadSet'+str(eLoad.tag) surfaceLoadSet= preprocessor.getSets.defSet(setName) surfaceLoadSet.elementalLoad= eLoad elemTags= eLoad.elementTags resultant= eLoad.getResultant(geom.Pos3d(0,0,0),True) #Total force over the elements. totalForce= geom.Vector3d(resultant.x,resultant.y,resultant.z) totalForceModulus= totalForce.getModulus() numberOfLoadedElements= len(elemTags) if(numberOfLoadedElements>1): #It seems there is a bug in totalArea= 0.0 for tag in elemTags: elem= domain.getMesh.getElement(tag) totalArea+= elem.getArea(True) if(elem): surfaceLoadSet.elements.append(elem) else: lmsg.error('element: '+ str(tag) + ' not found.') surfaceLoadSet.value= totalForceModulus surfaceLoadSet.vDir= [totalForce.x/surfaceLoadSet.value,totalForce.y/surfaceLoadSet.value,totalForce.z/surfaceLoadSet.value] surfaceLoadSet.value/= totalArea elementContours= surfaceLoadSet.elements.getContours(0.0) if(len(elementContours)>1): lmsg.error('surface load set: '+ setName + ' has more than one contour. Contours others than first are ignored.') surfaceLoadSet.polygon= elementContours[0] surfaceLoadSet.polygon.simplify(.01) #Deletes unnecesary vertices. loadSets.append(surfaceLoadSet) else: elem= domain.getMesh.getElement(elemTags[0]) #Only one element... pLoad= nld.PointForceRecord(destLoadCase, self.pointLoadCounter,elem.getPosCentroid(True),1.0) pLoad.value= totalForceModulus pLoad.vDir= [totalForce.x/pLoad.value,totalForce.y/pLoad.value,totalForce.z/pLoad.value] pLoad.tag= eLoad.tag destLoadCase.loads.punctualLoads.append(pLoad) self.pointLoadCounter+=1 eLoad= eLoadIter.next() for s in loadSets: sLoad= nld.SurfaceLoadRecord(destLoadCase, self.surfaceLoadCounter) sLoad.polygon= s.polygon sLoad.value= s.value sLoad.vDir= s.vDir destLoadCase.loads.surfaceLoads.append(sLoad) self.surfaceLoadCounter+=1
def getNeutralFormat(self, counter, typ, mapLoadCases): retval= dict() for key in self: comb= self[key] retval[comb.name]= nld.LoadComb(counter,comb.name,self.description,typ,nld.getComponentsFromStr(comb.expr,mapLoadCases)) counter+=1 return retval
def dumpSurfaceLoads(self, lp, destLoadCase): '''Dump loads over elements.''' eLoadIter = lp.loads.getElementalLoadIter el = eLoadIter.next() while el: eLoad = nld.ElementLoadRecord(destLoadCase, self.surfaceLoadCounter, 1.0) if (hasattr(el, "getGlobalForces")): eLoad.mode = False # Referred to local coordinate system. lf = xc.Vector([el.Wx, el.Wy, el.Wz]) eLoad.value = lf.Norm() if (eLoad.value >= 1e-3): eLoad.vDir = [ el.Wx / eLoad.value, el.Wy / eLoad.value, el.Wz / eLoad.value ] tags = el.elementTags for i in range(0, len(tags)): eLoad.tags.append(tags[i]) destLoadCase.loads.surfaceLoads.append(eLoad) else: lmsg.warning( 'loads2Neutral: vDir vector very small: "+ str(vDir) + " load ignored.' ) self.surfaceLoadCounter += 1 el = eLoadIter.next()
def dumpPointLoads(self, lp, destLoadCase): '''Dump loads over nodes.''' lIter= lp.loads.getNodalLoadIter nl= lIter.next() while nl: pLoad= nld.NodalLoadRecord(destLoadCase, self.pointLoadCounter,None,1.0) force= nl.getForce pLoad.value= force.Norm() pLoad.vDir= force.Normalized() pLoad.tag= nl.getNodeTag destLoadCase.loads.punctualLoads.append(pLoad) self.pointLoadCounter+=1 nl= lIter.next()
def dumpPointLoads(self, lp, destLoadCase): '''Dump loads over nodes as free punctual loads.''' lIter = lp.getNodalLoadIter nl = lIter.next() while nl: node = nl.getNode pLoad = nld.PointForceRecord(destLoadCase, self.pointLoadCounter, node.getInitialPos3d, 1.0) force = nl.getForce pLoad.value = force.Norm() pLoad.vDir = force.Normalized() pLoad.tag = nl.getNodeTag destLoadCase.loads.punctualLoads.append(pLoad) self.pointLoadCounter += 1 nl = lIter.next()