Fortran Coder

标题: 冲击响应の合成时域波形 [打印本页]
作者: hanshikai    时间: 2014-10-12 19:52
标题: 冲击响应の合成时域波形
程序的功能:
根据冲击响应谱试验的国家规范,来修正并合成时域波形(此结果不唯一)

也就是说在这个程序中,冲击响应是事先已知的(由规范给出,并且给出了误差限),利用小波合成的方法来得到时域波形(在这个冲击波的作用下,系统响应符合冲击响应规范)

程序运行过程:
1. 根据相关理论由规范确定小波合成[一定频率的基波序列]时候用到的各个参数(已校验多次,无错,A除外,A是程序中需要修正的变量)
2. 根据相关公式(已校验多次,公式无错)计算冲击响应谱
3. 修正时域波形直至满足误差(因为计算量不大,在误差判断之前先进行一次强行修正)
4. 输出时域频域数据供绘图
5. 完毕

问题:
1. 计算的冲击响应谱结果出奇的好,经过修正之后几乎与规范中的标准值重合;
2. 对于未修正的响应谱,只有低频段落在误差限外,其余均非常贴近规范标准.

对于这样的结果,我是持怀疑态度的这个"出奇的好"让我难以置信.

程序稍微长了点,270行,分了三个模块,保留了必要的注释.

附图为冲击响应谱的图,包括3条规范曲线(因为带有误差),以及修正前后的计算响应谱



下图从文献截取



WAVSYN.f90

8.05 KB, 下载次数: 12
作者: aliouying    时间: 2014-10-13 10:20
我对你的理论不是很清楚，是否需要添加误差，然后再进行修正？

看了下您的程序，貌似只是对理论的数据进行修正。
作者: hanshikai    时间: 2014-10-13 12:54
aliouying 发表于 2014-10-13 10:20
我对你的理论不是很清楚，是否需要添加误差，然后再进行修正？

看了下您的程序，貌似只是对理论的数据进行 ...

感谢您的关注!
1. 误差已经在module NORM_中给出即实常量Tol=3;
2. 本程序所要做的就是对幅值A进行修正,A(:)这个数组中各个元素是不发生联系的;
3. 时域波形是通过一系列特定频率的基波合成的,或者说把时域波形离散成一些列基波的和,在每个时间点上都要做这样的合成.程序中WavMat这个数组存储了所有的时域信息,他的行代表时间,列代表基波序列,按照列来代入有关计算公式(即subroutine Recs_Filt,一次仅代入一列)计算,得到是时间上的结果,挑选最大值就是该基频对应的响应,依次代入WavMat各列,得到各频率对应的响应.附图即是频率和响应关系图(横纵坐标均取常用对数);
4. 上面得到的响应不可避免的会带有误差,SelectCase(nUpDat)中的Case(2)就是修正数据的过程,在检查误差之前有一次强行修正.

再次感谢,不知道还有什么地方没讲清楚.
作者: hanshikai    时间: 2014-10-13 18:56
解决一些问题,并继续提出疑问

1. 关于未修正的响应只有少部分落在误差外
[如果将时间分的粗糙一些(iTimeSegs减小),例如100份,将会得到不太光滑的连线,并且高频部分也将有超出误差的,对此的解释是:采样粗糙丢失某些峰值,作用以后响应减小]

2. 修正后得到的响应太贴近标准值
[如果不进行强行修正,可以看到只有超出误差的基波的A被改善,这时在很大范围内最后的响应与标准值有较大的差别]

3.关于修正后的数据仍有些疑问
[凡是经过修正的均与标准非常贴近,而且调试时发现,各处修正仅需一次就使得误差非常小]

附图为iTimeSegs=100,没有强行修正的结果.



作者: pasuka    时间: 2014-10-14 08:44
1、冲击谱相关标准本来就是先有实验数据再有拟合经验公式，若是船舶的冲击谱，国外不知道炸多少条船做实验，汽车行业就更不用提了；
2、lz的问题和编程看不出有任何联系，理论问题还是请教原文作者或者翻阅标准和相关文献
作者: hanshikai    时间: 2014-10-14 08:55
pasuka 发表于 2014-10-14 08:44
1、冲击谱相关标准本来就是先有实验数据再有拟合经验公式，若是船舶的冲击谱，国外不知道炸多少条船做实验 ...

看来是内行啊!
我的问题和编程真的没有半毛钱关系?请指示
作者: pasuka    时间: 2014-10-14 09:14
hanshikai 发表于 2014-10-14 08:55
看来是内行啊!
我的问题和编程真的没有半毛钱关系?请指示

国防项目的话，通过单位或组织去联系中国工程物理研究院、海军装备研究院、中国船舶科技中心等部门
汽车、建筑的不了解
作者: hanshikai    时间: 2014-10-14 09:18
pasuka 发表于 2014-10-14 09:14
国防项目的话，通过单位或组织去联系中国工程物理研究院、海军装备研究院、中国船舶科技中心等部门
汽车 ...

非常感谢.




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2