ANSYS压电分析命令流详解

合集下载

ANSYS一段命令流解释

ANSYS命令流解释/PREP7 ! 进入前处理器ET,1,BEAM4 ! 定义单元beam4,是个三维梁单元KEYOPT,1,2,0 ! 定义单元的关键选项，如果后面是0，代表默认的，可以先不用理解，KEYOPT,1,6,0KEYOPT,1,7,1 ! 这里定义第七个关键选项，定义为编号1，也就是计算陀螺阻尼矩阵方程，要求! IYY 等于IZZ，也就是两个转动惯量要相等，这两个量要在实常数中定义! 也就是下面的命令RKEYOPT,1,9,0KEYOPT,1,10,0*SET,p,acos(-1) ! 定义三个参数，分别是派（4.1315）、第一个半径R1，第二个半径r2*SET,R1,5*SET,R2,60 ! 半径的单位! 定义单元的实常数，有两个，因为有两个半径，分别就是下面的R,1! 和r,2! r命令的定义中需要根据使用单元beam4来一一对应，不同的单元R命令! 定义的意思是不一样的,具体每个意思，下图看，一定要一一对应! 这里就先是截面积，Z向转动惯量，Y像转动惯量，是一样的，从上面的! 关键选项定义中可以看的出来，KEYOPT,1,7,1 R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1, ,RMORE, ,p*R1**4/2, , ,2175, , ! 这个也是定义实常数，因为命令只能定义6个数，从第七个! 就要使用这个命令，编号可以从下图中对应下R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2, , ! 定义第二个实常数R，2 RMORE, ,p*R2**4/2, , ,2175, ,MPTEMP,,,,,,,, ! 这里来定义材料属性，目前楼主提供的是GUI操作以后的! ，这个需要相应简化,因为GUi操作中一般包括了温度的考虑，所以使用了几个温度的命令! 目前看来只是有一个温度，所以不用考虑温度，简化后的，在后面的注释中，注释开始是!（感叹号）MPTEMP,1,0MPDATA,EX,1,,2e5MPDATA,PRXY,1,,.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,1e-10! 简化后的材料属性! mp,ex,1,2e5 ! 分别定义弹性模量，泊松比，密度，其中的1代表第一个材料属性! mp,PRXY,1,0.3! mp,DENS,1,1e-10MPTEMP,,,,,,,, ! 这里就是定义第二个材料属性了！同样包括弹性模量，泊松比，密度MPTEMP,1,0MPDATA,EX,2,,2E5MPDATA,PRXY,2,,.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8E-8K, ,,,, ! 定义两个关键点，第一个是0，第二个坐标是：(100,0,0) K, ,100,,,TYPE, 1 ! 选择第一个单元，也就是beam4MAT, 1 ! 选择第一个材料属性REAL, 1 ! 选择第一个实常数ESYS, 0 ! 选择默认的单元坐标系统LSTR, 1, 2 ! 用刚才生成的两个关键点，建立一条直线LESIZE,ALL, , ,200, ,1, , ,1, ! 把建立的这条直线划分为200个LMESH, 1 ! 根据划分的线的段数，网格这条线，从而生成有限元模型D,1,UX ! 对于第一个节点（也就是原点位置的那个节点）进行约束，三个平移方向都约束，但是转动方向不约束，! 因为单元beam4是有六个自由度的D,1,UYD,1,UZD,102,UY ! 对于第102个节点，这样来算他的位置，总共长度是100，划分为200分，那么102是那个位置，楼主应该知道了吧! 选择了Y方向和Z方向的平移约束D,102,UZFLST,2,1,2,ORDE,1 ! 后面的两个命令结合起来就是选择了第200个节点，但是没有具体操作，! 只是有了拾取选择这个点的操作，这两个命令一般是GUI直接生成的，写命令流的一般不使用这两个操作! 所以这两个命令其实什么任务都没有做目前FITEM,2,200对上面的命令一一般是需要简化的，简化后的就是：/PREP7ET,1,BEAM4KEYOPT,1,7,1*SET,p,acos(-1)*SET,R1,5*SET,R2,60R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1, ,RMORE, ,p*R1**4/2, , ,2175, ,R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2, ,RMORE, ,p*R2**4/2, , ,2175, ,mp,ex,1,2e5mp,PRXY,1,0.3mp,DENS,1,1e-10mp,ex,2,2e5mp,PRXY,2,0.3mp,DENS,2,8e-8K,1,,,,K,2,100,,,MAT, 1REAL, 1LSTR, 1, 2LESIZE,ALL, , ,200LMESH, 1D,1,UXD,1,UYD,1,UZD,102,UYD,102,UZFLST,2,1,2,ORDE,1FITEM,2,200简化后的命令楼主结合本人说的，可以理解下这段话因为只是在建模，所以没有执行什么计算的，思路其实我已经说的很明确了：首先定义单元，指定一个关键选项包括了陀螺阻尼效果，然后再定义两个实常数，然后定义了两个材料属性，然后建立一个一个线，使用第一个实常数，和第一个材料属性对这个线进行了网格！其次再对这个线进行了约束的施加，原点位置三个平移方向全约束，对102个节点位置，进行了两个方向的约束。

ANSYS压电分析命令流详解

