找回密码
 注册
Simdroid-非首页
查看: 220|回复: 5

[接触] 接触不收敛,改刚度又出现渗透问题

[复制链接]
发表于 2012-11-18 21:36:37 | 显示全部楼层 |阅读模式 来自 北京
悬赏50仿真币未解决
本帖最后由 lucymomo 于 2012-11-19 20:48 编辑

接触问题太高深了,求大牛指点
唉材料非线性,又有接触,怎么都收敛不了,调整了接触刚度到e-5次方,收敛了,可是渗透又很严重。不知道咋搞了
传不上去啊,才50k的文件都传不上去
模型中主要是变形体的钢管,外面是刚体面,钢管中加内压膨胀,收外面的刚体面约束,接触很不容易收敛。
接触设置的切向罚函数,slip tolerance=0.25 ,摩擦系数0.3.接触控制中设置接触刚度为1e-5才能收敛,但是渗透率接近7mm,太大了,但是又不能直接调整绝对渗透量,调整了就又不收敛
任何附件图片都传不上来,我贴py吧,版本是6.11的

# -*- coding: mbcs -*-
from part import *
from material import *
from section 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 *
import math
l=30000*1.84                         #管的长度
lo=10000                          #管取的长度
angle=30*1.84                        #弯曲的角度
angle2=20                           #曲面的角度
a=math.sin(angle*pi/180)        #角度的正弦
b=math.cos(angle*pi/180)        #角度的余弦
a2=math.sin(angle2*pi/180)        #角度的正弦
b2=math.cos(angle2*pi/180)        #角度的余弦
R=180*l/(pi*angle)
wf=R*(1-b)
height=R*a
r2=180*wf/(pi*angle2)
width=r2*a2
angle3=math.atan(width/height)
angle4=2*angle3
a3=math.sin(angle4*pi/180)
b3=math.cos(angle4*pi/180)
r=0.5*(width*width+height*height)/width              #弯曲半径
h1=0.5*width                         #高度余量
depth=12000                      #墩的厚度
R1=165*0.5                     #圆孔半径

#stone part 建模
mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=2000.0)
mdb.models['Model-1'].sketches['__profile__'].rectangle(point1=(0.0, 0.0),
    point2=(width*2, height+h1))
mdb.models['Model-1'].Part(dimensionality=THREE_D, name='Part-1', type=
    DISCRETE_RIGID_SURFACE)
mdb.models['Model-1'].parts['Part-1'].BaseSolidExtrude(depth=depth, sketch=
    mdb.models['Model-1'].sketches['__profile__'])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=52.36, name='__profile__',
    sheetSize=2094.44, transform=
    mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
    sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[1],
    sketchPlaneSide=SIDE1,
    sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[5],
    sketchOrientation=RIGHT, origin=(width, height+h1, 0.5*depth)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    0.0, -width*b2/a2), point1=(0.0, width/a2-width*b2/a2))
mdb.models['Model-1'].parts['Part-1'].PartitionFaceBySketch(faces=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(('[#2 ]',  
    ), ), sketch=mdb.models['Model-1'].sketches['__profile__'], sketchUpEdge=
    mdb.models['Model-1'].parts['Part-1'].edges[5])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=3688.76, name='__profile__'
    , sheetSize=147550.64, transform=
    mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
    sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[5],
    sketchPlaneSide=SIDE1,
    sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[11],
    sketchOrientation=RIGHT, origin=(0.0, 0.0, depth)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    r+width, height+h1), point1=(width, height+h1))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[5], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[3])
mdb.models['Model-1'].sketches['__profile__'].EqualDistanceConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[1], entity2=
    mdb.models['Model-1'].sketches['__profile__'].vertices[2], midpoint=
    mdb.models['Model-1'].sketches['__profile__'].vertices[5])
mdb.models['Model-1'].parts['Part-1'].PartitionFaceBySketch(faces=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(('[#20 ]',
    ), ), sketch=mdb.models['Model-1'].sketches['__profile__'], sketchUpEdge=
    mdb.models['Model-1'].parts['Part-1'].edges[11])
del mdb.models['Model-1'].sketches['__profile__']

mdb.models['Model-1'].parts['Part-1'].PartitionCellByExtrudeEdge(cells=
    mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#1 ]',
    ), ), edges=(mdb.models['Model-1'].parts['Part-1'].edges[5],
    mdb.models['Model-1'].parts['Part-1'].edges[6]), line=
    mdb.models['Model-1'].parts['Part-1'].edges[15], sense=FORWARD)
