频率响应分析¶
本分析使用 tutorial/17_freq_beam 中的数据。
分析对象¶
分析对象为悬臂梁，几何形状如图4.17.1所示，网格数据如图4.17.2所示。
| 项目 | 内容 | 备注 | 参考 |
|---|---|---|---|
| 分析类型 | 频率响应分析 | !SOLUTION,TYPE=EIGEN !SOLUTION,TYPE=DYNAMIC | |
| 节点数 | 55 | ||
| 单元数 | 126 | ||
| 单元类型 | 4节点一次四面体单元 | !ELEMENT,TYPE=341 | |
| 材料属性名 | Material-1 | !MATERIAL,NAME=Material-1 | |
| 边界条件 | 约束、集中载荷、特征值 | !EIGENREAD | |
| 矩阵求解方法 | CG/SSOR | !SOLVER,METHOD=CG,PRECOND=1 |
图4.17.1 悬臂梁的几何形状
图4.17.2 悬臂梁的网格数据
分析内容¶
完全约束分析对象悬臂梁的一端，并在另一端的两个节点上施加集中载荷，进行频率响应分析。
步骤是先在相同边界条件下进行到第10阶的模态分析，然后使用到第5阶的特征值和特征向量进行分析。频率响应分析用的分析控制数据如下所示。
分析控制数据 beam_eigen.cnt¶
# Control File for FISTR
!VERSION
3
!WRITE,RESULT
!WRITE,VISUAL
!SOLUTION, TYPE=EIGEN
!EIGEN
10, 1.0E-8, 60
!BOUNDARY
_PickedSet4, 1, 3, 0.0
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=NO,TIMELOG=YES
10000, 1
1.0e-8, 1.0, 0.0
!VISUAL,metod=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
分析控制数据 beam_freq.cnt¶
# Control File for FISTR
!VERSION
3
!WRITE,RESULT
!WRITE,VISUAL
!SOLUTION, TYPE=DYNAMIC
!DYNAMIC
11 , 2
14000, 16000, 20, 15000.0
0.0, 6.6e-5
1, 1, 0.0, 7.2E-7
10, 2, 1
1, 1, 1, 1, 1, 1
!EIGENREAD
eigen_0.log
1, 5
!BOUNDARY
_PickedSet4, 1, 3, 0.0
!FLOAD, LOAD CASE=2
_PickedSet5, 2, 1.
!FLOAD, LOAD CASE=2
_PickedSet6, 2, 1.
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=NO,TIMELOG=YES
10000, 1
1.0e-8, 1.0, 0.0
!VISUAL,metod=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
分析步骤¶
首先，将模态分析用的全局控制数据 hecmw_ctrl_eigen.dat 改名为 hecmw_ctrl.dat，并进行模态分析。
接着，将频率响应分析用的全局控制数据 hecmw_ctrl_freq.dat 改名为 hecmw_ctrl.dat，并将模态分析的分析结果日志文件 0.log 改名为 eigen_0.log（在频率响应分析用的分析控制数据中指定），然后进行频率响应分析。
cp hecmw_ctrl_eigen.dat hecmw_ctrl.dat
fistr1 -t 4
mv 0.log eigen_0.log
cp hecmw_ctrl_freq.dat hecmw_ctrl.dat
fistr1 -t 4
分析结果¶
使用 Microsoft Excel 绘制分析控制数据中指定的监测节点（节点编号1）的频率与位移幅值关系，并如图4.17.3所示。 此外，作为分析结果的数值数据，下面给出分析结果日志文件的一部分。
图4.17.3 监测节点的频率与位移幅值关系
分析结果日志 0.log（频率响应分析）¶
fstr_setup: OK
Rayleigh alpha: 0.0000000000000000
Rayleigh beta: 7.1999999999999999E-007
read from=eigen_0.log
start mode= 1
end mode= 5
start frequency: 14000.000000000000
end frequency: 16000.000000000000
number of the sampling points 20
monitor nodeid= 1
14100.000000000000 [Hz] : 8.3935554529723885E-002
14100.000000000000 [Hz] : 1 .res
14200.000000000000 [Hz] : 9.1211083509607632E-002
14200.000000000000 [Hz] : 2 .res
14300.000000000000 [Hz] : 9.9579777896922961E-002
14300.000000000000 [Hz] : 3 .res
14400.000000000000 [Hz] : 0.10914967594967491
14400.000000000000 [Hz] : 4 .res
14500.000000000000 [Hz] : 0.11992223203326918
14500.000000000000 [Hz] : 5 .res
14600.000000000000 [Hz] : 0.13164981801806747
14600.000000000000 [Hz] : 6 .res
14700.000000000000 [Hz] : 0.14360931008440975
14700.000000000000 [Hz] : 7 .res
14800.000000000000 [Hz] : 0.15436500205940235
14800.000000000000 [Hz] : 8 .res
14900.000000000000 [Hz] : 0.16180768408076251
14900.000000000000 [Hz] : 9 .res
15000.000000000000 [Hz] : 0.16388019610373711
15000.000000000000 [Hz] : 10 .res
15100.000000000000 [Hz] : 0.15982110598747551
15100.000000000000 [Hz] : 11 .res
15200.000000000000 [Hz] : 0.15074650286398145
15200.000000000000 [Hz] : 12 .res
15300.000000000000 [Hz] : 0.13885370598993371
15300.000000000000 [Hz] : 13 .res
15400.000000000000 [Hz] : 0.12618976409021948
15400.000000000000 [Hz] : 14 .res
15500.000000000000 [Hz] : 0.11405716994112736
15500.000000000000 [Hz] : 15 .res
15600.000000000000 [Hz] : 0.10306231010139058
15600.000000000000 [Hz] : 16 .res
15700.000000000000 [Hz] : 9.3374567545990342E-002
15700.000000000000 [Hz] : 17 .res
15800.000000000000 [Hz] : 8.4945897112663621E-002
15800.000000000000 [Hz] : 18 .res
15900.000000000000 [Hz] : 7.7641947016103510E-002
15900.000000000000 [Hz] : 19 .res
16000.000000000000 [Hz] : 7.1307642422355627E-002
16000.000000000000 [Hz] : 20 .res
start time: 0.0000000000000000
end time: 6.6000000000000005E-005
frequency: 15000.000000000000
node id: 1
num disp: 10
time= 0.0000000000000000 : 1 .res
time= 0.0000000000000000 : 1 .vis
time= 6.6000000000000003E-006 : 2 .res
time= 6.6000000000000003E-006 : 2 .vis
time= 1.3200000000000001E-005 : 3 .res
time= 1.3200000000000001E-005 : 3 .vis
time= 1.9800000000000000E-005 : 4 .res
time= 1.9800000000000000E-005 : 4 .vis
time= 2.6400000000000001E-005 : 5 .res
time= 2.6400000000000001E-005 : 5 .vis
time= 3.3000000000000003E-005 : 6 .res
time= 3.3000000000000003E-005 : 6 .vis
time= 3.9600000000000000E-005 : 7 .res
time= 3.9600000000000000E-005 : 7 .vis
time= 4.6200000000000005E-005 : 8 .res
time= 4.6200000000000005E-005 : 8 .vis
time= 5.2800000000000003E-005 : 9 .res
time= 5.2800000000000003E-005 : 9 .vis
time= 5.9400000000000000E-005 : 10 .res
time= 5.9400000000000000E-005 : 10 .vis