ANSYS压电分析命令流详解FINISH/CLEAR,START/PREP7!-------------------------------------------定义压电单元及参数!单元定义ET,1,SOLID5KEYOPT,1,1,3KEYOPT,1,3,0KEYOPT,1,5,0!定义材料MPTEMP,1,0MPDE,DENS,1MPDATA,DENS,1,,7700MPTEMP,1,0MPDE,PERX,1MPDATA,PERX,1,,826MPDE,PERY,1MPDATA,PERY,1,,919MPDE,PERZ,1MPDATA,PERZ,1,,919-------------------------------------------d矩阵,球坐标系下，r方向为3方向TBDE,PIEZ,1TB,PIEZ,1,,,1TBMODIF,1,1,3.56E-010TBMODIF,1,2,0TBMODIF,1,3,TBMODIF,2,1,-1.48E-010TBMODIF,2,2,0TBMODIF,2,3,TBMODIF,3,1,-1.48E-010TBMODIF,3,2,0TBMODIF,3,3,TBMODIF,4,1,0TBMODIF,4,2,4.32E-010TBMODIF,4,3,TBMODIF,5,1,0TBMODIF,5,2,TBMODIF,5,3,0TBMODIF,6,1,TBMODIF,6,2,0TBMODIF,6,3,4.32E-010-------------------------------------------d矩阵,球左边系下，r方向为3方向-------------------------------------------e矩阵,球左边系下，r方向为3方向TBDE,ANEL,1,,,TB,ANEL,1,1,21,1TBTEMP,0TBDATA,,11.087e-12,7.509e-12,7.509e-12,0,0,0TBDATA,,12.035e-12,7.5179e-12,0,0,0,12.035e-12TBDATA,,0,0,0,2.1053e-12,0,0TBDATA,,2.1053e-12,0,2.2584e-12,,,-------------------------------------------e矩阵,球左边系下，r方向为3方向!-------------------------------------------定义压电单元及参数!-------------------------------------------建模!-------------------------------------------r1内径!-------------------------------------------r2外径r1=7e-3r2=10e-3!-------------------------------------------球cyl4,,,r1,,r2,90VROTAT,all, , , , , ,3,2,90, ,!-------------------------------------------球!-------------------------------------------定义球坐标系，选择材料1，并弄将压电矩阵转置到球坐标系下CSWPLA,11,2,1,1,TYPE,1MAT,1REAL,ESYS,11SECNUM,!-------------------------------------------定义球坐标系，选择材料1，并弄将压电矩阵转置到球坐标系下!-------------------------------------------设置压电球网格密度CSYS,2 !球坐标系，X轴为r方向p1=5 !=============压电球厚度方向单元数LSEL,S,LOC,X,1.1*r1,0.9*r2LESIZE,all, , ,p1, , , , ,1LSEL,S,LOC,X,r1LESIZE,all, , ,4*p1, , , , ,1 !=============压电球周向单元数4*p1=============!-------------------------------------------设置压电球网格密度VMESH,all !画网格allsel!-------------------------------------------1/8球体位移边界条件ASEL,S,LOC,X,0APLOTNSLA,R,1D,ALL, ,0, , , ,UX, , , , ,ALLSELASEL,S,LOC,Y,0APLOTNSLA,R,1D,ALL, ,0, , , ,UY, , , , ,ALLSELASEL,S,LOC,Z,0APLOTNSLA,R,1D,ALL, ,0, , , ,UZ, , , , ,ALLSEL!-------------------------------------------进入求解器/SOLUALLSELANTYPE,TRANS !指定分析类型TRNOPT,FULL !指定瞬态分析方法NLGEOM,OFF !-----关闭大变形选项TINTP,,0.25,0.5,0.5 !压电瞬态分析中，瞬态积分参数EQSLV,SPAR !算法选择，稀疏矩阵算法OUTPR,ALL,ALL !控制结果的输出OUTRES,ALL,2 !保存所有的数据到结果文件中，2步输出一个结果文件TIMINT,ON,ALL !使用瞬态效应!-------------------------------------------施加电压1,v1=0CSYS,2ASEL,S,LOC,X,r1APLOTNSLA,R,1*SET,V1,0CP,1,VOLT,ALL !耦合电压*GET,N1,NODE,,NUM,MIND,N1,VOLT,V1 !施加电压荷载ALLSEL,ALL!-------------------------------------------施加电压2 ASEL,S,LOC,X,r2NSLA,R,1*SET,V2,200CP,2,VOLT,ALL !耦合电压*GET,N2,NODE,,NUM,MIN!D,N2,VOLT,V2 !施加电压荷载ALLSEL,ALLCSYS,0!----------------------------------施正弦加电压2，并求解f=100000 !激励电压频率T=1/fst=25 !每一个周期求解步数DT=T/ststepp=5*st !5个周期*DO,J,1,steppt1=f*DT*J !5个正弦波D,N2,VOLT,V2*SIN(2*3.1415926*t1) !施加电压荷载TIME,DT*J !载荷时间步NSUBST,1,,,1SOLVE*ENDDOTIME,2*stepp*DT !再算5周期NSUBST,50,,,1 !5周期，计算50步KBC,0SOLVEFINISH/POST26 !时间后处理器NSEL,S,LOC,X,r2 NSEL,r,LOC,y,NSEL,r,LOC,z,*GET,N3,NODE,,NUM,MINNSOL,3,N3,U,X, UX_3, PLVAR,3,。

ANSYS命令解释

ANSYS命令解释分享⼀下！BEGINNER'S GUIDE TO ANSYS COMMANDSThe symbol '*' corresponds to the following:* --> k, l, a, v, e, n, cm, et, mp, r where ==>k --> Keypointsl --> Linesa --> Areav --> Volumese --> Elementsn --> Nodescm --> componentet --> element typemp --> material propertyr --> real constant$ --> d, f, sf, bf, ic, where ==>d --> DOF constraint (ux... in Structural, Temp in thermal,f --> Force Load ( Heat in thermal)sf --> Surface load on nodesbf --> Body Force on NodesMore Commands can be generated by sensible combinations of " $* " family of commands. See the following list of $* possible options$* --> dk --> DOF constraints on KP (Vx,Vy,Pres... in CFD)dl --> DOF constraints on Linesda --> DOF constraints on Areasfk --> Force on Keypointssfl --> Surface load on Linessfa --> Surface load on Areassfe --> Surface load on element faces中国热模⽹⾸发bfk --> Body Force on Keypointsbfl --> Body Force on Linesbfa --> Body Force on Areabfv --> Body Force on Volumesbfe --> Body Force on Elementsic --> Initial Conditions ",p" --> If ",p" was issued at the end of the Command(in Input Window) the GUI based picking menu will be activated. Useful for listing, plotting, meshing, deleting, etc..**********************************************************1. Listing of picked Entities:COMMAND SYNTAX: *LIS,p & $*LIS,pA few Combinations of this command are:klis,p --> List KPllis,p --> Lists Linesalis,p --> Lists Areasvlis,p --> Lists Volumeselis,p --> Lists Elementsnlis,p --> Lists Nodescmlis,p --> Lists componentscslis,p --> Lists user created local co-ordinate systemsdlis,p --> Lists DOF constraints specfied on nodesdalis,p --> Lists DOF constraints applied on Areasflis,p --> Lists force on Nodes中国热模⽹⾸发sfllis,p --> Lists Surface Load on linesbfalis,p --> Lists body force load applied on Areasiclis,p --> Lists Initial condition on NodesIf ",p" was not issued, all entites currently selected will be listed.For certain commands ",p" cannot be issued. See the below mentioned commandsetlis --> Lists the different element types definedmplis --> Lists whatever Material propertiesrlis --> Lists whatever real constantscslis --> Lists all co-ordinate systemscmlis --> Lists all components*********************************************************2. Plotting of Entities: COMMAND SYNTAX: *plo KPLO / LPLO / APLO / VPLO / EPLO / NPLO / CMPLO / **********************************************************3. deleting of Entities:COMMAND SYNTAX: *DEL,p & $*DEL,pKDEL,p / LDEL,p / ADEL,p / VDEL,p / EDEL,p / NDEL,p / CMDEL,p / DDEL,p /DKDEL,p / DADEL,p / FDEL,p / SFDEL,p / SFEDEL,p / SFADEL,p / SFLDEL,p / 热点模具⽹BFADEL,p / ......The syntax for this command is very similar to LISTING command.**********************************************************4. distance between two entities:COMMAND SYNTAX: *DIS,pndis,p --> Distance between two nodeskdis,p --> Distance between two KPs**********************************************************5. Meshing of geometries:COMMAND SYNTAX: *MES,pKMES,p / LMES,p / AMES,p / VMES,p**********************************************************6. Size settings for Lines and Areas before meshing :COMMAND SYNTAX : *size,,p Lesiz,p / Aesize,p*********************************************************7. Clearing Meshes of already meshed geometries:COMMAND SYNTAX: *CLE,p KCLE,p / LCLE,p / ACLE,p / VCLE,p**********************************************************8. BOOLEAN Operations: IntersectCOMMAND SYNTAX : *IN* AINA,p / VINV,p / LINL,p / AINV,p / LINV,p / LINA,p 热点模具⽹**********************************************************9. BOOLEAN Operations: GLUECOMMAND SYNTAX : *GLUE VGLUE,p / AGLUE,p / LGLUE,p**********************************************************10. Boolean Operations: SUBTRACT/DIVIDE:COMMAND SYNTAX: *sb*,p See the following examples to understand how this works:asba,p --> Subtract Area from Areaasbl,p --> Divide Area by linevsba,p --> Divide volume by Arealsbw,p --> Divide line by Workplanevsbw,p --> Divide volume by Workplaneasbw,p --> Divide area by Workplanevsbv,p --> subtract Volume by another volumeMore combinations exist. The user needs to explore them for themselves --> forms a part of learning**********************************************************11. Boolean Operations: Overlap:COMMAND SYNTAX: *OVLAP,p AOVLAP,p / VOVLAP,p**********************************************************12. Concatenation of Lines / Areas --> for map meshing热点模具⽹COMMAND SYNTAX : *ccat,pLCCAT,p --> Concatenation of Lines for Map meshing AreaACCAT,p --> Concatenation of Areas for Map meshing Volume*********************************************************13. Dragging operationCOMMAND SYNTAX : *drag,pvdrag,p --> Drag areas along a line to create a new volumeadrag,p --> Drag line along a line to create a new arealdrag,p --> Drag KP along a line to create a new line**********************************************************14. Copy Geomtric entitiesCOMMAND SYNTAX : *GEN,,pKGEN,,p / LGEN,,p / AGEN,,p / VGEN,,pPlease note that *GEN commands are also used for MOVE operations. The difference lies in the value specified in the 10th field of these *GEN commands. By default it is 0 --> which does the COPY operation. If specfied as 1 --> it does the MOVE operation**********************************************************15. Bottom -to- Top modeling commands:COMMAND SYNTAX : *,p & **,p 热点模具⽹k,p ---> Allows user to pick KP in the Workplanel,p ---> Create lines from existing KPak,p ---> Create area from KPal,p ---> Create areafrom linesv,p ---> Create Volume from KPva,p ---> Create Volume from Arease,p ---> Create Elem from existing nodesen,p ---> Create Elem from nodes**********************************************************16. To apply common Boundary Conditions such as DOF constraint, Forces, Surface Loads, Bodyforce Loads and Initial conditions* --> is meant for the KLAVE entities only (KLAVEN stands for KP, Lines, Area, Volumes & ELem )16a. DOF constraint :COMMAND SYNTAX : $*,p ( Please Note: NOT all * are valid)See the valid combinations below:D,p --> To apply DOF on nodesDK,p --> To apply DOF on KeypointsDL,p --> Apply DOF on LinesDA,p --> Apply DOF on Areas ( symmetry or Anti-symmetry will be prompted)****************16b. FORCE Loading:COMMAND SYNTAX : $*,p 中国热模⽹⾸发See the valid combinations below:f,p --> Forces on nodesfk,p --> Force on Keypoints(fa,p or FV,p or FL,p ----> Since force cannot be applied on Lines or Area & volumes... this command does not exist.)****************16c. Surface Loads:COMMAND SYNTAX : $*,pSee the valid combinations below:sf,p --> Surface Load on a set of Nodessfl,p --> Surface Load on Linessfa,p --> Surface Load on Areasfe,p --> Surface Load on Element(SFk,p and SFV,p do not exist since pressure cannot be applied on a single Kp and neither can it be applied on a volume) ****************16d. BodyForce Load: COMMAND SYNTAX : bf*,pSee the valid combinations below:bf,p --> Bodyforce Load on a set of Nodesbfk,p --> Bodyforce Load on KPbfl,p --> Bodyforce Load on Linesbfa,p --> Bodyforce Load on Areasbfv,p --> Bodyforce Load on Volumesbfe,p --> Bodyforce Load on Element****************16e. Initial conditions:热点模具⽹ic,p --> Initial Conditions on Nodes(P.S: Initial Conditions can be applied only to nodes. )***********************************************************17. To refine a mesh :COMMAND SYNTAX : *ref,pkREF,p / kREF,p / aREF,p / eREF,p / nREF,p***********************************************************18. To TRANsfer loads from the Solid model to the FE model:COMMAND SYNTAX : $TRANdtran / ftran / sftran / bftran & SBCTRAN(SBCTRAN --> Transfers all solid model loads to FE model)***********************************************************19. Writing / Reading information to a file (ASCII)COMMAND SYNTAX : *read, & *write,NWRITE / MPWRITE / ETWRITE / RWRITE / EWRITE / CDWRITENREAD / MPREAD / ETREAD / RREAD / EREAD / CDREAD / LDREAD(Some of these commands ETWRITE/ETREAD , RWRITE/RREAD are undocumented. But they do work) The Commands CDWRITE and CDREAD are used to write/read all FE model related info (w or w/o geometry to ASCII files) Its recommended the user read the online help on these two commands before using them 热点模具⽹The LDREAD commands are used to read loads (LD) from other analysis types. For example: Tempfrom Thermal results file (*.rth) are applied onto Structural elements.好，我来补充⼀下楼上师兄的命令。