mdb.models['Model-1'].parts['Part-1'].PartitionCellByExtrudeEdge(cells=
    mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#1 ]',
    ), ), edges=(mdb.models['Model-1'].parts['Part-1'].edges[14], ), line=
    mdb.models['Model-1'].parts['Part-1'].edges[11], sense=REVERSE)

mdb.models['Model-1'].ConstrainedSketch(gridSpacing=3857.4, name='__profile__',
    sheetSize=154296.34, transform=(1.0, -1.22464679914735e-16, 0.0,
    9.18338008663387e-49, 7.49879891330929e-33, -1.0, 1.22464679914735e-16,
    1.0, 7.49879891330929e-33, width, height+h1,
    0.5*depth-r2+r2*b2))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -77148.17, 0.0), point2=(77148.17, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -77148.17), point2=(0.0, 77148.17))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    0.0, 0.0), point1=(0.0, R1))
mdb.models['Model-1'].parts['Part-1'].CutSweep(flipSweepDirection=ON, path=
    mdb.models['Model-1'].parts['Part-1'].edges.getSequenceFromMask(('[#1 ]',
    ), ), profile=mdb.models['Model-1'].sketches['__profile__'],
    sketchOrientation=RIGHT, sketchUpEdge=
    mdb.models['Model-1'].parts['Part-1'].edges[8])
del mdb.models['Model-1'].sketches['__profile__']

mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2679.59, name='__profile__'
    , sheetSize=107183.84, transform=
    mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
    sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[10],
    sketchPlaneSide=SIDE1,
    sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[44],
    sketchOrientation=RIGHT, origin=(2*width, height+h1, depth)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].rectangle(point1=(0.0,
    0), point2=(depth, -height-h1))
mdb.models['Model-1'].parts['Part-1'].CutExtrude(depth=2000.0,
    flipExtrudeDirection=OFF, sketch=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=RIGHT,
    sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[10],
    sketchPlaneSide=SIDE1, sketchUpEdge=
    mdb.models['Model-1'].parts['Part-1'].edges[44])
del mdb.models['Model-1'].sketches['__profile__']

mdb.models['Model-1'].parts['Part-1'].RemoveCells(cellList=
    mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(mask=(
    '[#7 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].checkGeometry()

#pipe part 建模 (分段前)

mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=200.0)
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    0.0, 29.665), point1=(0.0, 6.0))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, 29.665),
    point2=(0.0, 6.0))
mdb.models['Model-1'].sketches['__profile__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__profile__'].geometry[3])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, 29.665),
    point2=(22.9138604342438, 23.75))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[2], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[2])
mdb.models['Model-1'].sketches['__profile__'].autoDimension(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[3],
    mdb.models['Model-1'].sketches['__profile__'].geometry[4]))
mdb.models['Model-1'].sketches['__profile__'].dimensions[1].setValues(value=70)
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(22.2378258708986,
    21.571093308198), point2=(26.2156225995849, 32.5))
mdb.models['Model-1'].sketches['__profile__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[4], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[5])
mdb.models['Model-1'].sketches['__profile__'].FixedConstraint(entity=
    mdb.models['Model-1'].sketches['__profile__'].vertices[2])
mdb.models['Model-1'].sketches['__profile__'].autoDimension(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[5], ))
mdb.models['Model-1'].sketches['__profile__'].dimensions[2].setValues(value=
    17.881)
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(28.3534880537049,
    38.3737370604709), point2=(40.8809422007762, 33.8141166397836))
mdb.models['Model-1'].sketches['__profile__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[5], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[6])
mdb.models['Model-1'].sketches['__profile__'].FixedConstraint(entity=
    mdb.models['Model-1'].sketches['__profile__'].vertices[3])
mdb.models['Model-1'].sketches['__profile__'].autoDimension(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[6], ))
mdb.models['Model-1'].sketches['__profile__'].dimensions[3].setValues(value=
    11.5)
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    39.1599531927428, 34.4405054122257), point1=(28.3534880537049,
    38.3737370604709))
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
    0.0, 0.0), point1=(63.65, 0.0))
mdb.models['Model-1'].sketches['__profile__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].dimensions[0],
    mdb.models['Model-1'].sketches['__profile__'].dimensions[1]))
mdb.models['Model-1'].sketches['__profile__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].dimensions[2],
    mdb.models['Model-1'].sketches['__profile__'].dimensions[3]))
