一起学习板屈曲分析.doc

一起学习板屈曲分析simwe会员名：LXCAD一，材料为铁，适合做模板分析，动态显性,物质的非线形指南10.5 blast loading on a stiffened plate*Amplitude, name=Blast, smooth=0.250., 0., 0.001, 700000., 0.01, 700000., 0.02, 0.0.05, 0.* * MATERIALS*材料明确为铁 *Material, name=Steel*Density7800.,*Elastic 2.1e+11, 0.3* Plastic关键字规定使用Misess屈服准则*Plastic 3e+08, 0. 3.5e+08, 0.025 3.75e+08, 0.1 3.94e+08, 0.2 4e+08, 0.35* * BOUNDARY CONDITIONS* * Name: Fix edges Type: Symmetry/Antisymmetry/Encastre*BoundaryEdge, ENCASTRE* -* * STEP: Blast* *Step, name=BlastApply blast loading*显性动态分析*Dynamic, Explicit, 0.05*Bulk Viscosity!修改屈曲粘性参数，均使用缺省值0.06, 1.2* * LOADS* * Name: Pressure load Type: Pressure*Dsload, amplitude=Blast_PickedSurf6, P, 1.* 二，箱型结构使用动态求特征值例子手册2.2.4中，Eigenvalue analysis of a structure using the parallel Lanczos eigensolver*PARAMETER* modal input parametersnumModes = 200lowFreq = 0.0highFreq = 133.0*SHELL GENERAL SECTION, ELSET=ALL, MATERIAL=STEEL0.01*MATERIAL, NAME=STEEL*ELASTIC200.E9, 0.3*DENSITY7800.0* Step definition*STEPFrequency extraction*FREQUENCY, NUMBER INTERVAL=4,*BOUNDARYBOTLEFT, ENCASTREBOTRIGHT, ENCASTRE*output, field, freq=0*END STEP三，静力分析(屈曲崩塌分析 关键字*STATIC,STABILIZE的使用)例字手册1.2.5Unstable static problem: reinforced plate under compressive loads*ELEMENT,TYPE=SPRING1,ELSET=EPONT 50000,340 50001,2881*弹簧可模拟拉条*SPRING,ELSET=EPONT 3, 5.3E+08,*BEAM SECTION,MATERIAL=ACIER,SECTION=RECT,ELSET=E0000001 1.500000E-02, 1.500000E-02 0.000000E+00, 6.557236E-03,-9.999785E-01*MATERIAL,NAME=ACIER *ELASTIC 2.100E+11, 3.000E-01*PLASTIC 235.259E+06,0.0 471.9E+06,0.18837*DENSITY 7.850E+03, *BEAM SECTION,MATERIAL=ACIER,SECTION=RECT,ELSET=E0000002 1.500000E-02, 1.500000E-02 0.000000E+00, 8.695324E-03,-9.999622E-01*SHELL SECTION,ELSET=E0000003,MATERIAL=ACIER 7.000E-03, 3*SHELL SECTION,ELSET=E0000004,MATERIAL=ACIER 1.000E-02, 3*SHELL SECTION,ELSET=TOLE,MATERIAL=ACIER 5.000E-03, 3*SHELL SECTION,ELSET=E0000006,MATERIAL=ACIER 6.000E-03, 3*SHELL SECTION,ELSET=E0000007,MATERIAL=ACIER 7.000E-03, 3*SHELL SECTION,ELSET=E0000008,MATERIAL=ACIER 1.000E-02, 3*SHELL SECTION,ELSET=E0000009,MATERIAL=ACIER 1.000E-02, 3*SHELL SECTION,ELSET=E0000010,MATERIAL=ACIER 1.000E-02, 3*NSET,NSET=NM1 329, *ELSET,ELSET=HILOIRE E0000004,E0000009,E0000010*ELSET,ELSET=BARROTS E0000003,E0000007,E0000008*ELSET,ELSET=TOT TOLE,HILOIRE,BARROTS,E0000001,E0000002,E0000006*MPCBEAM,COTE1,NM1*BOUNDARYNM1,3,3NM1,4,4NM1,6,6COTE2,1,1COTE2,3,3COTE2,4,4COTE2,6,68896,2,2COTE3,2,2COTE3,4,4COTE3,6,6COTE4,2,2COTE4,4,4COTE4,6,6*RESTART,WRITE,F=30* *STEP,NLGEOM,INC=1GRAVITY LOAD*STATIC1.,1.*DLOADTOT,GRAV,9.81,0,0,-1.*NODE FILE,NSET=NM1,FREQ=1U,CF*EL FILE,ELSET=HILOIRE,FREQ=200S,E,PE *ENERGY FILE *OUTPUT,FIELD,FREQ=30*NODE OUTPUT U,*OUTPUT,HISTORY*NODE OUTPUT,NSET=NM1U,CF*ENERGY OUTPUT,VARIABLE=ALL*EL PRINT,FREQ=0*NODE PRINT,FREQ=0*END STEP*STEP,NLGEOM,INC=1000BUCKLING LOAD*STATIC,STABILIZE0.1,1.,*CLOADNM1,1,646.E+04*END STEP四，静力分析，焊点与接触的使用例子1.2.3Buckling of a column with spot welds*SHELL SECTION,ELSET=WPIL,MATERIAL=STEEL, SECTION INTEGRA=GAUSS 0.002,*SHELL SECTION,ELSET=BPIL,MATERIAL=STEEL, SECTION INTEGRA=GAUSS 0.003,*MATERIAL,NAME=STEEL*DENSITY 7850.,*ELASTIC 206.8E9,0.3*PLASTIC170.0E6, 0.0000000E+00180.0E6, 1.7205942E-03190.0E6, 3.8296832E-03200.0E6, 6.3897874E-03210.0E6, 9.4694765E-03220.0E6, 1.3143660E-02230.0E6, 1.7493792E-02240.0E6, 2.2608092E-02250.0E6, 2.8581845E-02260.0E6, 3.5517555E-02270.0E6, 4.3525275E-02280.0E6, 5.2722659E-02290.0E6, 6.3235357E-02300.0E6, 7.5197279E-02310.0E6, 8.8750519E-02320.0E6, 0.1040458330.0E6, 0.1212430340.0E6, 0.1405106350.0E6, 0.1620263360.0E6, 0.1859779370.0E6, 0.2125620380.0E6, 0.2419857390.0E6, 0.2744660400.0E6, 0.3102303410.0E6, 0.3495160420.0E6, 0.3925720430.0E6, 0.4396578440.0E6, 0.4910434450.0E6, 0.5470111460.0E6, 0.6078544470.0E6, 0.6738777480.0E6, 0.7453985490.0E6, 0.8227461500.0E6, 0.9062610510.0E6, 0.9962980*NSET,NSET=BASE,GEN 101,129,1 200,214,1*BOUNDARY BASE,ENCASTRE*SURFACE,TYPE=ELEMENT,NAME=BPIL BPIL1, *SURFACE,TYPE=ELEMENT,NAME=WPIL WPIL1, *RIGID BODY,ELSET=ROOF,REF=60000*NSET,NSET=WELDREFNODES,GEN 5203,45203,10000 5211,45211,10000*ELEMENT, TYPE=CONN3D2, ELSET=FASTCONN1000000,5203,1052031000001,15203,10152031000002,25203,10252031000003,35203,10352031000004,45203,10452032000000,5211,1052112000001,15211,10152112000002,25211,10252112000003,35211,10352112000004,45211,1045211*ORIENTATION, NAME=GLOBAL1.,0.,0.,0.,1.,0.*CONNECTOR SECTION, ELSET=FASTCONN, BEHAVIOR=BEHV,ELIMINATION=NOCARTESIAN,CARDANGLOBAL,*CONNECTOR BEHAVIOR, NAME=BEHV*CONNECTOR DERIVED COMPONENT, NAME=normal3,1.0*CONNECTOR DERIVED COMPONENT, NAME=normal4,5.3333,.3333*CONNECTOR DERIVED COMPONENT, NAME=shear6,.5,*CONNECTOR DERIVED COMPONENT, NAME=shear1,21.0,1.0*CONNECTOR ELASTICITY, COMPONENT=12.E+11,*CONNECTOR ELASTICITY, COMPONENT=22.E+11,*CONNECTOR ELASTICITY, COMPONENT=32.E+11,*CONNECTOR ELASTICITY, COMPONENT=42.E+11,*CONNECTOR ELASTICITY, COMPONENT=52.E+11,*CONNECTOR ELASTICITY, COMPONENT=62.E+11,*CONNECTOR DAMAGE INITIATION,CRITERION=FORCE,2.E+5*CONNECTOR POTENTIALnormal,.02,MACAULEYshear,.05,ABS*CONNECTOR DAMAGE EVOLUTION,TYPE=MOTION,AFFECTED COMPONENTS,SOFTENING=TABULAR1,2,3,4,5,60.0,0.01.0,1.E-7*CONNECTOR POTENTIAL1,2,3,4,5,6,*FASTENER, INTERACTION NAME=SPOTWELDS, ELSET=FASTCONN,PROPERTY=WELDPROPS BPIL, WPIL*FASTENER PROPERTY, NAME=WELDPROPS 0.001,*STEP*DYNAMIC,EXPLICIT ,10.E-3*CONTACT PAIR BPIL,*CONTACT PAIR,INTER=WELDS BPIL, WPIL*SURFACE INTERACTION,NAME=WELDS*FRICTION .3,*BOUNDARY,TYPE=VELOCITY 60000, 1,1, 0. 60000, 2,2, -25. 60000, 3,3, -2. 60000, 4,4, -7. 60000, 5,5, 0. 60000, 6,6, -78.5*OUTPUT,FIELD,TIME MARKS=YES,NUMBER INTER=50*NODE OUTPUT U, *OUTPUT,FIELD, NUMBER=2*NODE OUTPUT, NSET=QA_TEST U,*OUTPUT,HISTORY*ELEMENT OUTPUT, ELSET=FASTCONNCTF,CDIF,CDMG*ENERGY OUTPUT, VAR=PRESELECT*END STEP五，基准手册中1.2.6Buckling of a column with general contactBASEDHEADINGBUCKLING OF X SHAPED CROSS SECTIONWITH GENERAL CONTACT*PREPRINT,MODEL=YES,HISTORY=YES,CONTACT=YES*NODE 601,0.,-.2,0.0 611,0., .2,0.0 106,0., 0.,-.21106,0., 0., .2*10601,0.05,-.2,0.010611,0.05, .2,0.010106,0.05, 0.,-.211106,0.05, 0., .2*NGEN,NSET=A601,611,1106,1106,100*NGEN,NSET=A110601,10611,110106,11106,100*NODE200601,1.,-.2,0.0200611,1., .2,0.0200106,1., 0.,-.2201106,1., 0., .2*190601,0.95,-.2,0.0190611,0.95, .2,0.0190106,0.95, 0.,-.2191106,0.95, 0., .2*NGEN,NSET=B200601,200611,1200106,201106,100*NGEN,NSET=B1190601,190611,1190106,191106,100*NFILL,NSET=ALL1A1,B1,18,10000*ELEMENT,TYPE=S4R 1, 601,602,10602,10601101, 606,607,10607,10606201, 106,206,10206,10106301, 606,706,10706,10606*ELGEN,ELSET=FIN1 1, 5,1,1, 20,10000,5*ELGEN,ELSET=FIN2101, 5,1,1, 20,10000,5*ELGEN,ELSET=FIN3201, 5,100,1, 20,10000,5*ELGEN,ELSET=FIN4301, 5,100,1, 20,10000,5*ELSET,ELSET=SHELLSFIN1,FIN2,FIN3,FIN4*SHELL SECTION, SECTION INT=GAUSS,ELSET=SHELLS,MATERIAL=STEEL.003,*NODE500001,0.,-.45,-.45500002,0., .0,-.45500003,0., .45,-.45500004,0.,-.45,0.500005,0., .0,0.500006,0., .45,0.500007,0.,-.45,.45500008,0., .0,.45500009,0., .45, .45599995,0.,0.,0.*ELEMENT,TYPE=R3D4,ELSET=WALL1500001, 500001,500002,500005,500004500002, 500002,500003,500006,500005500003, 500004,500005,500008,500007500004, 500005,500006,500009,500008*NODE600001,1.,-.45,-.45600002,1., .0,-.45600003,1., .45,-.45600004,1.,-.45,0.600005,1., .0,0.600006,1., .45,0.600007,1.,-.45,.45600008,1., .0,.45600009,1., .45, .45699995,1.,0.,0.*ELEMENT,TYPE=R3D4,ELSET=WALL2600001, 600001,600002,600005,600004600002, 600002,600003,600006,600005600003, 600004,600005,600008,600007600004, 600005,600006,600009,600008*MATERIAL,NAME=STEEL*DENSITY7850., *ELASTIC200.E9,.3*PLASTIC250.E6,0.700.E6,1.*BOUNDARY599995,2,3599995,4,4599995,6,6699995,1,6*INITIAL CONDITIONS,TYPE=VELOCITY 599995, 1,1, 50. 599995, 5,5, 78.54*RESTART,TIMEMARKS=NO,WRITE,NUM=10*SURFACE,TYPE=ELEMENT, NAME=WALL1WALL1,SPOS*SURFACE,TYPE=ELEMENT, NAME=WALL2WALL2,SNEG*SURFACE,TYPE=NODE,NAME=A A,*SURFACE,TYPE=NODE,NAME=B B,*RIGID BODY, ELSET=WALL1, REF NODE=599995*RIGID BODY, ELSET=WALL2, REF NODE=699995*NSET,NSET=NRIGID 599995,699995*TIE, NAME=TIE1,POSITION TOLERANCE=1.E-4A, WALL1*TIE, NAME=TIE2,POSITION TOLERANCE=1.E-4B, WALL2*STEP*DYNAMIC,EXPLICIT,0.7E-2*BOUNDARY,TYPE=VEL* compress 35 percent axially and rotate 31 degrees during the step599995, 1,1, 50.599995, 5,5, 78.54*CONTACT*CONTACT INCLUSIONS,ALL ELEMENT BASED* Output requests for postprocessing*OUTPUT,FIELD,VAR=PRESELECT*OUTPUT,HISTORY,TIME INTERVAL=2.E-5*NODE OUTPUT, NSET=NRIGIDU, RF*ENERGY OUTPUTALLIE,ALLKE,ALLWK* Output requests for qa testing*ELSET,ELSET=QA_TESTELFIN4,*NSET,NSET=QA_TESTN,ELSET=FIN4*NSET,NSET=QA_TESTRBNRIGID,*OUTPUT,FIELD,NUMBER INTERVAL=1*ELEMENT OUTPUT,ELSET=QA_TESTELPEEQ,*NODE OUTPUT,NSET=QA_TESTNU, *NODE OUTPUT, NSET=QA_TESTRBU, RF*OUTPUT,HISTORY,TIME INTERVAL=1.75E-3*ENERGY OUTPUT,VAR=PRESELECT*FILE OUTPUT,NUMBER=4,TIMEMARKS=YES*ENERGY FILE*NODE FILEU,RF*END STEP六，基准手册中1.2.5Lateral buckling of an L-bracket*HEADING L-shaped bracket 17 element mesh, eigenvalue buckling analysis*NODE 101, 0., 225., 165, 225., 225., 173, 255., 225., 901, 0., 255., 965, 225., 255., 973, 255., 255., 1101, 225., 0., 1165, 225., 225., 1901, 255., 0., 1965, 255., 225.,*NGEN,NSET=EDGE1 101, 165 165, 173*NGEN,NSET=EDGE2 901, 965 965, 973*NGEN,NSET=EDGE3 1101, 1165*NGEN,NSET=EDGE4 1901, 1965*NFILL EDGE1, EDGE2, 8, 100*NFILL EDGE3, EDGE4, 8, 100*ELEMENT,TYPE=S4 101, 101, 109, 909, 901 1101,