ansys命令流最全详细介绍二

三生成关键点和线部分1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例:K,1,0,0,02.在激活坐标系生成直线LSTR,关键点P1,关键点P2例LSTR,1,23.在两个关键点之间连线L,关键点P1,关键点P2例L,1,2注:此命令会随当前的激活坐标系不同而生成直线或弧线4.由三个关键点生成弧线LARC,关键点P1,关键点P2,关键点PC,半径RAD例LARC,1,3,2,0.05注:关键点PC是用来控制弧线的凹向5.通过圆心半径生成圆弧CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG例:CIRCLE,1,0.05,,,,46.通过关键点生成样条线BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6例:BSPLIN,1,2,3,4,5,67.生成倒角线LFILLT,线NL1,线NL2,倒角半径RAD例LFILLT,1,2,0.0058.通过关键点生成面A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8...例:A,1,2,3,49.通过线生成面AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10例:AL,5,6,7,810.通过线的滑移生成面ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9例:ASKIN,1,4,5,6,7,8注:线1为滑移的导向线四目标:掌握常用的实体-面的生成生成矩形面1.通过矩形角上定位点生成面BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH例:BLC4,0,0,5,3,02.通过矩形中心定位点生成面BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH注:与上条命令的不同就在于矩形的定位点不一样例:BLC5,2.5,1.5,5,3,03.通过在工作平面定义矩形X.Y坐标生成面RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y坐标Y1,矩形上边界Y坐标Y2例:RECTNG,0,5,0,3生成圆面4.通过中心定位点生成实心圆面CYL4,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,圆面的内半径RAD1,内圆面旋转角度THETA1,圆面的外半径RAD2,外圆面旋转角度THETA2,圆面的深度DEPTH注:如要实心的圆面则不用RAD2,THETA2,DEPTH例:CYL4,0,0,5,3605.生成扇形圆面命令介绍如上例1实心扇形:CYL4,0,0,5,60例2扇形圆环:CYL4,0,0,5,60,10,60例3整的圆环:CYL4,0,0,5,360,10,360注:同时可通过定义圆面的深度以生成柱体6.通过在工作平面定义起始点生成圆面CYL5,开始点X坐标XEDGE1,开始点Y坐标YEDGE1,结束点X坐标XEDGE2,结束点Y坐标YEDGE2,圆面深度DEPTH例:CYL5,0,0,2,2,7.通过在工作平面定义内外半径和起始角度来生成圆面PCIRC,内半径RAD1,外半径RAD2,起始角度THETA1,结束角度THETA2例LCIRC,2,5,30,1808.生成面与面的倒角AFILLT,面1的编号NA1,面2的编号NA2,倒角半径RAD例:AFILLT,2,5,2下一讲:多边形面的生成。

常用-ansys命令流中文注释与讲解

ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。

用ANSYS软件分析压电换能器入门

用ANSYS软件分析压电换能器入门A：分析过程基本步骤一：问题描述（草稿纸上完成）1：画出换能器几何模型，包括尺寸2：选定材料3：查材料手册确定材料参数二：建立模型1：根据对称性确定待建模型的维数2：根据画出的几何模型确定关键点坐标，给关键点编好号码3：建立一个文件夹用于当前分析4：启动ANSYS软件，指定路径到建立的文件夹，5：定义单元类型压电换能器分析使用的单元类型：solid5：8个节点3D六面体耦合场单元（也可缩减为三角柱形单元或四面体单元）。

无实常数。

plane13：4个节点2D四边形耦合场单元（也可缩减为三角形单元）。

无实常数。

solid98：10个节点3D四面体耦合场单元。

无实常数。

Fluid30：8个节点3D六面体声学流体单元（也可缩减为三角柱形单元或四面体单元）。

