Create a gist now

Instantly share code, notes, and snippets.

anonymous /heat_aba.py
Created Jul 16, 2017

What would you like to do?
Download ZIP
Simple heat analysis in abaqus
Raw
heat_aba.py
# -*- coding: mbcs -*-
from part import *
from material import *
from section import *
from optimization import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=200.0)
mdb.models['Model-1'].sketches['__profile__'].rectangle(point1=(-5.0, 5.0),
point2=(5.0, -5.0))
mdb.models['Model-1'].sketches['__profile__'].ObliqueDimension(textPoint=(
-0.536808013916016, 8.91543197631836), value=10.0, vertex1=
mdb.models['Model-1'].sketches['__profile__'].vertices.findAt((5.0, 5.0), )
, vertex2=mdb.models['Model-1'].sketches['__profile__'].vertices.findAt((
-5.0, 5.0), ))
mdb.models['Model-1'].sketches['__profile__'].ObliqueDimension(textPoint=(
9.66257476806641, 1.50380706787109), value=10.0, vertex1=
mdb.models['Model-1'].sketches['__profile__'].vertices.findAt((5.0, -5.0),
), vertex2=mdb.models['Model-1'].sketches['__profile__'].vertices.findAt((
5.0, 5.0), ))
mdb.models['Model-1'].Part(dimensionality=THREE_D, name='Part-1', type=
DEFORMABLE_BODY)
mdb.models['Model-1'].parts['Part-1'].BaseSolidExtrude(depth=100.0, sketch=
mdb.models['Model-1'].sketches['__profile__'])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].Material(name='Copper')
mdb.models['Model-1'].materials['Copper'].Conductivity(table=((387.0, ), ))
mdb.models['Model-1'].materials['Copper'].Density(table=((8.9e-09, ), ))
mdb.models['Model-1'].materials['Copper'].SpecificHeat(table=((390000000.0, ),
))
mdb.models['Model-1'].materials['Copper'].Expansion(table=((1.67e-05, ), ))
mdb.models['Model-1'].HomogeneousSolidSection(material='Copper', name='Copper',
thickness=None)
mdb.models['Model-1'].parts['Part-1'].SectionAssignment(offset=0.0,
offsetField='', offsetType=MIDDLE_SURFACE, region=Region(
cells=mdb.models['Model-1'].parts['Part-1'].cells.findAt(((-5.0, -1.666667,
66.666667), ), )), sectionName='Copper', thicknessAssignment=FROM_SECTION)
mdb.models['Model-1'].rootAssembly.DatumCsysByDefault(CARTESIAN)
mdb.models['Model-1'].rootAssembly.Instance(dependent=ON, name='Part-1-1',
part=mdb.models['Model-1'].parts['Part-1'])
mdb.models['Model-1'].rootAssembly.Instance(dependent=ON, name='Part-1-2',
part=mdb.models['Model-1'].parts['Part-1'])
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].translate(vector=(
11.0, 0.0, 0.0))
mdb.models['Model-1'].rootAssembly.translate(instanceList=('Part-1-2', ),
vector=(59.0, 0.0, 0.0))
mdb.models['Model-1'].rootAssembly.translate(instanceList=('Part-1-2', ),
vector=(-10.0, 0.0, 0.0))
mdb.models['Model-1'].HeatTransferStep(amplitude=RAMP, name='Step-1', previous=
'Initial', response=STEADY_STATE)
mdb.models['Model-1'].Temperature(createStepName='Initial',
crossSectionDistribution=CONSTANT_THROUGH_THICKNESS, distributionType=
UNIFORM, magnitudes=(23.0, ), name='Predefined Field-1', region=Region(
cells=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].cells.findAt(
((-5.0, -1.666667, 66.666667), ), )+\
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].cells.findAt(((
55.0, -1.666667, 66.666667), ), ),
faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.findAt(
((5.0, 1.666667, 66.666667), ), ((1.666667, -5.0, 66.666667), ), ((-5.0,
-1.666667, 66.666667), ), ((-1.666667, 5.0, 66.666667), ), ((1.666667,
1.666667, 100.0), ), ((-1.666667, 1.666667, 0.0), ), )+\
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].faces.findAt(((
65.0, 1.666667, 66.666667), ), ((61.666667, -5.0, 66.666667), ), ((55.0,
-1.666667, 66.666667), ), ((58.333333, 5.0, 66.666667), ), ((61.666667,
1.666667, 100.0), ), ((58.333333, 1.666667, 0.0), ), ),
edges=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].edges.findAt(
((5.0, 2.5, 100.0), ), ((5.0, -5.0, 25.0), ), ((5.0, -2.5, 0.0), ), ((5.0,
5.0, 25.0), ), ((2.5, -5.0, 100.0), ), ((-5.0, -5.0, 25.0), ), ((-2.5,
-5.0, 0.0), ), ((-5.0, -2.5, 100.0), ), ((-5.0, 5.0, 25.0), ), ((-5.0, 2.5,
0.0), ), ((-2.5, 5.0, 100.0), ), ((2.5, 5.0, 0.0), ), )+\
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].edges.findAt(((
65.0, 2.5, 100.0), ), ((65.0, -5.0, 25.0), ), ((65.0, -2.5, 0.0), ), ((
65.0, 5.0, 25.0), ), ((62.5, -5.0, 100.0), ), ((55.0, -5.0, 25.0), ), ((
57.5, -5.0, 0.0), ), ((55.0, -2.5, 100.0), ), ((55.0, 5.0, 25.0), ), ((
55.0, 2.5, 0.0), ), ((57.5, 5.0, 100.0), ), ((62.5, 5.0, 0.0), ), ),
vertices=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].vertices.findAt(
((5.0, 5.0, 100.0), ), ((5.0, -5.0, 100.0), ), ((5.0, -5.0, 0.0), ), ((5.0,
5.0, 0.0), ), ((-5.0, -5.0, 100.0), ), ((-5.0, -5.0, 0.0), ), ((-5.0, 5.0,
100.0), ), ((-5.0, 5.0, 0.0), ), )+\
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].vertices.findAt(((
65.0, 5.0, 100.0), ), ((65.0, -5.0, 100.0), ), ((65.0, -5.0, 0.0), ), ((
65.0, 5.0, 0.0), ), ((55.0, -5.0, 100.0), ), ((55.0, -5.0, 0.0), ), ((55.0,
5.0, 100.0), ), ((55.0, 5.0, 0.0), ), )))
mdb.models['Model-1'].FilmConditionProp(dependencies=0, name='IntProp-1',
property=((0.003, ), ), temperatureDependency=OFF)
mdb.models['Model-1'].FilmCondition(createStepName='Step-1', definition=
PROPERTY_REF, interactionProperty='IntProp-1', name='Int-1', sinkAmplitude=
'', sinkDistributionType=UNIFORM, sinkFieldName='', sinkTemperature=40.0,
surface=Region(
side1Faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.findAt(
((5.0, 1.666667, 66.666667), ), ((1.666667, -5.0, 66.666667), ), ((-5.0,
-1.666667, 66.666667), ), ((-1.666667, 5.0, 66.666667), ), )+\
mdb.models['Model-1'].rootAssembly.instances['Part-1-2'].faces.findAt(((
65.0, 1.666667, 66.666667), ), ((61.666667, -5.0, 66.666667), ), ((55.0,
-1.666667, 66.666667), ), ((58.333333, 5.0, 66.666667), ), )))
mdb.models['Model-1'].predefinedFields['Predefined Field-1'].setValues(
magnitudes=(60.0, ))
mdb.models['Model-1'].parts['Part-1'].seedPart(deviationFactor=0.1,
minSizeFactor=0.1, size=10.0)
mdb.models['Model-1'].parts['Part-1'].setElementType(elemTypes=(ElemType(
elemCode=DC3D8, elemLibrary=STANDARD), ElemType(elemCode=DC3D6,
elemLibrary=STANDARD), ElemType(elemCode=DC3D4, elemLibrary=STANDARD)),
regions=(mdb.models['Model-1'].parts['Part-1'].cells.findAt(((-5.0,
-1.666667, 66.666667), ), ), ))
mdb.models['Model-1'].parts['Part-1'].generateMesh()
mdb.models['Model-1'].rootAssembly.regenerate()
# Save by soroush on 2017_07_10-11.43.22; build 2016.HF2 2016_01_20-04.05.00 126859
from part import *
from material import *
from section import *
from optimization import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
mdb.Job(atTime=None, contactPrint=OFF, description='', echoPrint=OFF,
explicitPrecision=SINGLE, getMemoryFromAnalysis=True, historyPrint=OFF,
memory=90, memoryUnits=PERCENTAGE, model='Model-1', modelPrint=OFF,
multiprocessingMode=DEFAULT, name='Job-Heat-Transfer-1',
nodalOutputPrecision=SINGLE, numCpus=1, numGPUs=0, queue=None,
resultsFormat=ODB, scratch='', type=ANALYSIS, userSubroutine='', waitHours=
0, waitMinutes=0)
mdb.jobs['Job-Heat-Transfer-1'].submit(consistencyChecking=OFF)
# Save by soroush on 2017_07_10-11.43.23; build 2016.HF2 2016_01_20-04.05.00 126859
from part import *
from material import *
from section import *
from optimization import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
mdb.jobs['Job-Heat-Transfer-1']._Message(STARTED, {'phase': BATCHPRE_PHASE,
'clientHost': 'Soroush', 'handle': 0, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(ODB_FILE, {'phase': BATCHPRE_PHASE,
'file': 'C:\\ali_project\\temperature\\Job-Heat-Transfer-1.odb',
'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(COMPLETED, {'phase': BATCHPRE_PHASE,
'message': 'Analysis phase complete', 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(STARTED, {'phase': STANDARD_PHASE,
'clientHost': 'Soroush', 'handle': 8004, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(STEP, {'phase': STANDARD_PHASE,
'stepId': 1, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(ODB_FRAME, {'phase': STANDARD_PHASE,
'step': 0, 'frame': 0, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(WARNING, {'phase': STANDARD_PHASE,
'message': 'THERE ARE 2 UNCONNECTED REGIONS IN THE MODEL.',
'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(MEMORY_ESTIMATE, {
'phase': STANDARD_PHASE, 'jobName': 'Job-Heat-Transfer-1', 'memory': 24.0})
mdb.jobs['Job-Heat-Transfer-1']._Message(WARNING, {'phase': STANDARD_PHASE,
'message': 'There is zero HEAT FLUX everywhere in the model based on the default criterion. please check the value of the average HEAT FLUX during the current iteration to verify that the HEAT FLUX is small enough to be treated as zero. if not, please use the solution controls to reset the criterion for zero HEAT FLUX.',
'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(ODB_FRAME, {'phase': STANDARD_PHASE,
'step': 0, 'frame': 1, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(STATUS, {'totalTime': 1.0,
'attempts': 1, 'timeIncrement': 1.0, 'increment': 1, 'stepTime': 1.0,
'step': 1, 'jobName': 'Job-Heat-Transfer-1', 'severe': 0, 'iterations': 1,
'phase': STANDARD_PHASE, 'equilibrium': 1})
mdb.jobs['Job-Heat-Transfer-1']._Message(END_STEP, {'phase': STANDARD_PHASE,
'stepId': 1, 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(COMPLETED, {'phase': STANDARD_PHASE,
'message': 'Analysis phase complete', 'jobName': 'Job-Heat-Transfer-1'})
mdb.jobs['Job-Heat-Transfer-1']._Message(JOB_COMPLETED, {
'time': 'Mon Jul 10 11:43:33 2017', 'jobName': 'Job-Heat-Transfer-1'})
# Save by soroush on 2017_07_10-11.44.02; build 2016.HF2 2016_01_20-04.05.00 126859
# Save by soroush on 2017_07_10-11.44.07; build 2016.HF2 2016_01_20-04.05.00 126859
from part import *
from material import *
from section import *
from optimization import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
mdb.saveAs(pathName='C:/ali_project/temperature/test of simple Heat transfer Analysis-Final.cae')
Raw
input.txt
Node Number X Y Z Temp
1 -5.00 -5.00 100.00 0.086
2 -5.00 5.00 100.00 0.0857
3 -5.00 -5.00 90.00 0.0862
4 -5.00 5.00 90.00 0.0852
5 -5.00 -5.00 80.00 0.0860
6 -5.00 5.00 80.00 0.0852
7 -5.00 -5.00 70.00 0.0862
8 -5.00 5.00 70.00 0.0849
9 -5.00 -5.00 60.00 0.0845
10 -5.00 5.00 60.00 0.0858
11 -5.00 -5.00 50.00 0.0851
12 -5.00 5.00 50.00 0.0860
13 -5.00 -5.00 40.00 0.0858
14 -5.00 5.00 40.00 0.0853
15 -5.00 -5.00 30.00 0.0859
16 -5.00 5.00 30.00 0.0843
17 -5.00 -5.00 20.00 0.0861
18 -5.00 5.00 20.00 0.0850
19 -5.00 -5.00 10.00 0.0862
20 -5.00 5.00 10.00 0.0850
21 -5.00 -5.00 0.00 0.0862
22 -5.00 5.00 0.00 0.0849
23 5.00 -5.00 100.00 0.0867
24 5.00 5.00 100.00 0.0869
25 5.00 -5.00 90.00 0.0873
26 5.00 5.00 90.00 0.0867
27 5.00 -5.00 80.00 0.0877
28 5.00 5.00 80.00 0.0862
29 5.00 -5.00 70.00 0.0875
30 5.00 5.00 70.00 0.0862
31 5.00 -5.00 60.00 0.0891
32 5.00 5.00 60.00 0.0871
33 5.00 -5.00 50.00 0.0885
34 5.00 5.00 50.00 0.0874
35 5.00 -5.00 40.00 0.0871
36 5.00 5.00 40.00 0.0862
37 5.00 -5.00 30.00 0.0874
38 5.00 5.00 30.00 0.0875
39 5.00 -5.00 20.00 0.0876
40 5.00 5.00 20.00 0.0862
41 5.00 -5.00 10.00 0.0875
42 5.00 5.00 10.00 0.0862
43 5.00 -5.00 0.00 0.0875
44 5.00 5.00 0.00 0.0862
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment