有关谐响应谱随机振动的总结文档资料下载.pdf

1、F矩阵和 u矩阵是简谐的，频率为:谐响应分析的运动方程：运动方程 FuKuCuM=+titiititiieuiueeuueFiFeeFF）（）（21max21max+=+=）（）（21212FiFuiuKCiM+=+运动方程 Fmax=载荷幅值 I =-1 =载荷函数的相位角 F1 =实部,Fmaxcos F2 =虚部,Fmaxsin umax=位移幅值 f =载荷函数的相位角 u1 =实部,umaxcos u2 =虚部,umaxsin 谐波载荷的本性 在已知频率下正弦变化；相角允许不同相的多个载荷同时作用，缺省值为零；施加的全部载荷都假设是简谐的，包括温度和重力。实部 虚部 复位移 在下列情

2、况下计算出的位移将是复数 具有阻尼 施加载荷是复数载荷（例如：虚部为非零的载荷）复位移滞后一个相位角（相对于某一个基准而言）可以用实部和虚部或振幅和相角的形式来查看 谐响应分析-术语和概念 求解方法 求解简谐运动方程的三种方法：完整法 为缺省方法，是最容易的方法；使用完整的结构矩阵，且允许非对称矩阵（例如：声学矩阵）。缩减法*使用缩减矩阵，比完整法更快；需要选择主自由度，据主自由度得到近似的 M矩阵和C矩阵。模态叠加法*从前面的模态分析中得到各模态；再求乘以系数的各模态之和；所有求解方法中最快的。求解方法 完整法完整法缩减法缩减法模态叠加法模态叠加法相对求解时间慢较快最快相对的使用容易程度最容

3、易较容易难允许元素载荷（例如压强）吗？允许不允许允许（一个载荷向量）允许非零位移载荷吗？允许允许不允许允许模态阻尼吗？不允许不允许允许能处理预应力吗？不能能能能进行“Restart“吗？能能不能允许非对称矩阵吗？允许不允许不允许需要为了求解而选择模态吗？不需要不需要需要需要选择主自由度吗？不需要需要需要（如果选用缩减法）步骤 四个主要步骤：建模 选择分析类型和选项 施加谐波载荷并求解 观看结果 建模 模型 只能用于线性单元和材料，忽略各种非线性；记住要输入密度；注意：如果ALPX（热膨胀系数）和T均不为零，就有可能不经意地包含了简谐热载荷。为了避免这种事情发生，请将ALPX设置为零.如果参考温

4、度 TREF与均匀节点温度 TUNIF不一致,那么T为非零值；建模 典型命令流/PREP7 ET,.MP,EX,.MP,DENS,!建立几何模型 !划分网格.选择分析类型和选项 选择分析类型和选项 进入求解器，选择谐响应分析；设置分析选项 1求解方法 2自由度输出格式 3是否使用集中质量逼近（用于结构的 一个方向的尺寸远小于另两个 方向的尺寸的情况中。例如：细长梁与薄壳。）典型命令:/SOLU ANTYPE,HARMIC,NEW HROPT,FULL HROUT,ON LUMPM,0 选择分析类型和选项 设置阻尼参数 从-阻尼、-阻尼和阻尼比中选 取 阻尼比最常用 典型命令：ALPHAD,0,

5、BETAD,0,DMPRAT,0.2,GUI:MainMenu SolutionLoad Step OptsTime/Frequenc Damping 施加谐波载荷并求解 施加谐波载荷并求解 所有施加的载荷以规定的频率（或频率范 围）简谐地变化“载荷”包括：位移约束-零或非零的 作用力 压强 注意：如果要施加重力和热载荷，它 们也被当作简谐变化的载荷来考虑！典型命令:DK,!或 D或DSYM DA,.DL,*AFUN,DEG FK,F,SFA,SFL,SFE,SF,施加谐波载荷并求解 规定谐波载荷时要包括：振幅和相角 频率 不同频率载荷具有不同的幅值时的处理 方法 载荷值（大小）代表振幅 Fm

6、ax 相角 是在两个或两个以上谐波 载荷间的相位差，单一载荷不需 要相角。实部 虚部 F2max 施加谐波载荷并求解 振幅和相角 ANSYS 不能直接输入振幅和相 角，而是规定实部和虚部分量。幅值/相角和实部/虚部的关系：可以使用APDL语言计算，但要确 保角度单位为度（缺省为弧度）。sincosampimagamprealFFFF=施加谐波载荷并求解 谐波载荷的频率：通过频率范围和在频率范围内的子步数量来规定每秒的循环次数（赫兹）；例如，在0-50频率范围内有10个子步时将给出在5，10，15.45和50Hz等频率上的解；而同一频率范围只有一个子步时，则只给出50Hz频率上的解。GUI:Ma

7、inMenu SolutionLoad Step OptsTime/Frequenc Freq and Substps 典型命令：HARFRQ,0,50,NSUBST,10,KBC,1 施加谐波载荷并求解 不同频率载荷具有不同的幅值时的处理方法 在施加谐波载荷后，下一步就是开始求解了，通常采用一个载荷 步，但是可以采用若干子步，且每个子步具有不同的频率范围。当 不同频率的载荷具有不同的幅值时，可以分多个载荷步施加。查看结果 1.绘制结构上的特殊点处的位移-频率曲线 2.确定各临界频率和相应的相角 3.观看整个结构在各临界频率和相角时的位移和应力 典型命令:/POST26 NSOL,PLVAR,

8、.查看结果 确定各临界频率 和相角 用图形显示最高振幅发生时的频率；由于位移与施加的载荷不同步（如果存在阻尼的话），需要确定出现最大振幅时的相角；要进行上述工作，首先要选择振幅+相位选项。查看结果 然后，列表显示变量。找出 响应最大值所对应的频率和 相角，据此确定最大振幅出 现的载荷步和子步号。进入POST1处理器，观察该 载荷子步下的结果，查看结 构的应力和变形情况。HRCPLX,PLDISP,2 PLNSOL,FINISH 谐响应分析实例 悬索拱桥的谐响应分析 输入文件：bridgeharmonic.cmds 305 307 谐响应分析实例 节点305、307在y向的 位移响应 紫色：30

9、7 蓝色：305 预应力谐响应分析 分为四个步骤 1.定义单元类型、材料属性、建模 2.进入求解器，指定分析类型静力分析,施加约束和预应力，定义预应力分析，求解。3.重新进入求解器（必须显式地，退出再进入），定义谐响应分析，删除预应力载荷，施加谐波载荷，求解。4.查看结果 注意：第2、3步中均开启预应力效应（PSTRES,ON）。预应力谐响应分析实例 张紧的吉他弦的谐响应分析 输入文件:presharmonic.cmds 126N预紧力 Y向谐波激励 预应力谐响应分析实例 预应力对响应（节点16，uy）的影响 无预应力 有预应力 谱分析 什么是谱分析？是模态分析的扩展，用于计算结构对地震及其它

10、随机激励的响应 在进行下述设计时要用到谱分析：建筑物框架及桥梁 太空船部件 飞机部件 承受地震或其它不稳定载荷的结构或部件 谱分析 谱分析的一种代替方法是进行瞬态分析，但是：瞬态分析很难应用于例如地震等随时间无规律变化载荷的分析;在瞬态分析中，为了捕捉载荷，时间步长必须取得很小，因而费时且昂贵.谱分析 主题包括：频谱的定义 响应谱如何用于计算结构对激励的响应：参与系数 模态系数 模态综合 谱分析 什么是频谱？用来描述理想化系统对激励响应的曲线，此响应可以是加速度、速度、位移和力;考虑安装于振动台上的四个单自由度弹簧质量系统它们的频率分别是f1，f2，f3及f4，而且f1f2f3 File Re

11、ad Input from.读取jobname.mcom文件 查看变形后的形状 应力和应变的图表显示 典型命令：/POST1/INPUT，mcom PLDISP，PLNSOL，FINISH 单点响应谱和多点响应谱 单点响应谱在模型的激励点（集）上指定一条响应谱曲线（族）。多点响应谱在不同的激励点（集）上指定不同的谱曲线（族）。谱分析实例 在这个实例中，将确定框架结构对于频谱激励的响应 输入文件：framespectrum.cmds 841 谱分析实例 在底座输入的激励谱 位移反应谱 节点841的位移响应 模态综合结果 模态综合后得到的应力和位移响应 随机振动分析 什么是随机振动分析?它是一种采

