Bmodes手册

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

Bmodes的计算基于:

1\有限元方法

2\线性化解析

3\科里奥利和离心效应

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

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

主要输入文件：

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

sec_loc  str_tw  tw_iner  mass_den flp_iner  edge_iner  flp_stff    edge_stff     tor_stff    axial_stff  cg_offst  sc_offst tc_offst
(-)      (deg)    (deg)   (kg/m)    (kg-m)    (kg-m)     (Nm^2)      (Nm^2)        (Nm^2)        (N)         (m)       (m)     (m)
0.00000  11.100   11.100 1447.607  646.044   646.044  7681.46E+06  7681.46E+06  2655.23E+06  17152.7E+06    0.000     0.000    0.0
0.02105  11.100   11.100  180.333   80.480    80.480  1169.87E+06  1169.87E+06   408.80E+06   2640.8E+06    0.000     0.000    0.0
0.05263  11.100   11.100  181.672   68.241    80.113  1020.62E+06  1092.28E+06   343.81E+06   2611.3E+06    0.032    -0.005    0.0
0.10526  11.100   11.100  183.905   47.842    79.502   771.88E+06   962.97E+06   235.50E+06   2562.1E+06    0.086    -0.014    0.0
0.15789  11.100   11.100  186.138   27.444    78.892   523.14E+06   833.66E+06   127.19E+06   2512.9E+06    0.140    -0.023    0.0
0.21053  11.100   11.100  188.370    7.045    78.281   274.40E+06   704.35E+06    18.87E+06   2463.6E+06    0.194    -0.032    0.0
0.26316   9.500    9.500  178.321    5.963    68.302   234.57E+06   614.65E+06    16.80E+06   2332.8E+06    0.188    -0.020    0.0
0.31579   7.900    7.900  168.271    4.881    58.323   194.74E+06   524.96E+06    14.72E+06   2202.0E+06    0.182    -0.007    0.0
0.36842   6.300    6.300  158.222    3.799    48.344   154.90E+06   435.26E+06    12.64E+06   2071.2E+06    0.176     0.005    0.0
0.42105   4.700    4.700  148.172    2.717    38.366   115.07E+06   345.57E+06    10.56E+06   1940.4E+06    0.170     0.018    0.0
0.47368   3.100    3.100  138.123    1.635    28.387    75.23E+06   255.87E+06     8.48E+06   1809.6E+06    0.164     0.030    0.0
0.52632   2.600    2.600  122.896    1.367    24.050    62.49E+06   217.87E+06     7.12E+06   1605.3E+06    0.168     0.038    0.0
0.57895   2.100    2.100  107.669    1.099    19.714    49.75E+06   179.86E+06     5.76E+06   1401.1E+06    0.172     0.047    0.0
0.63158   1.600    1.600   92.442    0.831    15.377    37.01E+06   141.86E+06     4.40E+06   1196.8E+06    0.176     0.055    0.0
0.68421   1.100    1.100   77.215    0.564    11.041    24.27E+06   103.85E+06     3.04E+06    992.6E+06    0.179     0.063    0.0
0.73684   0.600    0.600   61.988    0.296     6.704    11.53E+06    65.85E+06     1.68E+06    788.3E+06    0.183     0.071    0.0
0.78947   0.480    0.480   51.861    0.240     5.513     9.27E+06    54.25E+06     1.38E+06    654.3E+06    0.190     0.077    0.0
0.84211   0.360    0.360   41.734    0.185     4.322     7.01E+06    42.66E+06     1.08E+06    520.4E+06    0.198     0.082    0.0
0.89474   0.240    0.240   31.607    0.130     3.130     4.75E+06    31.06E+06     0.78E+06    386.4E+06    0.205     0.087    0.0
0.94737   0.120    0.120   21.480    0.074     1.939     2.49E+06    19.47E+06     0.48E+06    252.4E+06    0.212     0.092    0.0
1.00000   0.000    0.000   11.353    0.019     0.747     0.23E+06     7.87E+06     0.18E+06    118.5E+06    0.220     0.098    0.0