在ANSYS操作过程或条件语句中,常常需要知道有关模型的许多参数值,如选择集中的单元数、节点数,最大节点号等。此时,一般可通过*GET命令来获得这些参数。现在,对于此类问题,我们有了一个更为方便的选择,那就是查询函数 — Inquiry Function。
Inquiry Function类似于ANSYS的 *GET 命令,它访问ANSYS数据库并返回要查询的数值,方便后续使用。ANSYS每执行一次查询函数,便查询一次数据库,并用查询值替代该查询函数。
假如你想获得当前所选择的单元数,并把它作为*DO循环的上界。传统的方法是使用*GET命令来获得所选择的单元数并把它赋给一个变量,则此变量可以作为*DO循环的上界来确定循环的次数 *get, ELMAX,elem,,count *do, I, 1, ELMAX „ „
*enddo
现在你可以使用查询函数来完成这件事,把查询函数直接放在*DO循环内,它就可以提供所选择的单元数 *do, I, ELMIQR(0,13) „ „
*enddo
这里的ELMIQR并不是一个数组,而是一个查询函数,它返回的是现在所选择的单元数。括弧内的数是用来确定查询函数的返回值的。第一个数是用来标识你所想查询的特定实体(如单元、节点、线、面号等等),括弧内的第二个数是用来确定查询函数返回值的类型的(如选择状态、实体数量等)。
同本例一样,通常查询函数有两个变量,但也有一些查询函数只有一个变量,而有的却有三个变量。 查询函数的种类和数量很多,下面是一些常用、方便而快速快捷的查询函数 1 AREA—arinqr(areaid,key)
areaid—查询的面,对于key=12,13,14可取为0; key—标识关于areaidr的返回信息 =1, 选择状态 =12,定义的数目 =13,选择的数目 =14,定义的最大数 =-1,材料号 =-2,单元类型 =-3,实常数 =-4,节点数 =-6,单元数 „
arinqr(areaid,key)的返回值 对于key=1
=0, areaid未定义 =-1,areaid未被选择 =1, areaid被选择 „
2 KEYPOINTS—kpinqr(kpid,key)
kpid—查询的关键点,对于key=12,13,14为0 key —标识关于kpid的返回信息 =1,选择状态 =12,定义的数目 =13,选择的数目 =14,定义的最大数目 =-1,数料号 =-2,单元类型 =-3,实常数
=-4,节点数,如果已分网 =-7,单元数,如果已分网 kpinqr(kpid,key)的返回值 对于key=1 =-1,未选择 =0,未定义 =1, 选择
3 LINE—lsinqr(lsid,key)
lsid—查询的线段,对于key=12,13,14为0 key—标识关于lsid的返回信息 =1, 选择状态 =2, 长度
=12,定义的数目 =13,选择的数目 =14,定义的最大数 =-1,材料号 =-2,单元类型 =-3,实常数 =-4,节点数 =-6,单元数 „
4 NODE—ndinqr(node,key)
node—节点号,对于key=12,13,14为0 key—标识关于node的返回信息 =1, 选择状态 =12,定义的数目 =13,选择的数目 =14,定义的最大数 =-2,超单元标记 =-3,主自由度 =-4,激活的自由度 =-5,附着的实体模型
ndinqr(node,key)的返回值 对于key=1
=-1,未选择 =0,未定义 =1, 选择
5 VOLUMES—vlinqr(vnmi,key)
vnmi—查询的体,对于key=12,13,14为0 key—标识关于vnmi的返回信息 =1,选择状态 =12,定义的数目 =13,选择的数目 =14,定义的最大数目 =-1,数料号 =-2,单元类型 =-3,实常数 =-4,节点数 =-6,单元数 =-8,单元形状 =-9,中节点单元 =-10,单元坐标系
vlinqr(vnmi,key)的返回值 对于key=1 =-1,未选择 =0,未定义 =1, 选择
ANSYS能实现直接流-固耦合分析吗?
ANSYS流固耦合分析有三种形式,可以实现全直接或半直接耦合分析:
一: ANSYS/Mechanical模块或含该模块的软件包中的流固耦合分析功能,但此处的流体是非流动的流体,而是静流体,它计算流体由于重力、惯性力、波动压力等引起的分布压力载荷与结构的相互作用。 二:ANSYS/FLOTRAN/Structural模块或含该模块的软件包中的流固耦合分析功能,在建立好流体分析环境和结构分析环境的前提下,利用ANSYS5.7版本新增的FSSOLV命令,可自动实现流固耦合迭代计算,并可对迭代容差、流场收敛精度、流场网格变形等进行控制。
三:ANSYS/LS-DYNA模块或含该模块的软件包中的ALE任意欧拉-拉格朗日流固耦合分析功能,分析模型的流体部分用欧拉单元、直接求解流动方程,可通过速度、加速度、罚函数三种方式直接与结构单元耦合计算。
ANSYS中的地震响应分析
ANSYS中地震响应分析有二种基本的方法:时间历程分析和响应谱分析。
单自由度系统固有频率变化时对某一次地震地面运动的响应构成一条响应曲线,工程上的响应谱曲线实际
上是某一地震级别时各种地面运动响应曲线的包络线,这根包络线与单自由度的阻尼有关,因此在使用响应谱曲线时有二个关键因素:地震级别(对应最大加速度)和临界阻尼率.从理论上说,响应谱分析得到的数值应大于时间历程分析得到的最大响应,若时间历程分析得到的最大响应比响应谱分析得到的数值大,原因可能是阻尼的处理不当,如果无阻尼,由于某些频率的共振会使响应变大。
另外时间点的选取要看结构的固有频率,如果关心最大的结构固有频率为10Hz,则每秒钟至少有60-80个计算点
在ANSYS中建模后,能否把结构整体的质量、惯性矩等计算出来?
问:在ANSYS中建模后,能否把结构整体的质量、惯性矩等计算出来? 答:有两种办法。
1. Preprocessor>Operate>Calc Geom Items,可以求出面积、体积、形心、惯性矩等,如果给定了密度,也可以知道质量。 2. 在进行求解是你也可以看到程序会在OUTPUT窗口列出这些量。
在谱分析完成后,在组合模态之前如何获得结构每阶模态的响应
在谱分析完成后,进入POST1并利用*GET命令提取模态N的模态参与 系数,然后利用载荷工况(Load Case)功能创建载荷工况,并将提取 的模态参与系数指定为载荷工况的缩放系数,然后读入该载荷工况, 利用后处理功能处理。同理,可以获得其他模态的独立响应结果。 下面是获得一阶模态响应的命令流:
*GET,PF1,MODE,1,PFACT !* 提取1阶模态的参与系数PF1 /POST1
SET,FIRST !* 读入1阶模态结果
LCDEF,1,1,1, !* 将1阶模态结果定义载荷工况1 LCFACT,1,PF1, !* 载荷工况1的缩放系数为PF1 LCASE,1, !* 读入载荷工况1 /EFACE,1 !* 后处理操作 AVPRIN,0, , PLNSOL,U,SUM,0,1
在进行多点谱分析的时候,如何正确施加不同的谱曲线
在多点谱分析中,施加后面的谱曲线时,应该删除上次定义的谱曲线,如要对结构分别施加两个方向的激励,分别是,y和x方向。应如下处理: ...
!施加 y方向 PSDUNIT,1,ACCG,9.81 PSDFRQ,1,..... !频率值 PSDVAL,1,..... !谱值
D,BASENODE,uy,1.0 !施加谱曲线 PFACT,1,BASE
D,BASENODE,UZ,1.0 !删除上一次定义的谱
!施加 X 方向 PSDUNIT,2,ACCG,9.81 PSDFRQ,1,..... !频率值 PSDVAL,1,..... !谱值 D,BASENODE,UX,1.0
在非线性屈曲分析中,如何在结构上产生初始几何缺陷
有两种方式:
1、先进行特征值屈曲分析获得屈曲行为的理论解,这也可作为非线性屈曲分析的起点;
UPCOORD命令只是提供结构的初始位移,非常接近于真实缺陷情况。通过UPCOORD,,ON命令,就可以让我们从一个有轻微缺陷而没有任何初始应变的结构开始进行非线性屈曲分析。
2、还可以在结构上施加一个小的干扰力以引起结构的屈曲,这个方法简单而直接,也通常能得到合理结果,如果干扰力太大则会明显影响计算精度。 /PREP7
/TITLE, BUCKLING OF A BAR WITH HINGED ET,1,BEAM3 ! BEAM ELEMENT
R,1,.25,52083E-7,.5 ! AREA,IZZ,HEIGHT MP,EX,1,30E6 N,1 N,11,,100 FILL E,1,2 EGEN,10,1,1 FINISH /SOLU
ANTYPE,STATIC ! STATIC ANALYSIS
PSTRES,ON ! CALCULATE PRESTRESS EFFECTS M,2,UX,11,1 ! SELECT 10 UX DOF AS MASTERS D,1,ALL ! FIX SYMMETRY END
F,11,FY,-1 ! UNIT LOAD AT FREE END OUTPR,,1 SOLVE FINISH /SOLU
ANTYPE,BUCKLE ! BUCKLING ANALYSIS
BUCOPT,REDUC,1 ! USE HOUSEHOLDER SOLUTION METHOD, EXTRACT 1 MODE MXPAND,1 ! EXPAND 1 MODE SHAPE SOLVE /POST1 SET,LAST
UPCOORD,0.1,ON ! USE L/1000 FOR STARTING FACTOR
SAVE,PERT,DB ! SAVE TO NEW FILE TO PRESERVE ORIGINAL MODEL FINISH RESU,PERT,DB NROPT,FULL,,OFF /SOLU
ANTYPE,STATIC NLGEOM,ON AUTOTS,ON KBC,0 F,11,FY,-100 TIME,1
DELTIME,0.01,1E-4,0.5 SOLVE FINI
ANSYS后处理中如何显示三维实体模型表面结果云图和等值线
(1)将需要显示表面结果的三维实体模型的某些表面上的节点选出(Utility>Select) (2)将显示方式POWRGRPH设置为OFF(命令:/GRAPHICS,FULL)
(3)在/Post1下,绘制结果云图;或者在PlotCtrls>Device Options>Vector Mode„ON(命令/device,vector),可以绘制等值线。
在交互方式下如何施加任意矢量方向的表面载荷?
问:在交互方式下如何施加任意矢量方向的表面载荷?
答:若需在实体表面上施加任意方向的表面载荷,可通过在实体表面生成表面效应单元(比如SURF154单元)的方法来完成。
施加面载荷时,可施加在表面效应单元上,这样可以任意控制面力的方向。
加载过程中,选定表面效应单元,对话框中LKEY取值不同,则所加表面载荷的方向不同。(请仔细看一看surf154的单元手册)。
比如:LKEY=1(缺省),载荷垂直于表面;LKEY=2,载荷为+X切向;LKEY=3,载荷为+Y切向;LKEY=4,载荷垂直于表面;LKEY=5,则可输入任意矢量方向的载荷。 特别地:
LKEY=5,VALUE 项为均布压力值
VAL2、VAL3、VAL4 三项的值确定矢量的方向。
如何在ANSYS中存储动画?
Q: ANSYS中有好几种动画模式,但并非同时可用。那么有哪几种类型,何时可用呢?
A: ANSYS中有四种动画模式:Bitmap,AVI,Display List,Pixmap。该种模式是否可用取决于运行环境是UINX还是PC,以及是使用二维还是三维显示设备。可以通过选择菜单Utility Menu>PlotCtrls>Device Options(或键入适当的命令)来制作动画。下面是该命令选项的小结: 二维显示设备的PC:
Bitmap(命令-/device,anim,bmp)
AVI(命令-/device,anim,avi) 三维显示设备的PC:
Bitmap(命令-/dv3d,anim,1) AVI(命令-/dv3d,anim,2)
Display List(命令-/dv3d,anim,0) 二维显示设备的UNIX:
Pixmap(缺省,无相应命令) 三维显示设备的UNIX:
Pixmap(命令-/dv3d,anim,1)
Display List(命令-/dv3d,anim,0)
Q: 这些动画模式的区别是什么?
A: Display List选项(三维设备可用)与其他方式的不同在于:
Display List允许在播放动画的过程中对模型进行动态操作(如放大缩小等),而其它模式不能; 播放器不同:Bitmap,Display List,Pixmap动画模式在ANSYS图形窗口进行播放,可以通过ANSYS动画控制器面板控制动画;AVI(只在PC上)用WINDOWS媒体播放器播放。
需注意在ANSYS5.5中缺省的动画模式为AVI,播放器为媒体播放器。而在ANSYS5.6中二维设备的缺省动画模式为Bitmap,三维设备的缺省动画模式为Display List,动画在ANSYS图形窗口播放。
Q: 如何存储在ANSYS中创建的动画? A: 对于PC环境,Bitmap或AVI模式的动画会自动存储为jobname.avi(jobname是你指定的分析名称)。Display List动画在PC上不会自动存储,必须通过菜单Utility Menu> PlotCtrls >Animate>Save Animation指定存储,典型格式为jobname.anim。
对于UNIX环境,Pixmap及Display List都只能按ANSYS指定的格式进行存储,方法同上。
Q: A用Bitmap模式存储的.avi文件与用AVI模式存储的.avi文件有何不同? A: Bitmap模式的文件比AVI模式的文件要小,因为采用Bitmap模式存储的 .avi文件不包括回放的画面,它在ANSYS动画控制器中播放,可以向前或向后显示。而WINDOWS媒体播放器没有forward/backward选项,回放的画面必须包括在AVI动画文件中,因此文件较大。
Q: UNIX上制作的动画能否在PC上播放? A: 可以。可以使用ANSYS ANIMATE程序(ANIMATE.exe)在PC上播放.anim或.avi文件。还可以将.anim文件转变为.avi文件。所有的ANSYS产品都提供 ANIMATE程序,该程序位于ANSYS安装盘的根目录中。
Q: ANIMATE是否能读入所有的动画文件?
A: 不能。ANIMATE程序接受UNIX上的Pixmap模式的.anim文件或PC上的Bitmap或AVI模式的.avi文件。
如何提取模态质量
模态分析过程中打开振型型则化开关命令的Nrmkey设置为ON),ANSYS程序将自动将每阶模态的最大位移
单位化,就可以提取模态质量。计算方法如下:
1、利用SSUM对ETABLE动能数据求和获得结构总动能( );
2、 将结构总动能除以 得到 ,其中 是系统的角频率。下面是《ANSYS Verification Manual》VM89.DAT稍加修改后提取模态质量的例子: /PREP7
/TITLE, VM89, NATURAL FREQUENCIES OF A TWO-MASS-SPRING SYSTEM
C*** VIBRATION THEORY AND APPLICATIONS, THOMSON, 2ND PRINTING, PAGE 163,EX 6.2-2 ET,1,COMBIN14,,,2 ET,2,MASS21,,,4
R,1,200 ! SPRING CONSTANT = 200 R,2,800 ! SPRING CONSTANT = 800 R,3,.5 ! MASS = .5 R,4,1 ! MASS = 1 N,1 N,4,1 FILL
E,1,2 ! SPRING ELEMENT (TYPE,1) AND K = 200 (REAL,1) TYPE,2 REAL,3
E,2 ! MASS ELEMENT (TYPE,2) AND MASS = .5 (REAL,3) TYPE,1 REAL,2
E,2,3 ! SPRING ELEMENT (TYPE,1) AND K = 800 (REAL,2) TYPE,2 REAL,4
E,3 ! MASS ELEMENT (TYPE,2) AND MASS = 1 (REAL,4) TYPE,1 REAL,1
E,3,4 ! SPRING ELEMENT (TYPE,1) AND K = 200 (REAL,1) M,2,UX,3
OUTPR,BASIC,1 D,1,UY,,,4 D,1,UX,,,4,3 FINISH
/SOLU
ANTYPE,MODAL
MODOPT,subspa,2,,,2,ON MXPAND,2,,,YES SOLVE
FINISH /post1
set,1,1
etabl,kene,kene ssum
*get,keneval1,ssum,,item,kene *get,freqval1,mode,1,freq
eigen1=(2*3.14159*freqval1)**2 pmass1=2*keneval1/eigen1
set,1,2
etabl,kene,kene ssum
*get,keneval2,ssum,,item,kene *get,freqval2,mode,2,freq
eigen2=(2*3.14159*freqval2)**2 pmass2=2*keneval2/eigen2
finish
如何考虑结构分析中的重力
在结构分析中,如何模拟结构自重和设备重量是一个经常遇到的问题,对于结构自重有两点要注意: 1. 在材料性质中输入密度,如果不输入密度,则将不会产生重力效果。
2. 因为ANSYS将重力以惯性力的方式施加,所以在输入加速度时,其方向应与实际的方向相反。 对于结构上的设备重量可以用MASS21单元来模拟,该单元为一个空间“点”单元。设备重量可通过单元实常数来输入。下面附上一个小例子(设重力方向向下)。 /prep7 et,1,42 et,2,21
r,2,10,10,10 mp,ex,1,2e5 mp,nuxy,1,0.3 mp,dens,1,1 rect,,10,,1 esize,.5 amesh,all type,2
real,2
e,node(5,1,0) fini /solu dk,1,all dk,2,uy, acel,,10 solve fini /post1
plnsol,u,sum,2
/SOLU
ANTYPE,MODAL
MODOPT,subspa,2,,,2,ON MXPAND,2,,,YES SOLVE FINISH /post1
set,1,1
etabl,kene,kene ssum
*get,keneval1,ssum,,item,kene *get,freqval1,mode,1,freq
eigen1=(2*3.14159*freqval1)**2 pmass1=2*keneval1/eigen1
set,1,2
etabl,kene,kene ssum
*get,keneval2,ssum,,item,kene *get,freqval2,mode,2,freq
eigen2=(2*3.14159*freqval2)**2 pmass2=2*keneval2/eigen2
finish
声学中的阻尼等相关概念
明确ANSYS中的阻尼,声吸收,阻抗的含义:
阻尼是指动力学问题相关的能量损失,可以在瞬态或谐波声学中包括。声的吸收和阻抗指压力自由度相关的损失。ANSYS中的阻抗用来标识声表面可以吸收能量的开关,MU指能量在指定声表面被吸收的数量。这个用途对ANSYS是特殊的,意义比广义声学中更为严格。
通常的一个误解是约束的边界是吸收边界。实际上这种边界反射压力脉冲并将其反号。各种边界条件总结如下:
MU值 DOF(自由度约束) 结果边界条件 u=0 未约束 无压力反号
Mu=1 未约束 吸收边界(仿佛另一侧有相同材料) Mu=∞ 未约束 压力反向的反射边界 Mu=any 约束 压力反向的反射边界
Mu=0 模拟刚性壁条件:无吸收,100%反射声能。Mu<1表示(至少是典型如此)声波从低密度流体进入高密度流体。例如声波在空气中传播碰到空气/水界面就像遇到刚性墙壁,因此Mu会很小,为0.05。在谱的另一端,MU=∞相应于压力释放(P=0)边界。声在水中传播遇到空气/水界面就如同是p=0边界。这样大的MU值可以用于模拟声在水中传播的空气/水边界。如果要模拟声从高密度媒质到低密度媒质,设定的MU值应大于1。
下面例子示意了阻尼和声吸收的使用。这个问题是声学管,类似于管弦乐和弦,施加到一端的压力向另一端传递在尽头反射。问题包括压力波的几次反复,表明在管封闭端的吸收。包括了不同的阻尼值(对阻尼矩阵)和MU(吸声端)。阻抗值对全反射边界为0,有吸收的为1。 /show,acous,f33 *dim,dval,,5 *dim,mval,,5
dval(1)=.01 $ mval(1)=0 dval(2)=.05 $ mval(2)=.25 dval(3)=.1 $ mval(3)=.50 dval(4)=.2 $ mval(4)=.75 dval(5)=.3 $ mval(5)=1
*do,ii,1,5 ! loop on damping *do,jj,1,5 ! loop on absorption
parsav,all ! save parameters for /clear /clear, ! start new problem parres
/TITLE, TRAVELING ACOUSTIC WAVE, DAMP = %dval(ii)%, MU = %mval(jj)% /prep7
et,1,29,,1 ! 2d acoustic fluid n,1 ! define nodes
n,101,100 ! make tube 100 long fill
n,201,,1 n,301,100,1
fill
e,1,2,202,201 ! define elements *repeat,100,1,1,1,1 dens,1,9.633e-5
sonc,1,10000 ! define sonic for 100 hz freq
damp,1,dval(ii) ! transient problem, need some damping mu,1,mval(jj) ! specify sound absorption at end FINISH /solu
antype,tran outres,all,all
nsel,s,,,101,301,200 ! select nodes at end sf,all,impd,1 ! enable sound absorption nsel,all !
! apply pressure at x = 0 time,.0005
d,1,pres,1,,201,200 ! apply pressure at start deltim,.0001 ! delta time for solutions solve
d,1,pres,0,,201,200 time,.001 solve
time,.07 ! at least two traverses solve finish /post26
nsol,2,21,pres,,.2*LENG ! store pressure at some locations nsol,3,41,pres,,.4*LENG nsol,4,61,pres,,.6*LENG nsol,5,81,pres,,.8*LENG nsol,6,101,pres,,END
plvar,2,3,4,5,6 ! plot all finish *enddo *enddo
如何定制Beam188/189单元的用户化截面
ANSYS提供了几种通用截面供用户选用,但有时不能满足用户的特殊需求。为此,ANSYS提供了用户创建截面(库)的方法。如果你需要创建一个非通用横截面,必须创建一个用户网格文件。具体方法是,首先创建一个2-D实体模型,然后利用SECWRITE命令将其保存(Main Menu>Preprocessor>Sections> -Beam-Write Sec Mesh)。该过程的细节如下:
创建截面的几何模型(二维面模型)。
对所有线设置单元份数或者单元最大尺寸 (Main Menu>Preprocessor> -Meshing-Size Cntrls>-Lines-Picked Lines或使用MeshTool)。记住:保证模型中的区格(cell)数目不能超过250个。 选择菜单Main Menu>Preprocessor>Sections>-Beam-Write Sec Mesh,弹出一个拾取窗口,单击Pick All拾取包含区格的所有面。
ANSYS自动在所有面上创建区格。在划分网格时,ANSYS可能显示单元形状差的消息,也可以被忽略不显示,但是总能看到一条消息“Unable to mesh area....”。如果已经完成上述工作,你清除所有面上的单元(Main Menu>Preprocessor>-Meshing-Clear> Areas),并重复第2、3、4步,即必须重新控制网格密度。
如果第4步成功则弹出Write Section Library File对话框,File Name域填入一个未用过SECT文件名,Drives域指定一个截面文件存放驱动器,Directories域指定一个截面文件存放目录,然后单击按钮OK,完成用户截面文件建立。如果在第3步中给线指定太多单元份数,区格和节点的数目可能超过限制范围,必须清除所有面上的单元2~4步,直到获得合适数量的区格和节点。
一旦完成上述工作,用户可以在以后分析中读取用户网格文件 (Main Menu > Preprocessor > Sections > -Beam-Read Sect Mesh),定义成适当的截面ID号,就与通用梁完全一致。要想检查用户截面,只要绘制截面(网格)图(Main Menu>Preprocessor>Sections>Plot Section)或者列表截面属性(Main Menu>Preprocessor>Sections>List Sections)。
耦合及约束方程讲座 一、耦合
当需要迫使两个或多个自由度取得相同(但未知)值,可以将这些自由度耦合在一起。耦合自由度集包含一个主自由度和一个或多个其它自由度。 典型的耦合自由度应用包括: 模型部分包含对称;
在两重复节点间形成销钉、铰链、万向节和滑动连接; 迫使模型的一部分表现为刚体。 如何生成耦合自由度集
在给定节点处生成并修改耦合自由度集 命令:CP
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs
在生成一个耦合节点集之后,通过执行一个另外的耦合操作(保证用相同的参考编号集)将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集(即增、删节点或改变自由度标记)可用CPNGEN命令。(不能由GUI直接得到CPNBGEN命令)。
耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点(诸如一条缝处)尤为有用。 命令:CPINTF
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes 除耦合重复节点外,还可用下列替换方法迫使节点有相同的表现方式:
如果对重复节点所有自由度都要进行耦合,常用NUMMRG命令(GUI:Main Menu>Preprocessor>Numbering Ctrls>Merge Items)合并节点。
可用EINTF命令(GUI:Main Menu> Preprocessor>Create> Elements >At Coincid Nd)通在重复节点对之间生成2节点单元来连接它们。
用CEINTF命令(GUI:Main Menu>Preprocessor> Coupling/Ceqn >Adjacent Regions)将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。
生成更多的耦合集
一旦有了一个或多个耦合集,可用这些方法生成另外的耦合集:
用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。 命令:CPLGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same Nodes
用下列方法生成与已有耦合集不同(均匀增加的)节点编号但有相同的自由度标记的新的耦合集: 命令:CPSGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same DOF 使用耦合注意事项
每个耦合的节点都在节点坐标系下进行耦合操作。通常应当保持节点坐标系的一致性。 自由度是在一个集内耦合而不是集之间的耦合。不允许一个自由度出现在多于一个耦合集中。 由D或共它约束命令指定的自由度值不能包括在耦合集中。
在减缩自由度分析中,如果主自由度要从耦合自由度集中选取,只有主节点的自由度才能被指定为主自由度。
在结构分析中,耦合自由度以生成一刚体区域有时会引起明显的平衡破坏。不重复的或不与耦合位移方向一致的一个耦合节点集会产生外加力矩但不出现在反力中。
耦合及约束方程讲座
二、约束方程
约束方程提供了比耦合更通用的联系自由度的方法。有如下形式:
这里U(I)是自由度,N是方程中项的编号。
如何生成约束方程
1. 直接生成约束方程
o
直接生成约束方程: 命令:CE
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Constraint Eqn
下面为一个典型的约束方程应用的例子,力矩的传递是由BEAM3单元与PLANE42单元(PLANE42单元无平面转动自由度)的连接来完成的:
图12-1建立旋转和平移自由度的关系
如果不用约束方程则节点2处表现为一个铰链。下述方法可在梁和平面应力单元之间传递力矩,自由度之间满足下面的约束方程: ROTZ2 = (UY3 - UY1)/10 0 = UY3 - UY1 - 10*ROTZ2 相应的ANSYS命令为:
CE,1,0,3,UY,1,1,UY,-1,2,ROTZ,-10
o
修改约束方程
在PREP7或SOLUTION中修改约束方程中的常数项: 命令:CECMOD
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Modify ConstrEqn Main Menu>Preprocessor>Loads>Other>Modify ConstrEqn Main Menu>Solution>Other>Modify ConstrEqn
如果要修改约束方程中的其它项,必须在求解前在PREP7中用使CE命令(或相应GUI途径)。
2. 自动生成约束方程
o
生成刚性区域
CERIG命令通过写约束方程定义一个刚性区域。通过连接一主节点到许多从节点来定义刚性区。(此操作中的主自由度与减缩自由度分析的主自由度是不同的) 命令:CERIG
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Rigid Region
将CERIG命令的Ldof设置为ALL(缺省),此操作将为每对二维空间的约束节点生成三个方程。这三个方程在总体笛卡尔空间确定三个刚体运动(UX、UY、ROTZ)。为在二维模型上生成一个刚性区域,必须保证X─Y平面为刚性平面,并且在每个约束节点有UX、UY和ROTZ三个自由度。类似地,此操作也可在三维空间为每对约束节点生成六个方程,在每个约束节点上必须有(UX、UY、UZ、ROTX、ROY和ROTZ)六个自由度。
输入其它标记的Ldof域将有不同的作用。如果此区域设置为UXYZ,程序在二维(X,Y)空间将写两个约束方程,而在三维空间(X、Y、Z)将写三个约束方程。这些方程将写成从节点的平移自由度和主节点的平移和转动自由度。类似地,RXYZ标记允许生成忽略从节点的平移自由度的部分方程。其它标记的Ldof将生成其它类型的约束方程。
总之,从节点只需要由Ldof标记的自由度,但主节点必须有所有的平移和转动自由度(即二维的UX、UY和ROTZ;三维的UX、UY、UZ、ROTX、ROTY、ROTZ)。对由没有转动自由度单元组成的模型,应当考虑增加一个虚拟的梁单元以在主节点上提供旋转自由度。
o
将疏密不同的已划分网格区域连在一起
可将一个区域(网格较密)的已选节点与另一个区域(网格较稀)的已选单元用CEINTF命令(菜单途径Main Menu>Preprocessor>Coupling / Ceqn>Adjacent Regions)连起来生成约束方程。
这项操作将不相容网格形式的区域“系”在一起。在两区域的交界处,从网格稠密的区域选择节点A,从网格粗糙区域选择单元B,用区域B单元的形函数,在相关的区域A和B界面的节点处写约束方程。ANSYS允许这些节点位置使用两公差准则。节点在单元之外超过第一公差就认为节点不在界面上。节点贴近单元表面的距离小于第二公差则将节点移到表面上,见下图。
对CEINTF命令有些限制:应力或热通量可能会不连续地穿过界面。界面区域的节点不能指定位移。可用每节点有六个自由度的单元接合6自由度实体。
o
从已有约束方程集生成约束方程集
可用CESGEN命令从已有约束方程集生成约束方程。那么已有约束方程集内的节点编号将增加以生成另外的约束方程集。另外约束方程集的标记和系数保持与原集的一致。 命令:CESGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn >Gen w/same DOF
使用约束方程的注意事项
所有的约束方程都以小转动理论为基础。因此,它应用在大转动分析中〔NLGEOM〕应当限制在约束方程所包含的自由度方向无重大变化的情况。
约束方程的出现将产生不可预料的反力和节点力结果。
由于相邻区域网格疏密不同,边界上的相容性仍然存在。但是当网格越密,这种不相容的危害就越小。
ANSYS后处理查看某一剖面的数据云图的方法 无所谓 @ 2007-08-24 21:46
要观察某一剖面的数据云图,可以利用“工作平面”来实现,方法如下:
1、旋转、移动坐标系,使“工作面”(即:XOY面)和我们想观察的剖面重合。
2. 在PlotCtrls>Style>HiddenLineOptions里作两点设置: A. 设Type of Plot为Section;
B. 设Cutting plan is为Working plane.
这样,就可以在后处理中观察到该剖面的数据云图了,够简单吧!
求助:如何查看模型中内部某个零件的振型图和变形云图
方法一:(自己总结)
1. 先 select/entities/volumes(by number/pick),选中想要查看结果的零件或部件 2. 然后 select/entities/volumes (select below),再选中想要查看结果的零件或部件 3. 进入后处理,查看结果时,只有选择的零部件 方法二:
选中那个体,再用以下命令
ALLSEL,BELOW,VOLU 再进入后处理 方法三:
(1)先选择希望执行绘图操作的体,例如:vsel,s,,,volu-num等等; (2)再利用plnsol或者plesol命令执行绘图操作;
(3)练习以下,还可以用/trlcy命令对于部件和体等对象透明处理。
网格划分的问题
一个整体由于不规则,为了划分网格方便分成多个部分,用devide命令将其分割。
devide有(blank)和 sepo 两个选项,如果是共用面(线)(blank),分割后在划分网格前就不需要glue命令进行粘合操作了。
如果是先后建立不同的实体,应该是先GLUE再划分网格,不然的话在计算的时候,各个部分是分开脱离的,应力、应变不能从一部
分传递到另外一部分。
进行GLUE操作,是为了使同位置的节点耦合到一起,使他们不是彼此无联系的部分,否则计算结果肯定是不对的。
先GLUE再划分网格,但为了划分后合并、压缩重复节点和单元,最好进行以下操作:
1.Main Menu>Preprocessor>Numbering Ctrls>Merge Items
在弹出的对话框中,Label选项选择ALL,单击OK完成编号的合并;
2.Main Menu>Preprocessor>Numbering Ctrls>Compress Numers
在弹出的对话框中,Label选项选择ALL,单击OK完成编号的重排列。
这个操作还是很有必要的,也算是实现耦合吧,就是同位置的节点单元的合并,这里的压缩就是从新排列接点号,对于统计节点、单
元数量是必须的操作。
如何为不同的零件或部分指定不同的单元类型、材料
在划分网格之前,先设置 preprocessor/meshing/mesh attributes/default attribs
选择单元编号element number 或者材料编号material number ,,点击ok然后退出,进行网格划分,所划分的网格、材料就是指定的类型。
如何扩展对称结构的计算结果
结构对称的构建进行经理分析时,可以只分析1/2或者1/4,然后再对结果进行扩展,扩展操作如下:
plotctrls/style/symmetry expansion///然后可以指定轴对称、周期性对称、或者循环性对称,1/2, 1/4, 1/8等等。
我需要做一个结构瞬态动力学分析,不知道变化载荷怎么加,求教各位了 6楼:
将你的分析选择瞬态分析模块Transient后,就开始对模型进行加载:1)首先对模型进行约束;2)加载载荷,3)对载荷步进行设置Main Menu>Preprocessor>Load>Time/Frequenc>Time-Time Step;4)读入载荷步Main Menu>Preprocessor>Load>Write LS File。即完成一个载荷步的加载,如需加载第二个载荷步就重复步骤2—4的过程。5)求解。
注意:载荷步设置时,定义在本载荷步结束时的载荷时间。读入载荷步时,一定按照顺序1、2、3读入,这样方便计算时候设置。
7楼:
6楼说的极是,我没这样做过,学了一招,呵呵
弱问:云图中intensity和von mises的不同
stress intensity(应力强度),是由第三强度理论得到的当量应力,其值为第一主应力减去第三主应力。 von mises应力是由第四强度理论得到的当量应力。 材料力学里有关于四个强度理论的叙述。
第三强度理论认为最大剪应力是引起流动破坏的主要原因
如 低碳钢拉伸时在与轴线成45度的截面上发生最大剪应力,材料沿着这个平面发生滑移,出现滑移线。 这一理论比较好的解释了塑性材料出现塑性变形的现象。形式简单,但结果偏于安全。 第四强度理论认为形状改变比能是引起材料流动破坏的主要原因。结果更符合实际。 一般材料在外力作用下产生塑性变形,以流动形式破坏时,应该采用第三或第四强度理论。 如果你所研究的结构破坏是因为剪应力,就应该以intensity作为结构破坏准测。 多数情况时采用von mises作为破坏准测。
当然还有许多其他的破坏准测,具体采用何种要根据实际情况
内容:
-模态分析得到的结果是不是某个方向上的各阶频率啊我要得到各个方向的一阶频率能做到吗
-承是用来支撑要研究的对象把轴承等价为弹簧将其刚度加在对象的节点上这样就要定义刚度的大小在加是否还要定义单元类型然后才能定义其刚度
-处理动力问题如土层地震响应或基础振动反应可否从边界施加加速度荷载, 若可以如何瞬态动力反应分析时间子步为>200?
-板单元的后处理的显示在单元坐标还是在总体坐标? 比如位移,应力
-我现在需分析一个板梁结构板已用SHELL63 单元划分好梁我是用板上的一条线划分单元并添加截面而生成的但现在运算时发现板和梁是分开的它们之间互不影响请教各位高手怎样将板和梁合并为一个整体
-用SHELL63 单元划分的板结构想加垂直于板面的匀布载荷?
-ANSYS 可以自己定义输出格式吗?
-作模态分析时就是集中质量单元和刚度加法?
-我在原有的静态分析的基础上想继续作一下结构的模态分析应该采取那些步骤啊?
-1.在ANSYS 中如何设定小数的有效位数比如我输入的应力值为0.0004495568 而且我的一组数据都是最后几位差别但是ANSYSY 自做聪明把它四舍五入了结果弄的都一样了如何更改这个有效位数呢 2.ANSYS 在输出时如何修改输出的数字显示模式比如用科学记数法或者将数字缩放多少倍那样显示在什么地方能够修改呢
3.ANSYS5.7.1 中什么地方能够出现那种材料库表就是说不同ANSYS 的什么material models 对话框生成
而是通过table 输入材料属性诸如应力应变曲线之类的东西?
-我是个ANSYS 初学者现在碰到下面几个问题请各位高手帮忙看一下
1 是什么原因使得相邻ELEMENT 的共用节点解得的力不一样应该是方法步骤的原因而不是误差吧请简单介绍一下
2 假如我想用ANSYS 模拟陶瓷的烧结过程是不是有实际*作的可能性如果有其基本的原理和模型怎样 3 ANSYS 可以用来模拟热循环下的应力分配不均造成的焊点63%Sn&37%Pb 失效的过程 我看到有一篇文献介绍用在Sn 球里随机分布Pb 相和其它共晶相的方法建模的那么这个 随机过程在ANSYS 中怎么实现的而且如果要生成足够多的晶相小颗粒是不是只能手动一 个一个的建立
-在进行模拟计算结构多阶段施工过程计算时,在ansys 里面说可以用死活单元的方法但是死活单元每次只允许选择一个,我的模型比较大,如果采用这个方法不太现实,各位同志有什么比较好的方法没有请赐教
-如在beam4 和shell63 中显示的输出结果有应力和变形以及节点等效荷载请问如何才能显示截面弯矩扭矩和剪力呢
-最近在计算一个模型管状在地下分析中要施加地层给予的支撑我想简化为弹性支撑可是发现在ansys 中找不到这样的支撑情况
-用弹簧可是问题很多因为是一个管体整个管壁都有支撑简化为弹簧出现的问题很多了比如用多少个合适隔多长合适线性支撑和非线性支撑能等同吗我最近看校验手册发现V135 这个例子可以借鉴可是还是有许
多麻烦哪位大侠如果有这方面经验请指点指点
-如何梁单元同时施加两个方向分布荷载竖向侧向
-请问如何打印图形结果
-ansys5.6.2 计算一个压力容器,划分网格9 万余个.运行环境win2000,cpu:1G,内存640M,硬盘40G,模型建立,网格划分均无问题.但是运行到占用硬盘大约9.2G 的时候,就出错.提示内存无法read.请问该如何解决?是内存的问题?还是ansys 本身问题?还是系统问题?硬件问题?模型建立实在无法简化.
-我在求解是出现如下的提示
1. real constant 1 referenced by at least element types 1 and 2. 2. 2.coefficient ratio exceeds 1.0e8.
我全部忽略上述提示得出了结果结果能可信吗我有错在哪里呢模型我都检查了数遍了还是和初始文件名有关我是在一个初始名下面建立若干相关的模型
由于我的单元只是板梁的结合所以我重新检查了一下删除了一个影响不大的扁梁很扁尺寸不大系统可能认为是板了运行就通过了上述的提示可能是因为系统对某些单元有默认的尺寸限制
-slab on girder highway bridge 桥板为SOLID 单元由于girder 上部flange 嵌贴在板的下表面所以建模时我将桥板定义为SOLID 单元桥板下表面AREA 定义为SHELL63 单元但SOLID 和SHELL 单元的自由
度不同这样建模肯定有不对的地方运算结果与实验结果也相差了40% 有什么好的解决办法
-Q ansys 有没有评价一个网格划分好坏的参数或命令请告知
-Q (1)在gui 下运行宏时,重新定义原来已有的数组,ansys 就会有个提示,想不要这个提示,怎么办?(/uis,msgpop,3 这个行不通,我试过了)
(2)ansys 中数组的定义为a*b*c a<10040 b<255 c<7 而我实际上要定义更大的有什么好办法,难道非要写到文件里再读进来?
-Q 我现在做一个板壳的屈曲分析想要得到该结构在荷载作用下的极限承载力利用非线性屈曲分析请问极限承载力的荷载比例因子在什么命令下可以查到
-Q 如果材料一定为各向异性我用的材料为钢板该如何选用单元类型 ANSYS使用问答精华 Q为问A为回复
------------------------------------------------------------------------------- Q:模态分析得到的结果是不是某个方向上的各阶频率啊我要得到各个方向的一阶频率能做到吗
A模态分析得到的结果是你所选择的自由度内的振型如果需要得到某一个方向内的振型用RUDECED法选择你所希望的主自由度如UX将所求的频率值设置为一即可
Q不过大型复杂结构的振型一般都不是一个方向的应该是一个方向为主其它方向为次如果把所要求的方向
设为主自由度会不会跟真是实际结果有出入
Q对于复杂的大型实体其主自由度无法给出此时就无法应用reduce法可不可以改约束还是就认为一次得到前3阶频率就是x,y,z方向的一阶频率啊
A正是因为复杂结构的主自由度不好确定所以我还是倾向于用SUBSPACE和BLOCK LANCZOS法约束不要随便改它应该描述结构真实的状态不能说前3阶频率即为XYZ方向的第一阶频率你可以好好想想所用方程中[K]和[M]的意义
--------------------------------------------------------------------------------
Q轴承是用来支撑要研究的对象把轴承等价为弹簧将其刚度加在对象的节点上这样就要定义刚度的大小在加是否还要定义单元类型然后才能定义其刚度
A我曾经用LINK8单元模拟过根据K=EA/L确定刚度我当时取E=210e9,L=0.1,然后A=KL/E经过检验是可行的
A1用link单元模拟其刚度 k=EA/L通过假定E,A,L的值使之和k相等其密度要02用combin14模拟其刚度阻尼0
--------------------------------------------------------------------------------
QANSYS处理动力问题如土层地震响应或基础振动反应可否从边界施加加速度荷载若可以如何瞬态动力反应分析时间子步为>200?
A如果你所希望求的最高频率的周期为f则ANSYS取时间步长为1/20f为默认值 Q板单元的后处理的显示在单元坐标还是在总体坐标? 比如位移,应力 A应力位移在总体坐标弯矩不一定
A应力位移是显示在结果坐标系的顺便说明一下节点信息显示在显示坐标系
不过结果坐标系默认为全局坐标系而显示坐标系也是默认为全局坐标系 你可以这样修改结果坐标系
Main Menu:General Postprocessor>Options for Output... 你可以这样修改显示坐标系
Utility Menu:Workplane>Change Display CS to
A你可以用ETABLE命令先存储对应于各单元的位移等等再用ANSYS的柱状图或是别的画图软件PLOT出图形你认为如何
Q我现在需分析一个板梁结构板已用SHELL63单元划分好梁我是用板上的一条线划分单元并添加截面而生成的但现在运算时发现板和梁是分开的它们之间互不影响请教各位高手怎样将板和梁合并为一个整体 APreprocessor>Numbering Ctrls>Merge Items里element and node 合并 Q: 用SHELL63单元划分的板结构想加垂直于板面的匀布载荷?
A: 找到如下位置ANSYS HELP->ANSYS ELEMENTS REFERENCE->ELEMENT LIBRARY->SHELL63 ELASTIC SHELL 在INPUT SUMMARY 中有SURFACE LOADS PRESSURES 1
其中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我在原有的静态分析的基础上想继续作一下结构的模态分析应该采取那些步骤啊 我用solution>new 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阶!
Q1.在ANSYS中如何设定小数的有效位数比如我输入的应力值为0.0004495568而且我的一组数据都是最后几位差别但是ANSYSY自做聪明把它四舍五入了结果弄的都一样了如何更改这个有效位数呢
2.ANSYS在输出时如何修改输出的数字显示模式比如用科学记数法或者将数字缩放多少倍那样显示在什么地方能够修改呢
3.ANSYS5.7.1中什么地方能够出现那种材料库表就是说不同ANSYS的什么material models对话框生成而是通过table输入材料属性诸如应力应变曲线之类的东西? A用a.bEn来表示如2.45678e-9材料特性用tb
-------------------------------------------------------------------------------- Q我是个ANSYS初学者现在碰到下面几个问题请各位高手帮忙看一下
1是什么原因使得相邻ELEMENT的共用节点解得的力不一样应该是方法步骤的原因而不是误差吧请简单介绍一下
2假如我想用ANSYS模拟陶瓷的烧结过程是不是有实际*作的可能性如果有其基本的 2
原理和模型怎样
3ANSYS可以用来模拟热循环下的应力分配不均造成的焊点63%Sn&37%Pb失效的过程我看到有一篇文献介绍用在Sn球里随机分布Pb相和其它共晶相的方法建模的那么这个随机过程在ANSYS中怎么实现的而且如果要生成足够多的晶相小颗粒是不是只能手动一个一个的建立
A可能是你提取数据的时候思路有问题比如说当你提取的是有多个连接单元的节点数据而实际上你提取的是各个单元的
A相邻ELEMENT的共用节点你MERGE了吗ANSYS是通过单元的节点传递数值的从而得出结果的 Q在进行模拟计算结构多阶段施工过程计算时,在ansys里面说可以用死活单元的方法但是死活单元每次只允许选择一个,我的模型比较大,如果采用这个方法不太现实,各位同志有什么比较好的方法没有请赐教 A涉及到两个方面的问题一以TABLE形式定义参数二DO-LOOP循环你可以参考ANSYS GENERAL GUIDE中SOLUTION中对多重载荷步的叙述很简单
Q比如在beam4和shell63中显示的输出结果有应力和变形以及节点等效荷载请问如何才能显示截面弯矩扭矩和剪力呢
A在ELEMENT TABLE 中首先定义变量注意一你要知道BEAM设置了几个积分点注意二要知道SEQUENCE NUMBER值所以你需要参考ELEMENT REFERENCE 中的单元介绍
Q最近在计算一个模型管状在地下分析中要施加地层给予的支撑我想简化为弹性支撑可是发现在ansys中找不到这样的支撑情况
A首先你要知道弹性支撑的弹性系数具体的计算方法你应该知道的 然后用LINK8单元模拟要注意的是刚度k=EA/L由此确定LINK单元的面积A来体现单元的弹性刚度我对此作过验证结论是可行的
Q用弹簧可是问题很多因为是一个管体整个管壁都有支撑简化为弹簧出现的问题很多了比如用多少个合适隔多长合适线性支撑和非线性支撑能等同吗我最近看校验手册发现V135这个例子可以借鉴可是还是有许多麻烦哪位大侠如果有这方面经验请指点指点
-------------------------------------------------------------------------------- Q如何梁单元同时施加两个方向分布荷载竖向侧向
A同时施加两个方向的均布荷载相当与分别施加不同的均布荷载即分别加
主要是要选择好均布荷载的ID参考单元参考手册中的说明例如BEAM4 INPUT SUMMARY 中
Surface Loads Pressure: 有 face2(i-j)(-Y normal direction) 这里2就是LOAD ID以单元坐标的-Y轴为正你试一试看 Q请问如何打印图形结果
A用ANSYA中的功能就可以拉GUI菜单命令PLOTCTRLSCAPTURE IMAGES或RESTORE IMAGES A我是用plotctrl>hard copy>to file,然后在剪贴板中就可以看见了 A我都是用抓图工具抓下来然后再打印的
A/ui,copy,save是一个截取屏幕图形的命令,现在总算又实现了gui肯定可以通过相应的命令方式来实现的观点
AAnsys 5.6已经做到了把图形窗口COPY到文件如JPG/BMP...等格式PLOTCTRL/HARD COPY/TO FILE很好用
-------------------------------------------------------------------------------- 3
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的内存1GOS windowsNT,Total Workingspace是1400M,再高就控制不了Database 是56,我能解的有效节点数是140000,只是因为磁盘空间不够工作的分区是0G,当ansys探测到磁盘空间已满就自动退出尽量大的增大节点数并不是好办法因为运行速度太慢了利用对称性减小求解区域利用Lesize,esize等控制网格都是好的办法
A选用对称性的时候要注意边界条件但的确是个好办法
A以前在使用ALGOR时也遇到过类似情况当时有10万个结点对WIN2000来说已不存在WIN9X中文件最大只能2G的瓶颈但在处理超大虚拟内存需求时仍然会出问题处理方法有二
1使用PQMAGIC划分出一个20G的逻辑分区一定要跟WIN2000系统在不同的分区再在控制面板-系统中将虚拟内存文件指定到这个大分区上以尽量避免发生系统堆栈错误虚拟文件大小仍交WINDOWS控制 2将至少两个分区的硬盘升级至NTFS由于NTFS文件可以跨分区存储且要小10%左右比FAT32要好用得多 Q我在求解是出现如下的提示
1. real constant 1 referenced by at least element types 1 and 2. 2. 2.coefficient ratio exceeds 1.0e8.
我全部忽略上述提示得出了结果结果能可信吗我有错在哪里呢模型我都检查了数遍了还是和初始文件名有关我是在一个初始名下面建立若干相关的模型
由于我的单元只是板梁的结合所以我重新检查了一下删除了一个影响不大的扁梁很扁尺寸不大系统可能认为是板了运行就通过了上述的提示可能是因为系统对某些单元有默认的尺寸限制
A有些单元类型在默认值状态下不需要输入REAL CONSTANT但是当你改变了某项选项你就需要重新输入REAL CONSTANT你可以用UTILITY>SELECT>ENTITIES命令中的SELECT ELEMENT BY ATTRIBUTE来查看单元材料REAL CONSTANT等等是否定义明确然后用APPLY REPLOT命令来显示每个不同材料的单元
A我试了一下按照你所说的步骤出现了如下提示a default command range is not provided for the ESEL 4
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第二个问题经常出现.第一个问题有点问题,你看看单元定义对没有如果单元都没有划分对,结果当然不可靠啦
--------------------------------------------------------------------------------
Qslab 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对应点的UXUYUZ全部耦合起来让这些NODE的XYZ三个方向的位移相同但是结果仍然不理想我想约束方程的功能也进乎如此有没有更进一步的建议十分感谢
-------------------------------------------------------------------------------- Qansys有没有评价一个网格划分好坏的参数或命令请告知菜单*作
APREPROCESSOR>MESHING>CHECK MESH>两个命令中一个用来CHECK网格划分另一个用来CHECK连贯性只是CHECK完了别哭就好 APLOT BAD ELEMENT!
-------------------------------------------------------------------------------- Q(1)在gui下运行宏时,重新定义原来已有的数组,ansys就会有个提示,想不要这个提示,怎么办? (/uis,msgpop,3这个行不通,我试过了)
(2)ansys中数组的定义为a*b*c a<10040 b<255 c<7 而我实际上要定义更大的有什么好办法,难道非要写到文件里再读进来?
Aq1:有两种办法1在ultility menu / menuctrls/ message controls下设 (2)将 /uis,msgpop,3 加入到start.ans文件中去以后每一次运行ansys时都会默认这个设置
--------------------------------------------------------------------------------
Q我现在做一个板壳的屈曲分析想要得到该结构在荷载作用下的极限承载力利用非线性屈曲分析请问极限承载力的荷载比例因子在什么命令下可以查到
A非线性屈曲分析与特征屈曲的不同之处是它一开始就挠屈变形随着荷载的增大变形日益增大我们认为结构最危险点的变形时间历程曲线中变形若从某点开始发生较大的转折与先前完全不一样就发生了屈曲此时的荷载就是屈曲荷载故不存在荷载比例因子即特征值问题
A要进行非线性屈曲,必须加到可以使结构发生屈曲的荷载,即荷载要比结构的实际承载力要 5
大,在进行非线性屈曲分析(极限承载力分析)时荷载逐步增加到破坏为止.计算时一定要输出每一步的计算结果,这样在post1中就会有一系列的荷载系数,而最后收敛的那个系数就是相对于极限承载力的系数,极限承载力=系数和总荷载的乘积
-------------------------------------------------------------------------------- Q如果材料一定为各向异性我用的材料为钢板该如何选用单元类型 A只要材料参数中有Ex,Ey,Ez单独输入项就能考虑各向异性
A波纹板的EXEY可以输入但是EZ该如何确定呢我自己练习了一下如果不输入EZ则计算就不能进行下去提示我材料的弹性模量要输入正值这是为什么
A如果是钢板一般认为是各向同性的厚板因轧制次数少晶体不密实存在着Z向厚度方向弱化对厚板一般地Ex=EyEz较小通过试验确定对薄板<16mm认为各向同性
A不管是什么板材料的本构关系是微小单元之间应力与应变的关系因此波纹板的几何构成是用几何模型->有限元模型计算的而不是用比较大块的板试验得出的Ex,Ey来输入你说的波纹板是不是压型钢板 A是压型钢板EXEY是通过拉伸试验确定的EX与EY不相等但是既然是考虑板壳为什么还要输入EZ况且EZ并没有试验数据 A波纹板的Ex=Ey=Ez
--------------------------------------------------------------------------------
Q怎样使梁单元把弯矩传给实体单元实体单元只有应力使它与相连的梁单元保持连续性真是件很难的事情还请各位大侠指教
A关键不是实体单元只有应力而是实体单元每个点只有三个平动自由度能不能模拟成将梁插入到实体中一
小段用两个或多个实体单元节点固定梁单元根据Sevant原理忽略掉两种不同类型单元连接处的计算结果 A关键问题在于传递弯矩是为了保持变形的连续性实体单元中尽管没有弯矩但同一截面各点处的应力实际已构成了弯矩但这弯矩是多大怎样传递给相邻的梁恐怕不是简单的插入嵌固可以解决的从各结构杂志来看这一问题似乎很普遍还请各位一起深入探讨
A上述方法实质上就是施加了人工约束单元相邻近点的连线几乎是梁的切线这样就能保证转角和位移连续弯矩相等并不意味着位移和一阶导数相等这是因为两种单元的位移函数相差太大
-------------------------------------------------------------------------------- Q弹性地基梁中的弹簧(2维)在ANSYS 中应采用那一个单元?又如何*作? A很久以前我用过弹簧单元好象是COMBIN2D你试试看
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可能实施加重力时,单位不一致
-------------------------------------------------------------------------------- 6
QANSYS里椭圆有没有命令可以直接画 A你可以利用ansys的局部坐标系,非常方便
gui:workplane>local coordinate systems>creat local cs>at 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或乘于倾斜度再加上施加的压力值如水荷载倾斜度等于水的比重 Qansys默认面荷载是垂直于作用面而曲面承受的风荷载为恒定的方向那怎么办 Q如果想把上一次的分析结果作为下一次分析的初始条件该怎么做 A使用LDREAD命令 首先需要注意下面两个问题
1每一个ANSYS的实体模型的面或体都要定义对应的单元类型编号材料属性编号实常数编号单元坐标系编号这些参数在整个分析过程中保持不变而这些编号对应的属性在各个步骤中是不同的 2网格划分要满足所有步骤的要求单元类型必须兼容 步骤
1创建实体模型2 创建多个物理环境
设定一个物理环境中的单元类型材料属性实常数坐标系等 将这些参数的编号赋给实体模型的面或体 施加基本物理载荷和边界条件 设定求解选项
选择一个标题使用PHYSICS, WRITE命令将物理环境存入文件中 3清楚当前的物理环境命令是PHYSICS, CLEAR 4重复第二步准备下一个物理环境
例如可以按照下面的方法将磁场的计算结果施加到流场上 /solu 进行求解
physics,read,magnetics 读入磁场分析的物理环境 solve finish
/solve
physics,read,fluids 7
LDREAD,FORCE,,,,2,,rmg 读入洛伦兹力 solve
更多信息请参考ansys的耦合场分析指南和命令手册
-------------------------------------------------------------------------------- Q有一个形状很不规则的刚架结构外表面是玻璃建模的时候只建了刚架风荷载怎样传递到刚架上去啊 A建模时就应该把玻璃和刚架同时建好要注意两种材料的衔接面(或LINE)和材料特性定义风载加上去后就会由玻璃传递给刚架或者简化风载成为PRESSURE直接加到刚架上 如果手上有实验数据加以比较会避免很多失误
-------------------------------------------------------------------------------- Q管路穿过厚壁,管路采用梁单元,厚壁采用体单元,在两者交接的地方怎么办?
A建模时可以考虑先建管路再建厚壁体时厚壁体由外部的平面与管路外表曲面组成这样管路和厚壁SHARE同一个曲面荷载作用下力也可以传递材料属性分配的时候先PICK UP厚壁单元定义为所选取得体单元再PICK UP管路定义为梁单元其实很多梁单元也属于体单元没有什么冲突除非管路属于SHELL或面单元时要后PICK UP此面单元定义为所要的属性
至为重要的是在MESH这两个单元面体衔接的时候要注意衔接面体对应的ELEMENTS 和NODES的COMPATIBILITY否则将会出现大的错误REFINE时也不可忽视 QA1请详细解释一下
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
8
A使用SWEEP或CONCATENATE命令重新MESH桥板
可以在HELP文件中找到这两种命令的用法十分有用将桥面MESH成为HEX类的ELEMENT这样寻找加载位置会方便很多
A建议你将整个桥面划分为比较规则的四边形网格(六面体)当然网格的大小肯定与 加强筋的厚度等有关这样应该能够解决你所说的问题
-------------------------------------------------------------------------------- Q一个梁上有一滑块计重力和摩擦滑块在梁上滑动求梁的响应梁为1m*0.1*0.1,滑块为0.2*0.1*0.1ex=2e11,nuxy=0.3,dens(1)=7.8e3,dens(2)=2.7e3
A1命令CONCATENATE如果你要把一个面或者体MESH成为一个四方或HEX单元这个命令十分有用你可以CONCATENATE相邻的线或者面把所要MESH的面变成一个四边形或体变成一个六面体(例如5边形CONCATENATE两条相邻的线于是ANSYS认为此时的面是四边形你就可以进行3OR4SIDES的MESH)所以无论你的面由多少边 都可以进行MAPPED MESH
2SWEEP和上边的命令一样都是我所喜欢的命令只要你的体已经建成SWEEP好一个面以后整个体MESH可以一扫而过HEX的ELEMENT加载在什么地方或什么样的荷载都十分方便不需要在建模的时候就考虑加载位置 3至于选用单元具体问题具体分析但是衔接一定要注意特别是衔接面的ELEMENT要对应好以免出现错误 4自己的一点心得MESH十分重要别在这上面偷懒我曾经用FREE MESH与MAPPED MESH得到的值分别为19.2和25.9而理论值是26.3所以请小心
--------------------------------------------------------------------------------
Q我做的是一个三层的压力容器每层厚都是2mm请问我可不可以只定义一种单元类型没有shell91,99) 另外这三种材料的材料属性用什么命令分别定义给这三层 另外用轴对称分析的话应该怎么做
A划分网格前,用attribute>define可定义不同的材料属性
-------------------------------------------------------------------------------- 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我用一个静定结构进行分析结果输出后和实际的受力方向刚好相反受拉的成了受压但反力的方向是对的 -------------------------------------------------------------------------------- 9
Q有一个圆柱桶内装有70%的水如何将桶和水的重量考虑进去
Aload>gravity.. A把水另设为体单元
-------------------------------------------------------------------------------- Q四边固定的板用体单元有一条线用梁单元从板的任意一点穿过分析其固有频率 如果不将线体联为一体分析时各算各的问如何将线体联为一体 glue无法将线\\体连为一体
glue只能将同样的元素连为一体(如体和体,面和面,线和线)
-------------------------------------------------------------------------------- Q如何给模型加3g的冲击载荷
AIn short time(for example:0.0001 second),用一小time 步加载 A 3g指的是加速度冲击载荷例如在地震时的地基响应 例如
在固定端上用弹簧悬挂一根杆在固定端加一个3g的冲击载荷g指加速度求杆的响应 1之所以用弹簧悬挂是为了减振这个问题主要想知道如何加这种载荷 2类似于地震时的地基响应问题 3我在ANSYS的GUI中没有找到这个功能 4BTW在ANSYS中如何加弹簧元
--------------------------------------------------------------------------------
Q我想知道用ANSYS分析岩土洞室开挖效果好不好
我刚开始学ANSYS不知道ANSYS用于分析有裂隙岩体和松散土体方面效果好不好其精度如何望各位朋友能给予指教
Adyna做隧道开挖比较擅长而且看起来非常形象
A用ansys里的ls-dynadyna是绑在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检查单元prep7>meshing>check mesh若有不合格的单元用别的办法划单元还可以设limits有时候其实这种警告对结果没有多大的影响完全根据你的经验 A如果可以进行下去可忽略该warning
一般来说温度场分析对形状测试更严格一些一般不会对问题产生太大影响可通过剖分更密的网格来解决 -------------------------------------------------------------------------------- Q如何得到某一温度值的距离-时间关系曲线
在分析平板移动热源温度场时希望取得时间历程中熔化温度值在各时刻对应的平板上出现的位置或节点号即某一温度值的距离-时间关系曲线或者如何得到某一温度值在给定时刻出现的位置 A我想你可以把你需要的那些温度值的点选出来设置成路径在用路径来画图 在path operation下
A谢谢你提供的建议我不知道这些温度值点出现在什么位置因此路径方法不可取也许只能采用set设置时刻,sort排序各节点温度prnsol找出温度值对应的节点此法较为繁 10 琐
-------------------------------------------------------------------------------- Q5.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 33824. 3020 75691. 0.31964E+06 0.0000 18182. -9597.7 2285.4 3020 -0.24616E+06 -44015. 0.0000 -82203. 1491.2 33824.
我想弄明白为什么有的点(225)输出2行,有的(3020)输出4行.如果平均化了,那应该只有2行(Top & Bottom)?
A由于你的模型中有多种材料在材料交界处的节点结果跟交界的组数有关该节点共几组材料则有几种结果节点结果是由相邻单元的高斯点结果外推出来然后平均的在PowerGraphics ON 状态下ANSYS在材料交界处的节点自动按材料区分同一组材料的单元外推点参与平均不同材料的不参与平均所以节点交界处有几种材料该节点就有几组结果
-------------------------------------------------------------------------------- Q如何根据已创建的关键点划分网格使节点恰好在关键点上 A把关键点变成硬点即可或在关键点处生成硬点
--------------------------------------------------------------------------------
Q我在做一个saint vennent问题模型为平面薄板THK=1左边受x方向约束右边中点受一集中力作用现在全场的应力做出来了我想做沿力的作用线上的应力曲线请问如何做
ADEFINE PATH在POSTPROC中PATH OPERATION 命令定义PATH用NODE方便一点然后MAP你要所需的值ANSYS可以做出曲线
--------------------------------------------------------------------------------
Q如何获得应变最值出现的接点或单元的坐标
A用后处理里面的QUERY RESULTS 找到该接点或单元即可
A这样可以查到最值出现的单元或节点的序号但无法直接求出坐标还需另外查询不知又没有捷径此外我做的是一个圆环的应变分析想求出圆环内表面回转方向的应变怎11 么求
圆环只是整个结构的一部分因此转换为柱坐标系好像有问题不知如何解决
A用QUERY命令点住最值节点在拾取菜单Query subgrid result里就列出了坐标Global x= y= z=在后处理时把关心的圆环部分定义为局部柱坐标然后在列出结果即可
A问题是我不知道节点或单元的具体位置只知道其序号如何应用query命令 A用顶上菜单栏LIST>Picked Entities+,在ANSYS INPUT 对话框中输入序号即可 A可以用*GET命令请查COMMAND REFRENCE
-------------------------------------------------------------------------------- Q我在分析应变时需要知道最值区域究竟是处在弹性区还是已经进入塑性区有没有什么好的办法 A在后处理中 PLOT> NODAL SOLU...>STRAIN-PLASTIC即可
--------------------------------------------------------------------------------
Qansys 5.7 分析网架单元用link8 link10,beam结果用list line element可查但为何不能显示应力图形长度单位用mm,力用N,GRAVITY中自重方向单位是用9.8吗如何让ansys算自重 AANSYS中线单元的应力属于导出结果,需要定义单元表(ETABLE),然后绘制
--------------------------------------------------------------------------------
Q我正在做基坑方面的论文但是碰到了问题如下
1 在ansys5.7中如何在第二荷载步将第一荷载步计算的位移置零也就是位移恢复到第一荷载步未计算之前
2 什么单元既能使用d-p模型又能否考虑水压力的变化 3 在荷载步之间改变边界条件为何总是不收敛
A土的非线性用dp材料如想用其他需自己开发比较麻烦水压力可按输入荷载施加锚杆的模拟只能用杆单元按ansys模拟的不是很好不过ansys的dyna提供了较好的模拟
-------------------------------------------------------------------------------- Qansys的prucker-drager模型计算土石坝如何和邓肯模型的计算结果有何不同怎么算
A二者的比较不是ansys本身的问题而是两种计算方法本身的差别辅导老师也不知道采用dp材料只用在材料属性定义即可可详见非线性分析指南或帮助
-------------------------------------------------------------------------------- Q当对钢筋混凝土进行分析时候假如定义如下几项 tb,concr,1
tbdata,1,0.8,1.0,3.8,43.0 tb,dp,1
tbdata,1,2.9,32,28.0 (这里不讨论其中的具体数值及其含义)
按照指南培训手册帮助及理论手册中的说明用TB定义的选项几乎绝大多数同时就确定了屈服准则硬化准则及应力应变关系(但tb,hill除外)单独使用TB,MISO,1,就说明材料使用多线性的应力应变关系Von Mises屈服准则等向强化准则而单独用TB,DP,1,也就确定了理想弹塑性的应力应变关系D-P屈服准则(DP
材料不存在硬化问题)单独形如TB,CONCR,1,则确定了W-W屈服准则其它如理论手册中的那样(SOLID65部分)
而当如上述几条命令那样对一种材料同时定义了CONCR和DP则屈服准则强化准则及应力应变关系如何呢而当定义CONCR和MISO又如何呢 A不能这样定义只能定义一种非线性材料
-------------------------------------------------------------------------------- 12
Q我用坐标平面去切另一个面时为什么切不开呢
说unable to complete boolean operations.修改模型或改变tolerance值 A使用于切割的面成工作面然后DIVIDE
--------------------------------------------------------------------------------
Q我自己生成了大量的命令放在一个文件里面在INPUT窗口可以读入但是一次只能读入十几行太慢有没有更好的方法
AFile>Read Input from... A看一下/input命令
--------------------------------------------------------------------------------
Q在一个体中如何对不同的部分选择不同的划分模式比如说一个部分需要等间距划分而另一段则不等间距 A在发出ESIZE或LSIZE命令时对其中的域SPACE设置非0或1的数值
--------------------------------------------------------------------------------
Q在模拟路基填土时怎么实现分步加载我用time和nsubst命令结果他只给我算第一个有载加载我输入的是0加载-有载加载-0加载-然后用lssolve1n求解可是只有n=2时能算出来然后就不接着往下算了不知道时由于方程太多计算不来还是什么原因如果时方程太多有没有解决办法使他算下去 A你用的是荷载步文件最好别用这个求解方法一步步地算记得每步计算前重新生死单元
Q在ansys中main meun>preprocessor>sections>common sectns..当中对T型材计算所得到的Torsion Constant的值是不是指极惯性矩polar inertia但我从字面意思看好象是指扭转常数什么的问他到底是指什么
A即材料力学中的扭转常数Ip,ANSYS中缺省是Iyy+Izz,如用188/189则可以看到其是怎么计算出来的 A应该是材料力学中的关于z轴的惯性矩Izz对应ansys命令section显示截面图形的Iyy --------------------------------------------------------------------------------
Q有人用ansys做地基梁分析吗请问分析是采用的模型是什么 ansys中单元应该用什么是不是应该算做接触问题
A结构用梁模拟结构与地基的作用用弹簧单元模拟
-------------------------------------------------------------------------------- Q建立三维几何模型时能否用一个二维平面绕其对称轴旋转来生成如果可以具体怎么做 A2-D模型建好后用命令Pmodeling>operate>extrude>areas>about axis就可以了 A我在二维模型上已经划分了网格为什么生成三维体后看不到应该怎么划分{BANNED}格呢 --------------------------------------------------------------------------------
Q我要用ANSYS计算一个粘弹性材料夹心的夹层板的动力计算请问关于粘弹性材料属性的怎么定义应该选取什么样的单元对此夹层板进行建模
A不知道你看没看 ANSYS帮助中有一个类似问题的例子VM200这个例子如果没有的话我想他会对你有帮助的
--------------------------------------------------------------------------------
QA: 一直以来,ANSYS被人们认为的缺陷之一是没有UNDO命令最近了解到,可以通过一个未公开的命令实现UNDO.方法如下: 1 首先激活UNDO 13 /UNDO,ON
2 然后在任何时候需要进行UNDO时,发 /UNDO
这一命令在5.7版就有了,不过激活UNDO后,ANSYS的处理速度会有所下降.
-------------------------------------------------------------------------------- Q怎样在顺着梁的方向显示梁上的应力云图 A建立局部坐标系然后在局部坐标系下显示应力
--------------------------------------------------------------------------------
Q有一受弯扭的曲杆,用3Dsolid45单元(也不知道这样行不行?),杆件中的某一段受到轴承套的约束,仅能转动,不知如何施加边界位移约束?
A试试在柱坐标系中仅施加x向位移约束应该可解决该问题
--------------------------------------------------------------------------------
Q: 如何施加非法向的面力比如作用在斜杆上的竖向均布荷载 A在任意面施加任意方向任意变化的压力
在某些特殊的应用场合可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷当然这在一定程度上可以通过ANSYS表面效应单元实现如果利用ANSYS的参数化设计语言也可以非常完美地实现此功能下面通过一个小例子描述此方法
!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来 !!! finish /prep7 et,500,shell63 press=100e6 amesh,all esla,s nsla,s,1
! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将 ! 所有节点坐标系旋转到局部坐标系下. *get,enmax,elem,,num,max dofsel,s,fx,fy,fz
fcum,add !!!将力的施加方式设置为\"累加\而不是缺省的\"替代\"
*do,i,1,enmax *if,esel,eq,1,then
*get,ae,elem,i,area !此命令用单元真实面积如用投影面积请用下几条命令 ! *get,ae,elem,i,aproj,x !此命令用单元X投影面积如用真实面积请用上一条命令 ! *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文件. 你可以通过打开这个文件获得信息
我想你所说的批作业恐怕是指在GUI界面下直接输入命令它的信息都会保存在jobname.log这个文件中它也是一个ascii文件
至于你说的保存图形文件可以通过
plotctrls-->hard copy--->to file来保存 也可以在画图前
plotctrls-->redirect plots---->to .... file 19
你试一试吧
Q不过我不知道何谓ascii文件如何打开打开以后是否可以接着上次的进行分析还有仅用GUI模式不输入命令是否能完成一个三维结构的应力分析
Aascii文件通俗的说法就是文本文件至于你说的接着分析,我就有些不太明白了是因为某些原因ansys中断退出后的接着分析呢还是计算完了还想接着再来一遍这两种的做法可不一样
不过我建议你把ansys的帮助文件好好看一遍尤其是Basic analysis procedures guide绝大部分的问题都可以在里面找着答案
--------------------------------------------------------------------------------
Q我用ansys做了个项目需要在任意方向施加一个位移约束我定义了一个局部坐标系统并且把这个局部坐标系统设置为当前激活坐标系统然后施加位移约束可是不管我怎么加还是加在全局坐标系统上不知道何原因
A你可以把你要加约束地方的节点坐标旋转到你需要的约束方向然后再约束 preprocessor--->move/modify--->by angles
--------------------------------------------------------------------------------
Q现在先用block命令画了一块板然后我想在板的侧面上的一条线是侧面中间的一条上加线载荷以模拟预应力问题是在mesh以后想用
lsel,s,loc,x,-1000 lsel,a,loc,x,1000 lsel,r,loc,z,-25 sf,all,pres,0.03
加线载荷怎么老说没找到线怎么解决这个问题 A原因很简单因为你就根本没有先定义线
BLOCK所定义的只有一个体积六个面以及8条边线你可以用LIST看一下如果你想在侧面中线上加载的话只好如此
方法一用KEYPOINTS的方式建立模型定义是别忘了定义侧面中线的两个点在侧面生成两个面再以面建体积然后就可以在线上加荷载了
方法二正若DFENG所言在那条线上的NODE上加载但是要注意的是两端NODE的荷载只是中间NODES荷载的一半要计算清楚
A居我所知ANSYS不能在面上加线载荷所谓的线载其实指面载我建议在线上取点后均匀加上力 -------------------------------------------------------------------------------- Qopt运算器它的位置应该在main menu里面请问如何安装才能使用该运算器 A/opt或main menu>design opt,详细可参见ANSYS高级分析技术
Q做一个锥型工字梁的静力分析该梁由9个关键点分成8段还是一个整体在每个关键点处的高度缘条腹板的长宽均不同并非线性变化这么一个梁该具体如何处理,如果用BEAM44的话根据REAL属性用ESHAPE命令显示的梁是矩形梁,一个实际工程中梁的REAL 属性怎么得出,是否是对端点用 PLOT SECTION 方法?对于两端截面不等的I字梁,如用BEAM188(189)的话,只能与一个SECTION关联,能否对端点设置SECTION ?怎
么做?
A1 采用beam44截面变化时不能使用SECTION特性而使用SECTION关联就不能是渐变的截面 2 模拟BEAM44截面变化时你可以输入实际常数以模拟截面变化这就要你输入单元两端的特性数据(area1,iy1,iz1,h1...area2,iy2,iz2,h2...) 20
3 在非截面关联的情况下你所PLOT的单元截面形状当然只有矩形了
4 渐变的数据和参考VM34.dat中的BEAM44部分即是该例子中的单元两端的截面是渐变的 A不变截面的关联 /prep7 k,1 k,2,100 k,3,0,1 l,1,2 et,1,beam44 mp,ex,1,30e6 mp,prxy,1,0.3 sectype,1,beam,I secoffset,cent secdata,30,20,40,8,8,6
latt,1,,1,,3,,1 lesize,all,,,20,,1,,,1 lmesh,1 finish /solu /eshape,1
-------------------------------------------------------------------------------- Q我在ANSYS中分析一个受瞬态冲击荷载的悬臂梁
分两个时间步第一个时间步0.05秒受荷载一个子步第二个时间步4.95秒不受荷载99个子步 我想得到弹塑性情况下的梁末端的振动曲线但是梁的固端一旦进入塑性ANSYS的分析就中止输入结果为零实际上如果梁的固端只是部分进入塑性梁还是应该振动的请问这个问题应该如何解决 A塑性在完全瞬态动力学中是可以的而在其它动力学分析中是被忽略的仅保持最初的非线性特性 --------------------------------------------------------------------------------
Q我在计算一个渡槽的项目屈曲分析中结果输出只有几个屈曲模态的值不知道这几个值是什么意思还有屈曲临界应力怎么用这几个模态值来算
A每个屈曲模态都有一个fact系数则该系数与施加的计算荷载乘积即为该模态的屈曲荷载 --------------------------------------------------------------------------------
Q请教一个实体加载分网的问题实际上的问题是我在体上定义了关键点作加载点但是用meshtool分网时没有分到这几个点所以加载失败
A假如你必须在点上加载可
将你的那几个关键点第一为硬点但这样分网要注意以及此点存在应力集中问题或应力奇异 建议你用面加载模拟可能好些
-------------------------------------------------------------------------------- Qansys中对于温度应力的问题是怎么算的呢一个最简单的问题一根杆件两段固定在 21
温度变化下的应力怎么算呢
A这是一个耦合场计算问题单纯的热计算和结构运算都解决不了问题如果我没有记错的话你可以考虑plane13这个耦合场单元类型
-------------------------------------------------------------------------------- QSFA的用法哪位知道其用法最好有一个例子 ASFA的所有用法全在帮助文件里 SFA, AREA, LKEY, Lab, VALUE, VALUE2
1area是你要加载荷的面可以直接指定面的编号也可以先用asel命令选取然后用all参数即可 2lkey则是同element有关的载荷关键号不同的element会有不同的意义
3lab是用来指定载荷的意义比如说要加载压力就用PRES,要加载热流量就用HFLUX要加载对流就用对流热交换就用CONV.
4value就很简单了就是知载荷的大小罗针对某些载荷和某些单元而言一个value不足以表达得给出另外一个值来这时就用到了value2例如在加载热交换载荷CONV时往往还得给出环境温度或者是流体温度),
这个温度就用value2给出 SFA,3,CONV,14235,25 这个命令的意思就是
在第3个面上加上一个热交换载荷交换系数是14235环境温度是25度
-------------------------------------------------------------------------------- Q关于正交异性板有谁知道可以用什么单元模拟两个方向刚度不一样的板如何在ansys中实现 Amatetial props->orthotropic
A据我了解ansys的材料属性基本上都可以在xyz三个方向分别定义你仔细找一找应该没有问题 -------------------------------------------------------------------------------- Q我这里有两个问题向请教一下
1对于平面刚加结构来说比如说有一根杆是斜的我想对它施加x或是y方向的均部载荷而不是垂直于杆的均部载荷
2对于一个面如何施加剪应力
A1 梁元施加的荷载是针对单元的各个方向的例如垂直切向等而施加整体坐标下的X或Y方向的荷载只好分解后施加
2 查看单元的面力规定方向数字即可施加可查HELP
--------------------------------------------------------------------------------Q有一个二维模型大部分用原型法(solid modeling)建的模(存在keypoint,line,area等元素),模型上某个小部分用直接法(direct generation)建的模(即其中没有keypoint,line,area,只有node和element),请问我怎样才能overlap这两个部分
A曾经有人讲过这方面的事情(E文)但我忘记了至少有一点是肯定的即用原形法建立模型划分好单元后可以接着用直接法继续定义结点和单元是一体的FEA模型不需要任何*作
A如果实体模型已经剖分然后再在这个基础上定义节点element当然是没有问题的不过直接把没有剖分的实体模型和节点,element做布尔运算恐怕很不妥吧我个人认为布尔运算不能使用与节点和单元 --------------------------------------------------------------------------------Q怎么对三维实体加绕轴转动的约束就是只能绕轴转动 22
A如果在柱坐标下让所有节点的UX,UZ自由度为0是不是就可以了不过这里的自由度的方向是节点坐标系的方向所以你得先把节点得坐标系定义为柱坐标系
A你也可以试一试ROTXROTYROTZ这三个自由度说不定更方便
Aansys中还有一种惯性载荷其中可以定义角速度和角加速度命令是CGOMGA和DCGOMG你可以仔细琢磨琢磨
--------------------------------------------------------------------------------Q如果一个线既不是水平也不是垂直而是与水平线相交成一定的角度如何加与之平行的均布力 A查看单元图形不同受力的方向定义为不同数字选择切向就可以了
-------------------------------------------------------------------------------- Q怎么画一个圆柱面
A1建立圆柱体仅删除体留下的是圆柱面
2建立圆弧线沿一定的路径拖拉也可以建立圆柱面
--------------------------------------------------------------------------------Q请问ansys中的.full文件中如何看刚度矩阵和质量矩阵
Ayou can use the following steps to extract structural matrices: 1. run a modal analysis to get a file.full 2. generate a temp folder
3. copy binlib.dll, rdwrt.bat, rdfull.f, binlib.lib to temp folder (these files can be found in your ansys system)
4. rename rdfull.f to userprog.f
5. at command prompt, run: rdrwrt userprog 6. at command prompt, execute: userprog.exe 7. you can get MASS.MATRIX and STIFFNESS.MATRIX.
Ajobname.full是二进制文件它是以标准文件头作为第一条记录以后的若干条记录都有固定的模式在Interface Guide中对jobname.full的格式做了详细的介绍而且ANSYS提供了一系列的调用方法能访问各二进制文件不是很难的事
建议你仔细研读一下bintst.f好象还有几个重要的子程序是 binrd.f binwrt.f
--------------------------------------------------------------------------------Q我现在正在学习用ANSYS进行三维流场分析但在计算完毕后我不知道怎么看截面流场矢量图例如对称面上或者垂直于空气气流流向的截面上 A1将模型采用正视图或侧视
2用SELECT ELEMENT 菜单选box or circle,先用选取再用不选取就可得想看的截面
--------------------------------------------------------------------------------Q如何在屏幕上不显示ANSYS标记 A/plopts,vers,0
-------------------------------------------------------------------------------- Q请问如何在后处理图形中加坐标轴,以及注释
A加注释可采用utility menu>plotctrl>annotate..然后在屏上点填写即可假如是纯粹的注释性建议你形成bmp文件然后到画笔中爱干啥干啥
Q如何形成bmp文件,我自己也试过用utility menu>plotctrl>annotate>create 2-D annotation,但是写出来的字是倾斜的,该如何处理 23
AUtility menu>plotctrls>hardcopy>to file
在图上直接注解时先点一点然后拉出直线则字为正另外可以设置字体的大小字体颜色等. --------------------------------------------------------------------------------
Q已将预应力筋按LINK8单元划分30等分各单元的温降也已计算完毕现在想为这30个单元以读入数据文件方法加温可是总是只加第一个温度到所有单元我已经分别用了ARRAY和TABLE 都不行请指教 A你看看这个自己修改搞定如何 /prep7 et,1,link8 k k,,100
k,,200,50 l,1,2 l,2,3
lesize,all, , ,10 lmesh,all 以上是建模试验 -------------
!向线上施加温度的初步例题 !从数据文件读入数据并施加 !从线的一端到另一端
!选择线可采用其它方式数组大小可另定 *dim,tem,,20 !定义数组
*vread,tem(1),wd,txt !从文件wd.txt读入 (20f5.1) csys,0 lineno=2 allsel,all
lsel,s,,,lineno !选择线
nsll,r,1 !选择其中的结点
*get,k1,line,lineno,kp,1 !得到KP1 *get,ncount,node,0,count !得到结点总数 *get,nmax,node,0,num,max !得到最大结点号 *get,nmin,node,0,num,min !得到最小结点号
nnum=node(kx(k1),ky(k1),kz(k1)) !得到一端的结点号 *do,loopno,1,ncount,1 !循环结点总数 bf,nnum,temp,tem(loopno) !施加温度值 *if,nnum,eq,nmin,then !得到紧邻下一结点号 nnum=ndnext(nmin)+1 *elseif,nnum,eq,nmax,then nnum=ndnext(nmin) *else
nnum=ndnext(nnum) 24 *endif
*enddo !循环结束标志 allsel,all
bflis,all
!---------------------- 以下是wd.txt文件内容另存为文件
10.0 20.0 30.0 40.0 50.0 60.0 70.0 71.0 72.0 73.0 74.0 75.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 71.0 72.0 73.0 74.0 75.0
--------------------------------------------------------------------------------Q如何加密求得的等值线,我现在看到的太稀疏了 A设为等值线方式/device,vector,1 定义标记疏密/clabel,1,no1 设根数:/cont,1,no2
其中no1,no2是自己定义的试试疏密满意为止.
--------------------------------------------------------------------------------QA!/nopr !向线上施加荷载的宏 !需要四个参数 !lineno:线号 !p1:一端的荷载值 !p2:另端的荷载值 !prekey:荷载方向号
!------------------ !从对话框输入上述参数
multipro,'start',4 !对话框输入数值 *cset,1,3,lineno,'the line number' *cset,4,6,p1,'the pressure at one end' *cset,7,9,p2,'the pressure at other end' *cset,10,12,prekey,'the pressure key',1 multipro,'end' csys
!确定选择集(该线及其单元)并得到线长 lsel,s,,,lineno !选择线 esll,r !选择此线中的单元
*get,linecd,line,lineno,leng !得到线长 *get,qdh,line,lineno,kp,1 !得到KP1 !处理荷载(对选择集而言)
*get,ecount,elem,0,count !得到单元总数 *get,enum,elem,0,num,min !得到单元的最小序号 *do,loopno,1,ECOUNT,1 !循环单元总数
*get,i,elem,ENUM,node,1 !得到当前单元的i结点号 *get,j,elem,ENUM,node,2 !得到当前单元的j结点号 dx=nx-kx(qdh) !求I结点到KP1的距离 dy=ny-ky(qdh) 25
dz=nz-kz(qdh)
inloc=sqrt(dx*dx+dy*dy+dz*dz)
dx=nx(j)-kx(qdh) !求J结点到KP1的距离 dy=ny(j)-ky(qdh) dz=nz(j)-kz(qdh)
jnloc=sqrt(dx*dx+dy*dy+dz*dz)
ipress=(p2-p1)/linecd*inloc+p1 !计算I结点荷载 jpress=(p2-p1)/linecd*jnloc+p1 !计算j结点荷载 sfbeam,enum,prekey,pres,ipress,jpress !施加该单元荷载 *get,ENUM,elem,ENUM,nxth !从选择集中得到下一个单元号 *enddo !循环结束标志 ALLSEL,ALL !宏结束
-------------------------------------------------------------------------------- Q作高速旋转体的力分析除了加角速度载荷外还要加节点约束吗对于空心圆柱如何加
Q建立一个坝后背管的三维模型管道是弯的三维实体元8节点20节点都试过钢管内衬采用壳元48节点也都试过就是加管道的内水压力时老是加不上显示是双箭头 Q如何创建一对齿轮啮合的三维模型
Q在用ANSYS求解电磁轴承磁场磁力分析时不知如何建模和分析
Q请问那位前辈知道在钢筋混凝土中如何考虑钢筋的影响特别是密集配筋混凝土处于三向受压的情况 Q在隔振结构的动力学分析中两个平板之间安装有一个隔振器且已知隔振器3个方向的刚度和阻尼在ANSYS5.7中是否有可以同时表达3个方向的刚度和阻尼的单元一般而言一个单元只能设置一个方向的刚度和阻尼参数或者有其他方法可以同时考虑隔振器3个方向的刚度和阻尼的影响
Q质量弹簧单自由度系统在为正铉波的地基上振动激励为正铉位移怎么把位移约束加到系统中去使之达到减振效果
Q1弹性半无限边界的有限元简化模型2路面上作用车轮荷载其等代荷载是圆面或双圆面在ANSYS建型时如何体现3如何知道单元中IJKL等节点的编号顺序如何知道单元坐标系
Q我建了一个模型是箱形结构中间有隔板隔板间距不一我想将隔板与箱形glue起来但是只能glue原始隔板其他copy过来的隔板无法glue,我用overlap将其他隔板与箱形连接起来mesh完后加力加约束后无法求解我用的是shell单元
Q有哪位知道ANSYS中*MFOURI傅立叶变换的用法
Qansys中adams命令是用来生成.mnf文件的请问有谁知道其具体用法
Q请问能否在实体模型中的一条圆弧线上加载沿线分布的均布载荷载荷方向垂直于圆弧所在面
Q一个弹塑性问题想知道塑性变形是否要用瞬态分析先加载再卸载 Q请问怎么用命令实现File/List/Log File?
Q平板加热问题用热结构耦合算出结果后需要得到某一时刻平板上的等位移线我只知道有云图显示等后处理我想要得到一条一条的线不知怎么处理
Q我研究的是瞬态热分析在加载边界条件时有两个边界条件与模具的表面温度有关 26
即当模具的温度>104度时边界条件是热流密度当<100度时边界条件是对流换热系数模具温度是变化的我一直不知道如何区分加载步
Q我是一名结构工程专业的研究生在使用ANSYS5.7 版本中遇到了困难万望得到您的帮助现将具体问题描述如下模型是用来分析钢框架中的H型钢梁柱刚性节点H钢梁翼缘通过对接焊缝和H钢柱翼缘相连H钢梁腹板通过高强螺栓预拉力连到预先焊在柱翼缘的剪切板上模型采用SOLID45单元在梁腹板接触区定义了TARGET170和CONTA174接触单元模型共分24000个单元10000个节点约束条件为柱两端固定梁端侧向约束在梁端采用循环位移加载将梁端截面处的节点耦合在主节点加载PCG求解器大应变全牛顿-拉弗逊迭代法单元刚度矩阵是基于接触状态进行修正的在第三个循环求解时出现下列错误提示
\"PCG Solver error level 1. Possibly the model is unconstrained or additional iterations may be needed. Try running with a multiplier, MULT>1 in EQSLV command (3>MULT>1)\". 按下错误信息条的确定按钮以后程序二分荷载子步并进行继续求解但在第四个循环荷载步时程序中断ANSYS自动退出呈现了WINDOWS窗口在File.err文件中的错误信息为 \"****FATAL****
Random write error on unit =12 possible full disk I/O status error =28 ****************
The above error is non-recoverable by ANSYS ANSYS run terminated by the indicated error ***************\"
按照ANSYS 使用说明将数据库内存改小后仍出现上述错误
注本计算机硬盘为30 GANSYS 工作内存为360 MB数据库内存设为 80 MB Q请教高手用ansys如何进行预应力分析
Q请教两个问题1无限边界单元在所分析的区域远小于整个地域时要采用无限边界单元如对岩土路面等进行分析但ANSYS没有结构类无限单元
2应力强度因子理论上只有平面问题的计算公式ANSYS中对模型取积分路径实质上也是简化为平面问题处理若要真正求得空间开裂问题的应力强度因子该如何办 Q复杂结构网格化后能否的到任意截面的弯距?
Q在用壳单元SHELL63划分模型的时候单元的两面分别是紫色和蓝色未用颜色区分不同面, 即(PlotCtrls>Numbering...>Area Number Off请问 1这是为什么
2有时候划分完单元在模型的同一表面即同是模型内表面或外表面紫色和蓝色都有这是为什么会不会对结果有影响
因为我个人认为紫色和蓝色是用于区分壳单元的内外表面所以在同一面颜色应当一致不知对不对 Q问题是分析某一国防工程隧道的应力应变具体情况如下有一地下隧道在地下200m处隧道截面为直径4m的圆隧道上覆盖的岩石为新鲜混合岩岩性均匀无较大的断层通过具体参数为容重26.3KN/M3侧压力系数0.8围岩内摩擦系数为1.3凝聚力1.3MP隧道内承受内压1.9MP隧道内再做80CM后的钢筋混凝土衬砌在
衬砌内再做44MM厚的钢衬砌
现在要分析1围岩在毛洞情况下的应力应变2分析钢筋混凝土衬砌和钢衬砌在内压作用 27
下的应力应变
我用ANSYS对该问题作了初步分析具体过程如下 1计算岩体在自重作用下的应力
2把第一步计算出来的垂直方向的应力乘以侧压力系数得到水平方向的应力作为初始应力 3计算隧洞开挖形成的应力场具体做法为读入初始应力加节点等效荷栽杀死开挖掉的单元 4激活钢筋混凝土衬砌和钢衬砌并加内压 5完成计算
做完弹性分析我初步认为结果可靠继续做弹塑性分析在第三步时计算不收敛计算无法完成
现请问问题是出在初始应力场计算上还是岩石凝聚力太小并请问你在做初始应力场计算上是如何实现的 Qansys里有没有跟随力 Qansys如何输入地震波
Q在桥梁的抗震分析中建模时应如何处理橡胶支座呢
Q我建了一个实体上面再加一个弹簧单元请问弹簧单元与实体接触的地方怎么处理要考虑耦合吗 Q我在一个结构分析的课题中3d需要求出三维模型的某一条边上的各节点的应力分布曲线以及该曲线上的最值及相应节点坐标不知如何得到
Q在ansys中如何加上Biot固结理论是修改刚度距阵还是其它 Q在ansys中如何进行圆周阵列
Q请教OUTRES中如何控制CM选项我尝试了多次可每次的输出都不变 Q我用实体元建一个最简单的框架结构四根柱四根梁加一层板
1板梁作为block都已create好了怎样能够简单方便地在图形界面上把板作为block建起来不需输坐标值
2我已经求出了荷载下结构的反应能否得到梁中某一截面上的结果包括应力及内力
Q扶壁后的土体沉降及应变如何处理土体作为松散体选择哪种单元和收敛准则计算比较理想 Qansys中怎么模拟风荷载
Q请问除在谐响应分析中能加简谐载荷以其它分析中能不能加简谐载荷如何加 Q在ansys中如何仿真钢筋和混凝土的粘结滑移
Q有若干关键点我想用定义多义线的命令联闭合曲线但是不能闭合要求无接口并保留关键点请问如何处理 Q求解完一个多载荷步问题之后该选退出QUIT的哪一个选项 Q谁知道ansys中网面计算用什么单元
Q分析混凝土徐变的问题但是选用的徐变方程里的几个参数c1c2,c3,c4等不知道怎么定 Q混凝土带预制裂缝,裂纹扩展如何模拟
Q1松散土体是否可以用SOLID65单元模拟或有更好的方法 2土和混凝土接触该如何处理摩擦角等
Q本人正在做砌体裂缝开展的仿真的研究,使用ansys5.7,现遇到如下问题: 1. solid65带有混凝土的材
料特性,能不能赋予其它的材料特性?如何处理? 2. 其它的单元如solid95等,可不可以像solid65一样内涵混凝土材料?
Q斜拉桥的极限承载力中的索的弹性模量和混凝土开裂如何处理 28
Q如何得到3d模型上某条线上的最大应力节点坐标并将坐标值赋给自定义参数
Q在定义了钢材的非线性模型后对构件循环加载对构件施加很大的位移下程序仍能计算但实际构件早已破坏ANSYS不能确定构件在何时破坏即程序加载的终止点
Q桩锚支护结构中的锚杆怎样考虑水压力怎么计算土的非线性怎样考虑等等
Q请问疲劳分析中的SN曲线是什么单位fl命令中的应力集中系数如何定后处理结果如何看结果说明了什么后处理结果是最重要的问题
Qansys的子结构分析中的超单元可以包含裂纹和进行裂纹J积分计算否 Q如何在批处理模式下将某一个变形动画保存到文件中供以后参考
--------------------------------------------------------------------------------Q建立实体模型后计算完毕要配筋时有无方法直接将应力转换为内力计算配筋和裂缝我知道梁单元有输出那实体呢 A梁板才有MQ本来就是简化的概念在ANsys中用beamshell单元都可直接输出内力用所谓\"sequence number\"solid没有
-------------------------------------------------------------------------------- QANSYS里面好像只有STRESS 结果能不能察看内力值?
A是可以查看诸如轴力弯矩等结果的你要详细看一下单元的输出数据在帮助中 A当然可以了每种单元都有关于输出数据的定义
例如beam4想知道y方向的弯矩查到与之对应的是SMISC2SMISC8在general postprocelement
tabledefine tableadd
在弹出的对话框中选择by sequence num定义后在plot或者list element table中查看 -------------------------------------------------------------------------------- Q不知道杆单元有没有内力值
还有我要分析一个筏板基础用弹簧单元模拟土体不知道弹簧单元如何加到模型当中是先离散筏板还是先加弹簧单元如果先离散筏板如果单元不规则我如何计入弹簧单元的弹黄刚度弹黄刚度 KN/m=弹黄刚度KN/m3 * 单元面积m2
A杆单元有内力值, 查help看你所用的单元的element table分析筏板基础, 可选link10单元模拟土弹簧设定单元只能受压
A记住有一点查到的力和弯矩对应的是单元的局部坐标系
Q各位大虾我现在做桩土共同作用的分析想用ansys做数值模拟可怜我啥也不懂啥也不会又无老师指点希望各位大侠能指点一二不胜感激
A桩土共同作用 不知道你指的是以桩为主还是以土为主
如果是有限元同时分析土体和桩体土体用实体单元分析ANSYS可能难以很好的模拟因为它没有较好的土体本构关系目前岩土工程中较通用的土本构关系是Duncan-Chueng, Modified Cambridged model等 如果以分析桩即其他上部结构为主对桩体划分单元而土体用一系列非线性弹簧代替桩侧的摩阻力用剪切弹簧桩端阻力和水平受力桩的土阻力用压缩弹簧很多建筑结构和桥梁工程的共同作用问题就是这样做的 更一步简化将整个桩土系统单桩用在桩头处的三个方向的刚度来代替可以用于上部结构和地基及基础的共同作用分析中
可以参考计算土力学朱百里沈珠江上海科技出版社1990
Q由于现行的公路桥梁基础规范桩基的容许承载力是桩侧摩阻+桩端承载力安全系数2但是桩基的承载力
是桩侧摩阻先发挥到极限然后才是桩端承载力发挥到极限所以对两 29
者共同取2的安全系数是不科学的还有道桥桩基采用单参数的m法来计算桩侧土体的抗力重庆交院的吴恒力有个双参数法还有p-y曲线法我很想做个数值模拟看看土对桩的抗力分布如果有时间还想用随机有限元分别对桩侧摩阻和桩端承载力求取安全系数我硕士论文就准备做这个真希望peterman722等各位老师能多指导我不胜感激
若是ansys难以很好的模拟marc程序好不好一点呢
A你是指水平抗力吧 看来是要用有限元了而且是针对土体反映的研究要对土体划分单元m法等方法已经假定了土的抗力特征 Q请问
1ansys中的Drcker-Prager屈服准则用在土的弹塑性模拟中精度怎么样 2有没有人采用广西大学秦荣的QR及样条函数法进行过桩土共同作用的分析 3ansys有没有无限元单元和岩石节理元
AD-P模型虽然没有帽盖因为与Mohr-C模型的密切关系还是可以模拟土的很多特点但是国内用的人很少Duncan-Chueng, Modified Cambridged 等是在国内流行的
对于实际工程除了合适的模型外参数的获得是至关重要的参数的差异可能大于由于 模型带来的差异 你的意思好象是模型分析自己定参数对吗我认为限于定性研究D-P模型应该可以的但是精度我觉得这里没有精度概念因为你是理想化的研究
-------------------------------------------------------------------------------- 利用ANSYS随机振动分析功能实现随机疲劳分析
因篇幅问题不能全部显示,请点此查看更多更全内容