def getTangentStiffness(self, elemdat):
sData = getElemShapeData(elemdat.coords)
kin = Kinematics(3, 6)
elemdat.outlabel.append(["s11", "s22", "s33", "s23", "s13", "s12"])
elemdat.outdata = zeros(shape=(len(elemdat.nodes), 6))
for iData in sData:
b = self.getBmatrix(iData.dhdx)
kin.strain = dot(b, elemdat.state)
kin.dstrain = dot(b, elemdat.Dstate)
sigma, tang = self.mat.getStress(kin)
elemdat.stiff += dot(b.transpose(), dot(tang, b)) * iData.weight
elemdat.fint += dot(b.transpose(), sigma) * iData.weight
elemdat.outdata += outer(ones(len(elemdat.nodes)), sigma)
elemdat.outdata *= 1.0 / len(sData)
def getTangentStiffness(self, element_data):
sData = getElemShapeData(element_data.coords)
kinematic = Kinematics(2, 3)
element_data.outlabel.append("stresses")
element_data.outdata = zeros(shape=(len(element_data.nodes), 3))
for iData in sData:
B_matrix = self.getBmatrix(iData.dhdx)
kinematic.strain = dot(B_matrix, element_data.state)
kinematic.dstrain = dot(B_matrix, element_data.Dstate)
sigma, tangent_moduli = self.mat.getStress(kinematic)
element_data.stiff += dot(B_matrix.transpose(),
dot(tangent_moduli,
B_matrix)) * iData.weight
element_data.fint += dot(B_matrix.transpose(),
sigma) * iData.weight
element_data.outdata += outer(ones(len(element_data.nodes)), sigma)
element_data.outdata *= 1.0 / len(sData)
def getInternalForce(self, elemdat):
sData = getElemShapeData(elemdat.coords)
kin = Kinematics(3, 6)
elemdat.outlabel.append(["s11", "s22", "s33", "s23", "s13", "s12"])
elemdat.outdata = zeros(shape=(len(elemdat.nodes), 6))
for iData in sData:
b = self.getBmatrix(iData.dhdx)
kin.strain = dot(b, elemdat.state)
kin.dstrain = dot(b, elemdat.Dstate)
sigma, tang = self.mat.getStress(kin)
elemdat.fint += dot(b.transpose(), sigma) * iData.weight
self.appendNodalOutput(self.mat.outLabels(), self.mat.outData())
def getInternalForce ( self, elemdat ):
sData = getElemShapeData( elemdat.coords )
elemdat.outlabel.append("stresses")
elemdat.outdata = zeros( shape=(len(elemdat.nodes),3) )
kin = Kinematics(2,3)
for iData in sData:
b = self.getBmatrix( iData.dhdx )
kin.strain = dot ( b , elemdat.state )
kin.dstrain = dot ( b , elemdat.Dstate )
sigma,tang = self.mat.getStress( kin )
elemdat.fint += dot ( b.transpose() , sigma ) * iData.weight
elemdat.outdata += outer( ones(len(self)), sigma )
elemdat.outdata *= 1.0 / len(sData)
def getInternalForce(self, elemdat):
sData = getElemShapeData(elemdat.coords)
elemdat.outlabel.append("stresses")
elemdat.outdata = zeros(shape=(len(elemdat.nodes), 3))
kin = Kinematics(2, 3)
for iData in sData:
b = self.getBmatrix(iData.dhdx)
kin.strain = dot(b, elemdat.state)
kin.dstrain = dot(b, elemdat.Dstate)
sigma, tang = self.mat.getStress(kin)
elemdat.fint += dot(b.transpose(), sigma) * iData.weight
elemdat.outdata += outer(ones(len(self)), sigma)
elemdat.outdata *= 1.0 / len(sData)
def getInternalForce(self, element_data):
sData = getElemShapeData(element_data.coords)
element_data.outlabel.append("stresses")
element_data.outlabel.append("stress_elements")
element_data.stress_elements = zeros(shape=(len(element_data.nodes),
3),
dtype=float)
element_data.outdata = zeros(shape=(len(element_data.nodes), 3))
kinematic = Kinematics(2, 3)
for iData in sData:
B_matrix = self.getBmatrix(iData.dhdx)
kinematic.strain = dot(B_matrix, element_data.state)
kinematic.dstrain = dot(B_matrix, element_data.Dstate)
sigma, tangentModuli = self.mat.getStress(
kinematic) # Original code
#sigma = self.mat.mat.sigma # Minh's modification code
element_data.fint += dot(B_matrix.transpose(),
sigma) * iData.weight
element_data.outdata += outer(ones(len(self)), sigma)
element_data.stress_elements += outer(ones(len(self)), sigma)
element_data.outdata *= 1.0 / len(sData)
element_data.stress_elements *= 1.0 / len(sData)