12、用功率谱密度作为输入的谱分析 它是一种确定响应出现特定值的概率大小的分析方法 输入值：结构的自然频率及模态 功率谱密度曲线（将在后面讨论）输出值：以1s位移和应力表示最可能出现的结构响应 随机振动分析 什么是PSD?PSD 记录了响应的平均平方值随频率变化的函数关系：用于描述响应出现特定值的概率从算术上看，PSD曲线下的面积即为响应的协方差（标准偏差的平方）响应可以是位移、速度、加速度、力或压力 随机振动分析步骤 六个主要步骤：建模 获得模态解 转换成谱分析类型 定义和施加功率谱密度激励 求解 察看结果 随机振动分析步骤 模型：本步的注意事项与模态分析时的相同 只考虑线性单元及材料，忽略各种非

13、线性 注意输入密度！同时如果出现与材料相关的阻尼，必须在这一步中定义 随机振动分析步骤 获得模态解 模态提取：只有Block Lanczos法、子空间法或缩减法是有效的 提取足够的模态以包含频谱的频率范围 所有模态的扩展只有扩展的模态才能用于频谱的求解 典型命令：MODOPT，MXPAND，.随机振动分析步骤 载荷及边界条件：对于基础激励，一定要约束适当的自由度 对于压力PSD，在此步中施加压力于想要的表面上 文件：.mode文件含有特征向量而且是对于频谱求解所需要的 典型命令：DK，.!或 D 或 DSYM DL，DA，.!求解 SOLVE 随机振动分析步骤 转换成谱分析类型：退出并且重新进

14、入求解器 进行新的分析：谱分析 分析选项：将在后面讨论 阻尼：将在后面讨论 典型命令：FINISH /SOLU ANTYPE，SPECTR 随机振动分析步骤 分析选项：频谱类型：功率谱密度（PSD）模态数量：如果选项是 0或空缺，所有扩展的模态将被用于求解 单元计算：只有在模态分析中的单元计算处于开启（ON）状态的条件下，谱分析中的计算也才能处于开启（ON）典型命令：SPOPT，PSD，.随机振动分析步骤 阻尼：所有四种形式均可采用：-阻尼（质量阻尼）-阻尼（刚度阻尼）恒定阻尼比 与频率相关的阻尼比（模态阻尼）如果没有指明阻尼，ANSYS将以1%的恒定阻尼作为缺省值 典型命令：ALPHAD，B

15、ETAD，DMPRAT，MDAMP，随机振动分析步骤 定义并施加功率谱密度激励：指明功率谱密度设置 定义功率谱密度对频率的数据表 在想要的节点施加激励 随机振动分析步骤 功率谱密度设置：频谱类型（单位）：加速度（常规单位或 g2/Hz）速度 位移 力 压力 表号的缺省值为1，对于多条功率谱密度曲线，使用表号 典型命令：PSDUNIT，.随机振动分析步骤 功率谱密度对频率的数据表：指明表号（通常为1）然后一对一对地输入频率和功率谱密度值 典型命令：PSDFRQ，PSDVAL，.随机振动分析步骤 施加功率谱密度：施加步骤取决于功率谱密度类型 加速度、速度或位移功率谱密度：这些是基础上的激励并且只能

16、施加于以前约束过的节点上 在 UX，UY或UZ（激励方向）方向上作为约束来施加，其数值为1.0 典型命令：D，!值=1.0 随机振动分析步骤 作用力功率谱密度：节点激励 作为数值为1.0（或希望的比例因子）的作用力施加在 FX，FY或FZ（激励方向）上 压力功率谱密度：要求在模态分析中即施加压力 使用载荷矢量（在模态求解时计算出的）来施加压力功率谱密度激励 将数值设定为1.0或为希望的比例因子 典型命令：F，!对于力 PSD LVSCALE，!对于压力PSD（在模态分析步骤定义压力）随机振动分析步骤 求解：激活功率谱密度模态组合方法 指定计算的项目*计算参与系数*启动功率谱密度求解器*将在后面

