一个螺旋管道表面布满粒子的ABAQUS模型脚本

# -*- coding: utf-8 -*-
"""
Created on 10 mat 22 10:52:10 2022

@author: guozihao
"""

from abaqus import *
from abaqusConstants import *
import random
import math

myModel = mdb.models['Model-1']

def basic():
    s1=myModel.ConstrainedSketch(name='__profile__', sheetSize=20.0)
    s1.CircleByCenterPerimeter(center=(10, 0.0), point1=(20, 0.0))
    s1.Spot(point=(0.0, 0.0))
    p=myModel.Part(name='basic', dimensionality=THREE_D, type=DEFORMABLE_BODY)
    p.BaseSolidExtrude(sketch=s1, depth=100, pitch=200.0)
    del myModel.sketches['__profile__']

def sk(n):
    s1=myModel.ConstrainedSketch(name='__profile__', sheetSize=20.0)
    s1.ConstructionLine(point1=(0.0, -100.0), point2=(0.0, 100.0))
    s1.ArcByCenterEnds(cen