```python # -*- coding: mbcs -*- # # Abaqus/CAE Release 2020 replay file # Internal Version: 2019_09_14-01.49.31 163176 # Run by Administrator on Thu May 19 11:05:34 2022 # # from driverUtils import executeOnCaeGraphicsStartup # executeOnCaeGraphicsStartup() #: Executing "onCaeGraphicsStartup()" in the site directory ... from abaqus import * from abaqusConstants import * session.Viewport(name='Viewport: 1', origin=(0.0, 0.0), width=250.050506591797, height=145.686111450195) session.viewports['Viewport: 1'].makeCurrent() session.viewports['Viewport: 1'].maximize() from caeModules import * from driverUtils import executeOnCaeStartup executeOnCaeStartup() session.viewports['Viewport: 1'].partDisplay.geometryOptions.setValues( referenceRepresentation=ON) # RVE建立 - 模型建立 (基体) s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=200.0) g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=STANDALONE) s.rectangle(point1=(0.0, 0.0), point2=(10.0, 10.0)) p = mdb.models['Model-1'].Part(name='MATRIX', dimensionality=THREE_D, type=DEFORMABLE_BODY) p = mdb.models['Model-1'].parts['MATRIX'] p.BaseSolidExtrude(sketch=s, depth=10.0) s.unsetPrimaryObject() p = mdb.models['Model-1'].parts['MATRIX'] session.viewports['Viewport: 1'].setValues(displayedObject=p) del mdb.models['Model-1'].sketches['__profile__'] # RVE建立 - 模型建立 (纤维) p = mdb.models['Model-1'].parts['MATRIX'] f, e = p.faces, p.edges t = p.MakeSketchTransform(sketchPlane=f[4], sketchUpEdge=e[7], sketchPlaneSide=SIDE1, sketchOrientation=RIGHT, origin=(5.0, 5.0, 10.0)) s1 = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=34.64, gridSpacing=0.86, transform=t) g, v, d, c = s1.geometry, s1.vertices, s1.dimensions, s1.constraints s1.setPrimaryObject(option=SUPERIMPOSE) p = mdb.models['Model-1'].parts['MATRIX'] p.projectReferencesOntoSketch(sketch=s1, filter=COPLANAR_EDGES) s1.CircleByCenterPerimeter(center=(0.0, 0.0), point1=(1.72, 0.0)) p = mdb.models['Model-1'].parts['MATRIX'] f1, e1 = p.faces, p.edges p.SolidExtrude(sketchPlane=f1[4], sketchUpEdge=e1[7], sketchPlaneSide=SIDE1, sketchOrientation=RIGHT, sketch=s1, depth=10.0, flipExtrudeDirection=ON, keepInternalBoundaries=ON) s1.unsetPrimaryObject() del mdb.models['Model-1'].sketches['__profile__'] # 材料属性 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='MATRIX') mdb.models['Model-1'].materials['MATRIX'].Elastic(table=((7000.0, 0.38), )) mdb.models['Model-1'].Material(name='FIBER') mdb.models['Model-1'].materials['FIBER'].Elastic(table=((220000.0, 0.2), )) # 创建截面并赋予 mdb.models['Model-1'].HomogeneousSolidSection(name='FIBER-SEC', material='FIBER', thickness=None) mdb.models['Model-1'].HomogeneousSolidSection(name='MATRIX-SET', material='MATRIX', thickness=None) p = mdb.models['Model-1'].parts['MATRIX'] c = p.cells cells = c.getSequenceFromMask(mask=('[#1 ]', ), ) region = p.Set(cells=cells, name='Set-1') p = mdb.models['Model-1'].parts['MATRIX'] p.SectionAssignment(region=region, sectionName='FIBER-SEC', offset=0.0, offsetType=MIDDLE_SURFACE, offsetField='', thicknessAssignment=FROM_SECTION) p = mdb.models['Model-1'].parts['MATRIX'] c = p.cells cells = c.getSequenceFromMask(mask=('[#2 ]', ), ) region = p.Set(cells=cells, name='Set-2') p = mdb.models['Model-1'].parts['MATRIX'] p.SectionAssignment(region=region, sectionName='MATRIX-SET', offset=0.0, offsetType=MIDDLE_SURFACE, offsetField='', thicknessAssignment=FROM_SECTION) # 创建装配体 a = mdb.models['Model-1'].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=a) session.viewports['Viewport: 1'].assemblyDisplay.setValues( optimizationTasks=OFF, geometricRestrictions=OFF, stopConditions=OFF) a = mdb.models['Model-1'].rootAssembly a.DatumCsysByDefault(CARTESIAN) p = mdb.models['Model-1'].parts['MATRIX'] a.Instance(name='MATRIX-1', part=p, dependent=OFF) # Mesh session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=ON) session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues( meshTechnique=ON) a = mdb.models['Model-1'].rootAssembly c1 = a.instances['MATRIX-1'].cells pickedRegions = c1.getSequenceFromMask(mask=('[#3 ]', ), ) a.setMeshControls(regions=pickedRegions, elemShape=WEDGE) a = mdb.models['Model-1'].rootAssembly partInstances =(a.instances['MATRIX-1'], ) a.seedPartInstance(regions=partInstances, size=0.5, deviationFactor=0.1, minSizeFactor=0.1) a = mdb.models['Model-1'].rootAssembly partInstances =(a.instances['MATRIX-1'], ) a.generateMesh(regions=partInstances) session.viewports['Viewport: 1'].enableMultipleColors() session.viewports['Viewport: 1'].setColor(initialColor='#BDBDBD') cmap=session.viewports['Viewport: 1'].colorMappings['Material'] session.viewports['Viewport: 1'].setColor(colorMapping=cmap) session.viewports['Viewport: 1'].disableMultipleColors() # EasyPBC import sys sys.path.insert(8, r'e:/CAE/abaqus_plugins/EasyPBC V.1.4') import easypbc easypbc.feasypbc(part='Model-1', inst='MATRIX-1', meshsens=1E-07, CPU=1, E11=True, E22=True, E33=False, G12=True, G13=False, G23=False, onlyPBC=False, CTE=False, intemp=0, fntemp=100) ``` 以上代码是根据博客内容按操作逻辑顺序拼接而成的完整脚本。它遵循了Abaqus二次开发的标准流程：从初始化环境、创建几何部件（基体与纤维）、定义材料与截面属性、装配部件、进行网格划分，到最后调用`EasyPBC`插件进行周期性边界条件计算与均匀化弹性属性预测 [ref_1]。该脚本理论上可在Abaqus/CAE的命令行接口（CLI）或通过`abaqus cae noGUI=`方式运行，前提是用户已正确安装并配置了`EasyPBC`插件路径。需要注意的是，脚本中部分几何索引（如`f[4]`, `e[7]`）依赖于特定建模顺序下的对象ID，在完全空白的模型中首次运行时可能需要进行调整。