有限元分析_软件流程

六：Element tapes 不同的计算需要不同的网格类型，例如，四面体网格在
ABAQUS静强度计算中类型为C3D10M。 具体判断方法为：在分析软件CAE中设置计算类型，在
MESH模块查看网格类型。
1 对于ANASYS,要进行二次升级 2
七：Create 在完成以上工作以后，网格工作已经完成。下面创建 用于
S,E,SP,SINV *PRINT, FREQ=99, CONTACT=YES, MODEL CHANGE=YES, PLASTICITY=YES, RESIDUAL=YES, SOLVE=YES *CONTACT PRINT, FREQ=99 *CONTACT FILE, FREQ=99 *END STEP
提示：此处SURFAC 为螺杆中部任意截面，为保证截面平整，划分网格 时应特别处理。
2 压力
*DSLOAD side_load,P,0.625
压力大小，单位：MPa
热载荷
*Initial Conditions, type=TEMPERATURE node_all, 20.
*step
**调入温度场计算结果
文件后缀随意，一般为*.mesh
附：操作技巧
第三步：ABAQUS/ ANSYS
命令流方便、简洁、准确。计算设置主要是对命令流操作。 ABAQUS的命令流放在*.inp文件中 ANSYS的命令流放在*.log文件中
scratch="/cylhead_abq_orig_model/temperature_orig " 文件路径 pre_memory=100000000
注意： 特征线能保证网格与实体的确切对应，切忌随便Toggle； 有时还要创建特征线。
三：2V Mesh
二维网格划分有两种方法： 自动划分：复杂模型整体划分 手动调整：整体网格存在的前提下，局部调整； 在模型局部修改问题上应用较广；
Distribute node to edge Distribute node to edge
方程系数
约束方程
五：过盈
ABAQUS
*CONTACT PAIR, INTERACTION=SHL2-BP, SMALL, ADJUST=0.1
SHL2_BP, BP12 *SURFACE INTERACTION, NAME=SHL2-BP
接触建过盈
*CLEARANCE, MASTER=BP12, SLAVE=SHL2_BP, VALUE=-0.100
*BOUNDARY, OP=MOD 102680, 1,3, 0.00000E+00 SPC_Y_BS000001, 1,3, 0.00000E+00
节点编号
节点 Set 名字
单个节点 节点 Set
*BOUNDARY top,11,,32 firewall,11,,200
温度：32、200
热边界
八：载荷_力
! 进入前处理 ! 定义材料1的弹性模量 ! 定义材料1的泊松比 ! 定义材料1的热膨胀系数 ! 定义材料1的密度 ! 定义材料1的比热 ! 定义材料1的热传导系数
! 设置单元的材料编号 EMODIF,cyhead,MAT,1,
! 气缸盖为材料1
三：接触
ABAQUS
*CONTACT PAIR, INTERACTION=BL-BP_31, SMALL, ADJUST=0.1
1 螺栓预紧力
*NODE, NSET=BOLTPRET 999133, 1000., 1000., 1000. *PRE-TENSION SECTION,SURFACE=SURF_B33,NODE=999133 *CLOAD, OP=MOD 999133,1,182700. 预紧力大小，单位：N
预紧力施加方向
1.1200E-05, 0.0 *CONDUCTIVITY,TYPE=ISO
0.04479 ,0.0 *SPECIFIC HEAT
480.0 ,0.0 *DENSITY
7.3500E-06,0.0
定义材料－ANSYS
/PREP7
MP,EX,1,1.10e5 MP,PRXY,1,0.27 MP,ALPX,1,10.5e-6 MP,DENS,1,7.35e-6 MP,C,1,480 MP,KXX,1, 0.04479
1
2
点击
3
1
说明： 一 1处不选中，特征线保留以前的网格密度，选中则重新划分；
2处控制网格大小； 3处设置网格类型；
二 右图详细调节网格属性 。
Create elements Create elements Rule况。
Spline Ruled
BL31, BP31
*SURFACE INTERACTION, NAME=BL-BP_31
从面，主面
Tied
*EQUATION 2
节点编号
122318, 2, 1.00, 1224067, 2, -1.00
约束方程
两个节点 2方向自由度
ANSYS
节点编号 CE,1,0,15,UX,1,100208,UX,-1
首页
有限元分析
软件流程
第一步：Pro/E几何模型 依据实体几何尺寸，准确真实的划出Pro/E三维模型。
在Pro/E中保存副本, 文件格式为*.stp
第二步：Hyper Mesh网格划分 网格划分有以下几个步骤：
Import
Surface edit
2V mesh 3V mesh Order change Elem tapes Create Export
计算的Sets 、nodes 、surface。
确定Hyper mesh在ABAQUS 等计算软件环境下
Sets 、nodes依据需要创建
Create surface
注
3
意
命
名
规
则
2
4
1
对创建的面修改
注意
八： Export 为便于写*.inp文件，应将网格分块输出； 再分别导入，重新编号；各块之间编号留有余量
999999, 5, , 0.00
999999, 6, , 0.00
*EQUATION
三个节点
3
646966, 1, 1.00, 920795, 1, -1.00, 999999, 1, -1.00
约束方程
六：Step
*STEP, AMPLITUDE=RAMP Bolts Pre-Loads
*STATIC *HEAT TRANSFER, STEADY STATE
文件夹
文件名：网格、约束
ANSYS
*use,../mesh/cyhead.mesh
二：定义材料－ABAQUS
指定
弹性模量、迫松比 热膨胀系数 热传导系数 比热 密度
** -----------------------------------*SOLID SECTION, ELSET=SOLID_CRANKCASE ,
Break up into collectors Clear up Toggle
Distribute node to edge Create elements Edit elements
Sets 、nodes 、surface
一：Import 将Pro/E模型导入Hyper Mesh，方法如下图。
BL31, BP31
*SURFACE INTERACTION, NAME=BL-BP_31
*FRICTION, ELASTIC SLIP=0.002
0.15,
从面，主面：固定面、稀网格面作主面
ANSYS
四：MPC
ABAQUS
*CONTACT PAIR, INTERACTION=BL-BP_31, TIED, ADJUST=0.1
*temperature,file=CYL_ring *END STEP
提示：温度场可由 FIRE 等软件扫出，问题在于扫出的温度场是以节点给 出的，应该加在网格上。
九：输出
*STEP *NODE FILE, FREQ=99
U, *EL PRINT, FREQ=0
S, *EL FILE, FREQ=99
MATERIAL=SOLID_CRANKCASE *MATERIAL,NAME=SOLID_CRANKCASE ** -----------------------------------*ELASTIC,TYPE=ISOTROPIC
1.10000E+05, 2.70000E-01 *EXPANSION, ZERO = 20.0 , TYPE = ISO
Edit elements 优化网格交界处的连贯性，删除小、差网格，重建。
四：3V mesh 二维网格的好坏，决定三维网格质量。
利用此项，不会对二维网格进行调整，有利于保证一一对应
五：Order change 对于四面体网格，为计算精确，利用Order change 将其变成10节点。
中间节点的增加，将大大增加节点个数，对计算机要求更高。 在保证计算精度的前提下，应尽量增大网格密度
standard_memory=400000000
standard_buffer=5000000
ask_delete=OFF
#run_mode=INTERACTIVE
ABAQUS文件系统设置
odb_output_by_default=ON
cpus = 1
一：模型调入
ABAQUS
*HEADING
*INCLUDE, INPUT=../MESH_M/crankcase.mesh
十：提交计算
第四步：查看结果
Hyper View有很强的结果显示功能
将ABAQUS计算结果文件*.fil转化为Hyper View认可的格式
文件夹有此文件
二：Surface edit _Break up into collectors Break up into collectors： ➢将几何模型分块，放在各自的collector里, 注意命名规则。 ➢ collector的建立，有利于选中对象，便于操作。 ➢在操作时注意：每当创建新的物体时，提前修改Global， 确定新的物体被放在了当前collector。
*NODE, NSET=PROTRUSION 999999, 1000., 1000., 1000.
创建节点，坐标值随意
*BOUNDARY, OP=MOD
999999, 1, , 0.10
999999, 2, , 0.00
999999, 3, , 0.00 999999, 4, , 0.00
限定节点自由度偏移量
0.2, 1.
*BOUNDARY ..
*CLOAD ..
*NODE PRINT ..
*END STEP
Step 名称 静强度稳态 计算温度场
计算增量0.2
可以设置多个Step,各个Step设置不同的约束、载荷 和输出，有利于不同工况之间的叠加。
七：边界约束 ABAQUS
约束加在节点上，1,3， 0 表示约束1、2、3方向自由度
2
1
Clear up & Toggle
几何模型包含不规则交界面以及由特征线 组成的尖、小区域。 需要压缩、合并、重建特征线、点。
2
1
2
1
Clear up & Toggle_续
Toggle与否的断定方法：
两条特征线之间区域小于一个网格； 对相邻面划分二维网格，检测网格质量（眼测）， Toggle， 再划网格，确保网格质量；