BEGIN BULK DAREA,105,1,1,4.0E8 CONM2,3,2,,1.0 cROD,1,1,1,2 0.,0.,,2345 GRID,2,,1.,0.,0.,,23456 GRID,99,,-1.,0.,0 sPC,100,99,123456 CELAS2,5,4.0E8,99,1,1,1 MAT1,1,1.,,0 PARAM, USETPRT, 1 D,1,1,39.4784 TABLEDI, 8 /////,+TAB1 +TAB1,0.,0,,.025,0.,.05,0.,.075,0.,+TAB2 +TAB2,5.075,12.5,10.075;25.0,ENDT TLOAD1,105,105,,0,8 TSTEP,6,199,.025,1 ENDDATA (3) Lagrange乘子法 m=1.0 2 k=39.4784 B 1001 运动方程(静力部分)为
BEGIN BULK DAREA,105,1,1,4.0E8 CONM2,3,2,,1.0 CROD,1,1,1,2 GRID,1,,0.,0.,0.,,23456 GRID,2,,1.,0.,0.,,23456 GRID,99,,-1.,0.,0. SPC,100,99,123456 CELAS2,5,4.0E8,99,1,1,1 MAT1,1,1.,,0. PARAM,USETPRT,1 PROD,1,1,39.4784 TABLED1,8,,,,,,,,+TAB1 +TAB1,0.,0.,.025,0.,.05,0.,.075,0.,+TAB2 +TAB2,5.075,12.5,10.075,25.0,ENDT TLOAD1,105,105,,0,8 TSTEP,6,199,.025,1 ENDDATA （3）Lagrange 乘子法 运动方程（静力部分）为
kU= P 394784-39.47841.0 39478439.47840 2 0 1001 1001 1001 394784u1-394784u2+11001=0 -394784u1+394784u2 0 +x入1001=y1001 其中,x为一个小整数,用K2PP输入 Lagrange乘子 Lagrange Multiplier Method TIME 10 T工TTE=SDOE, GROUND MOT工oN SUBTI= LAGRANGE MULT DISPL K2PP=ST工F SET 1=1 DLOAD=105 TSTEP=6 DS卫I=ATT ACCEL=ALL BEG IN BULK DAREA,105,1001,0,1.0 coNM2,3,2,1.0 cRoD,1,1,1,2 DMIG, STIE,0,6.1. 0 DMIG,ST工F,1001,,,101,,-1.0E-10 DMIG,STIF,1001,,,1,1,,1 EPOINT 1001 GRTD,1,,0.,0.,0.,,23456 1.,0.,0.,23456 MA1,1,1.,,0 PARAM, US PROD, 1 1 4784 TABLEDl, 8,,,,,,,,+TAB1
其中，x 为一个小整数，用 K2PP 输入 Lagrange 乘子 Lagrange Multiplier Method SOL 109 TIME 10 CEND TITLE=SDOF, GROUND MOTION SUBTI= LAGRANGE MULT - DISPL K2PP=STIF SET 1=1 DLOAD=105 TSTEP=6 DISPL=ALL ACCEL=ALL OLOAD=ALL BEGIN BULK DAREA,105,1001,0,1.0 CONM2,3,2,,1.0 CROD,1,1,1,2 DMIG,STIF,0,6,1,0 DMIG,STIF,1001,,,101,,-1.0E-10 DMIG,STIF,1001,,,1,1, .1 EPOINT,1001 GRID,1,,0.,0.,0.,,23456 GRID,2,,1.,0.,0.,,23456 MAT1,1,1.,,0. PARAM,USETPRT,1 PROD,1,1,39.4784 TABLED1,8,,,,,,,,+TAB1