Import required packages and modules
# -*- coding: mbcs -*-
# Do not delete the following import lines
from abaqus import *
from abaqusConstants import *
import __main__
import section
import regionToolset
import displayGroupMdbToolset as dgm
import part
import material
import assembly
import step
import interaction
import load
import mesh
import optimization
import job
import sketch
import visualization
import xyPlot
import displayGroupOdbToolset as dgo
import connectorBehavior
Defines the range and interval of variables
specimen_length=0.18
specimen_width=0.12
specimen_hight=0.017
for bead_length in (0.05,0.055,0.06,0.065,0.07)
for bead_width in (0.0055,0.006,0.0065,0.007,0.0075)
for bead_hight in (0.0012,0.0014,0.0016,0.0018,0.002)
heat_input=int(1.1e+9)
Parametric modeling
s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=0.1)
g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints
s.sketchOptions.setValues(decimalPlaces=3)
s.setPrimaryObject(option=STANDALONE)
s.rectangle(point1=(specimen_length/2, specimen_width/2), point2=(-specimen_length/2, 0.0))
p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=THREE_D,
type=DEFORMABLE_BODY)
p = mdb.models['Model-1'].parts['Part-1']
p.BaseSolidExtrude(sketch=s, depth=specimen_hight)
s.unsetPrimaryObject()
p = mdb.models['Model-1'].parts['Part-1']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Model-1'].sketches['__profile__']
Define material properties, including density, elasticity, plasticity, specific heat, conductivity and expension rate.
session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=ON,
engineeringFeatures=ON)
session.viewports['Viewport: 1'].partDisplay.geometryOptions.setValues(
referenceRepresentation=OFF)
mdb.models['Model-1'].Material(name='AISI316L')
mdb.models['Model-1'].materials['AISI316L'].Density(table=((7966.0, ), ))
mdb.models['Model-1'].materials['AISI316L'].Elastic(temperatureDependency=ON,
table=((195600000000, 0.294, 20.0), (191200000000, 0.294, 100.0), (
185700000000, 0.294, 200.0), (179600000000, 0.294, 300.0), (
172600000000, 0.294, 400.0), (164500000000, 0.294, 500.0), (
155000000000, 0.294, 600.0), (144100000000, 0.294, 700.0), (
131400000000, 0.294, 800.0), (116800000000, 0.294, 900.0), (
100000000000, 0.294, 1000.0), (80000000000, 0.294, 1100.0), (
57000000000, 0.294, 1200.0), (30000000000, 0.294, 1300.0), (2000000000, 0.294, 1400.0)))
mdb.models['Model-1'].materials['AISI316L'].Plastic(temperatureDependency=ON,
table=((269119973.6, 0.0, 0.0), (248993663.6, 0.0, 100.0), (
228867353.5, 0.0, 200.0), (208741043.4, 0.0, 300.0), (188614733.4, 0.0,
400.0), (168488423.3, 0.0, 500.0), (148362113.3, 0.0, 600.0), (
128235803.2, 0.0, 700.0), (108109493.1, 0.0, 800.0), (87983183.08, 0.0,
900.0)))
mdb.models['Model-1'].materials['AISI316L'].SpecificHeat(
temperatureDependency=ON, table=((492, 20.0), (502,
100.0), (514, 200.0), (526, 300.0), (538,
400.0), (550, 500.0), (562, 600.0), (575,
700.0), (587, 800.0), (599, 900.0), (611,
1000.0), (623, 1100.0), (635, 1200), (647, 1300.0), (659, 1400.0)))
mdb.models['Model-1'].materials['AISI316L'].Conductivity(
temperatureDependency=ON, table=((14.12, 20.0), (15.26,
100.0), (16.69, 200.0), (18.11, 300.0), (19.54,
400.0), (20.96, 500.0), (22.38, 600.0), (23.81,
700.0), (25.23, 800.0), (26.66, 900.0), (28.08,
1000.0), (29.50, 1100.0), (30.93, 1200.0), (32.35, 1300.0), (33.78, 1400.0)))
mdb.models['Model-1'].materials['AISI316L'].Expansion(table=((1.456e-05, 20.0),
(1.539e-05, 100.0), (1.621e-05, 200.0), (1.686e-05, 300.0), (
1.737e-05, 400.0), (1.778e-05, 500.0), (1.812e-05, 600.0), (
1.843e-05, 700.0), (1.872e-05, 800.0), (1.899e-05, 900.0), (
1.927e-05, 1000.0), (1.953e-05, 1100.0), (1.979e-05, 1200.0), (
2.002e-05, 1300.0), (
155-2731-8020