ABAQUS脚本运行总是出错

WIN11系统

ABAQUS2023版本

运行结果

# -*- coding: utf-8 -*-
from abaqus import *
from abaqusConstants import *
import math

# ========================== 1. 基础参数定义 ==========================
# 工件尺寸 (mm)
WP_LENGTH = 500.0 # 长度
WP_HEIGHT = 200.0 # 高度
WP_THICK = 10.0 # 厚度

# 刀具参数
RAKE_ANGLE = 10.0 # 前角
CLEAR_ANGLE = 5.0 # 后角
TOOL_HEIGHT = 300.0 # 刀具高度
TOOL_LENGTH = 150.0 # 刀具宽度

# 切削参数
vc = 50.0 * 1000 / 60.0 # 50 m/min → 转换为 mm/s (Abaqus标准单位)
f = 0.15 # 进给量 mm/r

# 网格参数
MESH_SIZE = 5.0 # 网格大小

# ========================== 2. 创建模型 ==========================
modelName = "OrthogonalCutting_TC4"
mdb.Model(name=modelName, modelType=STANDARD_EXPLICIT)
mod = mdb.models[modelName]

# ========================== 3. 建立工件几何（500×200mm） ==========================
sketch_wp = mod.ConstrainedSketch(name="WP_Sketch", sheetSize=2000.0)
sketch_wp.rectangle(point1=(0.0, 0.0), point2=(WP_LENGTH, WP_HEIGHT))

part_wp = mod.Part(
name="Workpiece",
dimensionality=TWO_D_PLANAR,
type=DEFORMABLE_BODY
)
part_wp.BaseShell(sketch=sketch_wp)

# ========================== 4. 建立刀具几何（前角10°，后角5°） ==========================
sketch_tool = mod.ConstrainedSketch(name="Tool_Sketch", sheetSize=2000.0)
tan_r = math.tan(math.radians(RAKE_ANGLE))
tan_c = math.tan(math.radians(CLEAR_ANGLE))

# 刀具四点坐标
p1 = (0.0, 0.0)
p2 = (TOOL_LENGTH, 0.0)
p3 = (TOOL_LENGTH - tan_r * TOOL_HEIGHT, TOOL_HEIGHT)
p4 = (-tan_c * TOOL_HEIGHT, TOOL_HEIGHT)

sketch_tool.Line(point1=p1, point2=p2)
sketch_tool.Line(point1=p2, point2=p3)
sketch_tool.Line(point1=p3, point2=p4)
sketch_tool.Line(point1=p4, point2=p1)

# 修复：使用正确的常量 ANALYTIC_RIGID_SURFACE
part_tool = mod.Part(
name="Tool",
dimensionality=TWO_D_PLANAR,
type=ANALYTIC_RIGID_SURFACE
)
part_tool.AnalyticRigidSurf2DPlanar(sketch=sketch_tool)

# ========================== 5. 定义材料 Ti6Al4V (TC4) ==========================
mat = mod.Material(name="Ti6Al4V")
mat.Density(table=((4.43e-9,),)) # 密度 (t/mm³)
mat.Elastic(table=((110E3, 0.34),)) # 弹性参数 (MPa)

# Johnson-Cook 塑性本构（TC4标准参数）
mat.Plastic(
table=((860.0, 650.0, 450.0, 0.001, 1200.0, 0.31, 0.48),),
hardening=JOHNSON_COOK
)

# 热参数（温度场）
mat.Conductivity(table=((6.8,),))
mat.SpecificHeat(table=((560.0,),))
mat.Expansion(table=((9.0e-6,),))

# ========================== 6. 截面与网格（CPE4R） ==========================
# 截面属性
mod.HomogeneousSolidSection(
name="WP_Section",
material="Ti6Al4V",
thickness=WP_THICK
)
part_wp.Set(name="WP_Set", faces=part_wp.faces)
part_wp.SectionAssignment(region=part_wp.sets["WP_Set"], sectionName="WP_Section")

# 网格划分 CPE4R
part_wp.seedPart(size=MESH_SIZE)
part_wp.setMeshControls(regions=part_wp.faces, elemShape=QUAD, technique=STRUCTURED)
elemType = mesh.ElemType(elemCode=CPE4R, elemLibrary=EXPLICIT)
part_wp.setElementType(regions=part_wp.faces, elemTypes=(elemType,))
part_wp.generateMesh()

# ========================== 7. 装配模型 ==========================
assem = mod.rootAssembly
assem.DatumCsysByDefault(CARTESIAN)
instWP = assem.Instance(name="WP_Inst", part=part_wp, dependent=ON)
instTool = assem.Instance(name="Tool_Inst", part=part_tool, dependent=ON)

# 定位刀具到切削位置
assem.translate(instanceList=("Tool_Inst",), vector=(10.0, WP_HEIGHT, 0.0))

# ========================== 8. 显式分析步 ==========================
mod.ExplicitDynamicsStep(
name="Cutting_Step",
previous="Initial",
timePeriod=WP_LENGTH / vc + 0.01,
nlgeom=ON
)

# ========================== 9. 接触与摩擦（库仑摩擦） ==========================
# 定义接触属性
mod.ContactProperty("Friction")
mod.interactions["Friction"].Table(type=FRICTION, table=((0.15,),))

# 创建接触相互作用
mod.ContactStd(name="Tool_WP_Contact", createStepName="Initial")

# ========================== 10. 边界条件 ==========================
# 工件底部完全固定
bottom_edges = instWP.edges.findAt(((WP_LENGTH/2, 0.0, 0.0),))
assem.Set(edges=bottom_edges, name="Bottom_Edge")
mod.EncastreBC(name="Fix_Workpiece", createStepName="Initial", region=assem.sets["Bottom_Edge"])

# 刀具水平匀速切削
# 创建参考点并施加速度
tool_ref_point = assem.ReferencePoint(point=(10.0, WP_HEIGHT, 0.0))
tool_ref_id = tool_ref_point.id
assem.Set(referencePoints=(assem.referencePoints[tool_ref_id],), name="Tool_Ref_Point")
mod.VelocityBC(
name="Tool_Move",
createStepName="Cutting_Step",
region=assem.sets["Tool_Ref_Point"],
v1=-vc, v2=0.0
)

# ========================== 11. 输出结果 ==========================
mod.FieldOutputRequest(
name="F_Output",
createStepName="Cutting_Step",
variables=("S", "PE", "PEEQ", "TEMP", "HFL", "STATUS")
)

# ========================== 12. 提交作业 ==========================
job = mdb.Job(
name="Job_TC4_Cutting",
model=modelName,
description="2D Orthogonal Cutting Ti6Al4V"
)
job.submit()
job.waitForCompletion()

运行结果