17、讨论 典型命令：PSDCOM，随机振动分析步骤 将计算的项目：缺省的计算项目是计算相对于基础激励的位移解（包括应力和应变）相对于基础的或绝对的速度和加速度解也是可以获得的 随机振动分析步骤 计算参与因子：对于定义的每一个功率谱密度表必须计算 指定基础或节点激励 启动功率谱密度求解器：结果写入.rst文件中 典型命令：PFACT，SOLVE FINISH 随机振动分析步骤 察看结果：绘图和列表显示 1 时的各参数时的数值（POST1）生成响应功率谱密度（POST26）计算两个量之间的协方差（POST26）随机振动分析步骤 随机振动的结果为 1 时的各种参数的数值：1时的位移，1 时的应力等等 所

18、有的量都假定是平均值零的高斯（正态）分布 例如，最大位移 Umax=0.15表示Uimax 等为0.15 或更小值的机率为 68%（1）它同时也表示为：Umax等为0.15x2=0.3或更小的机率 为95%（2）Umax 等为 0.15x3=0.45 或更小的机率为98%（3）1 2 3 随机振动分析步骤 为了察看1s时的位移及应力：进入 POST1（常规后处理器）从载荷步3中读取结果，载荷步3是1s 结果存于结果文件中的位置 注意：如果需要的话，1s 速度和1s 加速度是分别存于载荷步4和载荷步5中 然后以图和表显示需要的物理量 随机振动分析步骤 结果文件最多可以有5个载荷步:载荷步 1：模

19、态结果，每阶模态一个子步 载荷步2：单位静态解，每个PSD表一个子步.如果不是基础激励则为空 载荷步3：1 位移解 载荷步4：1 速度解.如果没要求输出则为空.载荷步5：1 加速度解.如果没要求输出则为空.随机振动分析步骤 结果文件最多可以有5个载荷步:1 加速度解.如果没要求输出则为空.典型命令/POST1 SET，3 !1s 位移解 PLDISP，PLNSOL，FINISH 随机振动分析步骤 1 结果主要用于：第一步的主要计算：在给定的时间间隔内，某一自由度上的位移超过0.32个单位的概率是多少？疲劳计算：基于下述假设：应力水平在时间的68%以内为 1 ，在时间的27%（95-68）以内为

20、 2 ，在时间的3%（98-95）以内为 3 基于 Monte Carlo 仿真的可靠性（风险性）分析 随机振动分析步骤 响应功率谱密度：使工程师有某种响应量（如应力）是如何随频率变化的观念 结果文件包括 1s 时各量的数值，它是功率谱密度曲线下面积的平方根 POST26，时间-历程后处理器用于计算响应功率谱密度 随机振动分析步骤 为了计算响应功率谱密度：1.进入POST26，并首先存储频率矢量：在某个自然频率的两边各取1到10 个附加的数据点以绘制一条更光滑的频 率曲线，缺省值是各取5个附加的 数据点，变量1自动地赋给频率变量。典型命令：/POST26 STORE，PSD，.随机振动分析步骤

21、 2.确定结果中要对其进行功率谱密度响应计算的物理量：TimeHist Postpro Define Variables.可以是任何节点的或单元的结果项 典型命令：NSOL，ESOL，.随机振动分析步骤 3.计算并用图形显示功率谱密度响应：TimeHist Postpro Calc Resp PSD.TimeHist Postpro Graph Variables 典型命令：RPSD，PLVAR，.协方差：协方差表示的是两个量之间的相互关系 可在任何两个响应量之间进行计算，例如：计算模型中两个不同点之间的应力的协方差 用POST26时间-历程后处理器来计算协方差 为了计算协方差：1.重新设置或

22、退出POST26 并重新进入 POST26 2.确定对其要计算协方差的两个响应量 典型命令：FINISH/POST26 NSOL，ESOL，.3.计算并取出协方差：TimeHist Postpro Calc Covariance.使用*GET 命令来取出协方差：*GET，COVAR，VARI，#，EXTREM，CVAR -或-Utility Menu Parameters Get Scalar Data.典型命令：CVAR，*GET，FINISH 随机振动分析实例 在这个实例中，将确定电路板对于随机激励的响应 输入文件：psddemo.cmds 随机振动分析实例 有限元模型 随机振动分析实例 位移响应 1