请教ATV的计算
本算例利用ATV计算板的声辐射。总是到最后的 计算场点的时候,其他都正常。
:'(
提示的错误为:
............................
ANALYSIS NUMBER 32 / 37 FREQUENCY = 175.000 Hz
------------------------------------------------------------------------
CPU Stats --- Delta:0:00:00 ( 0.00)Total:0:06:39 ( 399.12)
LOADING POTENTIALS FROM DATABASE at Frequency 175.000 Hz
ERROR - CANNOT FIND RECORD OF TYPE POT1
ERROR - COULD NOT FIND POTENTIALS AT 175.00 HZ
...............................
附件为模型文件和命令流。
请大侠指教。:handshake
ATV.cmd
New Model 1 File ATV_structure.sdb Return
OptionFEM Frequency Structure Return
Import Mesh Format Free File plate.fre Return
Geometry Thickness 0.001
Elements All
Return
Material Plate
Name 'steel'
Young 2.1e11Poisson 0.3Rho 7800
Elements all
Return
Save Return
Set Name"Envelope" Envelope
Elements all
Return
Boundary
UX Real 0 Imag 0 UY Real 0 Imag 0 UZ Real 0 Imag 0
Nodes Set 1
Return
Save Return
Modes
Vector 50 Tolerance 1.0000e-6 Iteration 100
BLanczos
Return
Project Modes Structure Return
Save Return
Save Return
New Model 2 File ATV_fluid.sdb Return
OptionBEM Indirect Variational Uncoupled Unbaffled FrequencyReturn
Import Mesh Format Free File plate.fre Return
Set 11 Name"Envelope" Envelope
Elements All
Return
Boundary Jump PressureReal 0 Imag 0
Nodes Set 11
Return
Save Return
Activate Model 1 Return
Boundary
FZ Real 300 Imag 0
Nodes 50,60,73
Return
Save Return
Activate Model 2 Return
Point Plane0 0.5 1 To 1 0.5 1 Divide 5 To 0 0.5 0 Divide 5 Return
Compute ATV
Frequency 20 To 300 LinStep 10
Near 2
Far 5
Quadrature 2 2 1
Return
Save Return
Link
Model 1
Elements All
To
Model 2
Elements All
Behavior FLUID-STRUCTURE-WEAK
Return
Project ATV Return
Save Return
Parameter Model 2
MATV
FInterpolation Linear
Save Potentials none
Save Results Step 1
Store Results none
Return
Near 2
Far 5
Quadrature 2 2 1
Positive
Return
Parameter Model 1
Vector 50
Save Displacements Step 1
Store Displacements none
Return
Return
Solve
Frequency 20 To 200 LinStep 5
Return
Save Return
PostProcess
Points all
Frequency 20 To 200 LinStep 5
Near 2
Far 5
Quadrature 2 2 1
Save Results Step 1
Return
plate.fre
SYSNOISE MESH FILE
121 100 4
NODES
1 1 .00000000E+00 .00000000E+00 .00000000E+00
2 2 .10000000E+00 .00000000E+00 .00000000E+00
3 3 .20000000E+00 .00000000E+00 .00000000E+00
4 4 .30000000E+00 .00000000E+00 .00000000E+00
5 5 .40000000E+00 .00000000E+00 .00000000E+00
6 6 .50000000E+00 .00000000E+00 .00000000E+00
7 7 .60000000E+00 .00000000E+00 .00000000E+00
8 8 .70000000E+00 .00000000E+00 .00000000E+00
9 9 .80000000E+00 .00000000E+00 .00000000E+00
10 10 .90000000E+00 .00000000E+00 .00000000E+00
11 11 .10000000E+01 .00000000E+00 .00000000E+00
12 12 .00000000E+00 .10000000E+00 .00000000E+00
13 13 .10000000E+00 .10000000E+00 .00000000E+00
14 14 .20000000E+00 .10000000E+00 .00000000E+00
15 15 .30000000E+00 .10000000E+00 .00000000E+00
16 16 .40000000E+00 .10000000E+00 .00000000E+00
17 17 .50000000E+00 .10000000E+00 .00000000E+00
18 18 .60000000E+00 .10000000E+00 .00000000E+00
19 19 .70000000E+00 .10000000E+00 .00000000E+00
20 20 .80000000E+00 .10000000E+00 .00000000E+00
21 21 .90000000E+00 .10000000E+00 .00000000E+00
22 22 .10000000E+01 .10000000E+00 .00000000E+00
23 23 .00000000E+00 .20000000E+00 .00000000E+00
24 24 .10000000E+00 .20000000E+00 .00000000E+00
25 25 .20000000E+00 .20000000E+00 .00000000E+00
26 26 .30000000E+00 .20000000E+00 .00000000E+00
27 27 .40000000E+00 .20000000E+00 .00000000E+00
28 28 .50000000E+00 .20000000E+00 .00000000E+00
29 29 .60000000E+00 .20000000E+00 .00000000E+00
30 30 .70000000E+00 .20000000E+00 .00000000E+00
31 31 .80000000E+00 .20000000E+00 .00000000E+00
32 32 .90000000E+00 .20000000E+00 .00000000E+00
33 33 .10000000E+01 .20000000E+00 .00000000E+00
34 34 .00000000E+00 .30000000E+00 .00000000E+00
35 35 .10000000E+00 .30000000E+00 .00000000E+00
36 36 .20000000E+00 .30000000E+00 .00000000E+00
37 37 .30000000E+00 .30000000E+00 .00000000E+00
38 38 .40000000E+00 .30000000E+00 .00000000E+00
39 39 .50000000E+00 .30000000E+00 .00000000E+00
40 40 .60000000E+00 .30000000E+00 .00000000E+00
41 41 .70000000E+00 .30000000E+00 .00000000E+00
42 42 .80000000E+00 .30000000E+00 .00000000E+00
43 43 .90000000E+00 .30000000E+00 .00000000E+00
44 44 .10000000E+01 .30000000E+00 .00000000E+00
45 45 .00000000E+00 .40000000E+00 .00000000E+00
46 46 .10000000E+00 .40000000E+00 .00000000E+00
47 47 .20000000E+00 .40000000E+00 .00000000E+00
48 48 .30000000E+00 .40000000E+00 .00000000E+00
49 49 .40000000E+00 .40000000E+00 .00000000E+00
50 50 .50000000E+00 .40000000E+00 .00000000E+00
51 51 .60000000E+00 .40000000E+00 .00000000E+00
52 52 .70000000E+00 .40000000E+00 .00000000E+00
53 53 .80000000E+00 .40000000E+00 .00000000E+00
54 54 .90000000E+00 .40000000E+00 .00000000E+00
55 55 .10000000E+01 .40000000E+00 .00000000E+00
56 56 .00000000E+00 .50000000E+00 .00000000E+00
57 57 .10000000E+00 .50000000E+00 .00000000E+00
58 58 .20000000E+00 .50000000E+00 .00000000E+00
59 59 .30000000E+00 .50000000E+00 .00000000E+00
60 60 .40000000E+00 .50000000E+00 .00000000E+00
61 61 .50000000E+00 .50000000E+00 .00000000E+00
62 62 .60000000E+00 .50000000E+00 .00000000E+00
63 63 .70000000E+00 .50000000E+00 .00000000E+00
64 64 .80000000E+00 .50000000E+00 .00000000E+00
65 65 .90000000E+00 .50000000E+00 .00000000E+00
66 66 .10000000E+01 .50000000E+00 .00000000E+00
67 67 .00000000E+00 .60000000E+00 .00000000E+00
68 68 .10000000E+00 .60000000E+00 .00000000E+00
69 69 .20000000E+00 .60000000E+00 .00000000E+00
70 70 .30000000E+00 .60000000E+00 .00000000E+00
71 71 .40000000E+00 .60000000E+00 .00000000E+00
72 72 .50000000E+00 .60000000E+00 .00000000E+00
73 73 .60000000E+00 .60000000E+00 .00000000E+00
74 74 .70000000E+00 .60000000E+00 .00000000E+00
75 75 .80000000E+00 .60000000E+00 .00000000E+00
76 76 .90000000E+00 .60000000E+00 .00000000E+00
77 77 .10000000E+01 .60000000E+00 .00000000E+00
78 78 .00000000E+00 .70000000E+00 .00000000E+00
79 79 .10000000E+00 .70000000E+00 .00000000E+00
80 80 .20000000E+00 .70000000E+00 .00000000E+00
81 81 .30000000E+00 .70000000E+00 .00000000E+00
82 82 .40000000E+00 .70000000E+00 .00000000E+00
83 83 .50000000E+00 .70000000E+00 .00000000E+00
84 84 .60000000E+00 .70000000E+00 .00000000E+00
85 85 .70000000E+00 .70000000E+00 .00000000E+00
86 86 .80000000E+00 .70000000E+00 .00000000E+00
87 87 .90000000E+00 .70000000E+00 .00000000E+00
88 88 .10000000E+01 .70000000E+00 .00000000E+00
89 89 .00000000E+00 .80000000E+00 .00000000E+00
90 90 .10000000E+00 .80000000E+00 .00000000E+00
91 91 .20000000E+00 .80000000E+00 .00000000E+00
92 92 .30000000E+00 .80000000E+00 .00000000E+00
93 93 .40000000E+00 .80000000E+00 .00000000E+00
94 94 .50000000E+00 .80000000E+00 .00000000E+00
95 95 .60000000E+00 .80000000E+00 .00000000E+00
96 96 .70000000E+00 .80000000E+00 .00000000E+00
97 97 .80000000E+00 .80000000E+00 .00000000E+00
98 98 .90000000E+00 .80000000E+00 .00000000E+00
99 99 .10000000E+01 .80000000E+00 .00000000E+00
100 100 .00000000E+00 .90000000E+00 .00000000E+00
101 101 .10000000E+00 .90000000E+00 .00000000E+00
102 102 .20000000E+00 .90000000E+00 .00000000E+00
103 103 .30000000E+00 .90000000E+00 .00000000E+00
104 104 .40000000E+00 .90000000E+00 .00000000E+00
105 105 .50000000E+00 .90000000E+00 .00000000E+00
106 106 .60000000E+00 .90000000E+00 .00000000E+00
107 107 .70000000E+00 .90000000E+00 .00000000E+00
108 108 .80000000E+00 .90000000E+00 .00000000E+00
109 109 .90000000E+00 .90000000E+00 .00000000E+00
110 110 .10000000E+01 .90000000E+00 .00000000E+00
111 111 .00000000E+00 .10000000E+01 .00000000E+00
112 112 .10000000E+00 .10000000E+01 .00000000E+00
113 113 .20000000E+00 .10000000E+01 .00000000E+00
114 114 .30000000E+00 .10000000E+01 .00000000E+00
115 115 .40000000E+00 .10000000E+01 .00000000E+00
116 116 .50000000E+00 .10000000E+01 .00000000E+00
117 117 .60000000E+00 .10000000E+01 .00000000E+00
118 118 .70000000E+00 .10000000E+01 .00000000E+00
119 119 .80000000E+00 .10000000E+01 .00000000E+00
120 120 .90000000E+00 .10000000E+01 .00000000E+00
121 121 .10000000E+01 .10000000E+01 .00000000E+00
ELEMENTS
1 1 6 4 1 2 13 12
2 2 6 4 2 3 14 13
3 3 6 4 3 4 15 14
4 4 6 4 4 5 16 15
5 5 6 4 5 6 17 16
6 6 6 4 6 7 18 17
7 7 6 4 7 8 19 18
8 8 6 4 8 9 20 19
9 9 6 4 9 10 21 20
10 10 6 4 10 11 22 21
11 11 6 4 12 13 24 23
12 12 6 4 13 14 25 24
13 13 6 4 14 15 26 25
14 14 6 4 15 16 27 26
15 15 6 4 16 17 28 27
16 16 6 4 17 18 29 28
17 17 6 4 18 19 30 29
18 18 6 4 19 20 31 30
19 19 6 4 20 21 32 31
20 20 6 4 21 22 33 32
21 21 6 4 23 24 35 34
22 22 6 4 24 25 36 35
23 23 6 4 25 26 37 36
24 24 6 4 26 27 38 37
25 25 6 4 27 28 39 38
26 26 6 4 28 29 40 39
27 27 6 4 29 30 41 40
28 28 6 4 30 31 42 41
29 29 6 4 31 32 43 42
30 30 6 4 32 33 44 43
31 31 6 4 34 35 46 45
32 32 6 4 35 36 47 46
33 33 6 4 36 37 48 47
34 34 6 4 37 38 49 48
35 35 6 4 38 39 50 49
36 36 6 4 39 40 51 50
37 37 6 4 40 41 52 51
38 38 6 4 41 42 53 52
39 39 6 4 42 43 54 53
40 40 6 4 43 44 55 54
41 41 6 4 45 46 57 56
42 42 6 4 46 47 58 57
43 43 6 4 47 48 59 58
44 44 6 4 48 49 60 59
45 45 6 4 49 50 61 60
46 46 6 4 50 51 62 61
47 47 6 4 51 52 63 62
48 48 6 4 52 53 64 63
49 49 6 4 53 54 65 64
50 50 6 4 54 55 66 65
51 51 6 4 56 57 68 67
52 52 6 4 57 58 69 68
53 53 6 4 58 59 70 69
54 54 6 4 59 60 71 70
55 55 6 4 60 61 72 71
56 56 6 4 61 62 73 72
57 57 6 4 62 63 74 73
58 58 6 4 63 64 75 74
59 59 6 4 64 65 76 75
60 60 6 4 65 66 77 76
61 61 6 4 67 68 79 78
62 62 6 4 68 69 80 79
63 63 6 4 69 70 81 80
64 64 6 4 70 71 82 81
65 65 6 4 71 72 83 82
66 66 6 4 72 73 84 83
67 67 6 4 73 74 85 84
68 68 6 4 74 75 86 85
69 69 6 4 75 76 87 86
70 70 6 4 76 77 88 87
71 71 6 4 78 79 90 89
72 72 6 4 79 80 91 90
73 73 6 4 80 81 92 91
74 74 6 4 81 82 93 92
75 75 6 4 82 83 94 93
76 76 6 4 83 84 95 94
77 77 6 4 84 85 96 95
78 78 6 4 85 86 97 96
79 79 6 4 86 87 98 97
80 80 6 4 87 88 99 98
81 81 6 4 89 90 101 100
82 82 6 4 90 91 102 101
83 83 6 4 91 92 103 102
84 84 6 4 92 93 104 103
85 85 6 4 93 94 105 104
86 86 6 4 94 95 106 105
87 87 6 4 95 96 107 106
88 88 6 4 96 97 108 107
89 89 6 4 97 98 109 108
90 90 6 4 98 99 110 109
91 91 6 4 100 101 112 111
92 92 6 4 101 102 113 112
93 93 6 4 102 103 114 113
94 94 6 4 103 104 115 114
95 95 6 4 104 105 116 115
96 96 6 4 105 106 117 116
97 97 6 4 106 107 118 117
98 98 6 4 107 108 119 118
99 99 6 4 108 109 120 119
100 100 6 4 109 110 121 120 经过 confer指点已经解决了。这里特别感谢!
主要问题:
1.在定义耦合面时要将结构和流体分别定义成不同的SET。否则在耦合时分不清楚。
2.受到5.6详解这本书的误导(可能不适合5.5版本),其实不需要求解场点,因为在计算ATV时已经计算了。
只需要提取就可以了,具体用store命令。
附件是修改后的命令流。
New Model 1 File ATV_structure.sdb Return
OptionFEM Frequency Structure Return
Import Mesh Format Free File plate.fre Return
Geometry Thickness 0.001
Elements All
Return
Material Plate
Name 'steel'
Young 2.1e11Poisson 0.3Rho 7800
Elements all
Return
Save Return
Set Name"Envelope" Envelope
Elements all
Return
Set 20 Name"plate"
Elements All
Return
Boundary
UX Real 0 Imag 0 UY Real 0 Imag 0 UZ Real 0 Imag 0
Nodes Set 1
Return
Save Return
Modes
Vector 50 Tolerance 1.0000e-6 Iteration 100
BLanczos
Return
Project Modes Structure Return
Save Return
Save Return
New Model 2 File ATV_fluid.sdb Return
OptionBEM Indirect Variational Uncoupled Unbaffled FrequencyReturn
Import Mesh Format Free File plate.fre Return
Set 11 Name"Envelope" Envelope
Elements All
Return
Set 30 Name"fluid"
Elements All
Return
Boundary Jump PressureReal 0 Imag 0
Nodes Set 11
Return
Save Return
Activate Model 1 Return
Boundary
FZ Real 300 Imag 0
Nodes 50,60,73
Return
Save Return
Activate Model 2 Return
Point Plane0 0.5 1 To 1 0.5 1 Divide 5 To 0 0.5 0 Divide 5 Return
Compute ATV
Frequency 20 To 300 LinStep 10
Near 2
Far 5
Quadrature 2 2 1
Return
Save Return
Link
Model 1
Elements Set 20
To
Model 2
Elements Set 30
Behavior FLUID-STRUCTURE-WEAK
Return
Project ATV Return
Save Return
Parameter Model 2
MATV
FInterpolation Linear
Save Potentials none
Save Results Step 1
Store Results none
Return
Near 2
Far 5
Quadrature 2 2 1
Positive
Return
Parameter Model 1
Vector 50
Save Displacements Step 1
Store Displacements none
Return
Return
Solve
Frequency 20 To 200 LinStep 5
Return
Save Return
Store
Case Active
Results
Frequency 20 To 300 LinStep 10
Points all
Return
Return 我随后进行了 贡献量的分析 貌似有结果,但不知道是否正确。
而且在计算过程中同样出现了第一个帖子的提示,但有结果。
不知道有谁遇到过没有? 能否获得单元贡献量结果,即某个单元而不是节点对声场的贡献量? 我今天又算了一次,发现2个很奇怪的问题:
1.不管我选那个场点,最后还是全部都算出来了。
2.不管我选那个单元,最后算的只有那几个施加了速度边界条件的几个单元有。这点我不是很明白。
希望有人能指点一下。 会不会是以前的结果保留的? 我也觉得奇怪,我重新算过2次,不会是软件的问题吧! 我不太清楚你的问题,不过我用ATV分析后,再用单个节点施加法向振动(Parameters: Physical)可以验证加权矩阵中的对应元素。
回复 楼主 pengweicai 的帖子
真复杂。。。请问楼主,是不是要得到声压值,就得计算ATV啊?
我对sysnoise的步骤还是不怎么清楚。 不要,ATV只是其中的一种算法。。 ATV是计算某个/集合单元或某个/集合节点的声贡献量, 回复 5 # pengweicai 的帖子
是声学传递向量这种算法,只计算有振动速度的板之贡献,不振动,无贡献。所以只有加速度边界条件的单元有贡献结果
页:
[1]