Abaqus計算完成后，為輸出某部分模型的結(jié)果數(shù)據(jù)，需要用到region，并且配套使用的是getSubset函數(shù)。

示例一創(chuàng)建一個變量center，將一個叫做PUNCH的節(jié)點set賦值給這個變量，通過整個模型的位移場displacement得到該節(jié)點set的輸出數(shù)據(jù)，代碼如下：

center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']

centerDisplacement = displacement.getSubset(region=center)

centerValues = centerDisplacement.values

for v in centerValues:

    print v.nodeLabel, v.data

其中的displacement是已經(jīng)創(chuàng)建好的位移場數(shù)據(jù)，centerDisplacement變量得到該節(jié)點集的所有數(shù)據(jù)，centerValues變量得到該節(jié)點集的數(shù)值數(shù)據(jù)，通過循環(huán)輸出節(jié)點號和當前數(shù)據(jù)。

示例二如下：

#創(chuàng)建變量topCenter，得到單元集CENT

topCenter =odb.rootAssembly.instances['PART-1-1'].elementSets['CENT']

#創(chuàng)建變量stressField，得到位移場數(shù)據(jù)

stressField = odb.steps['Step-2'].frames[3].fieldOutputs['S']

#創(chuàng)建變量field，得到該單元集的積分點的數(shù)據(jù)

field = stressField.getSubset(region=topCenter,

    position=INTEGRATION_POINT, elementType = 'CAX4')

#創(chuàng)建變量fieldValues，得到field的數(shù)值數(shù)據(jù)

fieldValues = field.values

#循環(huán)輸出單元號，單元上的積分點號以及積分點的應力分量

for v in fieldValues:

    print 'Element label = ', v.elementLabel,    

    if v.integrationPoint:

        print 'Integration Point = ', v.integrationPoint

    else:

        print

# For each tensor component.

    for component in v.data:

# Print using a format. The comma at the end of the

# print statement suppresses the carriage return.

        print '%-10.5f' % component,

# After each tuple has printed, print a carriage return.

    Print

最終的輸出結(jié)果：

Element label =  1 Integration Point =  1

S : 0.01230    -0.05658   0.00892    -0.00015 

Element label =  1 Integration Point =  2

S : 0.01313    -0.05659   0.00892    -0.00106 

Element label =  1 Integration Point =  3

S : 0.00619    -0.05642   0.00892    -0.00023 

Element label =  1 Integration Point =  4

S : 0.00697    -0.05642   0.00892    -0.00108 

Element label =  11 Integration Point =  1

S : 0.01281    -0.05660   0.00897    -0.00146 

Element label =  11 Integration Point =  2

S : 0.01183    -0.05651   0.00897    -0.00257 

Element label =  11 Integration Point =  3 ...

另：對于getSubset的position有以下幾個可選：

INTEGRATION_POINT   ！積分點

NODAL                ！節(jié)點（平均）

ELEMENT_NODAL      ！單元節(jié)點（未平均）

 

CENTROID             ！單元質(zhì)心點

轉(zhuǎn)自公眾號——ABAQUS大世界

旨在分享，若侵即刪.