(完整版)ANSYS钢筋混凝土分离式建模
如何在ANSYS中模拟钢筋混凝土的计算模型
如何在ANSYS中模拟钢筋混凝土的计算模型最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。
就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。
就ANSYS而言，其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。
ANSYS 钢筋混凝土建模
ANSYS 钢筋混凝土建模一、简介钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。
Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。
希望大家一起讨论、批评指正(******************.cn)。
程序：ANSYS单元：SOLID65、BEAM188建模方式：分离暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。
二、单元选择以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。
最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。
只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。
建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。
暗支撑剪力墙数值模型看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。
但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计算提前发散。
LINK8+SOLID65的问题如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。
钢筋混凝土建模步骤
钢筋混凝土建模步骤在土木工程结构中，最为常用的一种结构形式就是钢筋混凝土结构，在各类房屋、水坝、桥梁、道路中都有广泛应用。
ANSYS软件提供了专门的钢筋混凝土单元和材料模型。
本算例将介绍ANSYS软件分析混凝土一些基本应用。
(1) 首先建立有限元模型，这里我们选用ANSYS软件自带的专门针对混凝土的单元类型Solid 65，进入ANSYS主菜单Preprocessor->ElementType->Add/Edit/Delete，选择添加Solid 65号混凝土单元。
(2) 点击Element types窗口中的Options，设定Stress relax after cracking 为Include，即考虑混凝土开裂后的应力软化行为，这样在很多时候都可以提高计算的收敛效率。
(3) 下面我们要通过实参数来设置Solid 65单元中的配筋情况。
进入ANSYS主菜单Preprocessor-> Real Constants->Add/Edit/Delete，添加实参数类型1与Solid 65单元相关，输入钢筋的材料属性为2号材料，但不输入钢筋面积，即这类实参数是素混凝土的配筋情况。
(4) 再添加第二个实参数，输入X方向配筋为0.05，即X方向的体积配筋率为5%。
(5) 下面输入混凝土的材料属性。
混凝土的材料属性比较复杂，其力学属性部分一般由以下3部分组成：基本属性，包括弹性模量和泊松比；本构关系，定义等效应力应变行为；破坏准则，定义开裂强度和压碎强度。
下面分别介绍如下。
(6) 首先进入ANSYS主菜单Preprocessor-> Material Props-> Material Models，在Define Material Model Behavior 窗口中选择Structural-> Linear -> Elastic-> Isotropic，输入弹性模量和泊松比分别为30e9和0.2(7) 下面输入混凝土的等效应力应变关系，这里我们选择von Mises屈服面，该屈服面对于二维受力的混凝土而言精度还是可以接受的。
(完整版)ansys钢筋混凝土梁的建模方法约束方程法
用约束方程法模拟钢筋混凝土梁结构问题描述建立钢筋线对钢筋线划分网格后形成钢筋单元bhP 位移载荷建立混凝土单元对钢筋线节点以及混凝土节点之间建立约束方程后施加约束以及位移载荷进入求解器进行求解；钢筋单元的受力云图混凝土的应力云图混凝土开裂fini/clear,nostart/config,nres,5000/filname,yue su fang cheng 5 jia mi hun nin tu /prep7/title,rc-beamb=150h=300a=30l=2000displacement=5!定义单元类型et,1,solid65et,2,beam188et,3,plane42!定义截面类型sectype,1,beam,csolid,,0secoffset,centsecdata,8,0,0,0,0,0,0,0,0,0sectype,2,beam,csolid,,0secoffset,centsecdata,4,0,0,0,0,0,0,0,0,0!定义材料属性，混凝土材料属性mp,ex,1,24000mp,prxy,1,0.2tb,conc,1,1,9tbdata,,0.4,1,3,-1!纵向受拉钢筋mp,ex,2,2e5mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350!横向箍筋，受压钢筋材料属性mp,ex,3,2e5mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!生成钢筋线k,,k,,bkgen,2,1,2,,,hk,,a,ak,,b-a,akgen,2,5,6,,,h-2*akgen,21,5,8,,,,-100 *do,i,5,84,1l,i,i+4*enddo*do,i,5,85,4l,i,i+1l,i,i+2*enddo*do,i,8,88,4l,i,i-1l,i,i-2*enddo!受拉钢筋lsel,s,loc,y,alsel,r,loc,x,alsel,a,loc,x,b-a lsel,r,loc,y,acm,longitudinal,line type,2mat,2secnum,1 lesize,all,50lmesh,allallscmsel,u,longitudinalcm,hooping reinforcement,line!箍筋，受压钢筋type,2mat,2secnum,2lesize,all,50lmesh,all/eshape,1!将钢筋节点建为一个集合cm,steel,node!生成面单元，以便拉伸成体单元a,1,2,4,3lsel,s,loc,y,0lsel,a,loc,y,hlesize,all,,,10lsel,alllsel,s,loc,x,0lsel,a,loc,x,blesize,all,,,20type,3amesh,all!拉伸成混凝土单元type,1real,3mat,1extopt,esize,30extopt,aclear,1vext,all,,,,,-lalls!建立约束方程cmsel,s,hooping reinforcement cmsel,a,longitudinalnsll,s,1ceintf,,ux,uy,uzallsel,all!边界条件约束nsel,s,loc,y,0nsel,r,loc,z,0d,all,uyd,all,uxnsel,s,loc,y,0nsel,r,loc,z,-ld,all,uyd,all,ux!施加外部荷载/solunsel,allnsel,s,loc,y,hnsel,r,loc,z,-1000d,all,uy,-displacement alls!求解nlgeom,on nsubst,200 outres,all,all neqit,100pred,oncnvtol,f,,0.05,2,0.5 allselsolvefinish/post1allselplcrack,0,1plcrack,0,2!时间历程后处理/post26nsel,s,loc,z,-l/2*get,Nmin,node,0,num,min nsol,2,nmin,u,yprod,3,2,,,,,,-1nsel,s,loc,y,0nsel,r,loc,z,0*get,Nnum,node,0,count *get,Nmin,node,0,num,min n0=Nminrforce,5,Nmin,f,y*do,i,2,ndinqr(1,13)ni=ndnext(n0)rforce,6,ni,f,yadd,5,5,6n0=ni*enddoprod,7,5,,,,,,1/1000/axlab,x,uy/axlab,y,p(kn) xvar,3 plvar,7。
ANSYS钢筋混凝土建模方法概述
利用大型通用有限元软件ANSYS进行钢筋混凝土结构的建模、计算分析、结果处理是目前针对钢筋混凝土进行数值模拟的重要步骤。
如何采用ANSYS进行钢筋混凝土建模，能否把握有限元模型的可行性、合理性是将有限元理论应用到实际工程中较为关键的一环。
按照目前在建模中对钢筋的处理方式，ANSYS钢筋混凝土建模方法主要分为三种：整体式、分离式以及组合式，每种方法都具有不同的建模特点，现略做总结如下。
一、整体式建模ANSYS采用Solid65单元来模拟混凝土,所谓整体式建模也即是在建模过程中，通过对65单元进行实常数的设置来考虑钢筋对混凝土结构的作用。
这种方法将钢筋弥散于整个单元中，并视单元为连续均匀材料。
与其他方法比较，整体式建模的单元刚度矩阵综合了钢筋和混凝土单元的刚度矩阵，并且是一次性求得综合的刚度矩阵。
因此，在采用整体建模方法时，在建模之前，应首先求得单元各个方向的配筋率，并设置实常数，一般适用于体量较大，配筋比较规整的钢筋混凝土结构。
整体式建模所得计算结果对比实验来讲，其计算的开裂荷载误差较小，但开裂荷载后的整体荷载位移曲线与实验相比误差较大。
但采用整体建模方法的主要好处是能有效避免因为单元细分导致的应力奇异问题，有利于提高整体计算的收敛性性能。
二、分离式建模与整体式建模方法不同，分离式建模是指在建模过程中，考虑钢筋与混凝土的相互作用，分别选用不同的单元来模拟钢筋和混凝土。
一般而言，钢筋采用线单元link8模拟，混凝土选用配筋率为0的素混凝土Solid65单元模拟。
由于采用不同单元建模，如果认为结构在受外部荷载作用时，钢筋与混凝土在相互约束情况下会产生相对滑移，这时可以在钢筋与混凝土之间添加粘结单元来模拟钢筋与混凝土之间的粘结与滑移，一般采用非线性弹簧conbin39。
如果认为两者之间连接紧密，不会出现滑移，可视为刚性连接，只需通过合并节点即可，也即是相当于两者节点耦合。
从上述表述可见，分离式模型可以揭示钢筋与混凝土之间相互作用的微观机理，而这也是整体式模型无法做到的。
ANSYS分析钢筋混凝土结构技巧及实例详解
0 前言利用ANSYS分析钢筋混凝土结构时，其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性，从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧，并以钢筋混凝土简支梁为例，采用分离式有限元模型，说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型)，MISO模型可包括20条不同温度曲线，每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型)，MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
ANSYS中混凝土的计算问题(分离式)解析
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。
就ANSYS而言，可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。
就ANSYS而言，其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。
而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。
钢筋混凝土简支梁模拟计算 ansys命令流
lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3பைடு நூலகம்00
lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0
ksel,r,loc,y,0
!方法：分离式；solid65和link8
!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!===================================================
lesize,all,,,20
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
finish
/solu
!施加约束
/config,nres,2000
/prep7
!定义单元及其材料特性等
rd0=20.0 !钢筋直径
et,1,solid65
et,2,link8
mp,ex,1,33e3
mp,prxy,1,0.20
r,1
hntra=28
hntrl=2.6
tb,concr,1
/post1
pldisp,1
用ANSYS建立钢筋混凝土梁模型
用ANSYS1立钢筋混凝土梁模型问题描述：钢筋混凝土梁在受到中间位移荷载的条件下的变形以及个组成部分的应力情况。
P=5mm位移L=2000mm图1钢筋混凝土结构尺寸图一、用合并节点的方法模拟钢筋混凝土梁1 .用solid65号单元以及beam188单元时材料特性钢材的应力应变关系混凝土的弹性模量采用线弹性TEMP建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元合并单元节点后施加约束以及位移载荷进入求解器进行求解钢筋单元的受力云图.4Q2F-0375.269 150,536 225.806 301 . 075 混凝土的应力云图混凝土开裂2使用单元solid45号单元与beam188钢筋的应力应变关系不变，而混凝土应力应变关系为：混凝土单元WK.355713 5.067 11.37S3,11119.644EPS3钢筋单元 力与位移曲线 .13257E弓・51611-0^^9uOS31^-652Q B 2L7 25-7S4~H ・ 793・0190363E a 9Q715&•号E 】 233.34311.113194 ・453272 ・227350p(kn)g105uy、用约束方程法模拟钢筋混凝土梁1 .用solid65号单元以及beam188单元时混凝土以及钢筋采用线弹性关系: 建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元对钢筋线节点以及混凝土节点之间建立约束方程WFOR.NMQMBFOR后施加约束以及位移载荷进入求解器进行求解；钢筋单元的受力云图MN011905 77.787 155.562 233.337 311.11238.899 116.675 194.45 272,225 350 混凝土的应力云图,1472166,969 13,79 20.612 27.43310.379 17 ・201 24 ・ DEE30.844混凝土开裂*11111112使用单元solid45号单元与beam188使用混凝土的本构关系曲线1255 7 9钢材的本构关系曲线钢筋的von mises应力.116491 77+86830.992155,62116.744 1.94,496233+372 311+124272-246 350混凝土的应力.1287895,695 11,662 17.429 23.195 3-012 B.779 14.545 20*312 26.079用在solid45号单元下，用合并节点法、约束方程法建立模中钢筋与混凝土之间的关系 的时候的一个力与位移全程曲线的比较。
钢筋混凝土分离式建模
一、简介钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。
Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。
希望大家一起讨论、批评指正(******************.cn)。
程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。
二、单元选择以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。
最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。
只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。
建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。
[center]暗支撑剪力墙数值模型[/center]看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。
但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计算提前发散。
[center]LINK8+SOLID65的问题[/center]如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
!跨中施加110KN的集中力
FINISH $/CLEAR $/PREP7!
AS0=380.1 $AS1=50.3 $A=30 $B=150!
H=300 $L=2650 $L0=125!
ET,1,SOLID65!
KEYOPT,1,1,0!
KEYOPT,1,5,1!
KEYOPT,1,6,3!
KEYOPT,1,7,1!
ET,2,LINK180!
ET,3,SOLID185,,3
R,1,AS0 $R,2,AS1 $R,3
MP,EX, 1,2.4E4 $MP,PRXY,1,0.2 $FC=25!
TB,CONCR,1,1,9!
TBDATA,,0.35,0.75,3.1125,-1!
TB,MISO,1,,15!
TBPT,,0.0002,4.8 $TBPT,,0.0004,9.375 $TBPT,,0.0006,13.51! TBPT,,0.0008,17.02 $TBPT,,0.001,19.83, $TBPT,,0.0012,21.95! TBPT,,0.0014,23.43 $TBPT,,0.0016,24.365 $TBPT,,0.0018,24.856! TBPT,,0.002,FC $TBPT,,0.0038,FC !
TBPLOT!
MP,EX,2,2E5 $MP,PRXY,2,0.25!
TB,BKIN,2 $TBDATA,,360!
MP,EX,3,2E5 $MP,PRXY,3,0.25!
TB,BKIN,3 $TBDATA,,210!
N,1,,B $N,9 $FILL,1,9!
NGEN,11,9,1,9,1,,,A!
NGEN,2,1000,1,99,1,75!
NGEN,3,1000,1001,1099,1,50!
NGEN,7,1000,3001,3099,1,75!
NGEN,4,1000,9001,9099,1,200/3!
NGEN,7,1000,12001,12099,1,75!
NGEN,2,1000,18001,18099,1,50!
/VIEW,1,-1,-1,1!
TYPE,2 $REAL,2 $MAT,3!
*DO,II,11,16,1 $E,II,II+1 $*ENDDO!
*DO,II,83,88,1 $E,II,II+1 $*ENDDO!
*DO,II,11,74,9 $E,II,II+9 $*ENDDO!
*DO,II,17,80,9 $E,II,II+9 $*ENDDO!
EGEN,20,1000,1,28,1!
*DO,II,83,18083,1000 $E,II,II+1000 $*ENDDO!
*DO,II,89,18089,1000 $E,II,II+1000 $*ENDDO!
TYPE,2 $REAL,1 $MAT,2
*DO,II,11,18011,1000 $E,II,II+1000 $*ENDDO!
*DO,II,17,18017,1000 $E,II,II+1000 $*ENDDO!
/ESHAPE,1 $EP!
BLC4,,,L/2,B,H $BLC4,75,,100,B,-40 $WPOFFS,,,H $BLC4,625,,200,B,40 $WPCSYS,-1! WPOFFS,75 $WPROTA,,,90 $VSBW,ALL!
WPOFFS,,,100 $VSBW,ALL $WPOFFS,,,450 $VSBW,ALL!
WPOFFS,,,200 $VSBW,ALL $WPOFFS,,,450 $VSBW,ALL!
WPCSYS,-1 $ALLSEL!
LSEL,S,LOC,Y,0 $LSEL,A,LOC,Y,150 $LSEL,R,LOC,X,0 $LESIZE,ALL,,,10!
LSEL,S,LOC,Z,0 $LSEL,A,LOC,Z,300 $LSEL,R,LOC,X,0 $LESIZE,ALL,,,8!
LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,0,75 $LESIZE,ALL,75!
LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,75,175 $LESIZE,ALL,25!
LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,175,625 $LESIZE,ALL,75/2!
LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,625,825 $LESIZE,ALL,100/3! LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,825,1275 $LESIZE,ALL,75/2! LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,1275,1325 $LESIZE,ALL,50!
LSEL,S,LOC,Z,340 $LSEL,R,LOC,X,625 $LESIZE,ALL,,,8!
LSEL,S,LOC,Y,0 $LSEL,R,LOC,X,625 $LSEL,R,LOC,Z,300,340 $LESIZE,ALL,,,1!
LSEL,S,LOC,Z,-40 $LSEL,R,LOC,X,75 $LESIZE,ALL,,,8!
LSEL,S,LOC,Y,0 $LSEL,R,LOC,Z,-40 $LESIZE,ALL,,,4!
LSEL,S,LOC,Y,0 $LSEL,R,LOC,X,75 $LSEL,R,LOC,Z,0,-40 $LESIZE,ALL,,,1!
VSEL,S,LOC,Z,0,H $VATT,1,3,1
MSHAPE,0,3D $MSHKEY,1
VMESH,ALL $ALLSEL
VSEL,S,LOC,Z,-40,0 $VSEL,A,LOC,Z,H,H+40 $VATT,2,3,3
MSHAPE,0,3D $MSHKEY,1 $VMESH,ALL
/VIEW,1,-0.2,-1,1 $EPLOT $ALLSEL
NUMMRG,ALL $NUMCMP,ALL $EPLOT
/SOLU!
NSEL,S,LOC,Z,-40 $NSEL,R,LOC,X,L0 $D,ALL,UY,,,,,UZ $ALLSEL!
ASEL,S,LOC,X,L/2 $DA,ALL,SYMM $ALLSEL!
NSEL,S,LOC,Z,H+40 $NSEL,R,LOC,X,725!
*GET,NODE1,NODE,,COUNT $F,ALL,FZ,-110000/NODE1 $ALLSEL!
ANTYPE,STATIC $NLGEOM,ON $NSUBST,80 $OUTRES,ALL,ALL $AUTOTS,1 $LNSRCH,1! CNVTOL,F,,0.05,2,0.5 $ALLSEL!
SOLVE $FINISH
/POST1 $SET,LAST $PRRSOL,FZ!
SET,LAST $PLDISP,1!
ESEL,S,TYPE,,2 $ETABLE,SAXL,LS,1 $PLLS,SAXL,SAXL!
ESEL,S,TYPE,,1 $/DEVICE,VECTOR,ON $PLCRACK!
/POST26 $NSOL,2,NODE(L/2,B/2,0),U,Z!
PROD,3,2,,,,,,-1 $PROD,4,1,,,LOAD,,,110!
/AXLAB,X,MID-UZ(MM) $/AXLAB,Y,P(KN)!
XVAR,3 $PLVAR,4。