mdb.models['Model-1'].sketches['__profile__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].constraints[17], ))
mdb.models['Model-1'].sketches['__profile__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].constraints[13], ))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[2], point1=(
    22.6671886444092, 35.9677124023438))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[4], point1=(
    16.4866161346436, 23.5132789611816))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[3], point1=(
    -0.311861038208008, 19.7055549621582))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[6], point1=(
    32.0172805786133, 36.760986328125))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[7], point1=(
    28.8477611541748, 29.0662097930908))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8], point1=(
    42.6351852416992, 47.3908843994141))
mdb.models['Model-1'].sketches['__profile__'].autoTrimCurve(curve1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[10], point1=(
    47.7762260437012, 42.0553894042969))
mdb.models['Model-1'].ConstrainedSketch(name='Sketch-1', objectToCopy=
    mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].offset(distance=7.5, objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[5],
    mdb.models['Model-1'].sketches['__profile__'].geometry[9],
    mdb.models['Model-1'].sketches['__profile__'].geometry[11],
    mdb.models['Model-1'].sketches['__profile__'].geometry[12]), side=LEFT)
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(angle=0.0,
    point1=(63.65, 0.0))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[9], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[17])
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[17])
mdb.models['Model-1'].sketches['__profile__'].copyMirror(mirrorLine=
    mdb.models['Model-1'].sketches['__profile__'].geometry[17], objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[5],
    mdb.models['Model-1'].sketches['__profile__'].geometry[9],
    mdb.models['Model-1'].sketches['__profile__'].geometry[11],
    mdb.models['Model-1'].sketches['__profile__'].geometry[12],
    mdb.models['Model-1'].sketches['__profile__'].geometry[13],
    mdb.models['Model-1'].sketches['__profile__'].geometry[14],
    mdb.models['Model-1'].sketches['__profile__'].geometry[15],
    mdb.models['Model-1'].sketches['__profile__'].geometry[16]))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, 13.5), point2=(
    0.0, 6.0))
mdb.models['Model-1'].sketches['__profile__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__profile__'].geometry[26])
mdb.models['Model-1'].sketches['__profile__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[13], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[26])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, -6.0), point2=(
    0.0, -13.5))
mdb.models['Model-1'].sketches['__profile__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__profile__'].geometry[27])
mdb.models['Model-1'].sketches['__profile__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].geometry[19], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[27])
mdb.models['Model-1'].Part(dimensionality=THREE_D, name='pipe', type=
    DEFORMABLE_BODY)
mdb.models['Model-1'].parts['pipe'].BaseSolidExtrude(depth=lo, sketch=
    mdb.models['Model-1'].sketches['__profile__'])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['pipe'].features['Solid extrude-1'].setValues(
    depth=50000.0)
mdb.models['Model-1'].parts['pipe'].regenerate()
mdb.models['Model-1'].ConstrainedSketch(name='__edit__', objectToCopy=
    mdb.models['Model-1'].parts['pipe'].features['Solid extrude-1'].sketch)
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__edit__'],
    upToFeature=
    mdb.models['Model-1'].parts['pipe'].features['Solid extrude-1'])
mdb.models['Model-1'].sketches['__edit__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__edit__'].geometry[18],
    mdb.models['Model-1'].sketches['__edit__'].geometry[19],
    mdb.models['Model-1'].sketches['__edit__'].geometry[20],
    mdb.models['Model-1'].sketches['__edit__'].geometry[21],
    mdb.models['Model-1'].sketches['__edit__'].geometry[22],
    mdb.models['Model-1'].sketches['__edit__'].geometry[23],
    mdb.models['Model-1'].sketches['__edit__'].geometry[24],
    mdb.models['Model-1'].sketches['__edit__'].geometry[25],
    mdb.models['Model-1'].sketches['__edit__'].geometry[26],
    mdb.models['Model-1'].sketches['__edit__'].geometry[27],
    mdb.models['Model-1'].sketches['__edit__'].constraints[83],
    mdb.models['Model-1'].sketches['__edit__'].constraints[85],
    mdb.models['Model-1'].sketches['__edit__'].constraints[86],
    mdb.models['Model-1'].sketches['__edit__'].constraints[88],
    mdb.models['Model-1'].sketches['__edit__'].constraints[90],
    mdb.models['Model-1'].sketches['__edit__'].constraints[93]))
mdb.models['Model-1'].sketches['__edit__'].ConstructionLine(angle=90.0, point1=
    (0.0, 0.0))
mdb.models['Model-1'].sketches['__edit__'].CoincidentConstraint(addUndoState=
    False, entity1=mdb.models['Model-1'].sketches['__edit__'].vertices[11],
    entity2=mdb.models['Model-1'].sketches['__edit__'].geometry[28])
