ABAQUS时程分析法计算地震反应的简单实例.docx

1、ABAQUS时程分析法计算地震反应的简单实例ABAQUS时程分析法计算地震反应的简单实例ABAQUS时程分析法计算地震反应的简单实例（在原反应谱模型上修改）问题描述:悬臂柱高12m，工字型截面（图 1），密度7800kg/m 3，EX=2.1e11 Pa，泊松比0.3，所有振型的阻尼比为2%，在3m高处有一集中质量 160kg，在6m、9m、12m 处分别有120kg的集中质量。反应谱按7度多遇地震，取地震影响系数为 0.08，第一组，山类场地,卓越周期Tg=0.45s。120k-第一部分：反应谱法几点说明:本例建模过程使用图1计算对象CAE ；*Spectrum 不可以在keyword ed

2、itor 中添加，keyword editor 不支持此关键词读入。ABAQUS的反应谱法计算过程以及后处理要比 ANSYS方便的多。操作过程为:(1)打开 ABAQUS/CAE，点击 create model database(2)进入Part模块，点击create part ，命名为column,3D、deformation 、wire。con ti nue Create FartHanw； cl untilModeling Sjace 3D 0 2T Planar O Ajii sjTnmetr i cTypftOptions(5)De forfli ableDiscrete ri韦id

3、 0血41亦诅tigi A O Euler i anNone availableBase FeatureOs皿4Q ShellO FointPl anarApprsxLUiteCflntimue.(3) Create lines，在E I Xj Pick a tar point for the enter Y：分别输入0, 0回车；0, 3回车；0, 6回车；0, 9回车；0, 12回车。mass den sity 7800 Hame： steelDensi tyDeleteGeneral He ctiani cal Thermal OtherDensityI I Us* I aif ara.

4、ture-dp nd4n t dataNwnlsr of field variables；Eat a S Dcnsilyi.jvaoaOKmecha ni cal-elasticity-elasticCanctl ,young s modulus : 2.1e11 , poisson sratio : 0.3. Edit lat erial甌terial Behavi oriDeiisi tyElasti ctime scale (for viscoelastic i ty) . LongtermI I Mo comprtssionI I teiisinData5Iodol as2 Iftll

5、PoiSDA St 耻jt.io03OKCancel name : Section-1 , category : beam , type : beam,(5) Create secti onCon ti nuecreate p rofile, n ame: Profile-1, sha pe:l.:邃昏Type： BeaonSe ction1 ntegr ation rGDuring analysiPrafilename：Frofiliships7p 讎匚ra盘te. . i H ame:Profile-1StiffnasEFluid litertiShapeV 匚raata. , Sect!

6、on Foisson z ratio： 0Temperature Tariitiort：（3）Line ST By gridi entsQ Interpfllted fror temperature pointsi qi! FireICircularRa 吐 ugolarIITrpsriiifilILITArbitraryIG-enar all i e J.I Cantiiuie. C蚯c&lOE（Cancel Edit ProfileFane- Frcil*-1在 profile name 选择 I, material name 选择 steel 。 Ok Edit Dea* Section

7、Ham： eutinTType： BaniStctifci iid电詐atioci； During ui辽yiis O E*f*re 畅aJyaiiFt血 ile st ape r I cSect ion foi ssor/ s ratio Z - CTwiprfcVurt vari *ti on：Li near ty 歹aii 电ntsQ Interpplt*d frm tmpratiur4 pointsCane* (6) Assign section ，选择全部，done，弹出的对话框选择 section : Section-1ok。 EditR gi OILE&gi fin.Sect!o

8、nMote； List cQjitaiiis onJy sectj pus apjlLcaljla to the selected regipnE.Hateriel sttel Caiigel (7)Assig n beam orie ntati onI，选择全部，默认值确定。(8)View-part dis play op ti ons，在弹出的对话框里勾选， ren der beamprofiles ,以可视化梁截面形状。 Part Display OptionsRender StyleO 爭i住frame O Hidden (V) ShadedGe time try叼 Eh口w lodg

9、es in 寸r ecider style 冈 SEow Eiltouel te tdgesOHighlight only visiblc antitieFace tii 纳ligMirt百：|5tij)pLiTig vCurve refinement： ICoarsaIdle： Th亡 refinement setting *111 be plied only to the current part. Show reference representationNashShow : i Hk terior ed-es |回 Show 屯in shided r eniei* styleII Hi

