abaqus经验总结.doc

例如要取缸套的主推力侧上一列节点的径向位移（缸套轴线平行于 Z 轴），然后按照 Z 坐标画出径向位移曲线总的方法是先选上要分析的节点，然后按空间位置顺序记录节点编号，根据编号生成一个 Path，再以这个 path 为横坐标，以要分析的量（如径向位移）为纵坐标画 X-Y 图1. 选出要分析区域节点：这列节点共有几十个，一个一个选太慢，打开要分析的 odb 文件，选择 Display Group 中的【Replace Selected】按钮，在目标类型中选择 【Nodes】，然后结合视图方向、框选方法（矩形、圆形、多边形框）和选择方式（individually 或者 by angle）选出需要的节点集合，例子中的一列节点只有 Z 坐标不同，因此把视图方向调整到 Z 轴垂直于屏幕，这时一列点在屏幕上变成了一个点，在这一个点的位置上框选就可以得到一列点了在individually 选择方式下，【shift+框选】是将选到的加入到已有选择中，【Ctrl+框选】是将选到的从已有选择中去除例如要在一个曲面上选一列节点，可以先用【by angle】方式将整个曲面选上，然后调整到合适的视图方向上切换到【individually】方式，用【Ctrl+ 框选】去掉多余的节点。

2. 按顺序记录节点编号：Path 对顺序很敏感，节点号顺序的调整会改变最后 plot 的数据点顺序要按照 Z 坐标大小顺序画出径向位移曲线，就必须按顺序记录节点编号上一步选好节点按中键确定后，屏幕上是空白，因为 ABAQUS不会显示单独的节点，打开显示节点编号开关（在【Common Plot Options->Labels】中），这样屏幕上会显示出要选的那些节点的编号，按顺序记录下来号码之间用逗号分隔，冒号表示连续和间隔，例如：1.2 （1 号和 2 号节点），1:10（表示 1,2,…,10），1:9:2（表示 1,3,5,…,9）3. 生成目标表格和曲线，两种方法：a）用上面排好的节点序列生成一个 Path，然后基于这个 path 生成一个 X-Y 图（Create X-Y data -> Path），在【X-Y Data Manager】中双击生成的图线可以得到相应的二维表格；b）菜单【Report->Field Output】，这种方法将二维表格数据写入一个文件，本身不能出图，但是可以对计算结果求和，因此可以用来求接触力、压力的合力（但要注意这里的求和是数值相加，不是矢量求和，因此只适用于各个力方向基本一致的情况下的合力估计）。

有时候历史变量输出的太多会出现“The number of history output requests (22466) in this analysis step has exceeded the maximum value of 10000 specified by the Abaqus environment variable 'max_history_requests.'”错误，解决办法是：在 abaqus_v6.env 中加一句“max_history_requests=0”即可abaqus_v6.env 文件在 C:\ABAQUS\6.7-1\site\下比如两个接触面之间有 0.1mm 的缝隙，在压力作用下发生接触，只要将Contact Controls 下的稳定因子设为 1 即可正常计算（Interaction 模块，主菜单 Interaction->Contact Controls->Create->Stabilization->Automatic Stabilization->factor 设为 1）；对于缝隙不均匀的情况，比如一端宽 0.1mm，一端宽 0.2mm，减小网格尺寸和稳定因子对计算结果几乎没有影响；[04/08] 关于 ABAQUS 历史输出变量 在帮助文档 Analysis User’s Manual 中 4.2.2 节 Nodal variables 部分定义了节点输出（用*Node Output 关键字表示）可以包括的变量，从中可以看到应力（S）是不能作为节点计算结果输出的，只能作为单元积分点计算结果输出。

