ANSYS地震响应分析讨论

地震响应分析地震响应分析 1 模态组合就是根据模态分析中的几阶振型 也可以少于这几阶 看你要求的精度 进行组合 类 似于结构最不利组合 从而求出地震响应地震响应的最大值 2 组合各振型反应的最大值 求得结构地震响应地震响应的最大值 这个问题在论坛上已经有很多人问过 也有各种各样的回答 但是至今没有令人满意的解答 我自 己试过很多种方法 加上论坛上其他人提到的方法 大致归类如下 1 先做静力恒载工况分析 打开预应力 pstres 开关 然后转到时程分析 结果 恒载对后面的时程计算不起作用 时程计算依然从 0 开始 2 直接在 antype trans 中考虑恒载 先把 timint off 加 acel 9 81 打开应力刚化 sstif on lswrite 1 然后 timint on 开始时程计算 结果 恒载 9 81 起作用了 但结果是错的 它被积分了 3 不用什么前处理 直接把 9 81 加在地震波上 acel 9 81 ac i 结果 同 2 9 81 带入了积分 这个 9 81 相当于阶跃荷载 而不是产生恒载 4 ansys 帮助中施加初始加速度的方法 篇幅限制请自己看帮助 结果 同 2 3 9 81 还是带进时间积分 5 这种是我受到别人的启发 通过结构受 ramp 荷载的特点施加的 可以近似的解决问题 即 1 求出结构的自振一阶频率 w 2 令 tr 1 w 3 定义 ramp 荷载为从 0 到 tr 加到 9 81 然后在整个时间积分中保持不变 4 antype trans 中分几个荷载步将荷载从 0 加到 9 81 5 在随后的荷载步中 acel 9 81 ac i 这种做法虽然也是将 9 81 加到地震波中 但是因为满足 TR 的要求 所以这个动力效应被削弱 到了静力效应 它作用在结构上就像静载一样 对于单自由度结构理论上跟静载是完全一样的 但 是多自由度会子静力效应上下很小的范围内波动 所以可以认为相当于静载的作用 这样我们就可 以达到考虑恒载的目的了 第 5 种是我至今为止考虑恒载的做法 我也很想知道还有没有更简单精确的方法 或者在前 4 种 方法中就有只是我使用不正确 希望大家能一起来讨论 彻底解决这个问题 谢谢 地震反应怎么考虑重力地震反应怎么考虑重力 SOLU ANTYPE TRANS TRNOPT FULL TIMINT OFF 先关闭时间积分效应 TIME 1E 8 设一个极短的积分时间 acel 9 8 NSUBST 2 有时候子步数要增大 KBC 1 LSWR 1 把这个写入第一步 TIMINT ON 然后再时间积分效应开关 以后就正常写载荷步了 这种方法应该是对的 ANSYS 帮助文件中也有提到 可是 有一个问题 由于是阶跃荷载 就会产生动力效应 整个结构的变形大于实际的情况吧 这样与实 际结构在重力下受到的变形就不一样了 我想 之所以用 KBC 1 的原因是产生阶跃荷载 第一个子步产生变形 然后后几个子步使其稳定 不产生 初始速度吧 自己作了一个例子 终于彻底搞明白啦 KBC 1 是产生阶跃荷载 不过 关键就在于此时关掉了时间积分效应 TIMINT OFF 那么 此时就相当于是静力分析了 当然就不会产生动力效应了 整个结构的变形就是实际变形 另外 KBC 1 的作用是不产生初始速度 总结下 非线性时程分析时 结构在地震波来之前 已经存在初始位移和初始内力 但是没有初始速度和 加速度 因此 一定要注意 KBC 1 同时 TIMINT OFF 设置 2 个以上的子步 很小的时间 一楼的总结有很多不妥的地方 你那个方法不知道正确不 而且太麻烦了 tmh cn 在第 12 楼总结的方法是对的 其他人说的方法都有问题 另外 不用打开应力刚化 fini config nres 20000 config nproc 2 lstep 1500 el DIM aei ARRAY lstep 1 1 CREATE ansuitmp VREAD aei 1 el txt 1 e16 6 END INPUT ansuitmp csys 0 solu antype trans trnopt full TIMINT OFF 先关闭时间积分效应 SSTIF On KBC 1 TIME 1e 3 设一个极短的积分时间 acel 8 7 NSUBST 2 有时候子步数要增大 solve timint on kbc 0 TINTP 0 25 0 5 1 0 5 1 ki 0 05 kj 0 05 oi 2 3 14159 0 09 oj 2 3 14159 0 16 ma 2 ki oj kj oi oi oj oj oj oi oi mb 2 kj oj ki oi oj oj oi oi ALPHAD ma BETAD mb allsel all do tm 1 1500 1 time tm 2 100 NSUBST 1 1 AUTOTS 1 acel aei tm 314 7 EI 波 只输入 X 向加速度 solve enddo fini 关于动力学分析模拟自重对结构的影响 这个问题 现在偶已经解决这个问题啦 几行命令 流就可以了 PREP7 ACEL 0 9800 0 先把重力加速度在预处理设好 以后就不在设了 FINISH SOLU ANTYPE TRANS TRNOPT FULL TIMINT OFF 先关闭时间积分效应 TIME 1E 8 设一个极短的积分时间 NSUBST 2 KBC 1 LSWR 1 把这个写入第一步 TIMINT ON 然后再时间积分效应开关 以后就正常写载荷步了 虽然 这只是个很简单的问题 但如果不注意 会使计算的动力响应值受重力影响而偏大 我当时 就很郁闷 怎么也不能把重力的影响去掉 现在终于搞定了 现与大家分享偶的一点心得 同时 也感谢各位大侠的关注 模态叠加法 时间可能取的太短了 求地震反应模态叠加法 时间可能取的太短了 求地震反应 步骤步骤 1 进行模态分析 进行模态分析 SOLU 1 MODE ANALYSIS AnType Modal ModOpt Lanb 20 EXTRACT 20 MODES Mxpand 20 1 SOLVE FINISH 时程分析命令流时程分析命令流 采用模态叠加法采用模态叠加法 set st 0 01 set endpoint 100 dim xforce array endpoint vread xforce 1 Seismic x txt 7F10 5 VOPER xforce 1 xforce 1 MULT 9 8 SOLU ANTYPE 4 TRNOPT MSUP 20 1 施加位移约束施加位移约束 NSEL S LOC Z 0 0 选取底层所有节点选取底层所有节点 D ALL ALL 约束所有自由度约束所有自由度 ALLSEL ALL dim shijian array endpoint vfill shijian 1 ramp st st Nt star 1 Nt end 100 Nt incr 1 DO Nt Nt star Nt end Nt incr TIME shijian Nt ACEL xforce Nt NSEL ALL LSWRITE Nt SOLVE ENDDO FINISH 模态扩展模态扩展 SOLU EXPASS ON NUMEXP all 0 1 OUTRES all all SOLVE FINISH 谱分析谱分析 我想 作谱分析之前应该先弄明白谱分析的工程意义 实际上 谱分析是一种线性化的动力分析思 路 简单点说 就是在获得结构自振特性以后 用一条符合本地场地土特性和设防烈度 地震响应地震响应 谱分析 的地震谱或者一条符合本地风场特性的风速谱 作为外部作用施加与结构上 以求得结构 在各个频率对应的谱值作用下的结构响应 他实际也是一种时程分析 只不过时间轴为频率值 所 以也称为频域分析 由于反应普的获得时基于结构的线性化假定 所以谱分析只能对线性结构进行 分析 先试贴 如有人需要给留言发送给你 Q 为问 A 为回复 Q 模态分析得到的结果是不是某个方向上的各阶频率啊 我要得到各个方向的一阶频率能做到吗 A 模态分析得到的结果是你所选择的自由度内的振型 如果需要得到某一个方向内的振型 用 RUDECED 法 选择你 所希望的主自由度 如 UX 将所求的频率值设置为一即可 Q 不过大型复杂结构的振型一般都不是一个方向的 应该是一个方向为主 其它方向为次 如果把所要求的方向设 为主自由度 会不会跟真是实际结果有出入 Q 对于复杂的大型实体 其主自由度无法给出 此时就无法应用 reduce 法 可不可以改约束 还是就认为一次得到 前 3 阶频率就是 x y z 方向的一阶频率啊 A 正是因为复杂结构的主自由度不好确定 所以我还是倾向于用 SUBSPACE 和 BLOCK LANCZOS 法 约束不要随便 改 它应该描述结构真实的状态 不能说前 3 阶频率即为 X Y Z 方向的第一阶频率 你可以好好想想所用方程中 K 和 M 的意义 Q 轴承是用来支撑要研究的对象 把轴承等价为弹簧 将其刚度加在对象的节点上 这样就要定义刚度的大小在加 是否还要定义单元类型 然后才能定义其刚度 A 我曾经用 LINK8 单元模拟过 根据 K EA L 确定刚度 我当时取 E 210e9 L 0 1 然后 A KL E 经过检验 是可 行的 A 1 用 link 单元模拟其刚度 k EA L 通过假定 E A L 的值使之和 k 相等 其密度要 0 2 用 combin14 模拟 其刚度 阻尼 0 Q ANSYS 处理动力问题 如土层地震响应地震响应 或基础振动反应 可否从边界施加加速度荷载 若可以 如何 瞬态 动力反应分析 时间子步为200 A 如果你所希望求的最高频率的周期为 f 则 ANSYS 取时间步长为 1 20f 为默认值 Q 板单元的后处理的显示在单元坐标还是在总体坐标 比如位移 应力 A 应力位移在总体坐标 弯矩不一定 A 应力 位移是显示在结果坐标系的 顺便说明一下 节点信息显示在显示坐标系 不过结果坐标系默认为全局坐标系 而显示坐标系也是默认为全局坐标系 你可以这样修改结果坐标系 Main Menu General PostprocessorOptions for Output 你可以这样修改显示坐标系 Utility Menu WorkplaneChange Display CS to A 你可以用 ETABLE 命令 先存储对应于各单元的位移 等等 再用 ANSYS 的柱状图 或是别的画图软件 PLOT 出 图形 你认为如何 Q 我现在需分析一个板梁结构 板已用 SHELL63 单元划分好 梁我是用板上的一条线划分单元并添加截面而生成的 但现在运算时 发现板和梁是分开的 它们之间互不影响 请教各位高手 怎样将板和梁合并为一个整体 A PreprocessorNumbering CtrlsMerge Items 里 element and node 合并 Q 用 SHELL63 单元划分的板结构 想加垂直于板面的匀布载荷 A 找到如下位置 ANSYS HELP ANSYS ELEMENTS REFERENCE ELEMENT LIBRARY SHELL63 ELASTIC SHELL 在 INPUT SUMMARY 中 有 SURFACE LOADS PRESSURES 其中 face 1 I J K L bottom in Z direction 的意思为 当 Load Key 1 时 pressure 加在由节点 IJKL 确定的面上 在单元坐标中以 Z 方向为正 也可以说垂直指向面 IJKL 时为正 至于均布 则将几个值都填成一样的就可以了 或者只填主要参数如 value at I 自己多看看帮助文件 做一个例子验 证一下 Q ANSYS 可以自己定义输出格式吗 A 可以 先打开一个文件 cfopen 再用 vwrite label 通过打开一个文件 读出数据 或以追加的方式写入数据 并可控制输入和输出的格式 cfopen vread f12 6 vwrite f12 6 cfclos Q 作模态分析时 就是集中质量单元和刚度加法 A 在选用求解方法时 例如 SUBSPACE 和 BLOCK LANCZOS 会有一个选项 LUMPM 默认为选用一致质量矩阵 你可以选为集中质量矩阵 A 用 mass21 通过 mesh keypoint 来把点划分成 mass21 单元即可 Q 我在原有的静态分析的基础上 想继续作一下结构的模态分析 应该采取那些步骤啊 我用 solutionnew analysis modal then current LS 出现如下的提示 1 the eigenvalue extraction method must be specified on MODOPT command SUBSP or LANB are usually recommanded 2 please expand modes if any downstream MODAL analysis based calculations need to be performed 这是怎么一回事 A 你最好先看看帮助文件中的 ANSYS STRUCTURE REFERENCE 中的 BUCKLING ANALYSIS 如果你对理论不清楚 我劝你不要盲目计算 不然结果也许并不可信 这是我的建议 A 系统默认为 0 阶 Q 1 在 ANSYS 中如何设定小数的有效位数 比如我输入的应力值为 0 而且我的一组数据都是最后几位差别 但是 ANSYSY 自做聪明把它四舍五入了 结果弄的都一样了 如何更改这个有效位数呢 2 ANSYS 在输出时如何修改输出的数字显示模式 比如用科学记数法 或者将数字缩放多少倍那样显示 在什么地方 能够修改呢 3 ANSYS5 7 1 中什么地方能够出现那种材料库表 就是说不同 ANSYS 的什么 material models 对话框生成 而是通过 table 输入材料属性诸如应力应变曲线之类的东西 A 用 a bEn 来表示 如 2 45678e 9 材料特性用 tb Q 我是个 ANSYS 初学者 现在碰到下面几个问题 请各位高手帮忙看一下 1 是什么原因使得相邻 ELEMENT 的共用节点解得的力不一样 应该是方法步骤的原因而不是误差吧 请简单介绍一 下 2 假如我想用 ANSYS 模拟陶瓷的烧结过程 是不是有实际 作的可能性 如果有 其基本的原理和模型怎样 3 ANSYS 可以用来模拟热循环下的应力分配不均造成的焊点 63 Snhard copyto file 然后在剪贴板中就可以看见了 A 我都是用抓图工具抓下来 然后再打印的 A ui copy save 是一个截取屏幕图形的命令 现在总算又实现了 gui 肯定可以通过相应的命令方式来实现的观点 A Ansys 5 6 已经做到了把图形窗口 COPY 到文件如 JPG BMP 等格式 PLOTCTRL HARD COPY TO FILE 很好用 Q 其实 APDL 语言是 ANSYS 的二次开发语言 是在 FORTRAN 编程语言的基础上形成的 你只要会 FORTRAN 语言 这就很容易理解 APDL 语言了 也有循环 选择等编程结构 把枯燥的 作用计算机自动实现 然后形成一个宏 作拉 Q 是的 APDL 语言确实很方便 而且和 FORTRAN 语言类似 很容易读懂 在计算电机电磁场中电机的反电势更加 方便 前几天刚算了一个需要循环运算 288 次的程序 使用 APDL 语言使的我可以让电脑自己运行程序 保存结果数 据 而我的任务就是睡一觉醒来 对结果数据进行处理 Q 我看见好多人在讨论用命令流处理程序 很不明白命令流和 gui 的形式有很大区别吗 在什么情况下更适用命令呢 A 很容易控制 Ansys 且结果文件小 功能很强大 可以进行循环控制 分支等 作 运算速度快 充分利用 APDL 语言 对 ansys 进一步开发 Q 用 ansys5 6 2 计算一个压力容器 划分网格 9 万余个 运行环境 win2000 cpu 1G 内存 640M 硬盘 40G 模型建立 网 格划分均无问题 但是运行到占用硬盘大约 9 2G 的时候 就出错 提示内存无法 read 请问该如何解决 是内存的问题 还 是 ansys 本身问题 还是系统问题 硬件问题 模型建立实在无法简化 A 我的硬盘是 30G 的 内存 1G OS windowsNT Total Workingspace 是 1400M 再高就控制不了 Database 是 56 我能解的有效节点数是 只是因为磁盘空间不够 工作的分区是 0G 当 ansys 探测到磁盘空间已满 就自动退出 尽量大的增大节点数 并不是好办法 因为运行速度太慢了 利用对称性减小求解区域 利用 Lesize esize 等控制网 格 都是好的办法 A 选用对称性的时候要注意边界条件 但的确是个好办法 A 以前在使用 ALGOR 时也遇到过类似情况 当时有 10 万个结点 对 WIN2000 来说 已不存在 WIN9X 中文件最大 只能 2G 的瓶颈 但在处理超大虚拟内存需求时仍然会出问题 处理方法有二 1 使用 PQMAGIC 划分出一个 20G 的逻辑分区 一定要跟 WIN2000 系统在不同的分区 再在控制面板 系统中将 虚拟内存文件指定到这个大分区上 以尽量避免发生系统堆栈错误 虚拟文件大小仍交 WINDOWS 控制 2 将至少两个分区的硬盘升级至 NTFS 由于 NTFS 文件可以跨分区存储 且要小 10 左右 比 FAT32 要好用得多 此帖已经发过多遍 请注意发帖规则此帖已经发过多遍 请注意发帖规则 日月星辰日月星辰 积分积分 28 帖子帖子 38 2 2004 6 17 20 42 Q 我在求解是出现如下的提示 1 real constant 1 referenced by at least element types 1 and 2 2 coefficient ratio exceeds 1 0e8 我全部忽略上述提示 得出了结果 结果能可信吗 我有错在哪里呢 模型我都检查了数遍了 还是和初始文件名有 关 我是在一个初始名下面建立若干相关的模型 由于我的单元只是板梁的结合 所以 我重新检查了一下 删除了一个影响不大的扁梁 很扁 尺寸不大 系统可 能认为是板了 运行就通过了 上述的提示可能是因为系统对某些单元有默认的尺寸限制 A 有些单元类型在默认值状态下不需要输入 REAL CONSTANT 但是当你改变了某项选项 你就需要重新输入 REAL CONSTANT 你可以用 UTILITYSELECTENTITIES 命令中的 SELECT ELEMENT BY ATTRIBUTE 来查看单元材 料 REAL CONSTANT 等等是否定义明确 然后用 APPLY REPLOT 命令来显示每个不同材料的单元 A 我试了一下 按照你所说的步骤 出现了如下提示 a default command range is not provided for the ESEL command the ESEL command is ignored 所以 没有任何结果可见 是不是需要对某些 default 的命令做些改动呢 A real constant 1 referenced by at least element types 1 and 2 此句意思好象是说 你为 1 号单元和 2 号单元或是 别的单元定义的 real constant 均为 1 你需要检查你的模型的单元类型和实常数部分 A 第二个问题经常出现 第一个问题有点问题 你看看单元定义对没有如果单元都没有划分对 结果当然不可靠啦 Q slab on girder highway bridge 桥板为 SOLID 单元 由于 girder 上部 flange 嵌贴在板的下表面 所以建模时 我 将桥板定义为 SOLID 单元 桥板下表面 AREA 定义为 SHELL63 单元 但 SOLID 和 SHELL 单元的自由度不同 这样建 模肯定有不对的地方 运算结果与实验结果也相差了 40 有什么好的解决办法 A 你可以根据需要为不同的单元增加自由度 DOF 顺便问一句 你的模型的两个部分的节点是否完全重合 提高 计算精度要用更细的网格划分 另外 要保证模型中每个单元的长宽比和理 是否有使用大应变的需要 在 ANSYS 非 线形手册中有影响精度的说明 A 本来 SOLID 体的下表面被我定义为 SHELL 这时面单元的节点同时也是体单元下表面的节点 但由于两种不同的 单元自由度不同 所以我放弃了这个模型 另建模型的 SHELL 单元和 SOLID 单元相贴面节点坐标完全相同 节点号不 同 这样才可能把它们耦合 但是结果仍不理想 可有什么方法 A 比较麻烦 可以通过编程输入命令流 QA 我已经把 SHELL 和 SOLID 对应点的 UX UY UZ 全部耦合起来 让这些 NODE 的 X Y Z 三个方向的位移相同 但是结果仍然不理想 我想约束方程的功能也进乎如此 有没有更进一步的建议 十分感谢 Q ansys 有没有评价一个网格划分好坏的参数或命令 请告知菜单 作 A PREPROCESSORMESHINGCHECK MESH两个命令中一个用来 CHECK 网格划分 另一个用来 CHECK 连 贯性 只是 CHECK 完了别哭就好 A PLOT BAD ELEMENT Q 1 在 gui 下运行宏时 重新定义原来已有的数组 ansys 就会有个提示 想不要这个提示 怎么办 uis msgpop 3 这个行不通 我试过了 2 ansys 中数组的定义为 a b c a10040 b255 c7 而我实际上要定义更大的有什么好办法 难道非要写到文 件里再读进来 A q1 有两种办法 1 在 ultility menu menuctrls message controls 下设 2 将 uis msgpop 3 加入到 start ans 文 件中去 以后每一次运行 ansys 时都会默认这个设置 Q 我现在做一个板壳的屈曲分析 想要得到该结构在荷载作用下的极限承载力 利用非线性屈曲分析 请问极限承载 力的荷载比例因子在什么命令下可以查到 A 非线性屈曲分析与特征屈曲的不同之处是 它一开始就挠屈变形 随着荷载的增大 变形日益增大 我们认为结构 最危险点的变形 时间历程曲线中变形若从某点开始发生较大的转折 与先前完全不一样 就发生了屈曲 此时的荷 载就是屈曲荷载 故不存在荷载比例因子 即特征值问题 A 要进行非线性屈曲 必须加到可以使结构发生屈曲的荷载 即荷载要比结构的实际承载力要大 在进行非线性屈曲分析 极限承载力分析 时荷载逐步增加到破坏为止 计算时一定要输出每一步的计算结果 这样在 post1 中就会有一系列的荷 载系数 而最后收敛的那个系数就是相对于极限承载力的系数 极限承载力 系数和总荷载的乘积 Q 如果材料一定为各向异性 我用的材料为钢板 该如何选用单元类型 A 只要材料参数中有 Ex Ey Ez 单独输入项 就能考虑各向异性 A 波纹板的 EX EY 可以输入 但是 EZ 该如何确定呢 我自己练习了一下 如果不输入 EZ 则计算就不能进行下去 提示我材料的弹性模量要输入正值 这是为什么 A 如果是钢板 一般认为是各向同性的 厚板因轧制次数少 晶体不密实存在着 Z 向 厚度方向 弱化 对厚板 一 般地 Ex Ey Ez 较小 通过试验确定 对薄板16mm 认为各向同性 A 不管是什么板 材料的本构关系是微小单元之间应力与应变的关系 因此波纹板的几何构成是用几何模型 有 限元模型计算的 而不是用比较大块的板试验得出的 Ex Ey 来输入 你说的波纹板是不是压型钢板 A 是压型钢板 EX EY 是通过拉伸试验确定的 EX 与 EY 不相等 但是既然是考虑板壳 为什么还要输入 EZ 况且 EZ 并没有试验数据 A 波纹板的 Ex Ey Ez Q 怎样使梁单元把弯矩传给实体单元 实体单元只有应力 使它与相连的梁单元保持连续性真是件很难的事情 还请 各位大侠指教 A 关键不是实体单元只有应力 而是实体单元每个点只有三个平动自由度 能不能模拟成将梁插入到实体中一小段 用两个或多个实体单元节点固定梁单元 根据 Sevant 原理 忽略掉两种不同类型单元连接处的计算结果 A 关键问题在于传递弯矩是为了保持变形的连续性 实体单元中尽管没有弯矩 但同一截面各点处的应力实际已构成 了弯矩 但这弯矩是多大 怎样传递给相邻的梁 恐怕不是简单的插入嵌固可以解决的 从各结构杂志来看 这一问 题似乎很普遍 还请各位一起深入探讨 A 上述方法实质上就是施加了人工约束 单元相邻近点的连线几乎是梁的切线 这样就能保证转角和位移连续 弯矩 相等并不意味着位移和一阶导数相等 这是因为两种单元的位移函数相差太大 Q 弹性地基梁中的弹簧 2 维 在 ANSYS 中应采用那一个单元 又如何 作 A 很久以前我用过弹簧单元 好象是 COMBIN 2D 你试试看 Q 我在平面应变的有限元分析中 要使用初始应力这个条件 该初始应力在没加重力的时候在各点的应力是相同的 于是 我在 WIN2000 中相当于 DOS 的平台上用 EDIT 生成了初始应力文件 在第一个载荷步中用 constsnt stress 这个选项读 入初始应力并逐步调整边界条件 使计算的结果与初始应力相差不大 然后我在第二个载荷步中施加重力 进一步形成迭 加应力场 在完成了上述步骤之后在第三载荷步中进行开挖 然而 在这一步中无论挖掉多少 要么结果十分不可靠或不收 敛 请众同行不惜赐教 指出错误之处 本人对您们的帮助万分感激 模型大概如下 模型边界为四方形 其边界取足够大 中间有圆孔需要在一定的初始应力条件下挖掉 这个初始应力条件为正应为 SX 12 7MPa SY 14 7MPa SZ 10 2MPa 剪应力 SXY 1 7MPa SYZ 2 1MPa SZX 3 0MPa 在下边界及左右边界法向约 束 上部边界施加稍大于 SZ 值的力 A 可能实施加重力时 单位不一致 Q ANSYS 里椭圆有没有命令可以直接画 A 你可以利用 ansys 的局部坐标系 非常方便 gui workplanelocal coordinate systemscreat local csat wp origin 选择 cylindrical 1 再用 part1 part2 设置长短轴 这种方法在工作平面原点产生局部椭圆坐标系 在局部椭圆坐标系下产生 2 keypoints 连点成线即可 Q 在椭圆坐标系下生成怎样两个点 A 每一种坐标系的坐标都用 x y z 来表示 但是当你设为不同的坐标系时 x y z 表示不同的含义 例如 在圆柱坐标下 x y z 非别代表 R z 输入相应的值即可 更简单 现在直角坐标系下生成两个点 然后在生成椭圆坐标系 可以参阅 ansys 基本过程手册关于坐标系的讲解 Q 一个圆锥形曲面刚架屋顶 怎样加风荷载 A 在 Ansys 中 施加面荷载时 可以设置倾斜度 首先用 SFGRAD 命令设置倾斜度和 0 点位置 然后选择面或节点施 加压力 0 点位置的压力等于施加的压力 其他点的压力值等于 由距 0 点的坐标差 笛卡尔坐标系 X 轴或 Y 轴或 Z 轴 球面坐标系 R 或 乘于倾斜度再加上施加的压力值 如水荷载 倾斜度等于水的比重 Q ansys 默认面荷载是垂直于作用面 而曲面承受的风荷载为恒定的方向 那怎么办 日月星辰日月星辰 积分积分 28 帖子帖子 38 3 2004 6 17 20 43 Q 有一个形状很不规则的刚架结构 外表面是玻璃 建模的时候只建了刚架 风荷载怎样传递到刚架上去啊 A 建模时就应该把玻璃和刚架同时建好 要注意两种材料的衔接面 或 LINE 和材料特性定义 风载加上去后就会由玻 璃传递给刚架 或者简化风载成为 PRESSURE 直接加到刚架上 如果手上有实验数据加以比较 会避免很多失误 Q 管路穿过厚壁 管路采用梁单元 厚壁采用体单元 在两者交接的地方怎么办 A 建模时可以考虑先建管路 再建厚壁体时厚壁体由外部的平面与管路外表曲面组成 这样管路和厚壁 SHARE 同一 个曲面 荷载作用下力也可以传递 材料属性分配的时候 先 PICK UP 厚壁单元 定义为所选取得体单元 再 PICK UP 管路 定义为梁单元 其实很多梁单元也属于体单元 没有什么冲突 除非管路属于 SHELL 或面单元时 要后 PICK UP 此面单元 定义为所要的属性 至为重要的是在 MESH 这两个单元面体衔接的时候 要注意衔接面体对应的 ELEMENTS 和 NODES 的 COMPATIBILITY 否则将会出现大的错误 REFINE 时也不可忽视 QA 1 请详细解释一下 1 再建厚壁体时厚壁体由外部的平面与管路外表曲面组成 这样管路和厚壁 SHARE 同一个曲面 2 至为重要的是在 MESH 这两个单元面体衔接的时候 要注意衔接面体对应的 ELEMENTS 和 NODES 的 COMPATIBILITY 否则将会出现大的错误 REFINE 时也不可忽视 并非英文不懂 2 我的情况是 已经建好了厚壁体 再建立管路 3 这样一个问题 该如何办 例如 一块厚板 x1 x2 y1 y2 z1 z2 0 1 0 1 0 0 1 一根管由 k1 x1 y1 z1 0 5 0 5 2 k2 x2 y2 z2 0 5 0 5 2 构成 求其固有频率 其他条件随便 管路为圆管 Q 模型是一个正交异性钢桥面板是个钢箱梁取了一部分来分析 由主梁 纵向加劲肋 横隔板 钢板还有沥青铺装层 组成横向 4800mm 纵向 9060mm 纵向加劲肋共 8 个 横隔板是 3m 一个 共四个钢板层厚度 14mm 沥青混凝土 50mm 除了上面的沥青铺装层其他都是钢的我原来考虑可以使用 shell63 单元但我要计算铺装层与钢板层的层间剪应力 似乎不好处理 所以全部采用 solid45 建模 在中间某一荷位作用一矩形均布荷载我现在建模完成了 在网格划分的时候如果按照 smartsize 的话可以划分 但是我 怎么施加荷载呢 所以好像必须在荷载作用区域有整数个 element 如果按照 smartsize 的话可以划分 但是我怎么施加 荷载呢 这样可以 press on element 但是我不会自己控制网格划分在中间某一荷位作用一矩形均布荷载建模型的时候就要预先 把那个部分设一个面 似乎很复杂这个面下面又有钢板又有加劲肋 另外我即使简化成一个集中荷载来计算 那也得中间有 node 呀 我怎么能保证我的网格在我需要的坐标值上恰好有个 node 我的模型并不是规则的啊 下面是梯形的加劲肋 如果按 map 的话它会提示非法拓扑形式即使我等效成集中力我在想如果定义 hard point 的话 可以加载到那个位置上但我好像没看到可以按坐标值定义 hard point A 使用 SWEEP 或 CONCATENATE 命令重新 MESH 桥板 可以在 HELP 文件中找到这两种命令的用法 十分有用 将桥面 MESH 成为 HEX 类的 ELEMENT 这样寻找加载位置会 方便很多 A 建议你将整个桥面划分为比较规则的四边形网格 六面体 当然网格的大小肯定与 加强筋的厚度等有关 这样应 该能够解决你所说的问题 Q 一个梁上有一滑块 计重力和摩擦 滑块在梁上滑动 求梁的响应 梁为 1m 0 1 0 1 滑块为 0 2 0 1 0 1 ex 2e11 nuxy 0 3 dens 1 7 8e3 dens 2 2 7e3 A 1 命令 CONCATENATE 如果你要把一个面或者体 MESH 成为一个四方或 HEX 单元 这个命令十分有用 你可以 CONCATENATE 相邻的线或者面 把所要 MESH 的面变成一个四边形或体变成一个六面体 例如 5 边形 CONCATENATE 两条相邻的线 于是 ANSYS 认为此时的面是四边形 你就可以进行 3OR4SIDES 的 MESH 所以无 论你的面由多少边 都可以进行 MAPPED MESH 2 SWEEP 和上边的命令一样 都是我所喜欢的命令 只要你的体已经建成 SWEEP 好一个面以后 整个体 MESH 可以一扫而过 HEX 的 ELEMENT 加载在什么地方或什么样的荷载都十分方便 不需要在建模的时候就考虑加载位置 3 至于选用单元 具体问题具体分析 但是衔接一定要注意 特别是衔接面的 ELEMENT 要对应好 以免出现错误 4 自己的一点心得 MESH 十分重要 别在这上面偷懒 我曾经用 FREE MESH 与 MAPPED MESH 得到的值分别为 19 2 和 25 9 而理论值是 26 3 所以 请小心 Q 我做的是一个三层的压力容器每层厚都是 2mm 请问我可不可以只定义一种单元类型 没有 shell91 99 另外这 三种材料的材料属性用什么命令分别定义给这三层 另外用轴对称分析的话应该怎么做 A 划分网格前 用 attributedefine 可定义不同的材料属性 Q 怎样进行螺纹的失效分析 我现在要做一个 Failure analysis of a thread joint 对于螺纹失效分析问题 我不知道 Ansys 里是如何处理的 在 creat 3D model 时 是否要建立完全体现 thread 的几何 model 然后直接用普通的 solid element 进行 mesh 再按照通常的分析过程分析 还是在建模时将 joint 简化为一普通轴 而采用特殊的 element 进行 mesh 然后再分析 对于前一方法 用 Ansys creat thread geometry 比较困难 至少我不知道怎么做 而且 mesh 时也可能会出错 thread geometry 比较复杂 类似 sharp corner or small area 对于后一方法 在 model 的时候比较简单 但不知道用什么样的 element 才能体现 tread 和 crew 的应力集中特性 Q 有几个很简单的问题 1 在 Ansys 中 如何在一个 2 节点的单元内 同时加多个集中荷载 我用了一个很土的办法 加了几个 LoadStep 2 如何使 Beam3 单元的内力以局部坐标输出 我用 List Results ElemSolution Nodal Force 后 结果仍和整 体坐标系下一样 3 我用一个静定结构进行分析 结果输出后和实际的受力方向刚好相反 受拉的成了受压 但反力的方向是对的 Q 有一个圆柱桶 内装有 70 的水 如何将桶和水的重量考虑进去 A loadgravity A 把水另设为体单元 Q 四边固定的板 用体单元 有一条线 用梁单元 从板的任意一点穿过 分析其固有频率 如果不将线 体联为一体 分析时各算各的 问 如何将线 体联为一体 glue 无法将线 体连为一体 glue 只能将同样的元素连为一体 如体和体 面和面 线和线 Q 如何给模型加 3g 的冲击载荷 A In short time for example 0 0001 second 用一小 time 步 加载 A 3g 指的是加速度冲击载荷 例如在地震时的地基响应 例如 在固定端上用弹簧悬挂一根杆 在固定端加一个 3g 的冲击载荷 g 指加速度 求杆的响应 1 之所以用弹簧悬挂是为了减振 这个问题主要想知道如何加这种载荷 2 类似于地震时的地基响应问题 3 我在 ANSYS 的 GUI 中没有找到这个功能 4 BTW 在 ANSYS 中如何加弹簧元 Q 我想知道用 ANSYS 分析岩土洞室开挖效果好不好 我刚开始学 ANSYS 不知道 ANSYS 用于分析有裂隙岩体和松散土体方面效果好不好 其精度如何 望各位朋友能给予 指教 A dyna 做隧道开挖比较擅长 而且看起来非常形象 A 用 ansys 里的 ls dyna dyna 是绑在 ansys 上的 你选用模块就可以 Q 请问这个错误是什么 划分网格时 出现 Shape testing revealed that 16 of the 30 new or modified elements violate shape warning limits To review test results please see the output file or issue the CHECK command 这个怎么改 A 检查单元 prep7meshingcheck mesh 若有不合格的单元 用别的办法划单元 还可以设 limits 有时候其 实这种警告对结果没有多大的影响 完全根据你的经验 A 如果可以进行下去可忽略该 warning 一般来说 温度场分析对形状测试更严格一些 一般不会对问题产生太大影响 可通过剖分更密的网格来解决 Q 如何得到某一温度值的距离 时间关系曲线 在分析平板移动热源温度场时 希望取得时间历程中 熔化温度值在各时刻对应的平板上出现的位置 或节点号 即某一温度值的距离 时间关系曲线 或者如何得到某一温度值在给定时刻出现的位置 A 我想你可以把你需要的那些温度值的点选出来 设置成路径 在用路径来画图 在 path operation 下 A 谢谢你提供的建议 我不知道这些温度值点出现在什么位置 因此路径方法不可取 也许只能采用 set 设置时刻 sort 排序各节点温度 prnsol 找出温度值对应的节点 此法较为繁琐 Q 5 4 版如何用 extrude 将 2D 单元拖拉出 3D 单元的 如何进行旋转体的结构分析 加 离心力 吗 怎么加 A 先用 2D 单元剖分面 再画出拖拉路径 直线 修改默认的设置 3D 单元类型 材料 实常数 extrude 即可 Q 一个后处理问题 我用 shell93 作的壳体 节点应力输出 PRINT S NODAL SOLUTION PER NODE POST1 NODAL STRESS LISTING PowerGraphics Is Currently Enabled LOAD STEP 1 SUBSTEP 1 TIME 1 0000 LOAD CASE 0 SHELL NODAL RESULTS ARE AT TOP BOTTOM FOR MATERIAL 1 THE FOLLOWING X Y Z valueS ARE IN ELEMENT COORDINATES NODE SX SY SZ SXY SYZ SXZ 225 0 51606E 06 0 49482E 06 0 0000 3464 6 10039 38600 225 0 50183E 06 0 53119E 06 0 48105E 06 7004 0 10031 38565 3020 61720 0 38298E 06 0 0000 3034 2 9606 2 2571 7 3020 0 21884E 06 10255 0 0000 81566 1491 2 338