用python计算每个单元体积与定义参数的乘积再加和: my code

wuweide127

用子程序做了个计算，单元结果中有一个自定义的参数SDV，希望每个单元的体积（这个用EVOL可以输出）与SDV相乘后，再把全部单元的乘积相加（就是EVOL1*SDV1+EVOL2*SDV2...+EVOLn*SDVn），请问用python如何实现？
我刚学python，水平很低，如果能有提供一个类似的例子指导一下也非常感谢！
谢谢大家！

-----------------------------------------------------------------------------------------------
后来自己根据abaqus中的帮助和例子，想了一个code，但是需要双重循环语句 （以下用红色注明了），运行起来速度比较慢，希望有高手可以帮我优化一下。谢谢！

wuweide127

my code:


from odbAccess import *

#
# get field
#

odb = openOdb(path='pytest.odb')
endSet = odb.rootAssembly.elementSets['ALLELE']
field_1 = odb.steps.values()[0].frames[-1].fieldOutputs['EVOL']
field_2 = odb.steps.values()[0].frames[-1].fieldOutputs['SDV11']
subField_1 = field_1.getSubset(region=endSet)
subField_2 = field_2.getSubset(region=endSet)

#
# sum values
#

sum = 0
vol = 0
for val_1 in subField_1.values:
   print 'ELENUMBER =', val_1.elementLabel
   print 'EVOL =', val_1.data
   for val_2 in subField_2.values:
      if val_1.elementLabel == val_2.elementLabel:
         print 'SDV11 =', val_2.data
         val = val_1.data * val_2.data
         break
   sum = sum + val
   vol = vol + val_1.data
   print 'val =', val
   print 'sum =', sum
ave = sum / vol

#
# print results
#

print 'All damage =', sum
print 'All volume =', vol
print 'Damage ratio =', ave

qdx天风海涛

你好，我找到一种方法能让计算快一点儿，但是我得出来的EVOL值却都是0，请大神给解答一下http://forum.simwe.com/thread-1116064-1-1.html

zjfswjtu

楼主，如果要提取11方向的平均应力怎么办？把S改成S11不行啊

ab小信

zjfswjtu 发表于 2017-2-10 14:49
楼主，如果要提取11方向的平均应力怎么办？把S改成S11不行啊

如果提取S是可以的，那么S11也一定是可以的

lyqday

楼主知道怎么把数组隔断