mdb.models['Model-1'].sketches['__edit__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__edit__'].geometry[28])
mdb.models['Model-1'].sketches['__edit__'].copyMirror(mirrorLine=
    mdb.models['Model-1'].sketches['__edit__'].geometry[28], objectList=(
    mdb.models['Model-1'].sketches['__edit__'].geometry[5],
    mdb.models['Model-1'].sketches['__edit__'].geometry[9],
    mdb.models['Model-1'].sketches['__edit__'].geometry[11],
    mdb.models['Model-1'].sketches['__edit__'].geometry[12],
    mdb.models['Model-1'].sketches['__edit__'].geometry[13],
    mdb.models['Model-1'].sketches['__edit__'].geometry[14],
    mdb.models['Model-1'].sketches['__edit__'].geometry[15],
    mdb.models['Model-1'].sketches['__edit__'].geometry[16]))
mdb.models['Model-1'].sketches['__edit__'].Line(point1=(-71.15, 0.0), point2=(
    -63.65, 0.0))
mdb.models['Model-1'].sketches['__edit__'].HorizontalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__edit__'].geometry[37])
mdb.models['Model-1'].sketches['__edit__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__edit__'].geometry[36], entity2=
    mdb.models['Model-1'].sketches['__edit__'].geometry[37])
mdb.models['Model-1'].sketches['__edit__'].Line(point1=(63.65, 0.0), point2=(
    71.15, 0.0))
mdb.models['Model-1'].sketches['__edit__'].HorizontalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__edit__'].geometry[38])
mdb.models['Model-1'].sketches['__edit__'].PerpendicularConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__edit__'].geometry[11], entity2=
    mdb.models['Model-1'].sketches['__edit__'].geometry[38])
mdb.models['Model-1'].parts['pipe'].features['Solid extrude-1'].setValues(
    sketch=mdb.models['Model-1'].sketches['__edit__'])
del mdb.models['Model-1'].sketches['__edit__']
mdb.models['Model-1'].parts['pipe'].regenerate()

#分段
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].vertices[28], vector=(0.0, 0.0,
    -9994.07))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].vertices[28], vector=(0.0, 0.0,
    -10005.93))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].vertices[30], vector=(0.0, 0.0,
    -9998.4))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].vertices[30], vector=(0.0, 0.0,
    -10001.6))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0,
    0.0, 67.4, 0.0, 25000.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-15005.93, 3.75),
    point2=(-14994.07, 3.75))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-14994.07, 3.75),
    point2=(-14998.4, -3.75))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-14998.4, -3.75),
    point2=(-15001.6, -3.75))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[14])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-15001.6, -3.75),
    point2=(-15005.93, 3.75))
mdb.models['Model-1'].parts['pipe'].CutSweep(flipSweepDirection=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#89240000 #49 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[12], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[39])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, -0.0, -1.0, 1.0, 0.0, -0.0, 0.0,
    -1.0, 0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(4995.95242185,
    3.7502707269287), point2=(5007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[20])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[20], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[21], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(5007.81242185,
    3.7502707269287), point2=(5003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[23], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(5003.48242185,
    -3.7497292730713), point2=(5000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[22])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[25], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(5000.28242185,
    -3.7497292730713), point2=(4995.95242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[27], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(1.0, 0.0, -0.0, 0.0, 0.0, 1.0, 0.0, -1.0,
    0.0, 71.15, 0.0, 50000.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-7.5, -9998.4),
    point2=(0.0, -9994.07))
mdb.models['Model-1'].sketches['__profile__'].ParallelConstraint(addUndoState=
    False, entity1=mdb.models['Model-1'].sketches['__profile__'].geometry[12],
    entity2=mdb.models['Model-1'].sketches['__profile__'].geometry[20])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[20], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[21], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, -9994.07),
    point2=(0.0, -10005.93))
mdb.models['Model-1'].sketches['__profile__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__profile__'].geometry[21])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[23], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0.0, -10005.93),
    point2=(-7.5, -10001.6))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[25], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(-7.5, -10001.6),
    point2=(-7.5, -9998.4))
mdb.models['Model-1'].sketches['__profile__'].VerticalConstraint(addUndoState=
    False, entity=mdb.models['Model-1'].sketches['__profile__'].geometry[23])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[27], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[12])
