ansys动力学瞬态分析详解
- 格式:doc
- 大小:445.50 KB
- 文档页数:73
ANSYS瞬态动力学分析完整教程瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。
可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。
载荷和时间的相关性使得惯性力和阻尼作用比较重要。
如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。
瞬态动力学的基本运动方程是:其中:[M] =质量矩阵[C] =阻尼矩阵[K] =刚度矩阵{ }=节点加速度向量{ }=节点速度向量{u} =节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力([M]{ })和阻尼力([C]{ })的静力学平衡方程。
ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。
两个连续时间点间的时间增量称为积分时间步长(integrationtime step)。
?3.2学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按―工程‖时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。
例如,可以做以下预备工作:1.首先分析一个较简单模型。
创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。
在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。
通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。
同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。
<<高级技术分指南>>中将讲述子结构。
?3.3三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。
完全法完全法采用完整的系统矩阵计算瞬态响应(没有矩阵缩减)。
它是三种方法中功能最强的,允许包括各类非线性特性(塑性、大变形、大应变等)。
注─如果并不想包括任何非线性,应当考虑使用另外两种方法中的一种。
这是因为完全法是三种方法中开销最大的一种。
完全法的优点是:·容易使用,不必关心选择主自由度或振型。
·允许各种类型的非线性特性。
·采用完整矩阵,不涉及质量矩阵近似。
·在一次分析就能得到所有的位移和应力。
·允许施加所有类型的载荷:节点力、外加的(非零)位移(不建议采用)和单元载荷(压力和温度),还允许通过TABLE数组参数指定表边界条件。
·允许在实体模型上施加的载荷。
完全法的主要缺点是它比其它方法开销大。
§3.4 完全法瞬态动力学分析首先,讲述完全法瞬态动力学分析过程,然后分别介绍模态叠加法和缩减法与完全法不相同的计算步骤。
完全法瞬态动力分析(在ANSYS/Multiphsics、ANSYS/Mechauioal及ANSYS/Structural中可用)由以下步骤组成:1.建造模型2.建立初始条件3.设置求解控制4.设置其他求解选项5.施加载荷6.存储当前载荷步的载荷设置7.重复步骤3-6定义其他每个载荷步8.备份数据库9.开始瞬态分析10.退出求解器11.观察结果§ 型在这一步中,首先要指定文件名和分析标题,然后用PREP7定义单元类型,单元实常数,材料性质及几何模型。
这些工作在大多数分析中是相似的。
<<ANSYS 建模与网格指南>>详细地说明了如何进行这些工作。
对于完全法瞬态动力学分析,注意下面两点:·可以用线性和非线性单元;·必须指定杨氏模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量)。
材料特性可以是线性的或非线性的、各向同性的或各向异性的、恒定的或和温度有关的。
§3.1瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。
可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。
载荷和时间的相关性使得惯性力和阻尼作用比较重要。
如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。
瞬态动力学的基本运动方程是:其中:[M] =质量矩阵[C] =阻尼矩阵[K] =刚度矩阵{}=节点加速度向量{}=节点速度向量{u} =节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力([M]{})和阻尼力([C]{})的静力学平衡方程。
ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。
两个连续时间点间的时间增量称为积分时间步长(integration time step)。
§3.2学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。
例如,可以做以下预备工作:1.首先分析一个较简单模型。
创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。
在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。
通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。
同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。
<<高级技术分指南>>中将讲述子结构。
§3.3三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。
§3.1瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。
可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。
载荷和时间的相关性使得惯性力和阻尼作用比较重要。
如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。
瞬态动力学的基本运动方程是:其中:[M] =质量矩阵[C] =阻尼矩阵[K] =刚度矩阵{}=节点加速度向量{}=节点速度向量{u} =节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力([M]{})和阻尼力([C]{})的静力学平衡方程。
ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。
两个连续时间点间的时间增量称为积分时间步长(integration time step)。
§3.2学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。
例如,可以做以下预备工作:1.首先分析一个较简单模型。
创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。
在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。
通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。
同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。
<<高级技术分指南>>中将讲述子结构。
§3.3三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。
•ANSYS瞬态动力学分析理论基础本文主要介绍了ansys软件进瞬态动力分析与计算的理论,通过介绍使读者可以更好的理解软件和操作软件以便进行相关的分析。
一假设和限制1、系统的初始条件已知,即速度和位移。
2、结构瞬态分析中当需要时可以考虑陀螺或科里奥力效应。
二结构和其他二阶系统分析对于线性结构的瞬态动力学平衡方程:(1)ANSYS里使用两种方法求解方程(1):向前差分时间积分和Newmark积分(包括改进后的算法称为HHT)。
向前差分方法适用于求解显示的瞬态分析。
Newmark和HHT方法使用隐式方法来求解瞬态问题。
Newmark方法使用有限差分法,在一个时间间隔内有,(2)(3)其中:α,δ:Newmark积分参数我们主要的目的就是计算下一时刻的位移u n+1,则在t n+1时刻的控制方程(1)为:(4)为了求解u n+1,可以把(2)和(3)重新排列,得(5)(6)其中:注意到(5)代入到(6)中,则,可以通过u n+1求出。
由(5)、(6)和(4)得(7)一旦求出u n+1,速度和加速度可以利用(5)和(6)求得。
对于初始施加于节点的速度或加速度可以利用位移约束并利用(3)计算得到。
根据Zienkiewicz的理论,利用(2)和(3)式得到的Newmark求解方法的无条件稳定必须满足:(8)Newmark参数根据下式输入:(9)其中:γ:振幅衰减因子通过观察(8)和(9)可以发现无条件稳定也可以表述为,并且γ≥0。
因此只要γ≥0,则求解就是稳定的。
对于压电分析参数设置为:α=0.25;δ=0.5并且θ=0.5。
通常情况下衰减因子γ=0.005。
当γ=0时即α=0.25,δ=0.5时Newmark方法为平均加速度法。
由于平均加速度法在位移幅值误差方面不产生任何数值阻尼。
如果其他方面也没有阻尼,缺乏数值阻尼在高频结构计算中会产生不可接受的数值噪声。
我们期望有一定水平的数值阻尼并且通过设置γ>0来实现。
ANSYS培训教程:瞬态动力学分析的基本步骤用不同的瞬态动力学方法进行分析时,进行瞬态动力学分析的过程不尽相同。
下面我们首先描述如何用完全法进行瞬态动力学分析的基本步骤,然后在列出用缩减法和模态叠加法时的不同地方。
完全法瞬态动力学分析过程由三个主要步骤组成:1.建模2.加载及求解3.结果后处理模型的建立建模过程和其它类型的分析类似,但应注意以下几点:1.可以用线性和非线性单元。
2.必须指定弹性模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量),材料特性可以是线性的或非线性的,各向同性的或各向异性的,恒定的或和温度有关的。
在划分网格时需要记住以下几点:1. 有限元网格需要足够精度以求解所关心的高阶模态;2. 感兴趣的应力应变区域的网格密度要比只关系位移的区域相对加密一些;3.如果想包含非线性,网格应当细到能够扑捉到非线性效果。
例如,对于塑性分析来说,它要求在较大塑性变形梯度的平面内有一定的积分点密度,所以网格必须加密;4.如果对波传播效果感兴趣,网格应当细到足以解算出波。
基本准则是沿波的传播方向每一波长至少有20个单元。
加载并求解在这一步中,要定义分析类型及选项,加载,指定载荷步选项,并开始有限元求解。
具体步骤如下:1.进ANSYS求解器命令:/SOLUGUI:Main Menu | Solution2.指定分析类型和分析选项(1)指定分析类型(ANTYPE)选择开始一个新的分析。
如果已经完成了静力学预应力或完全法瞬态动力学分析并准备对时间历程进行延伸,或者想重新启动一次失败的非线性分析,则可用Restart。
(Restart要求初始求解过程中生成的文件Jobname.EMAT,Jobname.ESAV及Jobname.DB存在。
新得到的解结果将被附加在初始结果文件Jobname.RST中)。
从弹出的对话框中选择瞬态动力学分析(Transient),并指定位完全法(Full)。
对于质量阵形成方法(Mass Matrix Formulation)建议在大多数分析应用中采用缺省的质量阵形成方式。
一、概述在工程领域中,研究物体相对移动的动力学行为具有重要意义。
在实际工程应用中,瞬态动力学分析是评估机械设计的重要手段之一。
本文将通过ANSYS软件进行瞬态动力学设置,研究两物体相对移动的问题。
二、瞬态动力学分析基本原理1. 瞬态动力学分析瞬态动力学分析是指在物体受到外部力或扭矩作用下,物体产生瞬时运动或者受到瞬时力的影响时的动力学分析方法。
该方法适用于应用于工程领域中需要考虑加速度、惯性力、阻尼等瞬态动力学因素的问题。
2. ANSYS软件ANSYS软件是一种用于工程仿真和设计的有限元分析软件。
它能够模拟和分析多种工程问题,包括结构分析、热分析、流体力学分析等。
在瞬态动力学分析中,ANSYS软件可以模拟物体的瞬时运动、应力分布等。
三、两物体相对移动问题分析1. 问题描述假设有两个物体A和B,它们之间通过一根弹簧相连。
当施加外力使得物体A移动时,弹簧会受到拉力,同时对物体B施加相等反作用力。
我们希望通过瞬态动力学分析,研究物体A和B在相对移动过程中的动力学行为。
2. ANSYS设置我们需要建立物体A和B的几何模型,并在ANSYS中导入。
根据物体的材料属性、外部力的施加情况等,设置瞬态动力学分析的条件和参数。
在设置过程中,需注意考虑物体的刚度、弹簧的刚度、阻尼等因素。
3. 模拟过程在模拟过程中,我们可以通过ANSYS软件对物体A施加外力,观察物体A和B在相对移动过程中的运动状态、应力分布等动力学行为。
通过分析模拟结果,可以得出两物体在相对移动过程中所受到的动力学影响。
四、模拟结果分析1. 动态响应通过模拟分析,我们可以观察到物体A受到外力作用后的瞬时加速度、速度和位移变化。
物体B也会在弹簧的作用下产生相对运动。
通过观察动态响应,我们可以得出两物体相对移动的动力学特性。
2. 应力分布在瞬态动力学分析中,我们还可以观察到物体A和B在相对移动过程中受到的应力分布情况。
弹簧受到的拉力、物体产生的应力等都可以得到清晰的分析和展示。
第10章 瞬态动力学分析
瞬态动力学分析(亦称时间历程分析)是用于确定承受任意随时间变化的载荷的结构动力学响应的一种方法。
利用瞬态动力学分析可以确定结构在静载荷、瞬态载荷和简谐载荷的随意组合下随时间变化产生的位移、应变、应力及力。
★ 了解瞬态动力学分析。
10.1 瞬态动力学分析概述
瞬态动力学分析(Transient Structural Analysis)给出的是结构关于时间载荷的响应,它不同于刚体动力学分析,在Workbench中瞬态动力学的模型可以是刚体,也可以是柔性体,而对于柔性体可以考虑材料的非线性特征,由此可得出柔性体的应力和应变值。
在进行瞬态动力学分析时,需要注意:
当惯性力和阻尼可以忽略时,采用线性或非线性的静态结构分析来代替瞬态动力学分析。
当载荷为正弦形式时,响应是线性的,采用谐响应分析更为有效。
当几何模型简化为刚体且主要关心的是系统的动能时,采用刚体动力学分析更为有效。
除上述三种情况外,其余情况均可采用瞬态动力学分析,但其所需的计算资源较其他方法要大。
10.2 瞬态动力学分析流程
在ANSYS Workbench左侧工具箱中Analysis
Systems下的Transient Structural上按住鼠标左键拖动到
项目管理区的A6栏,即可创建瞬态动力学分析项目,
如图10-1所示。
当进入Mechanical后,单击选中分析树中的
Analysis Settings即可进行分析参数的设置,如图10-2
图10-1 创建瞬态动力学分析项目。
§3.1瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。
可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。
载荷和时间的相关性使得惯性力和阻尼作用比较重要。
如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。
瞬态动力学的基本运动方程是:其中:[M] =质量矩阵[C] =阻尼矩阵[K] =刚度矩阵{}=节点加速度向量{}=节点速度向量{u} =节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力([M]{})和阻尼力([C]{})的静力学平衡方程。
ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。
两个连续时间点间的时间增量称为积分时间步长(integration time step)。
§3.2学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。
例如,可以做以下预备工作:1.首先分析一个较简单模型。
创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。
在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。
通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。
同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。
<<高级技术分指南>>中将讲述子结构。
§3.3三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。
ANSYS/Professional产品中只允许用模态叠加法。
在研究如何实现这些方法之前,让我们先探讨一下各种方法的优点和缺点。
§3.3.1完全法完全法采用完整的系统矩阵计算瞬态响应(没有矩阵缩减)。
它是三种方法中功能最强的,允许包括各类非线性特性(塑性、大变形、大应变等)。
注─如果并不想包括任何非线性,应当考虑使用另外两种方法中的一种。
这是因为完全法是三种方法中开销最大的一种。
完全法的优点是:·容易使用,不必关心选择主自由度或振型。
·允许各种类型的非线性特性。
·采用完整矩阵,不涉及质量矩阵近似。
·在一次分析就能得到所有的位移和应力。
·允许施加所有类型的载荷:节点力、外加的(非零)位移(不建议采用)和单元载荷(压力和温度),还允许通过TABLE数组参数指定表边界条件。
·允许在实体模型上施加的载荷。
完全法的主要缺点是它比其它方法开销大。
§3.3.2模态叠加法模态叠加法通过对模态分析得到的振型(特征值)乘上因子并求和来计算结构的响应。
此法是ANSYS/Professional程序中唯一可用的瞬态动力学分析法。
模态叠加法的优点是:·对于许多问题,它比缩减法或完全法更快开销更小;·只要模态分析不采用PowerDynamics方法,通过 LVSCALE 命令将模态分析中施加的单元载荷引入到瞬态分析中;·允许考虑模态阻尼(阻尼比作为振型号的函数)。
模态叠加法的缺点是:·整个瞬态分析过程中时间步长必须保持恒定,不允许采用自动时间步长;·唯一允许的非线性是简单的点点接触(间隙条件);·不能施加强制位移(非零)位移。
§3.3.3缩减法缩减法通过采用主自由度及缩减矩阵压缩问题规模。
在主自由度处的位移被计算出来后,ANSYS可将解扩展到原有的完整自由度集上。
(参见“模态分析”中的“矩阵缩减”部分对缩减过程的详细讨论。
)缩减法的优点是:·比完全法快且开销小。
缩减法的缺点是:·初始解只计算主自由度的位移,第二步进行扩展计算,得到完整空间上的位移、应力和力;·不能施加单元载荷(压力,温度等),但允许施加加速度。
·所有载荷必须加在用户定义的主自由度上(限制在实体模型上施加载荷)。
·整个瞬态分析过程中时间步长必须保持恒定,不允许用自动时间步长。
·唯一允许的非线性是简单的点—点接触(间隙条件)。
§3.4 完全法瞬态动力学分析首先,讲述完全法瞬态动力学分析过程,然后分别介绍模态叠加法和缩减法与完全法不相同的计算步骤。
完全法瞬态动力分析(在ANSYS/Multiphsics、ANSYS/Mechauioal及ANSYS/Structural中可用)由以下步骤组成:1.建造模型2.建立初始条件3.设置求解控制4.设置其他求解选项5.施加载荷6.存储当前载荷步的载荷设置7.重复步骤3-6定义其他每个载荷步8.备份数据库9.开始瞬态分析10.退出求解器11.观察结果§3.4.1建造模型在这一步中,首先要指定文件名和分析标题,然后用PREP7定义单元类型,单元实常数,材料性质及几何模型。
这些工作在大多数分析中是相似的。
<<ANSYS 建模与网格指南>>详细地说明了如何进行这些工作。
对于完全法瞬态动力学分析,注意下面两点:·可以用线性和非线性单元;·必须指定杨氏模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量)。
材料特性可以是线性的或非线性的、各向同性的或各向异性的、恒定的或和温度有关的。
划分合理的网格密度:·网格密度应当密到足以确定感兴趣的最高阶振型;·对应力或应变感兴趣的区域比只考察位移的区域的网格密度要细一些;·如果要包含非线性特性,网格密度应当密到足以捕捉到非线性效应。
例如,塑性分析要求在较大塑性变形梯度的区域有合理的积分点密度(即要求较密的网格);·如果对波传播效果感兴趣(例如,一根棒的末端准确落地),网格密度应当密到足以解算出波动效应。
基本准则是沿波的传播方向每一波长至少有20个单元。
§3.4.2建立初始条件在执行完全法瞬态动力学分析之前,用户需要正确理解建立初始条件和正确使用载荷步。
瞬态动力学分析顾名思义包含时间函数的载荷。
为了定义这样的载荷,用户需要将载荷—时间关系曲线划分成合适的载荷步。
载荷—时间曲线上的每个“拐角”对应一个载荷步,如图3.1所示。
图3.1载荷—时间关系曲线第一个载荷步通常被用来建立初始条件,然后为第二和后继瞬态载荷步施加载荷并设置载步选项。
对于每个载荷步,都要指定载荷值和时间值,同时指定其它的载荷步选项,如采用阶梯加载还是斜坡加载方式施加载荷以及是否使用自动时间步长等。
然后,将每个载荷步写入载荷步文件,最后一次性求解所有载荷步。
施加瞬态载荷的第一步是建立初始条件(即零时刻时的情况)。
瞬态动力学分析要求给定两种初始条件(因为要求解的方程是两阶的):初始位移()和初始速度()。
如果没有进行特意设置,和都被假定为0。
初始加速度()一般假定为0,但可以通过在一个小的时间间隔内施加合适的加速度载荷来指定非零的初始加速度。
下面的段落描述了如何施加不同组合形式的初始条件。
§3.4.2.1零初始位移和零初始速度这是缺省的初始条件,即如果= = 0,则不需要指定任何条件。
在第一个载荷步中可以加上对应于载荷/时间关系曲线的第一个拐角处的载荷。
§3.4.2.2非零初始位移及/或非零初始速度可以用IC命令设置这些初始条件。
命令:ICGUI:Main Menu>Solution>-Loads-Apply>Initial Condit’n>Define注意:不要定义矛盾的初始条件。
例如,在某单一自由度处定义了初始速度,则在所有其它自由度处的初始速度将为0.0,潜在地会产生冲突的初始条件。
在大多数情形下要在模型的每个未约束自由度处定义初始条件。
如果这些条件对各自由度是不同的,那么就可以较容易地明确指定初始条件,如下所述。
关于TIMINT和IC命令的说明参见<<ANSYS命令参考手册>>。
非零速度是通过对结构中需指定速度的部分加上小时间间隔上的小位移来实现的。
比如如果=0.25,可以通过在时间间隔0.004内加上0.001的位移来实现,命令流如下:...TIMINT,OFF! Time integration effects offD,ALL,UY,.001! Small UY displ. (assuming Y-direction velocity)TIME,.004! Initial velocity = 0.001/0.004 = 0.25LSWRITE! Write load data to load step file (Jobname.S01) DDEL,ALL,UY! Remove imposed displacementsTIMINT,ON! Time integration effects on...§3.4.2.4非零初始位移和非零初始速度和上面的情形相似,不过施加的位移是真实数值而非“小”数值。
比如,若= 1.0且= 2.5,则应当在时间间隔0.4内施加一个值为1.0的位移:...TIMINT,OFF! Time integration effects offD,ALL,UY,1.0! Initial displacement = 1.0TIME,.4! Initial velocity = 1.0/0.4 = 2.5LSWRITE! Write load data to load step file (Jobname.S01) DDELE,ALL,UY! Remove imposed displacementsTIMINT,ON! Time integration effects on...需要用两个子步[NSUBST,2]来实现,所加位移在两个子步间是阶跃变化的[KBC,1]。
如果位移不是阶跃变化的(或只用一个子步),所加位移将随时间变化,从而产生非零初速度。
下面的例子演示了如何施加初始条件= 1.0,= 0.0:...TIMINT,OFF! Time integration effects off for staticsolutionD,ALL,UY,1.0! Initial displacement = 1.0TIME,.001! Small time intervalNSUBST,2! Two substepsKBC,1! Stepped loadsLSWRITE! Write load data to load step file (Jobname.S01)!transient solutionTIMINT,ON! Time-integration effects on for transientsolutionTIME,...! Realistic time intervalDDELE,ALL,UY! Remove displacement constraintsKBC,0! Ramped loads (if appropriate)!Continue with normal transient solution procedures...§3.4.2.6非零初始加速度可以近似地通过在小的时间间隔内指定要加的加速度[ACEL]实现。