10、 畑igM only vi sibl tntili tsO shell ttuclmezsOK L 帥ly I Def包nils CaaceL mass/i nertia(9 )添加集中质量，Sp ecial-i nertia-create,name : massl , type : pointcontinue，选择（0, 3）位置点done , mass : 160 , ok。 Edit InertiaKeel on： (f LtkeiJIMiigni Eaftping |Mass :伯5Rotary Iiker+i.a.OSpecify off-diagoTial termsILU12Z:

11、133： I CSYS: (Global)跖 J J 丄，匚 1电直 I 直jlote： Values nfill Iv applied pr point.OEcreate , name : mass2 , type : point mass/inertia , continue ，选择 0 , 6; 0, 9;0, 12 位置点(按 shift 多选)，done , mass : 120 , ok , dismiss。ok。,name : step-1 ,procedure type 选 freqency , continue(11) Step-create step在basic选项卡中，e

12、igensolver 选择频率提取方法， 本例选用lanczos 法，number of eigenvalues request ，选 value，输入 10.ok 。 Edit StppKame: SteplTyph FrequtncyBasi c OtherDescription, frtquftncy analysis of columnNlcQm: OffELgezisalver LaneiQs O SuBspace 0 AMSNumber of ei geTivttlues requested： All in frequency rw Value： : 10Maxinufh numb

13、 ar of block Luiczok stpi： Dfult O Vlin：| Use Slfi-baied linear tiynami cs procedir esCancel jI1 Iikcludo r QS m4 J4-5OjJ再 create step , create step , name : step-2 , procedure type 选 response spectrumcontinue在 basic 选项卡中，excitati ons 选择单向 sin gle direct ion , sumati ons 选择 squareroot of the sum of

14、squares (SRSS) 法，use response spectrum : sp(反应谱的 name , 后面再inp中添加），方向余弦（0, 0, 1） , scale factor : 1.r Edit 5tepSame： 2讥旷Typ： Risponss ipclru(n0T:01:0.02.ok。Use response spectrum： sp.-econ-j d; rec . i or. , . - Direction CQSin?. K；,?hi rd d:re：t打匸 rScale factor：何K：Cancel j进入 damping 选项卡，阻尼使用直接模态( di

15、rect modal )，勾选 direct dampingdata , start mode : 1 , end mode : 8, critical damping fraction Edit StppMarne: Step-ZTyp; RtEponie ipctrivnI Basic Dampinf |fy dsiiipiik wvci i 血 of： = qiieiic i e iM Creat e BoimdaryKam*: ixdStep : Ini ti alO Otherdisplacement vielcc jty accclra.ti onContinue.匚ait住1 选

16、择0, 0点，done，勾选u1ur3 所有6个自由度。Ok。 Edit IQoundary ConditionN劭4 ；Type：D1 spL a c em e nt/E. tatSlap :Ini ti!Kegi ori：CFicke d)CSIE:(Global) |Edit J1 g0也 回協0旳0 imi回叱0邂Th ft digpLax 酗电 nt v&lue will la fh&int9.ine4 m subsesttjs.OK(13)进入mesh 模块，object 选择 part , 点 seed edge by number ，选择所有杆，B tii t Constrai

17、nts-El .sO 助 5de Sone O O 贏讪L直 mimberSizing CanLfolsHuiubtr of elements：2屮Cr e*t on1 Create wi th Ttane : Edge SeeJs-lK Apply Jef鱼alts 匚aiic&l 点 assign element type，选择全部杆，done，默认B31 , ok。点mesh part ,yes。文件。进入 job manager ，点击 write input ，在工作目录生成 demo-spc.inp(15)进入ABAQUS工作目录，使用 UltraEdit软件(或其他类似软件)打开

18、demo-spc.inp ，*Boundary 关键词的后面加如下根据问题叙述确定的反应谱:*Sp ectrum,t yp e=accelerati on,n ame=s p 0.1543,0.167,0 0.1915,0.25,0 0.2102,0.333,0 0.2241,0.444,0 0.25,0.5,0 0.3295,0.667,0 0.4843,1,0 0.5987,1.25,0 0.7868,1.667,0 1.0342,2.222,0 1.0342,10,0 0.3528,10000,0第一列为加速度，第二列为频率，第三列为阻尼比。玖今画 J已閘mi翼或ai社畛T-gj,赴a迢

19、用ijT ms号n屯保存。(16)进入 job 模块，create job ,name : spc, source 选择 input file ,nput file select Czeate JobSource： iTiput file vL ContinuE一亠. CaneeI 默认，ok，进入job manager ，选择 spc , submit，计算成功!Freque ncy must be in creas ing continu ously in a sp ectrum defi niti on士讲解的ansys结果3.1611cm 基本一致。可以查看工作目录下的 spc.dat

20、文件查看详细的频率和模态分析结果。第二部分：时程分析(1)进入 steP 模块，删除原 step1、step2。建立 step1 (static general )，用于施 加重力 Create St ep(2)将step1结果作为动态分析的初始状态，time period 设置为1e-10 (很短时间)。 Edit StepName： Step1BasicTypt： Static丿 GvittralIiuzrementaii on II Other IDescription：Time period： | le-lQ Qhis stttinc controls th inclufiQii of

21、 RQulinvur afftcti0 nf IhtEp di acAmPTht 幷nxi 7+ pp )c stabilization： NoneI I Includs aditibfttic hfalin專 affectsCancel 建立ste p2 (dyn amic imp licit )，进行动力时程分析X， Crsafte StepHusc： Step-Insert tew step aftrHni 11 alProced.Tir& type： GeneralCoupled t emp * di 3pl acementDirect cyclicjeostati cSoilsSta

22、tic, GeneralStatic, Bik: Visea大增量数量设置为 2000步，将初始时间增量设置为 0.02，最小增量设置为 1e-15，最大增量设置为 0.02 , half-step residual toleranee:100 （控制automatic 求解精度的值,在地震分析中应该设置多大为好？还没弄清楚！请大家赐教！）。Hame： St ep-2Typ#： Tlynftftic, ImplicitIncremei-t zi ze Z 0. H松cirvim incr电iiEt： O Anilysii appLic&ti&nJ) Specify： D. 02Hslf-inc

23、reneiit Etiidual Supir ass calculationVote： Hay h automatically suppressed application is not set to transient fidli tyO Analysis product defaultTolset-create（在tools中设置，用于观察顶端节点的反应情况），同样的方法，底端节点设置为set-2 Create Set: IbSBIType Geometry O Hodt O EL电niitVaf&iiLg： Falivs nod* gd Igmantwi11 invslidalsdif

24、the mtahI Contiivu电. Cancel在out put 中设置需要输出结果，在 edit history out put request将 domain 改为 set，选择set1，在displacement 里面选择 U。I Hifftory Output Requests onagerfuHn-ntput- 1Stej-1St*p-2Pr opagat&iSlip procedure : Static GiL*rilVari iblts： c+ei JefaiiltfiStatus； Create! in this step卜 Contact0 Energyo Failur

25、c/Fr acture卜 ThermalQ ElectricalO Volune/Thi ckneis/Coorlinates卜 I I Mass center kixtemalicsI I。叭pat rttr(Jutput at shell leam_, and layered section points：Use defiults OSpecify： |1 QElautal )Creat H-Output-2 ，选择 set-2，同上在材料中补充dam ping ，使用瑞利阻尼，质量系数alpha为0.15，刚度系数beta为0.01。Xl Edit KatorialName： steel

26、d. Bihayi ors1 Damping 1Densi tyElastic1 |P 心 tvDunpiikg0.01Covpp te；StructTiral:Cancel ，将 fix 在 step2 的 propagated(3)进入 load 模块，boundary condition manager改为 in active (点击 deactivate )create 一个新的边界条件（在 step2 ），取消z向位移约束（以在该方向施加加速度） Create Boundary ConditionVel0city/Angar velocityAccalerat 1911/Angula

27、r 4ccelaratioR Connection ii 电ntXFelodi tyutlCoimectOTCoiuiectirCancelConHinig. . . IJ Edit Boundary ConditionName;BC-2Ty严；D is pl a 匚 ement/RotationStep;Step-2 (SnarmiG Irnplidt)Region;(Picked) Edit Region,hMethod; Specify Constr日irrtsU2; U3;Notfr： The displaceme-rrt value will hemaintained in subsequent step5.再 create 一个边界条件 (在 step2 ) ,ype 为 acceleration/angular acceleration,continue Create Boundary Conditiono otherCoimectof lisplacement velocityChTiIiftCt OF iV Lt i Oik匚 ancelContiiHie.选择基底节点，勾选 A1，输入1 （加速度记录单位 m/s2 ）,可 Edit Boundary C