mdb.models['Model-1'].parts['pipe'].SolidSweep(keepInternalBoundaries=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:2 #1fe ]', ), ), profile=mdb.models['Model-1'].sketches['__profile__']
    , sketchOrientation=RIGHT, sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[78])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[2], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[3], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[4], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[5], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, -0.0, -1.0, 1.0, 0.0, -0.0, 0.0,
    -1.0, 0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(14995.95242185,
    3.7502707269287), point2=(15007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[24])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15007.81242185,
    3.7502707269287), point2=(15000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].delete(objectList=(
    mdb.models['Model-1'].sketches['__profile__'].geometry[25], ))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15007.81242185,
    3.7502707269287), point2=(15003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15003.48242185,
    -3.7497292730713), point2=(15000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[27])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15000.28242185,
    -3.7497292730713), point2=(14995.95242185, 3.7502707269287))
mdb.models['Model-1'].parts['pipe'].CutSweep(path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:3 #1fe0 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[44], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[113])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0,
    0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(14995.95242185,
    3.7502707269287), point2=(15007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[32])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[32], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[20])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[33], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15007.81242185,
    3.7502707269287), point2=(15003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[35], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15003.48242185,
    -3.7497292730713), point2=(15000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[34])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[37], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[20])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(15000.28242185,
    -3.7497292730713), point2=(14995.95242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[39], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[20])
mdb.models['Model-1'].parts['pipe'].SolidSweep(flipSweepDirection=ON,
    keepInternalBoundaries=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:3 #ab540 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[48], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[116])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[8], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[9], vector=(0.0, 0.0, -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[10], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[11], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0,
    0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(24995.95242185,
    3.7502707269287), point2=(25007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[36])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25007.81242185,
    3.7502707269287), point2=(25003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25003.48242185,
    -3.7497292730713), point2=(25000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[38])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25000.28242185,
    -3.7497292730713), point2=(24995.95242185, 3.7502707269287))
mdb.models['Model-1'].parts['pipe'].CutSweep(flipSweepDirection=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:4 #ab5400 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[66], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[152])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, -0.0, -1.0, 1.0, 0.0, -0.0, 0.0,
    -1.0, 0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(24995.95242185,
    3.7502707269287), point2=(25007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[44])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[44], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[32])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[45], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25007.81242185,
    3.7502707269287), point2=(25003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[47], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25003.48242185,
    -3.7497292730713), point2=(25000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[46])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[49], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[32])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(25000.28242185,
    -3.7497292730713), point2=(24995.95242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[51], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[32])
mdb.models['Model-1'].parts['pipe'].SolidSweep(keepInternalBoundaries=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:5 #7f80 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[76], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[179])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[14], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[15], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[16], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].parts['pipe'].DatumPointByOffset(point=
    mdb.models['Model-1'].parts['pipe'].datums[17], vector=(0.0, 0.0,
    -10000.0))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, -0.0, -1.0, 1.0, 0.0, -0.0, 0.0,
    -1.0, 0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(34995.95242185,
    3.7502707269287), point2=(35007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[48])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35007.81242185,
    3.7502707269287), point2=(35003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35003.48242185,
    -3.7497292730713), point2=(35000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[50])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35000.28242185,
    -3.7497292730713), point2=(34995.95242185, 3.7502707269287))
mdb.models['Model-1'].parts['pipe'].CutSweep(path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:6 #7f800 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[94], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[215])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=2500.01, name='__profile__'
    , sheetSize=100000.46, transform=(0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0,
    0.0, 67.3997292730713, 0.0, 45001.88242185))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(
    -50000.23, 0.0), point2=(50000.23, 0.0))
mdb.models['Model-1'].sketches['__profile__'].ConstructionLine(point1=(0.0,
    -50000.23), point2=(0.0, 50000.23))
mdb.models['Model-1'].parts['pipe'].projectReferencesOntoSketch(filter=
    COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(34995.95242185,
    3.7502707269287), point2=(35007.81242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[56])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[56], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[44])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[57], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35007.81242185,
    3.7502707269287), point2=(35003.48242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[59], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[8])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35003.48242185,
    -3.7497292730713), point2=(35000.28242185, -3.7497292730713))
mdb.models['Model-1'].sketches['__profile__'].HorizontalConstraint(
    addUndoState=False, entity=
    mdb.models['Model-1'].sketches['__profile__'].geometry[58])
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[61], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[44])
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(35000.28242185,
    -3.7497292730713), point2=(34995.95242185, 3.7502707269287))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
    addUndoState=False, entity1=
    mdb.models['Model-1'].sketches['__profile__'].vertices[63], entity2=
    mdb.models['Model-1'].sketches['__profile__'].geometry[44])
