! *get,ae,elem,i,aproj,y !此命令用单元Y投影面积 ! *get,ae,elem,i,aproj,z !此命令用单元Z投影面积 xe=centrx !单元中心X坐标(用于求解压力值) ye=centry !单元中心Y坐标(用于求解压力值) ze=centrz !单元中心Z坐标(用于求解压力值)
! 下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化. p_e=(xe-10)*press+(ye-5)*press f_tot=p_e*ae esel,s,elem,,i 14
nsle,s,corner *get,nn,node,,count f_n=f_tot/nn *do,j,1,nn
f,nelem(i,j),fx,f_n !压力的作用方向为X方向 ! f,nelem(i,j),fy,f_n !压力的作用方向为Y方向 ! f,nelem(i,j),fz,f_n !压力的作用方向为Z方向 *enddo *endif esla,s *enddo aclear,all
fcum,repl !!!将力的施加方式还原为缺省的\替代\dofsel,all allsel
-------------------------------------------------------------------------------- QA*DOWHILE,PARA
如果PARA为真,则进行循环*DOWHILE,1.0将无限循环下去 这条命令对使用APDL的朋友会有用的
-------------------------------------------------------------------------------- Q在ansys中如何施加与某一坐标轴成一角度的集中载荷 A可以通过旋转结点坐标系来达到目的
Ayou may decompose the multiple-force to the coordinate axis. -------------------------------------------------------------------------------- Q模态分析中输出的基本频率是Hz还是圆频率 A是圆频率omiga,又简称频率单位是rad/s;
--------------------------------------------------------------------------------
Q在ANSYS/DYNA中已经进行输出子步控制但所有的数据却只写入最后一子步不知为何 A可曾使用outres,all,all? 你应该是没有输出每一子步的结果使用上述命令看看 -------------------------------------------------------------------------------- Q请问索的弹性模量在ansys中如何处理 A非线性弹性
A考虑全桥用ernst公式修正即可 Q: 如何在ansys中实现?
--------------------------------------------------------------------------------
Q既然混凝土这方面问题多多是否请各位大虾建立一小模型然后将全命令流贴上来大家一起讨论和分析岂不快哉
A一受均布荷载的简支梁跨度为3000mm,截面为100x200mm配有两根钢筋承受的均布载荷为0.04N/mm*mm命令中采用了1/4模型材料参数详见命令中由于选择时采用的是实体号而不是坐标可能在有些系统上会受到影响各位注意 !钢筋混凝土简支梁分析 /COM, Structural
!----------定义单元及材料等--------------------- /PREP7 15
et,1,link8 !定义link8单元 et,2,solid65 !定义solid65单元 keyopt,2,7,1
r,1,314 !定义link8单元的面积 r,2 !定义solid65的实常数号
mp,ex,1,2e5 !定义link8单元的弹性模量 mp,prxy,1,0.3 !定义link8单元的泊松系数 mp,ex,2,4e4 !定义solid65单元的弹性模量 mp,prxy,2,0.3 !定义solid65单元的泊松系数 tb,concr,2 !定义2号为混凝土
tbdata,,0.9,1,1.8,50 !定义混凝土的c1,c2,Rl,Ra !----------建立几何模型------------------------- blc4, , ,50,200,1500 !定义梁体 /view,1,1,1,1 !定义ISO查看 /ang,1
vplot !绘制梁体
kwpave,5 !工作平面移动到关键点5 wpoff,25 !工作平面移动25mm wprot,0,0,90 !工作平面旋转 vsbw,1 !分割梁体
wpoff,0,40 !工作平面移动40mm wprot,0,90 !工作平面旋转 vsbw,all !分割梁体
wpoff,300 !再将梁体分割为三个区域 wprot,0,0,90 !(原为控制单元密度而设) vsbw,all wpoff,0,0,900 vsbw,all
wpstyl !关闭工作平面显示 nummrg,all !整理编号 numcmp,all
!------------划分单元网格--------------------------- lsel,s,,,41,57,16 !定义line41,line57为新的选择集
lsel,a,,,9 !定义line9也在选择集中 latt,1,1,1 !定义选择集的属性
lesize,9,,,20,,,,,1 !将line9划分为20个单元 lesize,41,,,20,,,,,1 !将line41划分为20个单元 lesize,57,,,60,,,,,1 !将line57划分为60个单元 lmesh,all !对线划分单元 !-------------
allsel,all !新的选择集为所有的实体 vsel,s,,,all !定义所有体为选择集 vatt,2,2,2 !定义选择集的属性
lesize,22,,,2,,,,,0 !定义体线的单元数以控制体单元 16 lesize,23,,,2,,,,,0 lesize,33,,,2,,,,,0 lesize,28,,,4,,,,,0
mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器 vmesh,all !对体进行划分单元 finish
!--------------定义约束及荷载------------------------- /solu
asel,s,,,2,22,20 asel,a,,,14,17,3
da,all,uy !对支撑端面施加Uy方向约束 asel,s,,,3,23,20 asel,a,,,15,19,4
da,all,uz !对正对称面施加Uz方向约束 asel,s,,,32,36,4 asel,a,,,50,52,2 asel,a,,,41,43,2
da,all,ux !对侧对称面施加Ux方向约束 asel,s,,,21,27,6 asel,a,,,42,46,4 asel,a,,,35,51,16
sfa,all,1,pres,0.04 !对顶面施加面荷载(N/mm*mm) allsel,all
!--------------定义非线性分析选项-------------------- cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE nsubst,25 !定义子步数
outres,all,all !输出每一子步的结果 autots,1 !打开自动时间步控制 lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出 neqit,50 !每一子步中方程的迭代次数限值 pred,on !打开预测器
time,50 !定义时间 solve !求解
!--------------进入后处理过程----------------------- finish /post1
pldisp,0 !绘制变形图(缺省为最后一步的) etable,sig,ls,1 !定义应力表 plls,sigi,sigi,1 !绘制上述应力 plnsol,s,z,0,1 !绘制Z方向的应力 /device,vector,on
A大家都提到了收敛困难的问题为加速收敛应该注意一下几个问题 17
1 收敛容差ANSYS缺省的收敛准则会根据单元的不同而检查不同的收敛力素和容差例如当采用solid65和link8时,缺省的要检查F和DISP两个力素其容差也是缺省的(Help中有)对于钢筋混凝土结构一般而言其位移比较小仅使用F力素收敛即可但其容差也同时放松一般采用5%即可(缺省是5)命令:cnvtol,f,,0.05,2 2 其它选项的设置
自动时间步打开此选择可以让程序决定子步间荷载增量的大小及其是增加或是减小收敛速度较快(命令autots,1)打开后似乎定义的子步数不起控制作用了 打开线性搜索可以帮助收敛的速度(命令:lnsrch,1) 打开预测器可以帮助收敛的速度(命令red,on)
平衡迭代次数在每一子步中的迭代次数缺省是25,将其增加例如改为50(命令: neqit,50) NSUBST此值不宜过小否则计算过程中老是调整影响计算速度
当然对于比较简单的算例或是分布模型可能不需要如此多的选项但对于复杂的模型是需要的各位可以试试
A1 结构有限元分析流行的几种模型
比较公认的有分离式模型(discrete model), 分布式模型(smeared model)(也有称整体式模型)组合式模型(也有称埋藏式模型)三种
分离式模型把钢筋和混凝土作为不同的单元处理例如混凝土采用solid65而钢筋采用link8等可以根据不同的单元类型进行组合在该模型中可以插入联结单元考虑粘结和滑移当然如果认为粘结很好也可以不考虑联结单元问题
分布式模型是把钢筋以一定的角度分布于整个单元中并认为二者粘结很好单元是连续均匀的材料 组合式模型也认为二者之间粘结很好没有相对滑移又分为分层组合式带钢筋膜的单元等方式与分布式模型比较二者的差异在于组合是分别求出各自的刚度矩阵进行组合而分布是在弹性矩阵中考虑各自对刚度的贡献
2 ANSYS可以考虑的模型
应该是两种即分离式模型(discrete model)和分布式模型(smeared model)分离式你可以选择不同的单元划分方式例如钢筋用体单元而非link8但是这种模型的收敛性常存在问题对于分布式模型solid65可以考虑三种钢筋材料应该也够用了这种模型比较容易得到收敛的解(这里的收敛是指尚没有达到破坏时候的解,当结构破坏时不能收敛是正常的)
--------------------------------------------------------------------------------
Q需要生成3d模型上某一条边的应力曲线必须先选择该边上的所有节点不知道用什么命令. Aselect nodes attached to lines
1. Utility menu -> select entities -> select the line from your 3d model 2. Utility menu -> select entities -> select nodes attached to lines
then all nodes on this line will be selected.
Q用选中的节点作一条应力曲线发现曲线上的应力最值与节点列表中的应力最值不一样在我的实例中列表中显示第55号节点s1为最值点110mpa曲线中包含了55节点但曲线的s1最值只有88mpa为何?
A我知道了原来是网格划分有问题
--------------------------------------------------------------------------------
Q钢筋混凝土在ANSYS中模态分析但其具有几何材料的非线性和两者之间的粘贴滑移 A模态分析中不能包含非线性若要考虑用直接求响应 18
AANSYS模态分析是线性分析任何非线性特性即使定义也将被忽略 --------------------------------------------------------------------------------
Q请问工字钢梁的几何特性怎么设置是否就把IAREAHEIGHT直接输进去就可以
A因所选单元不同输入的实常数不同如用BEAM3,则面积惯性矩梁高等就可以如采用实体则不必输入如采用梁截面则输入的信息又不同了
-------------------------------------------------------------------------------- Qcrack和crash有什么区别?
ACrack 指裂缝 crash指大面积的破坏 Qdamage和failure有什么区别和联系?
Adamage指损伤当材料或构件中产生微小或宏观裂纹时就表明存在损伤failure指失效即构件已经不能正常工作通常指工件断裂
--------------------------------------------------------------------------------
Q请问钢筋混凝土结构完全瞬态动力分析的阻尼矩阵如何输入谁知道钢筋混凝土框架结构的和值大概是多少隔震结构呢
AYou can use a formula in ansys help file to calculate an initial alfa and beita ! You can also determine those by curve fitting from the experimental data! Q材料性质中的damper和求解选项中的damper都要输入吗
A我想to calculate an initial alfa and beita 还需要知道结构试验中的各阶模态阻尼因子不知对不对
Q先进行模态分析就可求出初始的的好像是对吗
AIn order to calculate Alfa and Beita from formula, you need not modal experiment damping ratio.
In order find Alfa and Beita from experimental damping ratio by curve fitting, you need modal experimental damping ratio and modal frequency. Both are feasible from my experience.
--------------------------------------------------------------------------------
Q要得到一叶片的振型在ansys中利用其模态分析得到叶片的animate,但需要叶片振动的数学方程不知如何导出
AIf you 要得到一叶片的振型, you can *GET all modal Vector Data, you need not know the whole equation.
You can get the displacement in X, Y, Z direction for every node. --------------------------------------------------------------------------------
Q1.如用批处理形式输错的命令如何改这些命令输完后如何保存以便下次上机时进行调试 2.用交互模式进行结构分析后得出变形图如何保存下次怎么打开
A所谓批作业本身就是把一连串的命令放在一个文件中恐怕没有你说的保存问题它的所有输出信息都会放在一个文件中这个文件通常是jobname.out的ascii文件. 你可以通过打开这个文件获得信息