应用于近场水和远场水。

实常数为参考声压，可缺省。

Fluid130：4个节点面无穷吸收水声学流体单元（也可缩减为三角形面单元）。

实常数：半径，球心X，Y，Z坐标值。

6：定义材料参数对一般均匀各向同性材料要给出材料密度，杨氏模量，泊松系数。

（静态分析不用密度）对压电材料：一般使用的压电方程：e型压电方程，因此输入的常数为注意！一般顺序为：XX，YY，ZZ，YZ，XZ，XY。

在ANSYS中为XX，YY，ZZ，XY，YZ，XZ。

因此，前两矩后三行和后三列要做相应变化。

7：建立关键点8：把关键点连成线9：把线段围成面10：通过适当的方法生成体11：指定单元类型和材料参数12：划分线段13：划分体单元14：坐标转换，（转换到柱坐标系下）15：节点转换三：加载约束条件1：加载边界约束条件2：电极上加电压四：求解1：模态分析2：谐响应分析五：查看结果1：查看模态分析结果，计算导纳。

2：各模态的动态演示3：查看谐响应分析结果，计算导纳、发射与接收响应。

六：生成命令流文件1：给程序分块，添加适当的注释2：把相应参数具体值改成变量，同时给变量赋值B：空气中建模过程一：问题描述弯曲式换能器实体模型为轴对称结构。

ANSYS 命令流解释大全

,例：vsel,s,volu,,14
vsel,a,volu,,17,23,2
上面的命令选中了实体编号为 14，17，19，21，23的五个实体
u rforce, nvar, node, item, comp, name 指定待存储的节点力数据
nvar: 变量号
node: 节点号
item comp
F x, y.z
mu: 摩擦系数
dens: 质量密度
mat: 材料编号（缺省为当前材料号）
c 材料特性值，或材料之特性，温度曲线中的常数项
c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数
二、定义DP材料：
首先要定义EX和泊松比：MP，EX，MAT，……
MP，NUXY，MAT，……
定义DP材料单元表（这里不考虑温度）：TB，DP，MAT
!第二个（或后继）载荷步：
TIME,...
ESEL,...
EKILL,... !杀死选择的单元
ESEL,...
EALIVE,... !重新激活选择的单元
...
FDELE,... !删除不活动自由度的结点载荷
D,... !约束不活动自由度
...
F,... !在活动自由度上施活的自由度上的约束
Command:EKILL
GUI: Main Menu>Solution>-Load Step Opts-Other>Kill Elements
用下列命令重新激活单元：
Command: EALIVE
GUI: Main Menu>Solution>-Load Step Opts-Other>Activate Elem
SAVE

Ansys压电分析

ansys压电分析压电效应分析是一种结构－电场耦合分析。

当给石英和陶瓷等压电材料加电压时，它们会产生位移，反之若使之振动，则会产生电压。

压力传感器就是压电效应的一种典型的应用。

压电分析（ANSYS/Multiphysics或ANSYS /Mechanical软件包提供这种分析）类型可以是静力、模态、预应力模态、谐波、预应力谐波和瞬态分析。

压电分析只能用下列单元类型之一：PLANE13, KEYOPT(1)=7，耦合场四边形实体单元SOLID5, KEYOPT(1)=0或3，耦合场六面体单元SOLID98, KEYOPT(1)=0或3，耦合场四面体单元KEYOPT选项激活压电自由度：位移和电压。

对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。

注意：如果模型中激活了至少一个带有压电自由度（位移和VOLT）的单元，则需要用到VOLT自由度的所有单元必须是上面三种压电单元其中之一。

而且，所有的这些单元均需激活压电自由度。

如果不希望在这些单元中存在压电效应，则需给材料定义非常小的压电特性。

压电KEYOPT用NLGEOM，SSTIF，PSTRES命令可用大挠度和应力刚化作用（参见《ANSYS Commands Reference》对这些命令的更多信息，参见《ANSYS Structural Analysis Guide》及《ANSYS，Inc. Theory Reference》的第三章关于大挠度及应力刚化功能的更多信息）。

对PLANE13，通过设置KEYOPT（1）＝7可用大挠度及应力刚化功能。

对SOLID5及SOLID98通过设置KEYOPT（1）＝3可用大挠度及应力刚化功能。

而且小挠度及应力刚化选项可以通过KEYOPT（1）＝0使用。

注意－对压电分析不能使用自动求解控制。

SOLCONTROL缺省设置只能对纯结构或纯热分析使用。

对大挠度压电分析，必须用非线性求解命令定义有关设置。

ansys命令流中文说明

ansys命令流中文说明展开全文KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同，仅对梁或壳有效，对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反，仅对梁或壳有效，对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的，而不是选单元，如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点， Elem,单元，kp: 关键点（也合并线，面及点）mat: 材料，type: 单元类型，Real: 实常数cp：耦合项，CE：约束项，CE: 约束方程，All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留（low, high）注意：可以先选择一部分项目，再执行合并。

如果多次发生合并命令，一定要先合并节点，再合并关键点。

合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。

u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点（如面内线，关键点处的节点）u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组（？）Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目（缺省）VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。

Ansys命令流大全剖析

Ansys彭彭(沈阳化工大学机械工程学院，辽宁沈阳110142)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。

点要依次序输入，输入的顺序会决定面的法线方向。

如果超过四个点，则这些点必须在同一个平面上。

Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。

7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。

8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。

Ansys压电分析

ansys压电分析压电效应分析是一种结构－电场耦合分析。

当给石英和陶瓷等压电材料加电压时，它们会产生位移，反之若使之振动，则会产生电压。

压力传感器就是压电效应的一种典型的应用。

压电分析（ANSYS/Multiphysics或ANSYS /Mechanical软件包提供这种分析）类型可以是静力、模态、预应力模态、谐波、预应力谐波和瞬态分析。

对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。

注意：如果模型中激活了至少一个带有压电自由度（位移和VOLT）的单元，则需要用到VOLT自由度的所有单元必须是上面三种压电单元其中之一。

而且，所有的这些单元均需激活压电自由度。

如果不希望在这些单元中存在压电效应，则需给材料定义非常小的压电特性。

对PLANE13，通过设置KEYOPT（1）＝7可用大挠度及应力刚化功能。

对SOLID5及SOLID98通过设置KEYOPT（1）＝3可用大挠度及应力刚化功能。

而且小挠度及应力刚化选项可以通过KEYOPT（1）＝0使用。

注意－对压电分析不能使用自动求解控制。

SOLCONTROL缺省设置只能对纯结构或纯热分析使用。

对大挠度压电分析，必须用非线性求解命令定义有关设置。

(完整版)ANSYS命令流总结(全)

ANSYS结构解析单元功能与特征/POST1/可以构成一一些命令，一般是一种整体命令（ session),三十也有特别，比方是办理 ! 是说明说明符号，，与其余软件的说明是相同的， ansys 不作为命令读取，*此符号一般是 APDL 的表记符，也就是 ansys 的参数化语言，如 *do ,,,*enddo 等等NSEL 的意思是node select，即选择节点。

s 就是 select，选择。

DIM是定义数组的意思。

array 数组。

MP 命令用来定义资料参数。

K 是建立要点点命令。

K, 要点点编号 ,x 坐标 ,y 坐标， z 坐标。

K, NPT, X, Y , Z 是定义要点点， K 是命令， NPT 是要点点编号， XYZ 是坐标。

NUMMRG , keypoint 用这个命令，要保证要点点的地点完整相同，不过要点点号不一样样的才行。

这个命令关于重复的线面都可以用。

这个很简单，压缩要点。