mdb.models['Model-1'].parts['pipe'].SolidSweep(flipSweepDirection=ON,
    keepInternalBoundaries=ON, path=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#0:6 #2ad5000 ]', ), ), profile=
    mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
    sketchPlane=mdb.models['Model-1'].parts['pipe'].faces[98], sketchUpEdge=
    mdb.models['Model-1'].parts['pipe'].edges[218])
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['pipe'].Set(cells=
    mdb.models['Model-1'].parts['pipe'].cells.getSequenceFromMask(('[#ba ]', ),
    ), name='noweld')
mdb.models['Model-1'].parts['pipe'].Set(cells=
    mdb.models['Model-1'].parts['pipe'].cells.getSequenceFromMask(('[#145 ]',
    ), ), name='weld')
mdb.models['Model-1'].parts['pipe'].Mirror(keepInternalBoundaries=ON,
    keepOriginal=ON, mirrorPlane=mdb.models['Model-1'].parts['pipe'].faces[76])
mdb.models['Model-1'].parts['pipe'].Set(cells=
    mdb.models['Model-1'].parts['pipe'].cells.getSequenceFromMask(('[#174ba ]',
    ), ), name='noweld')
mdb.models['Model-1'].parts['pipe'].Set(cells=
    mdb.models['Model-1'].parts['pipe'].cells.getSequenceFromMask(('[#28b45 ]',
    ), ), name='weld')
   
#组合
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='pipe-1', part=
    mdb.models['Model-1'].parts['pipe'])
mdb.models['Model-1'].rootAssembly.rotate(angle=90.0, axisDirection=(1.0, 0.0,
    0.0), axisPoint=(0.0, 0.0, 0.0), instanceList=('pipe-1', ))
mdb.models['Model-1'].rootAssembly.translate(instanceList=('pipe-1', ), vector=
    (width, 50000+height+h1, 0.5*depth-r2+r2*b2))
mdb.models['Model-1'].rootAssembly.rotate(angle=90.0, axisDirection=(0.0,
    50000.0, 0.0), axisPoint=(width, 50000+height+h1, 0.5*depth-r2+r2*b2),
    instanceList=('pipe-1', ))
#时间步
mdb.models['Model-1'].StaticStep(initialInc=0.5, maxInc=0.5, maxNumInc=10000,
    minInc=1e-08, name='Step-1', nlgeom=ON, previous='Initial')
mdb.models['Model-1'].StaticStep(initialInc=0.001, maxInc=0.5, maxNumInc=10000,
    minInc=1e-08, name='Step-2', previous='Step-1')
mdb.models['Model-1'].StaticStep(initialInc=0.5, maxInc=0.5, maxNumInc=10000,
    minInc=1e-08, name='Step-3', previous='Step-2')
mdb.models['Model-1'].StaticStep(initialInc=0.001, maxInc=0.5, maxNumInc=10000,
    minInc=1e-08, name='Step-4', previous='Step-3')

#属性
mdb.models['Model-1'].rootAssembly.regenerate()
mdb.models['Model-1'].Material(name='pipe')
mdb.models['Model-1'].materials['pipe'].Elastic(table=((210000.0, 0.3), ))
mdb.models['Model-1'].materials['pipe'].Plastic(table=((369.41973, 0.0), (
    415.89808, 0.02663), (467.56252, 0.05128), (520.65527, 0.0869), (568.97666,
    0.14544), (603.72026, 0.19959)))
mdb.models['Model-1'].Material(name='weld')
mdb.models['Model-1'].materials['weld'].Elastic(table=((210000.0, 0.3), ))
mdb.models['Model-1'].materials['weld'].Plastic(table=((434.0, 0.0), (678,
    0.2839)))
mdb.models['Model-1'].HomogeneousSolidSection(material='pipe', name='pipe',
    thickness=None)
mdb.models['Model-1'].HomogeneousSolidSection(material='weld', name='weld',
    thickness=None)
mdb.models['Model-1'].parts['pipe'].SectionAssignment(offset=0.0, offsetField=
    '', offsetType=MIDDLE_SURFACE, region=
    mdb.models['Model-1'].parts['pipe'].sets['noweld'], sectionName='pipe',
    thicknessAssignment=FROM_SECTION)
mdb.models['Model-1'].parts['pipe'].SectionAssignment(offset=0.0, offsetField=
    '', offsetType=MIDDLE_SURFACE, region=
    mdb.models['Model-1'].parts['pipe'].sets['weld'], sectionName='weld',
    thicknessAssignment=FROM_SECTION)


#接触
mdb.models['Model-1'].rootAssembly.regenerate()
mdb.models['Model-1'].ContactProperty('IntProp-1')
mdb.models['Model-1'].StdContactControl(name='ContCtrl-1',
    stiffnessScaleFactor=1e-05)