原因是节点应力是根据单元积分点应力计算结果外插得到的，因此节点应力只能在场变量（而且是节点所在单元有输出）的情况下才能输出如果在 abaqus/cae 中强行定义节点应力输出，INP 文件中不会有相应的关键字，如果在 INP 文件中强行定义，求解时处理 INP 文件阶段会出错退出由于应力历史变量只能在单元积分点上输出，而在 abaqus/cae 中不能生成基于 element 型的 set，除非直接使用网格模型做 Part，但这样的模型往往不好操作，变通的方法是写由abaqus/cae 生成 INP 文件，再手工修改这个 INP 文件，在其中添加一个 element 型的set，然后让其输出应力历史变量具体步骤：1. 记住要输出的单元编号，例如为 27，在 INP 文件的*Assembly 部分定义一个element 型的 set：*Elset, elset=Set-Name, instance=Part-Name27,2. 在*End Step 前加上：*Output, history, frequency=10*Element Output, elset=Set-Elem-forMeshPartMISES,[03/29] ABAQUS 子程序 VDLOAD 学习笔记[折叠] 使用的 INP 文件（beam.inp）：*Headingload was added in cae** Job name: Job-CAE Model name: Explicit-NoMove*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=beam*Node1, 0., 0.2, 50., 0.3, 100., 0.4, 150., 0. 5, 200., 0.6, 250., 0.7, 300., 0.8, 350., 0.9, 400., 0.10, 450., 0.11, 500., 0.*Element, type=B211, 1, 22, 2, 33, 3, 44, 4, 55, 5, 66, 6, 77, 7, 88, 8, 99, 9, 1010, 10, 11*Elset, elset=_PickedSet2, internal, generate1, 10, 1*Elset, elset=_PickedSet3, internal, generate1, 10, 1** Section: Section-1 Profile: Profile-1*Beam Section, elset=_PickedSet2, material=Material-1, poisson = 0.3, temperature=GRADIENTS, section=RECT10., 5.0.,0.,-1.*End Part** **** ASSEMBLY***Assembly, name=Assembly** *Instance, name=beam-1, part=beam*End Instance** *Elset, elset=forLoad, instance=beam-16, 7*Nset, nset=_PickedSet12, internal, instance=beam-11,*Nset, nset=_PickedSet13, internal, instance=beam-111,*Surface, type=ELEMENT, name=forLoadSPos, internalforLoad, SPOS*End Assembly** ** MATERIALS** *Material, name=Material-1*Damping, alpha=5.*Density7.8e-09,*Elastic210000., 0.3** ** BOUNDARY CONDITIONS** ** Name: BC-1 Type: Displacement/Rotation*Boundary_PickedSet12, 1, 1_PickedSet12, 2, 2** Name: BC-2 Type: Displacement/Rotation*Boundary_PickedSet13, 2, 2** ----------------------------------------------------------------** ** STEP: Step-1** *Step, name=Step-1*Dynamic, Explicit, 0.5*Bulk Viscosity0.06, 1.2** ** LOADS** ** Name: Load-1 Type: Concentrated force*DsloadforLoadSPos, PNU, 1.** ** OUTPUT REQUESTS** *Restart, write, number interval=1, time marks=NO** ** FIELD OUTPUT: F-Output-1** *Output, field, time interval=0.005*Node OutputU, *Element Output, directions=YESS, ** ** HISTORY OUTPUT: H-Output-1** *Output, history, variable=PRESELECT*End Step使用的 for 文件（Gas.for）：subroutine vdload (1 nblock, ndim, stepTime, totalTime,1 amplitude, curCoords, velocity, dirCos, jltyp, sname,value )Cinclude 'vaba_param.inc'Cdimension curCoords(nblock,ndim), velocity(nblock,ndim),1 dirCos(nblock,ndim,ndim), value(nblock)character*80 snameCopen(unit=17,file='e:\aa.txt',status='old',position='append')write(17,*)'stepTime= ',stepTimewrite(17,*)'ndim= ',ndimwrite(17,*)'curCoords= ',curCoordswrite(17,*)'jltyp= ',jltypvalue = 1write(17,*)'sname= ',snamewrite(17,*)'value= ',valuewrite(17,*)' 'close(17)returnend命令行：abaqus job=beam user=gas interactive分析：ABAQUS/Explicit 每个增量步（Increment。