ABAQUS计算随机振动设置及输出

ABAQUS计算随机振动设置及输出

1.分析步设置

随机振动主要包括两个分析步:频率和随机振动
在这里插入图片描述

1.1 频率设置

频率这里需要注意的是最高频率最好是扫频范围的2-2.5倍
比如随机频率区间是[0-2000hz],最高频率应该大于4000Hz,才能保证精度
在这里插入图片描述

1.2 随机响应设置

这里需要设置的主要为基本信息和阻尼
基本信息:低频率为扫频区间下限,高频率为上限,缩放选择对数,Bias默认3即可,点数这里涉及到精度,默认是20,但是20的话计算过慢,减少这个值可以加快计算速度。
在这里插入图片描述
点数等于5的效果

Bias的影响

阻尼信息:可以选择根据模态阶数或者直接对频率设置,阻尼系数按照要求填写即可,一般为0.02。
在这里插入图片描述

2.分析步输出

由于随机振动关注的应力是Rmises,因此需要进行一些设置才能保证输出
其他的关于加速度之类的不再叙述

2.1 频率输出

频率分析这里场输出选择Mises,S(这里不选择的话,后续随即设置输出会报错)在这里插入图片描述

2.1 随机输出

历程输出和场输出勾选RMISES-等效均方根应力
在这里插入图片描述

3.载荷边界

主要包括边界:位移边界和加速度基运动
在这里插入图片描述

3.1 位移边界

不多叙述,固定6个自由度即可
在这里插入图片描述

3.2 加速度基运动

基本信息:自由度也就是施加激励方向,U3也就是整体坐标系Z向
在这里插入图片描述
Correlation:对于单方向激励,Approach选择不想关即可;
PSD谱是关键输入,后续的实部虚部默认1,0即可
在这里插入图片描述
PSD是在幅值中定义的,本模型采用的是M制,因此输出选择中重力,并且X9.8,mm制度X9800
关于PSD计算,可以参考另一篇博客:https://blog.csdn.net/weixin_44446998/article/details/108697419?spm=1001.2014.3001.5501

在这里插入图片描述

4.结论

一般而言,材料在随机振动中的许用应力以3σ为原则,其值为屈服应力的1/3(若以屈服作为失效判定),但由于分析无法模拟接触和冲击,得到的RMISES值偏低,故须对许用应力值进行修改,该修改值应以实际情况和以往试验结果为标准进行重新定义。
ABAQUS做随机分析速度慢,模型大点就计算不动,而用nastran,workbench计算特别快,速度差10倍以上,希望有专业人士看看能不能解决

### Abaqus 中进行随机振动分析并输出加速度的设置方法 #### 功率谱密度(PSD)作为输入 在Abaqus中执行随机振动分析时,输入通常采用功率谱密度(PSD)。PSD可以表示为位移谱、速度谱或加速度谱。为了获得加速度响应,在定义载荷时应选择合适的PSD形式[^1]。 ```python # 定义功率谱密度(PSD)作为加载条件的一部分 psd_load = model.PsdDefinition( name='Acceleration_PSD', psdData=((frequency_values, acceleration_psd), ), type=ACCELERATION, timeSpan=TOTAL ) ``` #### 创建模态分析步 随机振动分析依赖于先前完成的模态分析结果。因此,需先创建一个用于提取结构固有特性的模态分析步: ```python modal_step = model.StaticStep(name='Modal_Analysis', previous='Initial') model.setValuesInStep(stepName='Modal_Analysis', nlgeom=ON) model.FieldOutputRequest(name='F-Output-1', createStepName='Modal_Analysis', variables=('U','UR')) ``` #### 添加随机振动分析步 接着增加一个新的分析步专门处理随机振动问题。此步骤会利用之前求解得到的模式形状来进行后续计算: ```python random_vibration_step = model.RandomResponseStep( name="Random_Vibration", prev="Modal_Analysis" ) # 设定阻尼参数 model.DampingRatioByFrequencyRange( stepName="Random_Vibration", frequencyRanges=(min_freq,max_freq), dampingRatios=damping_ratio_value ) ``` #### 请求特定输出项 为了让程序能够输出所需的加速度数据,必须显式指定希望保存哪些变量。对于加速度而言,可以通过`*Node Output`命令请求节点上的加速度信息;而对于单元内部,则可通过`Element Output`来实现这一点。 ```python # 节点输出设定 node_output_request = model.HistoryOutputRequest( name='Nodal_Acceleration_Output', createStepName='Random_Vibration', variables=('A',) ) element_output_request = model.FieldOutputRequest( name='Elemental_Acceleration_Output', createStepName='Random_Vibration', variables=('EACC',) ) ``` 以上配置完成后提交作业运行,最终可以从ODB文件读取到所关心位置处随时间变化的加速度记录。
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值