mdb.models['Model-1'].interactionProperties['IntProp-1'].TangentialBehavior(
    dependencies=0, directionality=ISOTROPIC, elasticSlipStiffness=None,
    formulation=PENALTY, fraction=0.25, maximumElasticSlip=FRACTION,
    pressureDependency=OFF, shearStressLimit=None, slipRateDependency=OFF,
    table=((0.3, ), ), temperatureDependency=OFF)
mdb.models['Model-1'].SurfaceToSurfaceContactStd(adjustMethod=NONE,
    clearanceRegion=None, contactControls='ContCtrl-1', createStepName=
    'Initial', datumAxis=None, initialClearance=OMIT, interactionProperty=
    'IntProp-1', master=Region(
    side1Faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
    mask=('[#c18 ]', ), )), name='Int-1', slave=Region(
    side1Faces=mdb.models['Model-1'].rootAssembly.instances['pipe-1'].faces.getSequenceFromMask(
    mask=('[#7f00a5a6 #69800000 #1fc029 #5007f7f0 #fc00054b #cb4a01 #14b4a00',
    ' #3f8 #e00a5a50 #8fe0000f #94b5003f #fe0000 #14b4a ]', ), )), sliding=
    FINITE, thickness=ON)
# 载荷
mdb.models['Model-1'].Pressure(amplitude=UNSET, createStepName='Step-2',
    distributionType=UNIFORM, field='', magnitude=-100.0, name='Load-1',
    region=Region(
    side1Faces=mdb.models['Model-1'].rootAssembly.instances['pipe-1'].faces.getSequenceFromMask(
    mask=('[#0:3 #800 #0:5 #20000000 ]', ), )))
mdb.models['Model-1'].Pressure(amplitude=UNSET, createStepName='Step-4',
    distributionType=UNIFORM, field='', magnitude=30.0, name='Load-2', region=
    Region(
    side1Faces=mdb.models['Model-1'].rootAssembly.instances['pipe-1'].faces.getSequenceFromMask(
    mask=('[#80005a59 #967f803f #e0000016 #a800000f #3fab4 #ff34b400 #b4b100',
    ' #f003f800 #5a58f #1fc000 #4b4a8000 #7f #1fe34b4 ]', ), )))
mdb.models['Model-1'].loads['Load-1'].setValuesInStep(magnitude=0.0, stepName=
    'Step-3')
mdb.models['Model-1'].DisplacementBC(amplitude=UNSET, createStepName='Step-1',
    distributionType=UNIFORM, fieldName='', fixed=OFF, localCsys=None, name=
    'fix-node', region=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['pipe-1'].faces.getSequenceFromMask(
    mask=('[#0 #40 #0:5 #40000 ]', ), )), u1=0.0, u2=UNSET, u3=0.0, ur1=UNSET,
    ur2=UNSET, ur3=UNSET)
mdb.models['Model-1'].boundaryConditions['fix-node'].deactivate('Step-3')
mdb.models['Model-1'].DisplacementBC(amplitude=UNSET, createStepName='Initial',
    distributionType=UNIFORM, fieldName='', localCsys=None, name='fix-z',
    region=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
    mask=('[#140 ]', ), )), u1=UNSET, u2=UNSET, u3=SET, ur1=UNSET, ur2=UNSET,
    ur3=UNSET)
mdb.models['Model-1'].DisplacementBC(amplitude=UNSET, createStepName='Initial',
    distributionType=UNIFORM, fieldName='', localCsys=None, name='fix-x',
    region=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
    mask=('[#30000 ]', ), )), u1=SET, u2=UNSET, u3=UNSET, ur1=UNSET, ur2=UNSET,
    ur3=UNSET)
mdb.models['Model-1'].DisplacementBC(amplitude=UNSET, createStepName='Initial',
    distributionType=UNIFORM, fieldName='', localCsys=None, name='fix-y',
    region=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
    mask=('[#40080 ]', ), )), u1=UNSET, u2=SET, u3=UNSET, ur1=UNSET, ur2=UNSET,
    ur3=UNSET)
mdb.models['Model-1'].DisplacementBC(amplitude=UNSET, createStepName='Step-1',
    distributionType=UNIFORM, fieldName='', fixed=OFF, localCsys=None, name=
    'BC', region=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['pipe-1'].faces.getSequenceFromMask(
    mask=('[#0 #40 #0:5 #40000 ]', ), )), u1=UNSET, u2=-1000.0, u3=UNSET, ur1=
    UNSET, ur2=UNSET, ur3=UNSET)
