Zebulon 软件 相关命令

Zebulon 安装：

用户环境变量：

Path=%PATH%;%Z7PATH%\bin；  //定义Zebulon安装位置，将来inp文件要通过位置查找材料属性信息。

Z7PATH=C:\Zebulon_nt

Z7_MAX_NB_DOF=999999  //计算最大自由度数

win10 用户修改：

将Zebulon文件夹下renum.bat  more 修改为：type %Z7PATH%\bin\renumtd | Zrenum %1

 

Zebulon 命令：

1 ZBatch_depou  //%Z7PATH%\bbA\Batch_depou_NT.exe %1 生成网格图片

2 ZExport_def_mesh NissanD90_Dynamic 0 1 CC_CONE  NissanD90_Dynamic_def   //重新生成geof文件, 更改一些格式 0 1 前一个数字是 结果的步数

需要文件夹内有 .node .geof .ut文件  前一个数字是分析的increeasement数，后一个是squence数

3 copy /A NissanD90_Dynamic_def.geof  “NissanD90_Dynamic_kz-.geof” 复制一份用于Kz-计算

4 renum NissanD90_Dynamic_kz- //单元重新编号

5 Zrun NissanD90_Dynamic_kz-

6 Zideas  //unv文件转化为geof文件

7 zrun_ns  // zrun 模型较大用zrun_ns

8 Zebu_to_abaqus // 参数 geof文件名，将geof网格文件转化为abaqus的inp文件（geom文件）

 

关键字：

9 swaging

**radial

Nset_name origin direction value[table]

mobio 0. 0. 0. d2 -1.0 tabU0

//origin 轴向量经过的点 swaging方向必须和global一致

//d1 d2 d3指向轴心方向为Z，则为d3

 

10 ***eigen lanczos

20 2000  //20 Number of Eigen frequencies  , 2000 Max Frequency

**frequency_scanning

0 2000. 20. //0~2000Hz 间隔20

11 rotation

2d

**rotation [node]

Nset_name ((origin)|num_node)

Angle [table]

3d

**rotation [node]  //[node] 旋转中心，abaqus中参考点

Nset_name  [num_node](axe)  node节点编号

[(origin)]  angle [table]  //origin 旋转方向

Example：

**rotation

MOBI  (0.0 0.0 1.0) (0.0 0.0 0.0)  15. tabU3

Nset  旋转轴方向向量  旋转轴通过点 角度

 

12 contact

// Impactor (nset slave)  target(faset master)

***contact soft

**zone

*impactor contact

*target contact

*friction 1.

*warning_distance 0.8

*gap 1.e-1

**conv 1.e-5 30

**soft_param 1.e3 1.

13 关键字

****calcul none/ eigen / harmonic_response_by_modal_superposition/ harmonic_direct_response

普通分析/模态提取/模态叠加法扫频，直接法扫频

动刚度扫频计算要设置：**elset CC small_deformation

橡胶计算，橡胶模态提取： **elset CC_CONE total_lagrangian_mixte_u_ps

14

***eigen lanczos
100 2000.             求模态

***harmonic_response    //谐响应
**eigens_file  ******.rst
**frequency_scanning     //扫频
0. 2000. 20.

***harmonic_response   //直接法
**frequency_scanning
0. 2000. 20.
**damping 0.2

***resolution newton

15

***material
**elset CC_CONE
*file rubber.txt

材料文件格式：

普通材料： 密度模量泊松比

***coefficient
masvol 1.1e-9
***behavior linear_elastic
**elasticity
young 1.84
poisson 0.49
***return

橡胶材料（静态分析）：

***coefficient
masvol 1.1e-9
***behavior hyper_elastic_mixte_mooney
**mooney
moon1 0.926
moon2 0.000
**thermal_strain
alpha 1.4e-4
****return

橡胶材料（动态分析）：模量已乘系数1.4

%=============================================
%Fichier materiau genere pour un c/c de fleche 2.00
% Avec une rigidification de 1.4
%module de cisaillement G = 2.02
%=============================================
***coefficient
masvol 1.1e-9
***behavior hyper_elastic_mixte_mooney
**mooney
moon1 1.009
moon2 0.000
**thermal_strain
alpha 1.4e-4
****return

 

16

***output

**component //定义输出区域

component1 2

**fequency //定义输出结果的频率

*dtime

*increment

****return

默认输出 和文件结尾

****post  //后处理文件设置

***curve 生成曲线

**factor 2   //曲线系数，半模型 2

**x_component
*name frequency
*scale 0 2000 20
**y_component
*name mod_char3
*position MOBI     //曲线设置

**save_image_curve D:\kevin.sun\ANT2019\NissanWZ2A\CAE\Z\nissan_wz2a_dx.png
**save_curve_value D:\kevin.sun\ANT2019\NissanWZ2A\CAE\Z\nissan_wz2a_dx.pnt

//保存曲线数据和曲线图片

****return

17 定义曲线

**sequence 2  //分析步数 sequences or blocks of time steps

*time 1. 2.    //时间 用于定义time 和 value

*ratio 1.      //精度

*increment 4 4   //每block 分多少步计算，默认为1

*iteration 20   //  最大迭代步数

*algorithm p1p2p3   //Newton-Raphson algorithm where the matrix is always updated
during iterations 橡胶一般用牛顿拉夫森法

*limit_dof  // for stability with bulging

U1 1. 1.

U2 1. 1.

U3 1. 1.

**name tabU1

*time    0. 1. 2.

*value   0.   0.0000e+000   0.0000e+000

18 释放MOBI U2自由度 ，计算降温位移

***equation
**mpc1 MOBI U2

***equation
**mpc4_wd IMPA2 fa_CIBLE2 0.5 //将两面tie到一起，搜索范围0.5   nset  linset(faset)

***equation    //从1.001s开始 FULIDE面填充物质，计算k1
**nul_div_u FLUIDE     #[correction dir]
*time_start 1.001

***curve
**factor 2
**x_component   曲线X轴
*name U3           x轴数据
*position MOBI
*scale -9 0 9       坐标范围
**y_component
*name char3
*position MOBI
**compute_rate      计算斜率
*x_value_sup -8
**save_image_curve D:\kevin.sun\ANT2018\Mazda_D-typeA\MazdaTypeA_v14-2_k1.png
**save_curve_value D:\kevin.sun\ANT2018\Mazda_D-typeA\MazdaTypeA_v14-2_k1.pnt