Ngen 复制节点e,节点号码：这个命令式经过节点来形成单元NUMCMP,ALL ：压缩所有编号，这样你所有的线都会挨次次重新编号 ~你若是需要固定的线固定的标号NSUBST,100,500,50 ：经过指定子步数来设置载荷步的子步LNSRCH 线性搜寻是求解非线性代数方程组的一种技巧，此法会在一段区间内，以必定的步长逐渐搜寻根，对比常用的牛顿迭代法所要耗费的计算量大得多，但它可以防备在一些状况下牛顿迭代法出现的跳跃现象。

LNSRCH激活线性搜寻PRED 激活自由度求解展望NEQIT 指定一个荷载步中的最大子步数AUTOTS自动求解控制打开自动时间步长.KBC -指定阶段状也许用跳板装载里面一个负荷步骤。

SPLINE :P1， P2， P3，P4， P5， P6， XV1 ， YV1 ， ZV1 ， XV6 ，YV6 ， ZV6 （生成分段样条曲线）*DIM ， Par，Type ，IMAX ，JMAX ， KMAX ， Var1，Var2， Var3（定义载荷数组的名称）【注】 Par: 数组名Type： array 数组，仿佛fortran, 下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8 个字符）tableIMAX ， JMAX ， KMAX各维的最大下标号Var1， Var2，Var3 各维变量名，缺省为row,column,plane( 当 type 为 table 时 )/config 是设置 ansys 配置参数的命令格式为 /CONFIG, Lab, V ALUELab 为参数名称value 为参数值比方： /config ， MXEL ，10000 的意思是最大单元数为10000杆单元 : LINK1、 8、 10、 11、 180梁单元： BEAM3、 4、 23、 24，44， 54， 188， 189管单元 : PIPE16， 17， 18， 20， 59， 602D实体元 : PLANE2， 25， 42， 82， 83， 145，146， 182， 1833D实体元 : SOLID45， 46， 64，65， 72， 73，92， 95， 147，148， 185， 186，187， 191壳单元 : SHELL28， 41， 43， 51， 61， 63， 91， 93， 99， 143， 150， 181，208， 209弹簧单元 : COMBIN7， 14， 37，39， 40质量单元 : MASS21接触单元 : CONTAC12， 52， TARGE169， 170， CONTA171， 172， 173， 174， 175， 178矩阵单元 : MATRIX27， 50表面效应元 : SURF153， 154粘弹实体元 : VISCO88， 89， 106， 107， 108，超弹实体元 : HYPER56， 58， 74， 84， 86， 158耦合场单元 : SOLID5， PLANE13， FLUID29， 30，38， SOLID62， FLUID79， FLUID80，81，SOLID98， FLUID129， INFIN110 ， 111， FLUID116，130界面单元 : INTER192， 193， 194， 195显式动力解析单元 : LINK160， BEAM161， PLANE162， SHELL163， SOLID164， COMBI16杆单元单元名称简称节点数节点自由度特征备注LINK12D杆2Ux,Uy EPCSDGB常用杆元LINK83D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉EDGB模拟缆索的废弛及或仅受压杆缝隙LINK113D线性调理EGB模拟液压缸和大转器动LINK1803D有限应变杆EPCDFGB另可考虑粘弹塑性E- 弹性 (Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度deflection)， F- 大应变 (Large strain)或有限应变(Finite strain)，B-单元存亡(Birth and dead),G-化 (Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应降落(Adaptive descent)等。

ANSYS命令流解释大全

A N S Y S命令流解释大全集团企业公司编码：（LL3698-KKI1269-TM2483-LUI12689-ITT289-一、定义材料号及特性mp,lab,mat,co,c1,…….c4lab:待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex:弹性模量nuxy:小泊松比alpx:热膨胀系数reft:参考温度reft:参考温度prxy:主泊松比gxy:剪切模量mu:摩擦系数dens:质量密度mat:材料编号（缺省为当前材料号）c材料特性值，或材料之特性，温度曲线中的常数项c1-c4:材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,...!设定时间值（静力分析选项）NLGEOM,ON!打开大位移效果NROPT,FULL!设定牛顿－拉夫森选项ESTIF,...!设定非缺省缩减因子（可选）ESEL,...!选择在本载荷步中将不激活的单元EKILL,...!不激活选择的单元ESEL,S,LIVE!选择所有活动单元NSLE,S!选择所有活动结点NSEL,INVE!选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0!约束所有不活动的结点自由度（可选）NSEL,ALL!选择所有结点ESEL,ALL!选择所有单元D,...!施加合适的约束F,...!施加合适的活动结点自由度载荷SF,...!施加合适的单元载荷BF,...!施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF 命令得到更详细的解释。

ANSYS命令流详解（超全）

ANSYS命令流详解（超全）一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。

(完整版)ANSYS最常用命令流+中文注释(超级大全)

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。

Ansys压电分析

ansys压电分析压电效应分析是一种结构－电场耦合分析。

当给石英和陶瓷等压电材料加电压时，它们会产生位移，反之若使之振动，则会产生电压。

压力传感器就是压电效应的一种典型的应用。

压电分析（ANSYS/Multiphysics或ANSYS /Mechanical软件包提供这种分析）类型可以是静力、模态、预应力模态、谐波、预应力谐波和瞬态分析。

对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。

注意：如果模型中激活了至少一个带有压电自由度（位移和VOLT）的单元，则需要用到VOLT自由度的所有单元必须是上面三种压电单元其中之一。

而且，所有的这些单元均需激活压电自由度。

如果不希望在这些单元中存在压电效应，则需给材料定义非常小的压电特性。

对PLANE13，通过设置KEYOPT（1）＝7可用大挠度及应力刚化功能。

对SOLID5及SOLID98通过设置KEYOPT（1）＝3可用大挠度及应力刚化功能。

而且小挠度及应力刚化选项可以通过KEYOPT（1）＝0使用。

注意－对压电分析不能使用自动求解控制。

SOLCONTROL缺省设置只能对纯结构或纯热分析使用。

对大挠度压电分析，必须用非线性求解命令定义有关设置。

ANSYS命令流学习(中文解释)

ANSYS命令流学习（中文解释）第一天目标:熟悉ansys基本关键字的含义sf-->surfaceforceonnodes表面载荷bf-->bodyforceonnodes体载荷ic-->initialconditions初始条件第二天目标:了解命令流的整体结构,掌握每个模块的标识!文件表明段/batch/titile,testanalysis!定义工作标题/filename,test!定义工作文件名/prep7!步入前处置模块标识!定义单元,材料属性,实常数段et,1,shell63!指定单元类型et,2,solid45!指定体单元mp,ex,1,2e8!指定弹性模量mp,prxy,1,0.3!输入泊松比mp,dens,1,7.8e3!输入材料密度r,1,0.001!选定壳单元实常数-厚度......!建立模型k,1,0,0,,!定义关键点k,2,50,0,,k,3,50,10,,k,4,10,10,,k,5,10,50,,k,6,0,50,,a,1,2,3,4,5,6,!由关键点生成面......!分割网格esize,1,0,amesh,1......finish!前处理结束标识/solu!步入解模块标识!施予约束和载荷dl,5,,allsfl,3,pres,1000sfl,2,pres,1000......solve!解标识finish!求解模块结束标识/post1!步入通用型后处理器标识....../post26!进入时间历程后处理器……/exit,save!选择退出并计算机上安装以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理log文件时有所帮助/angle!选定拖轴转动视图/dist!表明对视图展开翻转/device!设置图例的显示,如:风格,字体等/replot!重新显示当前图例/reset!恢复缺省的图形设置/view!设置观察方向/zoom!对图形表明窗口的某一区域展开翻转第三天生成关键点和线部分1.分解成关键点k,关键点编号,x坐标,y坐标,z坐标例:k,1,0,0,02.在转化成坐标系分解成直线lstr,关键点p1,关键点p2基准:lstr,1,23.在两个关键点之间连线l,关键点p1,关键点p2例:l,1,2备注:此命令可以随当前的转化成坐标系相同而分解成直线或弧线4.由三个关键点生成弧线larc,关键点p1,关键点p2,关键点pc,半径rad基准:larc,1,3,2,0.05注:关键点pc是用来控制弧线的凹向5.通过圆心半径分解成圆弧circle,关键点圆心,半径rad,,,,圆弧段数nseg例:circle,1,0.05,,,,46.通过关键点分解成样条线bsplin,关键点p1,关键点p2,关键点p3,关键点p4,关键点p5,关键点p6例:bsplin,1,2,3,4,5,67.分解成倒角线lfillt,线nl1,线nl2,倒角半径rad例:lfillt,1,2,0.005(如果不是圆角呢？)8.通过关键点分解成面a,关键点p1,关键点p2,关键点p3,关键点p4,关键点p5,关键点p6,p7,p8...例:a,1,2,3,4（关键点有没有顺序？）9.通过线分解成面al,线l1,线l2,线l3,线l4,线l5,线l6,线l7,线l8,线l9,线l10例:al,5,6,7,810.通过线的位移分解成面askin,线nl1,线nl2,线nl3,线nl4,线nl5,线nl6,线nl7,线nl8,线nl9例:askin,1,4,5,6,7,8注:线1为滑移的导向线第四天目标:掌控常用的实体-面的分解成生成矩形面1.通过矩形角上定位点分解成面blc4,定位点x方向坐标xcorner,定位点y方向坐标ycorner,矩形宽度width,矩形高度height,矩形深度depth例:blc4,0,0,5,3,02.通过矩形中心定位点分解成面blc5,定位点x方向坐标xcenter,定位点y方向坐标ycenter,矩形宽度width,矩形高度height,矩形深度depth注:与上条命令的不同就在于矩形的定位点不一样例:blc5,2.5,1.5,5,3,03.通过在工作平面定义矩形x.y座标分解成面rectng,矩形左边界x坐标x1,矩形右边界x坐标x2,矩形下边界y坐标y1,矩形上边界y坐标y2例:rectng,0,5,0,3分解成圆面4.通过中心定位点生成实心圆面cyl4,定位点x方向座标xcenter,定位点y方向座标ycenter,圆面的内半径rad1,内圆面转动角度theta1,圆面的外半径rad2,外圆面转动角度theta2,圆面的深度depth注:如要实心的圆面则不用rad2,theta2,depth例:cyl4,0,0,5,3605.分解成扇形圆面命令了解例如上例1实心扇形:cyl4,0,0,5,60基准2扇形圆环:cyl4,0,0,5,60,10,60基准3整的圆环:cyl4,0,0,5,360,10,360注:同时可通过定义圆面的深度以生成柱体6.通过在工作平面定义初始点分解成圆面cyl5,开始点x坐标xedge1,开始点y坐标yedge1,结束点x坐标xedge2,结束点y坐标yedge2,圆面深度depth例:cyl5,0,0,2,2,7.通过在工作平面定义内外半径和初始角度去分解成圆面pcirc,内半径rad1,外半径rad2,起始角度theta1,结束角度theta2例:pcirc,2,5,30,1808.分解成面与面的倒角afillt,面1的编号na1,面2的编号na2,倒角半径rad例:afillt,2,5,2下一谈:多边形面的分解成第五天目标:掌握多边形面和体的生成1.分解成多边形面命令:rpr4,多边形的边数nsides,中心定位点x坐标xcenter,中心定位点y坐标ycenter,中心定位点距各边顶点的距离radius,多边形旋转角度theta基准:rpr4,4,0,0,0.15,30注:这条命令可通过定义不同的nsides生成三边形,四边形,...,八边形2.分解成多边形体命令:rpr4,多边形的边数nsides,中心定位点x坐标xcenter,中心定位点y坐标ycenter,中心定位点距各边顶点的距离radius,多边形旋转角度theta,多边形的深度depth例:rpr4,4,0,0,0.15,30,0.1备注:多边形体和面命令唯一的相同就是深度depth的定义到此,关键点,线,面的生成讲解已结束,下一讲:体的生成第六天目标:掌控体的分解成命令1.通过关键点生成体命令:v,关键点p1,关键点p2,p3,p4,p5,p6,p7,p8基准:v,4,5,6,7,15,24,252.通过面生成体命令:va,面a1,面a2,a3,a4,a5,a6,a7,a8,a9,a10基准:va,3,4,5,8,103.通过长方形角上定位点生成体命令:blc4该命令前面在谈分解成面的时候已并作了解,唯一的相同是深度depth的定义.4.通过长方形中心定位点生成面命令:blc55.通过定义长方体初始边线分解成体命令:block,开始点x坐标x1,结束点x坐标x2,y1,y2,z1,z2例:block,2,5,0,2,1,36.分解成圆柱体基本命令通生成圆形面,不同在于depth的定义基本命令:cyl4基本命令:cyl5基本命令:cylind7.分解成棱柱基本命令通生成多边形,不同在于depth的定义基本命令:rpr48.通过球心半径分解成球体命令:sph4,球心x坐标xcenter,球心y坐标ycenter,半径rad1,半径rad2基准:sph4,1,1,2,59.通过直径上起始点坐标生成球体命令:sph5,起点x座标xedge1,起点y座标yedge1,完结点x座标xedge2,完结点y座标yedge2基准:sph5,2,5,7,610.在工作平面起点通过半径和转动角度生成球体命令:sphere,半径rad1,半径rad2,旋转角度theta1,旋转角度theta2基准:sphere,2,5,0,6011.生成圆锥体命令:cone,底面半径rbot,顶面半径rtop,底面低z1,顶面低z2,旋转角度theta1,旋转角度theta2基准:cone,10,20,0,50,0,180下一讲:布尔操作第七天目标:掌控常用的布尔操作命令1.沿法向延伸面生成体命令:voffst,面的编号narea,面弯曲的长度dist,关键点增量kinc基准:voffst,1,2,,2.通过坐标的增量延伸面生成体命令:vext,面1的编号na1,面2的编号na2,增量ninc,x方向的增量。

用ANSYS软件分析压电换能器入门

用ANSYS软件分析压电换能器入门压电换能器是一种将电能和机械能相互转换的装置，其原理是利用压电效应将机械应力转化为电荷或电位差，或将电场作用下产生的应变引起物体发生位移。

压电换能器在自动化控制、传感器、电能转换和储能装置等领域中有广泛的应用。

为了分析压电换能器的工作特性和性能，我选择使用ANSYS软件进行研究。

ANSYS是一种强大的有限元分析软件，可以对结构、固体力学、流体力学等领域进行模拟和分析，适用于各种复杂的物理现象和工程问题。

首先，我们需要建立压电换能器的几何模型。

在ANSYS中，可以通过几何建模工具创建三维模型，包括导入现有几何模型、绘制几何、创建参数化模型等方法。

我们可以按照实际的压电换能器尺寸和形状进行创建，或者通过参数化建模来进行优化设计。

建立几何模型后，需要将其进行网格划分，转化为有限元模型。

这可以通过ANSYS中的网格划分工具来完成，可以选择适合问题特性的网格划分方法，如结构网格、四面体网格、六面体网格等。

合适的网格划分可以更好地保证计算结果的准确性和稳定性。

接下来，需要定义材料的物性参数和边界条件。

在ANSYS中，可以设定材料的压电系数、弹性系数、密度等参数，并选择适合的材料模型。

例如，对于压电材料，可以选择线性压电模型或非线性压电模型。

同时，还需要设置边界条件，如施加的电压或电荷，施加的外力或位移等。

在定义完初始条件后，可以进行仿真计算并分析结果。

ANSYS提供了各种求解器和分析工具，可以得到压电换能器的电位差、电荷、应变、位移等物理量的变化规律。

可以通过改变输入条件或改变材料参数来研究其对压电换能器性能的影响，进而优化设计。

除了静态分析，ANSYS还可以进行动态分析，模拟压电换能器在不同频率和振幅下的工作特性。

可以通过ANSYS中的谐振分析、模态分析等工具来获取压电换能器的固有频率、模态形状及模态质量等信息。

最后，在完成仿真计算后，可以根据分析结果进行结果评估和后处理。

超全ANSYS命令流详细解析

1.A，P1，P2，…，P17，P18（以点定义面）2. AADD，NA1，NA2，…NA8，NA9（面相加）3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性）【注】ESYS为坐标系统号、SECN为截面类型号。

4. *ABBR，Abbr，String（定义一个缩略词）5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词）6. ABBSAVE，Lab，Fname，Ext（将当前定义的缩略词写入文件）7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值）【注】*************8. ACCAT，NA1，NA2（连接面）9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度）10. ACLEAR，NA1，NA2，NINC（清除面单元网格）11. ADAMS，NMODES，KSTRESS，KSHELL【注】*************12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC【注】*************13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算）14. ADELE，NA1，NA2，NINC，KSWP（删除面）【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。

15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面）16. AESIZE，ANUM，SIZE（指定面上划分单元大小）17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面）18. AFSURF，SAREA，TLINE（在既有面单元上生成重叠的表面单元）19. *AFUN, Lab（指定参数表达式中角度单位）20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE（复制面）21. AGLUE，NA1，NA2，…，NA8，NA9（面间相互粘接）22. AINA，NA1，NA2，…，NA8，NA9（被选面的交集）23. AINP，NA1，NA2，…，NA8，NA9（面集两两相交）24. AINV，NA，NV（面体相交）25. AL，L1，L2，…，L9，L10（以线定义面）26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息）【注】Lab=HPT时，显示面上硬点信息，默认为空。

压电驱动器ansys有限元分析命令流

压电驱动器ansys有限元分析命令流fini/CLEAR, START/CLEAR/VERIFY,VM175/PREP7ET,1,SOLID98MP,EX,4,1.79E5MP,PRXY,4,0.3MP,PERX,3,730 ! PERMITTIVITY (X AND Y DIRECTION) MP,PERY,3,730MP,PERZ,3,635 ! PERMITTIVITY (Z DIRECTION)TB,PIEZ,3 ! DEFINE PIEZ. TABLETBDATA,16,12.7 ! E61 PIEZOELECTRIC CONSTANT TBDATA,14,12.7 ! E52 PIEZOELECTRIC CONSTANT TBDATA,3,-5.2 ! E13 PIEZOELECTRIC CONSTANT TBDATA,6,-5.2 ! E23 PIEZOELECTRIC CONSTANT TBDATA,9,15.1 ! E33 PIEZOELECTRIC CONSTANT TB,ANEL,3 ! DEFINE STRUCTURAL TABLE TBDATA,1,13.9E4,7.78E4,7.43E4 ! INPUT [C] MATRIX TBDATA,7,13.9E4,7.43E4TBDATA,12,11.5E4TBDATA,16,2.56E4TBDATA,19,2.56E4TBDATA,21,3.06E4boptn,keep,offL=3500wa=200wb=40hp=2hsi=10aerfa=1.5volt1=43guding_kuan=200guding_chang=40r_yuanjiao=10PI=acos(-1)delta_di=L+2*r_yuanjiao/tan(aerfa/2/180*PI)delta_gao=delta_di/2*tan(aerfa/180*PI)zong_gao=delta_gao+wa/2+wb/cos(aerfa/180*PI)-guding_kuan/2*tan(aerfa/180*PI)+guding_chang k,1,-delta_di/2,-wa/2k,2,0,-wa/2k,3,0,-wa/2-delta_gaol,1,2l,1,3l,2,3lfillt,1,2,r_yuanjiaoal,1,2,3,4k,6,0,0k,7,0,-delta_gao-wa/2-wb/cos(aerfa/180*PI)k,8,0,-zong_gaok,9,-guding_kuan/2,-zong_gaok,10,-guding_kuan/2,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+guding_kuan/2*tan(aerfa/180*PI) k,11,-L/2-40-r_yuanjiao+10,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+(L/2+40)*tan(aerfa/180*PI) k,12,-L/2-40-r_yuanjiao+10,0l,6,8l,8,9l,9,10l,10,11l,11,12l,12,6al,5,6,7,8,9,10asba,2,1,,delete,deletevext,3,,,,,-hsivsymm,x,1vsymm,y,1vsymm,y,2vadd,1,2,3,4numcmp,volu !压缩体编号numcmp,area !压缩面编号!去掉多余的面asel,s,loc,z,0,0aadd,allallsel,allasel,s,loc,z,-hsi,-hsiaadd,allallsel,allasel,s,loc,y,-wa/2,-wa/2aadd,allallsel,allasel,s,loc,y,wa/2,wa/2aadd,allallsel,allasel,s,loc,y,-zong_gao,-zong_gaoaadd,allallsel,allasel,s,loc,y,zong_gao,zong_gaoaadd,allallsel,allasel,s,loc,x,-L/2-40-r_yuanjiao+10,-L/2-40-r_yuanjiao+10 aadd,allallsel,allasel,s,loc,x,L/2+40+r_yuanjiao-10,L/2+40+r_yuanjiao-10 aadd,allallsel,allnumcmp,area !压缩面编号!去掉z=0处多余的线asel,s,loc,z,0,0lsla,slsel,r,loc,y,-wa/2,-wa/2lcomb,all,,0allsel,allasel,s,loc,z,0,0lsla,slsel,r,loc,y,wa/2,wa/2lcomb,all,,0allsel,allasel,s,loc,z,0,0lsla,slsel,r,loc,y,-zong_gao,-zong_gaolcomb,all,,0allsel,allasel,s,loc,z,0,0lsla,slsel,r,loc,y,zong_gao,zong_gaolcomb,all,,0allsel,allasel,s,loc,z,0,0lsla,slsel,r,loc,x,-L/2-40-r_yuanjiao+10,-L/2-40-r_yuanjiao+10 lcomb,all,,0allsel,allasel,s,loc,z,0,0lsla,slsel,r,loc,x,L/2+40+r_yuanjiao-10,L/2+40+r_yuanjiao-10 lcomb,all,,0allsel,all!去掉z=-hsi处多余的线asel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,y,-wa/2,-wa/2lcomb,all,,0allsel,allasel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,y,wa/2,wa/2lcomb,all,,0allsel,allasel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,y,-zong_gao,-zong_gaolcomb,all,,0allsel,allasel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,y,zong_gao,zong_gaolcomb,all,,0allsel,allasel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,x,-L/2-40-r_yuanjiao+10,-L/2-40-r_yuanjiao+10 lcomb,all,,0allsel,allasel,s,loc,z,-hsi,-hsilsla,slsel,r,loc,x,L/2+40+r_yuanjiao-10,L/2+40+r_yuanjiao-10 lcomb,all,,0allsel,allnumcop,line !压缩线编号numcop,kp !压缩关键点编号k,,-L/2-40-r_yuanjiao+10,wa/2k,,-L/2-40-r_yuanjiao+10,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+(L/2+40)*tan(aerfa/180*PI) k,,-guding_kuan/2,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+guding_kuan/2*tan(aerfa/180*PI) k,,-guding_kuan/2,-zong_gaok,,guding_kuan/2,-zong_gaok,,guding_kuan/2,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+guding_kuan/2*tan(aerfa/180*PI) k,,L/2+40+r_yuanjiao-10,-delta_gao-wa/2-wb/cos(aerfa/180*PI)+(L/2+40)*tan(aerfa/180*PI) k,,L/2+40+r_yuanjiao-10,wa/2a,47,48,49,50,51,52,53,54k,,-delta_di/2,-wa/2k,,delta_di/2,-wa/2k,,0,-wa/2-delta_gaol,55,56l,55,57l,56,57lfillt,75,76,r_yuanjiaolfillt,75,77,r_yuanjiaoal,75,78,76,77,79asba,25,26,,delete,deletevext,27,,,,,hpvglue,1,2! vdele,3,6,1numcmp,volu !压缩体编号MAT,4VMESH,2MAT,3ESIZE,5VMESH,1FINISH/SOLUANTYP,STATICasel,s,loc,y,-zong_gao,-zong_gao DA,all,all,0allsel,allasel,s,loc,z,hp,hpnsla,snsel,r,loc,x,-L/2,L/2nsel,r,loc,y,-wa/2,wa/2D,all,VOLT,0allsel,allasel,s,loc,z,0,0nsla,snsel,r,loc,x,-L/2,L/2nsel,r,loc,y,-wa/2,wa/2D,all,VOLT,volt1allsel,allSOLVE/POST1/POST1/ESHAPE,1.0 !显示单元/POST1PLNSOL, U, SUM , 0,1.0 !显示总体的位移云图PLNSOL, U, X , 0,1.0 !显示X位移云图PLNSOL, U, Y , 0,1.0 !显示Y位移云图PLNSOL, U, Z , 0,1.0 !显示Z位移云图finiallselnsort,u,sum,0,0,all*get,max_u,sort,0,max*get,max_u,sort,0,minallselnsort,u,x,0,0,all*get,max_u,sort,0,max*get,max_u,sort,0,minallselnsort,u,y,0,0,all*get,max_u,sort,0,max*get,max_u,sort,0,minallselnsort,u,z,0,0,all*get,max_u,sort,0,max*get,max_u,sort,0,min。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

MPTEMP,1,0
MPDE,DENS,1
MPDATA,DENS,1,,7700
MPTEMP,1,0
MPDE,PERX,1
MPDATA,PERX,1,,826
MPDE,PERY,1
MPDATA,PERY,1,,919
MPDE,PERZ,1
MPDATA,PERZ,1,,919
!!!-------------------------------------------e矩阵,球左边系下，r方向为3方向
TBDE,ANEL,1,,,
TB,ANEL,1,1,21,1
TBTEMP,0
TBDATA,,11.087e-12,7.509e-12,7.509e-12,0,0,0
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,V1 !施加电压荷载
ALLSEL,ALL
!-------------------------------------------施加电压2
ASEL,S,LOC,X,r2
NSLA,R,1
*SET,V2,200
TBMODIF,4,3,
TBMODIF,5,1,0
TBMODIF,5,2,
TBMODIF,5,3,0
TBMODIFห้องสมุดไป่ตู้6,1,
TBMODIF,6,2,0
TBMODIF,6,3,4.32E-010
!!!-------------------------------------------d矩阵,球左边系下，r方向为3方向
!!!-------------------------------------------d矩阵,球坐标系下，r方向为3方向
TBDE,PIEZ,1
TB,PIEZ,1,,,1
TBMODIF,1,1,3.56E-010
TBMODIF,1,2,0
TBMODIF,1,3,
!-------------------------------------------r2外径
r1=7e-3
r2=10e-3
!-------------------------------------------球
cyl4,,,r1,,r2,90
VROTAT,all, , , , , ,3,2,90, ,
REAL,
ESYS,11
SECNUM,
!-------------------------------------------定义球坐标系，选择材料1，并弄将压电矩阵转置到球坐标系下
!-------------------------------------------设置压电球网格密度
!-------------------------------------------定义压电单元及参数
!-------------------------------------------建模
!-------------------------------------------r1内径
TBDATA,,12.035e-12,7.5179e-12,0,0,0,12.035e-12
TBDATA,,0,0,0,2.1053e-12,0,0
TBDATA,,2.1053e-12,0,2.2584e-12,,,
!!!-------------------------------------------e矩阵,球左边系下，r方向为3方向
CP,2,VOLT,ALL !耦合电压
*GET,N2,NODE,,NUM,MIN
!D,N2,VOLT,V2 !施加电压荷载
ALLSEL,ALL
CSYS,0
!----------------------------------施正弦加电压2，并求解
f=100000 !激励电压频率
TIMINT,ON,ALL !使用瞬态效应
!-------------------------------------------施加电压1,v1=0
CSYS,2
ASEL,S,LOC,X,r1
APLOT
NSLA,R,1
*SET,V1,0
CP,1,VOLT,ALL !耦合电压
SOLVE
*ENDDO
TIME,2*stepp*DT !再算5周期
NSUBST,50,,,1 !5周期，计算50步
KBC,0
SOLVE
FINISH
/POST26 !时间后处理器
NSEL,S,LOC,X,r2
NSEL,r,LOC,y,
NSEL,r,LOC,z,
!-------------------------------------------球
!-------------------------------------------定义球坐标系，选择材料1，并弄将压电矩阵转置到球坐标系下
CSWPLA,11,2,1,1,
TYPE,1
MAT,1
T=1/f
st=25 !每一个周期求解步数
DT=T/st
stepp=5*st !5个周期
*DO,J,1,stepp
t1=f*DT*J !5个正弦波
D,N2,VOLT,V2*SIN(2*3.1415926*t1) !施加电压荷载
TIME,DT*J !载荷时间步
NSUBST,1,,,1
*GET,N3,NODE,,NUM,MIN
NSOL,3,N3,U,X, UX_3,
PLVAR,3,
FINISH
/CLEAR,START
/PREP7
!-------------------------------------------定义压电单元及参数
!单元定义
ET,1,SOLID5
KEYOPT,1,1,3
KEYOPT,1,3,0
KEYOPT,1,5,0
!定义材料
CSYS,2 !球坐标系，X轴为r方向
p1=5 !=============压电球厚度方向单元数
LSEL,S,LOC,X,1.1*r1,0.9*r2
LESIZE,all, , ,p1, , , , ,1
LSEL,S,LOC,X,r1
LESIZE,all, , ,4*p1, , , , ,1 !=============压电球周向单元数4*p1=============
TBMODIF,2,1,-1.48E-010
TBMODIF,2,2,0
TBMODIF,2,3,
TBMODIF,3,1,-1.48E-010
TBMODIF,3,2,0
TBMODIF,3,3,
TBMODIF,4,1,0
TBMODIF,4,2,4.32E-010
NSLA,R,1
D,ALL, ,0, , , ,UZ, , , , ,
ALLSEL
!-------------------------------------------进入求解器
/SOLU
ALLSEL
ANTYPE,TRANS !指定分析类型
TRNOPT,FULL !指定瞬态分析方法
NLGEOM,OFF !-----关闭大变形选项
TINTP,,0.25,0.5,0.5 !压电瞬态分析中，瞬态积分参数
EQSLV,SPAR !算法选择，稀疏矩阵算法
OUTPR,ALL,ALL !控制结果的输出
OUTRES,ALL,2 !保存所有的数据到结果文件中，2步输出一个结果文件
!-------------------------------------------设置压电球网格密度
VMESH,all !画网格
allsel
!-------------------------------------------1/8球体位移边界条件
ASEL,S,LOC,X,0
APLOT
NSLA,R,1
D,ALL, ,0, , , ,UX, , , , ,
ALLSEL
ASEL,S,LOC,Y,0
APLOT
NSLA,R,1
D,ALL, ,0, , , ,UY, , , , ,
ALLSEL
ASEL,S,LOC,Z,0
APLOT