• 洲周州粥
    了解作者
  • Python
    开发工具
  • 2KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 0
    下载次数
  • 2020-03-15 19:02
    上传日期
可以建模一个立方体并提交作业展示作业结果
cube.rar
  • cube.py
    5.2KB
内容介绍
# -*- coding: mbcs -*- # Do not delete the following import lines # 单元体提交作业 #session.viewports视角更改 from abaqus import * from abaqusConstants import * import __main__ import random import section import regionToolset import displayGroupMdbToolset as dgm import part import material import assembly import step import interaction import load import mesh import optimization import job import sketch import visualization import xyPlot import displayGroupOdbToolset as dgo import connectorBehavior #建模 s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=30.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=(20.0, 20.0)) p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=THREE_D, type=DEFORMABLE_BODY) p = mdb.models['Model-1'].parts['Part-1'] p.BaseSolidExtrude(sketch=s, depth=20.0) s.unsetPrimaryObject() p = mdb.models['Model-1'].parts['Part-1'] session.viewports['Viewport: 1'].setValues(displayedObject=p) 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='steel')#材料命名 mdb.models['Model-1'].materials['steel'].Elastic(table=((random.randint(200000,210000), 0.3), ))#输入弹性数值 #创建截面 mdb.models['Model-1'].HomogeneousSolidSection(name='Section-1', material='steel', thickness=None) #指派截面 p = mdb.models['Model-1'].parts['Part-1'] c = p.cells cells = c.getSequenceFromMask(mask=('[#1 ]', ), ) region = p.Set(cells=cells, name='Set-1') p = mdb.models['Model-1'].parts['Part-1'] p.SectionAssignment(region=region, sectionName='Section-1', 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['Part-1'] a.Instance(name='Part-1-1', part=p, dependent=OFF) #创建分析步 session.viewports['Viewport: 1'].assemblyDisplay.setValues( adaptiveMeshConstraints=ON) mdb.models['Model-1'].StaticStep(name='Step-1', previous='Initial', maxNumInc=10000, initialInc=0.01) session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Step-1') session.viewports['Viewport: 1'].assemblyDisplay.setValues(loads=ON, bcs=ON, predefinedFields=ON, connectors=ON, adaptiveMeshConstraints=OFF) #创建载荷 a = mdb.models['Model-1'].rootAssembly s1 = a.instances['Part-1-1'].faces side1Faces1 = s1.getSequenceFromMask(mask=('[#2 ]', ), ) region = a.Surface(side1Faces=side1Faces1, name='Surf-1') mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-1', region=region, distributionType=UNIFORM, field='', magnitude=-100.0, amplitude=UNSET) session.viewports['Viewport: 1'].view.setValues(nearPlane=48.2148, farPlane=98.1267, width=50.305, height=22.6254, cameraPosition=( 52.6647, 51.406, 52.6647), cameraUpVector=(-0.584997, 0.577299, -0.569653)) session.viewports['Viewport: 1'].view.setValues(nearPlane=50.0063, farPlane=98.34, width=52.1742, height=23.4661, cameraPosition=(47.7133, -39.8411, 49.9411), cameraUpVector=(-0.263031, 0.768255, 0.583609), cameraTarget=(10.4196, 9.16085, 10.4196)) #边界条件 a = mdb.models['Model-1'].rootAssembly f1 = a.instances['Part-1-1'].faces faces1 = f1.getSequenceFromMask(mask=('[#8 ]', ), ) region = a.Set(faces=faces1, name='Set-1') mdb.models['Model-1'].PinnedBC(name='BC-1', createStepName='Step-1', region=region, localCsys=None) #划分网格 session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=ON, loads=OFF, bcs=OFF, predefinedFields=OFF, connectors=OFF) session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues( meshTechnique=ON) a = mdb.models['Model-1'].rootAssembly partInstances =(a.instances['Part-1-1'], ) a.seedPartInstance(regions=partInstances, size=1.0, deviationFactor=0.1, minSizeFactor=0.1) a = mdb.models['Model-1'].rootAssembly partInstances =(a.instances['Part-1-1'], ) a.generateMesh(regions=partInstances) session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=OFF) session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues( meshTechnique=OFF) #提交作业 mdb.Job(name='Job-1', model='Model-1', description='', type=ANALYSIS, atTime=None, waitMinutes=0, waitHours=0, queue=None, memory=90, memoryUnits=PERCENTAGE, getMemoryFromAnalysis=True, explicitPrecision=SINGLE, nodalOutputPrecision=SINGLE, echoPrint=OFF, modelPrint=OFF, contactPrint=OFF, historyPrint=OFF, userSubroutine='', scratch='', resultsFormat=ODB, multiprocessingMode=DEFAULT, numCpus=1, numGPUs=0) mdb.jobs['Job-1'].submit(consistencyChecking=OFF) session.mdbData.summary() o3 = session.openOdb(name='d:/Temp/Job-1.odb') session.viewports['Viewport: 1'].setValues(displayedObject=o3) session.viewports['Viewport: 1'].odbDisplay.display.setValues(plotState=( CONTOURS_ON_DEF, ))
评论
    相关推荐
    • 建模计算脚本.rar
      用于立方胞元的6面有限元计算脚本,计算了点阵结构的抗拉强度
    • GaussDB_100_1.0.1-DATABASE-REDHAT-64bit.tar.gz
      guassdb100在redhat上安装包,单机部署的包,安装步骤请看我的文中介绍,经过大量实验搭建总结出来的文档
    • SIM800C_MQTT.rar
      使用SIM800C模块,使用MQTT协议,连接中国移动onenet平台,能实现数据的订阅、发布、存储等
    • 卷积神经网络
      这是卷积神经网络的一个实际用例,已经调试好了,能够在matlab上成功运行,适合从事卷积神经网络(CNN)研究的人员学习使用。
    • android从bootloader到launcher启动流程整理
      讲述android 开机流程 从boot rom---bootloader---init--zygote---systemserver---ams 并附上自己整理的每个流程的流程图 ,清晰熟悉android 启动流程
    • 分数阶混沌系统:分数阶混沌系统的数值解。-matlab开发
      该工具箱包含可用于模拟一些著名的分数阶混沌系统的函数,例如: - 陈的系统, - Arneodo的系统, - Genesio-Tesi 的系统, - 洛伦兹系统, - 牛顿-莱普尼克系统, - 罗斯勒的系统, - Lotka-Volterra系统, - 达芬的系统, - 范德波尔的振荡器, -伏打的系统- 陆氏系统, - 刘的系统, - Chua的系统, - 金融系统, - 3 细胞 CNN。 这些函数以数值方式计算描述混沌系统的分数阶非线性微分方程的解。 每个函数返回总模拟时间的状态轨迹(吸引器)。 更多详情请看书: Ivo Petras,分数阶非线性系统:建模、分析和仿真,Springer,系列:非线性物理科学,2011,ISBN 978-3-642-18100-9。 http://www.springer.com/engineering/control/book/978-3-
    • matlabpam代码-dsp-library:用于光通信的MATLABDSP函数库
      matlab pam代码 DSP库 该项目包含几种不同的功能,可将DSP算法应用于光通信。 这些功能适用于相干和非相干(PAM,DMT)光通信。 大多数功能是相互独立的。 因此,该代码中的功能可以轻松使用,并与其他DSP功能结合使用。 用法 功能列表以及简短说明在文件中。 输入和输出参数的描述在每个函数的标题中。 讯号 通常,输入(和输出)信号在第一维度上具有时间(例如,列向量),而第二维度用于一次管理多个信号(例如,不同的极化,不同的参数等)。 之所以选择这种约定,是因为MATLAB通过将列保留在内存的连续部分中来存储矩阵,因此,这种约定比其他方法(行向量)要快。 参数 大多数功能使用参数结构作为输入参数。 函数中使用的参数的描述通常在标头中,而coherent-dsp函数的默认参数在文件中。 参考 执照 此代码在下发布。
    • 有关多目标跟踪的PHD滤波的一些资料
      一些关于多目标跟踪的新的文献,主要是有关概率假设密度(PHD)的
    • 基于python开发的全国新工商采集工具 v1.2版本
      这个软件是通过scrapy爬虫框架结合代理IP池再加上request模拟请求技术以及验证码识别技术,可以做到日更新采集全国新工商信息。采集的数据自动存储在mysql数据库表里,可下载全量1.8亿多企业工商基本信息和36维度的详细信息.支持sql和excel导出数据包格式。
    • Aerosim Blockset
      The AeroSim aeronautical simulation blockset provides a complete set of tools for the rapid development of nonlinear 6-degree of freedom aircraft dynamic models. In addition to the basic aircraft dynamics blocks, the library also includes complete aircraft models which can be customized through parameter files.