动力弹塑性分析直接积分法

动力弹塑性分析的直接积分法MIDASIT开发二组1

动力弹塑性分析功能中的恢复力特性恢复力特性•Kinematichardening(标准型)

•Isotropichardening(标准型) •Mixedhardening(标准型)•Trilinearpeak-oriented(最大点指向型)•Trilinearorigin-oriented(原点指向型)•Degradingbilinear•Degradingtrilinear•Trilinearelastic•Ramberg-Osgood(Ramberg-Osgood型)•Takeda(武田型) •Takedawithslip•Trilinearslip•Degradingstiffnessmodel-Clough&Johnston•逆行型 •ひび割れ域剛性低下型•JR總硏剛性低減RC型 •JR總硏剛性低減SRC型•マルチリニアの逆行型 •3次係數逆行型

2

User-definedHysteresisModel•Degradingstiffness&strength

•Pinching&slip•Fibermodel•Multi-axialspringmodel•Tri-axialspringmodel

AxialSpringModels3ComparisonofVariousAxialSpringModelsFiberModelMulti-axialSpringModelTri-axialSpringModelSectionDivisionConcreteSteelTakedaModelModels4ComparisonofTri-axialSpringModelwithTestResultsSpecimenLoadHistory5ComparisonofTri-axialSpringModel6EX-T1.FreeVibrationof2DFrame:Linear&NonlinearTimeHistoryAnalysisbyModalAnalysis&DirectIntegration7LinearAnalysisModelwithProportionalDamping2kk2kkc=2a1kc=a1k2mm2mmppppProportionalDamping:C=a1KNonlinearAnalysisModelwithNL-LinkLLLL8•SectionDimension1stStoryColumn:b=0.1(m),h=0.1(m)2ndStoryColumn:b=0.1(m),h=0.1(m)•ColumnLength: L=4(m)•Material: E=2E+11(N/m2)•NodalMass: m=10000(kgf/g)•Damping: C=a1·K a1=5.60171411E-2•Quasi-staticLoad: p=6150(N)ModelProperties9Quasi-StaticLoadandInitialDisplacementppx2(0)=2cmx1(0)=1cm10DisplacementResultQuasi-staticLoadingFreeVibration11Time(sec)TheoreticalLinearModalAnalysisNonlinearModalAnalysisLinearDirectIntegrationNonlinearDirectIntegration1.8-1.705E-02-1.708E-02(1.002)-1.708E-02(1.002)-1.708E-02(1.002)-1.708E-02(1.002)3.61.447E-021.450E-02(1.002)1.450E-02(1.002)1.451E-02(1.003)1.450E-02(1.002)5.4-1.222E-02-1.226E-02(1.003)-1.226E-02(1.003)-1.227E-02(1.003)-1.226E-02(1.003)7.21.028E-021.032E-02(1.004)1.031E-02(1.003)1.033E-02(1.005)1.032E-02(1.004)Theoretical&NumericalSolutionof2ndFloorDisplacement(m)12EX-T2.SimpleBeamsupportedbyViscousDamper:NonlinearTimeHistoryAnalysisbyModalAnalysis&DirectIntegration13cpTheoreticalAnalysisModelLLmmcE,IE,IBeam&ViscousDamperEquivalentSDFModel14NumericalAnalysisModelwith10BeamElements&NL-linkpStepLoadingFunction15•SectionDimension: b=0.2(m),h=0.2(m)•EntireBeamLength: L=10(m)•Material: E=2E+11(N/m2)•NodalMass: m=3000(kgf/g)•ViscousDampingCoefficientinDamper: cd=3.88297266E+4•StepLoad: p(t)=0(N)(t<0) 10000(N)(t≥0)ModelProperties16ComparisonofModalAnalysis&DirectIntegrationCenterDisplacementCenterAcceleration17Theoretical&NumericalSolutionofCenterDisplacement(m)Time(sec)TheoreticalNonlinearModalAnalysisNonlinearDirectIntegration(1)(2)(2)/(1)(3)(3)/(1)0.5-1.345E-02-1.347E-021.001-1.348E-021.0021.0-3.818E-03-3.817E-031.000-3.802E-030.9961.5-1.058E-02-1.060E-021.002-1.062E-021.0042.0-5.932E-03-5.929E-030.999-5.904E-030.9952.5-9.058E-03-9.076E-031.002-9.102E-031.0053.0-7.011E-03-7.009E-031.000-6.984E-030.99618EX-T3.Base-Isolated3DFrameStructure:NonlinearTimeHistoryAnalysisbyModalAnalysis&DirectIntegration19NL-linkPropertiesDampingProperties20GroundAccelerationX-DirectionY-Direction21RoofDispl.(X-Direction)RoofDispl.(Y-Direction)TimeHistoryComparison22TimeHistoryComparisonNL-link1HisteresisCurve(Shear-y)NL-link1HisteresisCurve(Shear-y)23MaximumDisp.ComparisonMaximumFloorDisplacement(X-direction)MaximumFloorDisplacement(Y-direction)FloorModalAnalysis(ft)DirectIntegration(ft)10.22010.220020.22700.226030.23170.2296FloorModalAnalysis(ft)DirectIntegration(ft)10.35070.352320.36950.370330.37850.379524ComputationTimeComparison(sec)(CPU:IntelPentium4)ModalAnalysisDirectIntegration8.00104.1525EX-T4.InelasticBeamElement:NonlinearTimeHistoryAnalysisbyDirectIntegration26HingeDefinitionofBeamElementFlexibilityMBendingMomentCurvatureEI1FlexibilityDistributionTri-linearSkeletonCurve27YieldSurfaceofBeamSectionMzMyPTwoSurfaceModel•CrackSurface•YieldSurface28InelasticAnalysisModelMzMyPMzMy•BeamProperty:

L=4mb=50cmh=70cmE=230000kgf/cm2SectionHingeProperty

•5hingesareused

•Yieldforce&momentPy=91840kgfMy=30612kgf-m