10-31-06-Python-Scripting-in-ABAQUS (2)

mdb.models['Model-1'].Part(dimensionality=TWO_D_PLANAR, name='Part-1', type= DEFORMABLE_BODY)

mdb.models['Model-1'].parts['Part-1'].BaseShell(sketch= sketch1)

del mdb.models['Model-1'].sketches['__profile__'] part1 = mdb.models['Model-1'].parts['Part-1']

#Create partition at corner

mdb.models['Model-1'].ConstrainedSketch(gridSpacing=0.03, name='__profile__', sheetSize=1.51, transform= part1.MakeSketchTransform(

sketchPlane=part1.faces.findAt((0,

0, 0.0), (0.0, 0.0, 1.0)), sketchPlaneSide=SIDE1, sketchOrientation=RIGHT, origin=(0, 0, 0.0)))

sketch1 = mdb.models['Model-1'].sketches['__profile__']

mdb.models['Model-1'].sketches['__profile__'].Line(point1=(r1, 0.0), point2=(r1,r1))

mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0, r1), point2=(r1,r1))

part1.PartitionFaceBySketch(faces= part1.faces.findAt(((0, 0,

0.0), (0.0, 0.0, 1.0)), ), sketch=

mdb.models['Model-1'].sketches['__profile__']) del mdb.models['Model-1'].sketches['__profile__']

#Sketch the partition to make the fiber

mdb.models['Model-1'].ConstrainedSketch(gridSpacing=0.35, name='__profile__', sheetSize=14.14, transform= part1.MakeSketchTransform(

sketchPlane=part1.faces.findAt((fiber_radius,

fiber_radius, 0.0), (0.0, 0.0, 1.0)), sketchPlaneSide=SIDE1, sketchOrientation=RIGHT, origin=(0, 0, 0.0)))

sketch1 = mdb.models['Model-1'].sketches['__profile__'] sketch1.sketchOptions.setValues( gridSpacing=0.35)

part1.projectReferencesOntoSketch(filter= COPLANAR_EDGES, sketch=sketch1) sketch1.ArcByCenterEnds(center=(0,

0), direction=COUNTERCLOCKWISE, point1=(fiber_radius,0), point2=(0,fiber_radius))

sketch1.CoincidentConstraint(entity1=

sketch1.vertices.findAt((fiber_radius, 0), ), entity2=sketch1.geometry.findAt(( fiber_radius, 0), ))

sketch1.CoincidentConstraint(entity1=

sketch1.vertices.findAt((0, fiber_radius),

), entity2=sketch1.geometry.findAt(( 0, fiber_radius), ))

part1.PartitionFaceBySketch(faces=

part1.faces.findAt(((fiber_radius, fiber_radius, 0.0), (0.0, 0.0, 1.0)), ), sketch= sketch1)

del mdb.models['Model-1'].sketches['__profile__']

#Create outer partition

mdb.models['Model-1'].ConstrainedSketch(gridSpacing=0.03, name='__profile__', sheetSize=1.51, transform= part1.MakeSketchTransform(

sketchPlane=part1.faces.findAt((length*c30,0,0),(0,0,1)), sketchPlaneSide=SIDE1, sketchOrientation=RIGHT, origin=(0, 0, 0.0))) part1.projectReferencesOntoSketch(filter=

COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__']) mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=( 0, 0), point1=(r2, 0))

part1.PartitionFaceBySketch(faces= part1.faces.findAt(((length*c30, 0, 0.0), (0.0, 0.0, 1.0)), ), sketch=

mdb.models['Model-1'].sketches['__profile__']) del mdb.models['Model-1'].sketches['__profile__']

#Instance the part

mdb.models['Model-1'].rootAssembly.DatumCsysByDefault(CARTESIAN)

mdb.models['Model-1'].rootAssembly.Instance(dependent=ON, name='Part-1-1', part=part1)

#Seed and mesh the part

part1.seedPart(deviationFactor=0.1, size=element_size) part1.seedEdgeByNumber(edges=

part1.edges.findAt(((r2*c30, r2*s30, 0.0), ), ((fiber_radius*c30, fiber_radius*s30, 0.0),) , ),

number=40,constraint=FINER) part1.seedEdgeByNumber(edges=

part1.edges.findAt(((c30*(length-length/3), length/2+length/3*s30, 0.0),

), ), number=edgeseed,constraint=FIXED)

part1.setMeshControls(regions=

part1.faces.findAt(((.9*fiber_radius*c30, .9*fiber_radius*s30, 0.0), (0.0, 0.0, 1.0)), ((.9*r2*c30, .9*r2*s30, 0.0), (0.0, 0.0, 1.0)), (( length*c30, 0, 0.0), (0.0, 0.0, 1.0)), ), technique=STRUCTURED) part1.generateMesh()