mdb.models['Model-1'].boundaryConditions['BC'].setValuesInStep(stepName=
    'Step-2', u2=-50000.0)
mdb.models['Model-1'].boundaryConditions['BC'].deactivate('Step-3')

#PIPE网格
mdb.models['Model-1'].parts['pipe'].seedPart(deviationFactor=0.1,
    minSizeFactor=0.1, size=500.0)
mdb.models['Model-1'].parts['pipe'].seedEdgeByNumber(constraint=FINER, edges=
    mdb.models['Model-1'].parts['pipe'].edges.getSequenceFromMask((
    '[#ffffffff #7ffffff #e0000000 #6abfffff #ffffffd5 #1ffffff #3fffffe',
    ' #a0000000 #fffff55a #ff807fff #ffff55aa #f807ffff #ffffffff #1ffffff',
    ' #ffe00000 #ff5aafff #ffffffff #ffffe07f #f #fffff5aa #aafe07ff',
    ' #fffffff5 #fe07 ]'), ), number=3)
mdb.models['Model-1'].parts['pipe'].generateMesh()

#刚体面修改
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#3200 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#140 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#4000 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#1800 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#7 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#108 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#4 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#20 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].RemoveFaces(deleteCells=False, faceList=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(mask=(
    '[#10 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].seedPart(deviationFactor=0.1,
    minSizeFactor=0.1, size=560.0)
mdb.models['Model-1'].parts['Part-1'].seedEdgeByNumber(constraint=FINER, edges=
    mdb.models['Model-1'].parts['Part-1'].edges.getSequenceFromMask(('[#511 ]',
    ), ), number=8)
mdb.models['Model-1'].parts['Part-1'].seedEdgeByNumber(constraint=FINER, edges=
    mdb.models['Model-1'].parts['Part-1'].edges.getSequenceFromMask(('[#a44 ]',
    ), ), number=8)
mdb.models['Model-1'].parts['Part-1'].generateMesh()
mdb.models['Model-1'].parts['Part-1'].deleteMesh(regions=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(('[#f ]',
    ), ))
mdb.models['Model-1'].parts['Part-1'].setMeshControls(regions=
    mdb.models['Model-1'].parts['Part-1'].faces.getSequenceFromMask(('[#f ]',
    ), ), technique=STRUCTURED)
mdb.models['Model-1'].parts['Part-1'].generateMesh()
mdb.models['Model-1'].rootAssembly.regenerate()
mdb.models['Model-1'].rootAssembly.ReferencePoint(point=(0.0, 0.0, 0.0))
mdb.models['Model-1'].RigidBody(bodyRegion=Region(
    faces=mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
    mask=('[#f ]', ), )), name='Constraint-1', refPointRegion=Region(
    referencePoints=(mdb.models['Model-1'].rootAssembly.referencePoints[15],
    )))
mdb.models['Model-1'].boundaryConditions['fix-x'].suppress()
mdb.models['Model-1'].boundaryConditions['fix-y'].suppress()
mdb.models['Model-1'].boundaryConditions['fix-z'].suppress()
mdb.models['Model-1'].boundaryConditions['fix-z'].resume()
mdb.models['Model-1'].boundaryConditions['fix-z'].suppress()
mdb.models['Model-1'].EncastreBC(createStepName='Initial', localCsys=None,
    name='BC-6', region=Region(referencePoints=(
    mdb.models['Model-1'].rootAssembly.referencePoints[15], )))
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='30-20-165', nodalOutputPrecision=SINGLE,
    numCpus=1, queue=None, scratch='', type=ANALYSIS, userSubroutine='',
    waitHours=0, waitMinutes=0)


发表于 2012-11-18 22:14:54 | 显示全部楼层 来自 北京
Simdroid开发平台
没有文件,也不知道是哪里出问题啊
回复

使用道具 举报

发表于 2012-11-18 22:22:54 | 显示全部楼层 来自 北京
能不能贴个图出来看看
回复

使用道具 举报

 楼主| 发表于 2012-11-19 10:01:26 | 显示全部楼层 来自 北京
飞翔的fresh 发表于 2012-11-18 22:22
能不能贴个图出来看看

谢谢关注,贴不上图,都上传不了,只能放py了
回复

使用道具 举报

 楼主| 发表于 2012-11-19 21:07:52 | 显示全部楼层 来自 北京
没有大牛关注啊
回复

使用道具 举报

 楼主| 发表于 2012-11-19 21:28:42 | 显示全部楼层 来自 北京
顶上去..
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-29 02:52 , Processed in 0.038661 second(s), 9 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表