Bmodes使用手册

Bmodes可以计算叶片或塔架的耦合模态.对于基于模态线形叠加的程序如Fast和Bladed程序,耦合模态的计算对气弹模拟就至关重要. 其主要是为fast而开发的。

 

Bmodes的计算基于:

1\有限元方法

2\线性化解析

3\科里奥利和离心效应

每个单元有两个边界点3个内部节点,共15个自由度.

需要定义结构属性分布、转速、风轮锥角、桨距角、叶片几何。fast当前只用到了非耦合的模态，而忽略了扭转自由度和剪切中心、拉伸中心、质心的偏置，而这将导致严重的动态耦合。fast利用bmodes计算的挥舞-摆振-扭转耦合模态可以通过隐含地考虑扭转自由度来克服限制。

 

主要输入文件：

--------- General parameters ---------------------------------------------------------------------
true     Echo        Echo input file contents to *.echo file if true.
1         beam_type   1: blade, 2: tower (-)
60.0      rot_rpm:    rotor speed, automatically set to zero for tower modal analysis (rpm)
1.0       rpm_mult:   rotor speed muliplicative factor (-)
35.0      radius:     rotor tip radius measured along coned blade axis OR tower height (m)
1.75      hub_rad:    hub radius measured along coned blade axis OR tower rigid-base height (m)
0.        precone:    built-in precone angle, automatically set to zero for a tower (deg)
0.        bl_thp:     blade pitch setting, automatically set to zero for a tower (deg)
1         hub_conn:   hub-to-blade or tower-base boundary condition [1: cantilevered; 2: free-free; 3: only axial and torsion constraints] (-)
20        modepr:     number of modes to be printed (-)
t         TabDelim    (true: tab-delimited output tables; false: space-delimited tables)
f         mid_node_tw  (true: output twist at mid-node of elements; false: no mid-node outputs)

--------- Blade-tip or tower-top mass properties --------------------------------------------
0.        tip_mass    blade-tip or tower-top mass (kg)
0.        cm_loc      tip-mass c.m. offset from the blade axis measured along the tip section y reference axis (m)
0.        cm_axial    tip-mass c.m. offset tower tip measures axially along the z axis (m)
0.        ixx_tip     blade lag mass moment of inertia about the tip-section x reference axis (kg-m^2)
0.        iyy_tip     blade flap mass moment of inertia about the tip-section y reference axis (kg-m^2)
0.        izz_tip     torsion mass moment of inertia about the tip-section z reference axis (kg-m^2)
0.        ixy_tip     cross product of inertia about x and y reference axes(kg-m^2)
0.        izx_tip     cross product of inertia about z and x reference axes(kg-m^2)
0.        iyz_tip     cross product of inertia about y and z reference axes(kg-m^2)

--------- Distributed-property identifiers --------------------------------------------------------
1         id_mat:     material_type [1: isotropic; non-isotropic composites option not yet available]
'blade_sec_props.dat' sec_props_file   name of beam section properties file (-)

Property scaling factors..............................
1.0       sec_mass_mult:   mass density multiplier (-)
1.0       flp_iner_mult:   blade flap or tower f-a inertia multiplier (-)
1.0       lag_iner_mult:   blade lag or tower s-s inertia multiplier (-)
1.0       flp_stff_mult:   blade flap or tower f-a bending stiffness multiplier (-)
1.0       edge_stff_mult:  blade lag or tower s-s bending stiffness multiplier (-)
1.0       tor_stff_mult:   torsion stiffness multiplier (-)
1.0       axial_stff_mult: axial stiffness multiplier (-)
1.0       cg_offst_mult:   cg offset multiplier (-)
1.0       sc_offst_mult:   shear center multiplier (-)
1.0       tc_offst_mult:   tension center multiplier (-)

--------- Finite element discretization --------------------------------------------------
12        nselt:     no of blade or tower elements (-)
Distance of element boundary nodes from blade or flexible-tower root (normalized wrt blade or tower length), el_loc()
0.  0.08  0.16  0.24  0.32  0.40  0.48  0.56  0.64  0.72  0.80  0.90   1.0

截面属性输入文件