ANSYS在斜拉桥的稳定分析中的应用
- 格式：pdf
- 大小：84.08 KB
- 文档页数：2
3. 5 高应变动力试桩法
在我国 ,高应变动力试桩法的研究是起自 20 世纪 80 年代中 后期 ,90 年代初期已有相关的软硬件问题 ,其实际应用效果已不 弱于国外 ,其后面向国内大量的灌注桩检测 ,已有单位在模型改 进 、拟合技巧 、参数选定等方面进行了大量工作 ,也有应用者在桩 如何才算被充分激发方面进行了研究 。值得一提的是 ,桩基动测 方面 ,国产仪器和软件业已达到国际先进水平 ,许多方面甚于更 具有中国特色 。
[ K ] DL 考虑了几何非线性 , 因为考虑几何非线性后 , 变形导 致坐标参考系的改变 ,因此构成刚度矩阵的几何参数变化受{Δδ} 影响 ,构成非线性关系 。
1. 2 考虑几何非线性和材料非线性法[2]
考虑几何非线性和材料非线性的斜拉桥结构增量平衡方程
和计算情形与几何非线性一致 ,但是切线刚度矩阵发生变化 : [ K ] T = [ K ] DL + [ K ] ep + [ K ]σ。
33 (19) :1032104. [ 3 ] 付建新. 桩基检测研析[J ]. 山西建筑 ,2007 ,33 (20) :1082109.
The construction and detection of highway bridge pile f oundation
L I Jian2zhong Abstract : Combined wit h t he application of cast in place bored pile and manual digging pile in highway bridge pile foundation , t he key matters of bot h file foundation construction were discussed. The main technology of present pile foundation detection was introduced so as to complete construction technology of highway bridge pile foundation and enhance pile foundation construction quality. Key words : manual digging pile , cast in place bored pile , construction , detection , quality
山 西 建 筑 第 35 200
卷第 9年
1 1
期 月
SHAN
XI
ARC HI T EC
TU R
E
Vol . Jan.
35 No . 2009
1
·329 ·
文章编号 :100926825 (2009) 0120329202
ANSYS 在斜拉桥的稳定分析中的应用
2. 1 特征值屈曲分析方法
ANSYS 中特征值屈曲分析的步骤如下 : 1) 建立模型 。这一步和通常的建模方法相同 ,进入前处理器 定义单元类型 、单元实常数 、材料性质 、模型几何形状 。但是建立
的模型有如下要求 :a. 只允许线性行为 ,如果定义了非线性单元 , 则将按线性对待 ; b. 材料的弹性模量或某种形式的刚度必须定 义 ,非线性的材料性质即使被定义也将被忽略 。
4) 展开解 。展开解就是根据需要选定前几阶特征值求对应 的特征向量 。考虑到只有最小特征值才有实际意义 ,本文仅取了 第一阶特征值 。
5) 观察结果 。进入一般后处理器查看各阶特征值对应的变 形图即为对应的失稳模态 。因为特征值分析的结果本身是计算 值的上限 ,且只允许线性行为 。
2. 2 特征值屈曲分析方法
3. 4 应力波反射法完整性检测
尽管近年来国内外对于这种方法的研究未见本质性的进展 , 但在实用和普及方面国内却有较大提高 ,这些不仅表现在国产桩 基动测仪和配套用传感已达到或接近国外先进仪器方面 ,也表现 在许多单位认真研究各个测试细小环节和分析环节方面 ,更主要 的是表现在许多管理部门已开始认真总结应力波反射法完整性 检测的得与失 ,开始使这种方法的应用回归到一种正常的位置 , 如广东省正在将这种方法定位于为后续检测提供前期技术准备 , 这种定位已与该方法解决问题的真实能力完全对应 。
3. 6 动静法
由于高应变动力试桩法力的作用时间过短 ,桩只能被视为弹 性体进行分析 ,国外有人提出了一种动静法 ,采用技术将力的作 用时间延长 ,使沿桩身传播的应力波波长大于实际桩长 ,进而将 桩视为刚体 ,回避了应力波的传播问题 。应该说这种方法既克服 了传统静载试验的笨重与费时 ,也克服了高应力方法的过分间接 性 ,是一种较好的方法 ,但由于该方法对锤的配重要求太高 ,具体 操作仍有较大难度 。
收稿日期 :2008209211 作者简介 :鄢生全 (19832 ) ,男 ,助理工程师 ,广东长大公路工程有限公司第一分公司 ,广东 广州 511431
柴海峰 (19832 ) ,男 ,助理工程师 ,中交二航局二公司 ,重庆 400074 田浩亮 (19822 ) ,男 ,助理工程师 ,中交第一公路勘察设计研究院有限责任公司 ,陕西 西安 710075
1 第一类稳定和第二类稳定计算的有限元实现
在有限元计算中 ,考虑了结构的几何刚度 ,该刚度反映了结 构的不稳定因素 。当外力增加 λ倍时 ,杆力和几何刚度矩阵也增 大了 λ倍 ,因而可以写出下式 :
( [ K ] D +λ[ K ]σ) {δ} =λ{ F} 。 如果 λ足够大 ,使得结构达到随遇平衡状态 , 节点位移矩阵 由{δ} 变为{δ} + {Δδ} ,上式平衡方程也能满足 ,则 :
其中 ,[ K ] ep为弹塑性刚度矩阵 。 第二类稳定的安全系数可以采用下式计算 :
λ=
Pcritical 。 Pexisting
其中 , Pcritical为结构的临界荷载 ; Pexisting为结构的实际荷载 。
2 斜拉桥稳定计算在 ANSYS 中的实现
ANSYS 程序提供两种结构屈曲载荷和屈曲模态的分析方法 为 :特征值屈曲分析和非线性屈曲分析 。
1. 1 几何非线性法[2 ]
斜拉桥结构几何非线性增量平衡方程为 : [ K ] T{Δδ} = {Δ F} 。
[ K ] T = [ K ] DL + [ K ]σ。
其中 ,[ K ] T 为结构切线刚度矩阵 ; [ K ] DL 为考虑结构非线性 的总体弹性刚度矩阵 ;{Δ F} 为外荷载增量 。
© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
·330 ·
第 35 200
卷第 9年
1 1
期 月
山
西
建ห้องสมุดไป่ตู้
筑
( [ K ] D +λ[ K ]σ) ({δ} + {Δδ} ) =λ{ F} 。
同时满足上面两式的条件是 : ( [ K ] D +λ[ K ]σ) {Δδ} = { 0} 。
仅在使用的方便程度上有了质的飞跃 ,而且在分析手段上也有了 很大提高 ,声失时判读已不再是唯一的选择 ,声幅和声频已开始 进入了分析判断领域 ,尤其令人欣慰的是 ,声波 CT 已步入实用阶 段 ,为声波透射法的后续研究提供了广阔的前景 。
ANSYS 中几何非线性屈曲分析的步骤如下 : 1) 建立模型 。 2) 求解 。 a. 进入求解器 ,施加约束 ; b. 特征值屈曲分析 ;c. 根据特征值 分析的结果施加荷载和初始缺陷 (通常是将单倍荷载乘以第一阶 特征值作为几何非线性分析的荷载上限来施加 ,根据特征值分析 的第一阶失稳模态给结构施加初始缺陷) ;d. 设定分析选项 ,主要 包括 打 开 大 变 形 开 关 ( NL GEOM , ON) , 牛 顿 —拉 普 森 选 项 (NROP T ,ON) 等 ;e. 指定荷载步选项 ,包括指定 TIM E 值和 AU2 TO T 自动荷载分步项等 ;f . 设定收敛准则 ; g. 确定最大平衡迭代 次数和求解终止准则 ;h. 求解发散 ,在后处理中检查构件内力 ,查 看各子步计算结果中单元受力状态 ,如果第 K 子步构件实际受 力超过了该构件的极限承载力 ,则记录该构件对应的单元 ,在以 后的计算中将这些单元杀死 ,重启动至第 K - 1 子步计算结束状 态并继续加载 ,一直循环计算到不再删除单元 ,计算结束 。 3) 进入一般后处理器或时间历程后处理器查看结果 。
其中 ,[ K ] D 为结构整体弹性刚度矩阵 ; [ K ]σ 为结构整体几 何刚度矩阵 ;λ为整体安全系数 。
这就是计算稳定安全系数的特征方程式 ,如果方程有 n 阶 , 那么理论上存在 n 个特征值λ1 ,λ2 , …,λn ,工程问题中只有最低 阶特征值才有意义 。
由于第一类稳定问题是分支点失稳 ,实际工程结构都存在初 始缺陷 ,是极值失稳问题即第二类稳定 。由于斜拉桥存在几何非 线性和材料非线性 ,因此 ,第二类稳定问题必须考虑线性影响 。
2) 获得静力解 。先定义约束及荷载 ,并求解 。在这一步里必 须打开预应力开关命令 ,以记录内力并求得应力刚度矩阵 。求解 完成后退出求解器 。
3) 获得特征值屈曲解 。重新进入求解器 :定义分析类型为特 征值屈曲分析 ;定义分析选项 ,分析选项包括特征值提取方式 、需 提取的特征值数目 、特征值计算的起始点 ;定义荷载步选项 ;求解 。
非线性稳定分析的实质是求全过程 P —Δ 曲线 ,其中 ,Δ 为 关心节点的位移 ; P 为加载系数 ,取值 0～1 ,其含义是迭代荷载与 预加极限荷载的比值 。非线性稳定包括几何非线性稳定和材料 非线性稳定 。
几何非线性屈曲分析可以直接利用特征值屈曲分析的模型 , 因为非线性的材料性质将被忽略 。
在我国 ,高应变动力试桩法的研究是起自 20 世纪 80 年代中 后期 ,90 年代初期已有相关的软硬件问题 ,其实际应用效果已不 弱于国外 ,其后面向国内大量的灌注桩检测 ,已有单位在模型改 进 、拟合技巧 、参数选定等方面进行了大量工作 ,也有应用者在桩 如何才算被充分激发方面进行了研究 。值得一提的是 ,桩基动测 方面 ,国产仪器和软件业已达到国际先进水平 ,许多方面甚于更 具有中国特色 。
[ K ] DL 考虑了几何非线性 , 因为考虑几何非线性后 , 变形导 致坐标参考系的改变 ,因此构成刚度矩阵的几何参数变化受{Δδ} 影响 ,构成非线性关系 。
1. 2 考虑几何非线性和材料非线性法[2]
考虑几何非线性和材料非线性的斜拉桥结构增量平衡方程
和计算情形与几何非线性一致 ,但是切线刚度矩阵发生变化 : [ K ] T = [ K ] DL + [ K ] ep + [ K ]σ。
33 (19) :1032104. [ 3 ] 付建新. 桩基检测研析[J ]. 山西建筑 ,2007 ,33 (20) :1082109.
The construction and detection of highway bridge pile f oundation
L I Jian2zhong Abstract : Combined wit h t he application of cast in place bored pile and manual digging pile in highway bridge pile foundation , t he key matters of bot h file foundation construction were discussed. The main technology of present pile foundation detection was introduced so as to complete construction technology of highway bridge pile foundation and enhance pile foundation construction quality. Key words : manual digging pile , cast in place bored pile , construction , detection , quality
山 西 建 筑 第 35 200
卷第 9年
1 1
期 月
SHAN
XI
ARC HI T EC
TU R
E
Vol . Jan.
35 No . 2009
1
·329 ·
文章编号 :100926825 (2009) 0120329202
ANSYS 在斜拉桥的稳定分析中的应用
2. 1 特征值屈曲分析方法
ANSYS 中特征值屈曲分析的步骤如下 : 1) 建立模型 。这一步和通常的建模方法相同 ,进入前处理器 定义单元类型 、单元实常数 、材料性质 、模型几何形状 。但是建立
的模型有如下要求 :a. 只允许线性行为 ,如果定义了非线性单元 , 则将按线性对待 ; b. 材料的弹性模量或某种形式的刚度必须定 义 ,非线性的材料性质即使被定义也将被忽略 。
4) 展开解 。展开解就是根据需要选定前几阶特征值求对应 的特征向量 。考虑到只有最小特征值才有实际意义 ,本文仅取了 第一阶特征值 。
5) 观察结果 。进入一般后处理器查看各阶特征值对应的变 形图即为对应的失稳模态 。因为特征值分析的结果本身是计算 值的上限 ,且只允许线性行为 。
2. 2 特征值屈曲分析方法
3. 4 应力波反射法完整性检测
尽管近年来国内外对于这种方法的研究未见本质性的进展 , 但在实用和普及方面国内却有较大提高 ,这些不仅表现在国产桩 基动测仪和配套用传感已达到或接近国外先进仪器方面 ,也表现 在许多单位认真研究各个测试细小环节和分析环节方面 ,更主要 的是表现在许多管理部门已开始认真总结应力波反射法完整性 检测的得与失 ,开始使这种方法的应用回归到一种正常的位置 , 如广东省正在将这种方法定位于为后续检测提供前期技术准备 , 这种定位已与该方法解决问题的真实能力完全对应 。
3. 6 动静法
由于高应变动力试桩法力的作用时间过短 ,桩只能被视为弹 性体进行分析 ,国外有人提出了一种动静法 ,采用技术将力的作 用时间延长 ,使沿桩身传播的应力波波长大于实际桩长 ,进而将 桩视为刚体 ,回避了应力波的传播问题 。应该说这种方法既克服 了传统静载试验的笨重与费时 ,也克服了高应力方法的过分间接 性 ,是一种较好的方法 ,但由于该方法对锤的配重要求太高 ,具体 操作仍有较大难度 。
收稿日期 :2008209211 作者简介 :鄢生全 (19832 ) ,男 ,助理工程师 ,广东长大公路工程有限公司第一分公司 ,广东 广州 511431
柴海峰 (19832 ) ,男 ,助理工程师 ,中交二航局二公司 ,重庆 400074 田浩亮 (19822 ) ,男 ,助理工程师 ,中交第一公路勘察设计研究院有限责任公司 ,陕西 西安 710075
1 第一类稳定和第二类稳定计算的有限元实现
在有限元计算中 ,考虑了结构的几何刚度 ,该刚度反映了结 构的不稳定因素 。当外力增加 λ倍时 ,杆力和几何刚度矩阵也增 大了 λ倍 ,因而可以写出下式 :
( [ K ] D +λ[ K ]σ) {δ} =λ{ F} 。 如果 λ足够大 ,使得结构达到随遇平衡状态 , 节点位移矩阵 由{δ} 变为{δ} + {Δδ} ,上式平衡方程也能满足 ,则 :
其中 ,[ K ] ep为弹塑性刚度矩阵 。 第二类稳定的安全系数可以采用下式计算 :
λ=
Pcritical 。 Pexisting
其中 , Pcritical为结构的临界荷载 ; Pexisting为结构的实际荷载 。
2 斜拉桥稳定计算在 ANSYS 中的实现
ANSYS 程序提供两种结构屈曲载荷和屈曲模态的分析方法 为 :特征值屈曲分析和非线性屈曲分析 。
1. 1 几何非线性法[2 ]
斜拉桥结构几何非线性增量平衡方程为 : [ K ] T{Δδ} = {Δ F} 。
[ K ] T = [ K ] DL + [ K ]σ。
其中 ,[ K ] T 为结构切线刚度矩阵 ; [ K ] DL 为考虑结构非线性 的总体弹性刚度矩阵 ;{Δ F} 为外荷载增量 。
© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
·330 ·
第 35 200
卷第 9年
1 1
期 月
山
西
建ห้องสมุดไป่ตู้
筑
( [ K ] D +λ[ K ]σ) ({δ} + {Δδ} ) =λ{ F} 。
同时满足上面两式的条件是 : ( [ K ] D +λ[ K ]σ) {Δδ} = { 0} 。
仅在使用的方便程度上有了质的飞跃 ,而且在分析手段上也有了 很大提高 ,声失时判读已不再是唯一的选择 ,声幅和声频已开始 进入了分析判断领域 ,尤其令人欣慰的是 ,声波 CT 已步入实用阶 段 ,为声波透射法的后续研究提供了广阔的前景 。
ANSYS 中几何非线性屈曲分析的步骤如下 : 1) 建立模型 。 2) 求解 。 a. 进入求解器 ,施加约束 ; b. 特征值屈曲分析 ;c. 根据特征值 分析的结果施加荷载和初始缺陷 (通常是将单倍荷载乘以第一阶 特征值作为几何非线性分析的荷载上限来施加 ,根据特征值分析 的第一阶失稳模态给结构施加初始缺陷) ;d. 设定分析选项 ,主要 包括 打 开 大 变 形 开 关 ( NL GEOM , ON) , 牛 顿 —拉 普 森 选 项 (NROP T ,ON) 等 ;e. 指定荷载步选项 ,包括指定 TIM E 值和 AU2 TO T 自动荷载分步项等 ;f . 设定收敛准则 ; g. 确定最大平衡迭代 次数和求解终止准则 ;h. 求解发散 ,在后处理中检查构件内力 ,查 看各子步计算结果中单元受力状态 ,如果第 K 子步构件实际受 力超过了该构件的极限承载力 ,则记录该构件对应的单元 ,在以 后的计算中将这些单元杀死 ,重启动至第 K - 1 子步计算结束状 态并继续加载 ,一直循环计算到不再删除单元 ,计算结束 。 3) 进入一般后处理器或时间历程后处理器查看结果 。
其中 ,[ K ] D 为结构整体弹性刚度矩阵 ; [ K ]σ 为结构整体几 何刚度矩阵 ;λ为整体安全系数 。
这就是计算稳定安全系数的特征方程式 ,如果方程有 n 阶 , 那么理论上存在 n 个特征值λ1 ,λ2 , …,λn ,工程问题中只有最低 阶特征值才有意义 。
由于第一类稳定问题是分支点失稳 ,实际工程结构都存在初 始缺陷 ,是极值失稳问题即第二类稳定 。由于斜拉桥存在几何非 线性和材料非线性 ,因此 ,第二类稳定问题必须考虑线性影响 。
2) 获得静力解 。先定义约束及荷载 ,并求解 。在这一步里必 须打开预应力开关命令 ,以记录内力并求得应力刚度矩阵 。求解 完成后退出求解器 。
3) 获得特征值屈曲解 。重新进入求解器 :定义分析类型为特 征值屈曲分析 ;定义分析选项 ,分析选项包括特征值提取方式 、需 提取的特征值数目 、特征值计算的起始点 ;定义荷载步选项 ;求解 。
非线性稳定分析的实质是求全过程 P —Δ 曲线 ,其中 ,Δ 为 关心节点的位移 ; P 为加载系数 ,取值 0～1 ,其含义是迭代荷载与 预加极限荷载的比值 。非线性稳定包括几何非线性稳定和材料 非线性稳定 。
几何非线性屈曲分析可以直接利用特征值屈曲分析的模型 , 因为非线性的材料性质将被忽略 。
