第九章强迫运动 9.1概述 (1)用于分析具有地基输入加速度、速度和位移的约束结构 (2)常见例子是地震(瞬态分析)、正弦扫描器测试仿真(频响分析) (3) Nastran中将作用在未约束结构上的力转化为等效约束结构的强迫运动 (4)使用的方法有:大质量法、大刚度法、 Lagrange乘子法 9.2瞬态分析中的强迫运动 (1)如果卡片 TLOAD中的第5字段选择,给定的加速度、速度、位移不会执行 (2) Nastran中仅能作用力{P(t)}在结构上,如果选择强迫运动,即假设用户在大质 量上施加了一个运动,因此作用在大质量上的力是与加速度成正比的 9.3瞬态分析中的大质量法 Large Mass Structure of M=(103-105m nterest Rigid Body (1)对强迫加速度ib P(M large mass 2)对作用在大质量上的力P可以得到期望的ib,质量M在卡片 CMASSI或 cONM2中输入;比例因子ML(对力)在 DAREA或 DLOAD中输入 (3)对瞬态响应中期望i,lb,用卡片 TLOAD1或 TLOAD2的第5字段进行区分 (4)大质量应该为结构质量的103~108倍,106倍效果最好 (5)根据卡片 TLOADI输入的强迫加速度、速度、位移得到加速度历程 a)定义
第九章 强迫运动 9.1 概述 （1） 用于分析具有地基输入加速度、速度和位移的约束结构 （2） 常见例子是地震（瞬态分析）、正弦扫描器测试仿真（频响分析） （3） Nastran中将作用在未约束结构上的力转化为等效约束结构的强迫运动 （4） 使用的方法有：大质量法、大刚度法、Lagrange乘子法 9.2 瞬态分析中的强迫运动 （1） 如果卡片TLOADi中的第5字段选择，给定的加速度、速度、位移不会执行 （2） Nastran中仅能作用力{P(t)}在结构上，如果选择强迫运动，即假设用户在大质 量上施加了一个运动，因此作用在大质量上的力是与加速度成正比的 9.3 瞬态分析中的大质量法 （1）对强迫加速度 b u （2） 对作用在大质量上的力 P 可以得到期望的 b u ,质量 ML 在卡片 CMASSi 或 CONM2 中输入；比例因子 ML (对力)在 DAREA 或 DLOAD 中输入 （3） 对瞬态响应中期望 ub ub , ，用卡片 TLOAD1 或 TLOAD2 的第 5 字段进行区分 （4） 大质量应该为结构质量的 103~108倍，106倍效果最好 （5） 根据卡片 TLOADi 输入的强迫加速度、速度、位移得到加速度历程 a) 定义
FN=TLOADi at time N(Atn) FN= computed force at time N(Atn) b)强迫加速度 F N c)强迫速度 N 0 d)强迫位移 NN N_1-「N-2 At ftN子0 △t 0 if tN=0 e)使用大质量法,,同样要乘比例因子 F= Ma Mass scaling for ↑ Computed fr Calculation of TLOADi Input Applied Force 94瞬态分析中的大刚度法 Structure of Interest Large spring
b) 强迫加速度 c) 强迫速度 d) 强迫位移 e) 使用大质量法， b u 同样要乘比例因子 9.4 瞬态分析中的大刚度法
(KL large spring) (1)对强迫位移u (2)作用力P于弹簧或结构表面得到期望的ab,弹簧K在卡片 CELAS1或 CELAS2中输入;比例因子KL在 DAREA或 DLOAD中输入 (3)大质量应该为结构质量的102~104倍O2,其中2为激励的截断频率 9.5瞬态分析中的 LAGRANGE乘子法 (1)通过DMG或模拟约束的瞬态函数增加动态矩阵的行和列 (2)每行表示如下形式的约束 u =Y= enforced displacement) (3)每列表示称为 Lagrange乘子的附加自由度λ,其值为第j个约束自由度的约束 (4)在直接瞬态分析中, EPOINTS, SPOINTS或GRID用于描述 Lagrange乘子,模 态响应分析中仅用 EPOINTS (5) Nastran输入中需要设置卡片 SDCOMP=32 9.6强迫位移例子 Rod Constra ed Unconstrained structure Structure for msc/ Nastran Enforced Motion Simulation Enforced Displacement Ramp
（1） 对强迫位移 b u ， （2） 作用力 P 于弹簧或结构表面得到期望的 b u ,弹簧 KL 在卡片 CELAS1 或 CELAS2 中输入；比例因子 KL 在 DAREA 或 DLOAD 中输入 （3） 大质量应该为结构质量的 102~104倍 2 c ，其中 2 c 为激励的截断频率 9.5 瞬态分析中的LAGRANGE乘子法 （1） 通过 DMIG 或模拟约束的瞬态函数增加动态矩阵的行和列 （2） 每行表示如下形式的约束 （3） 每列表示称为 Lagrange 乘子的附加自由度 j ，其值为第 j 个约束自由度的约束 力 （4） 在直接瞬态分析中，EPOINTs, SPOINTs 或 GRID 用于描述 Lagrange 乘子，模 态响应分析中仅用 EPOINTs （5） Nastran 输入中需要设置卡片 SDCOMP = 32 9.6 强迫位移例子
For this example m=1.0 k=(2)2 To get k, use rod with A=39.4784, E=1.0, L=1.0 Solution*: u(t)=d[t-1/o sin ot] u(t=do sin ot] where(=6.28318 d=2.5 (1)大质量法 1.0 m k=394784 M=1x10 在 TLOAD1的第5字段中,指定uB为强迫运动 (2)大刚度法 m=1.0 k=394784 kB=4.0E8 B 在 TLOAD1中指定PB为力,乘以强迫运动值kB,在 DAREA中输入 Large Mass Method
（1）大质量法 在TLOAD1的第5字段中，指定uB为强迫运动 （2）大刚度法 在 TLOAD1 中指定 PB为力，乘以强迫运动值 kB,在 DAREA 中输入 Large Mass Method
SOL 109 TIME 10 TITLE=S DOF, GROUND MOTION SUBTIE LARGE MAS 工ABEL=ⅩV工AⅩ SET 1=1 DLOAD=105 TSTEP=6 DISPLSALL ACCELSALL OLOAD=1 BEGIN BULK DAREA105,1,1,1.0E8 CONM2,3,2,1.0 CONM2,999,1,,1.0E8 cROD,1,1,1,2 GRID,1,,0.,0.,0.,,23456 GRID,2,,1.,0.,0.,,23456 MAT1,1,1.,,0 PARAM, USETPRT, 1 PROD,1,1,39.4784 TABLEDI,8,,/,+TABl +TAB1,0.,0,.025,0.,05,0,.075,0.,+TAB2 +TAB2,5.075,12.5,10.075,25.0,ENDT TLOAD1,105,105,,1,8 TSTEP,6,199,.025,1 ENDDAT Large stiffness Method SOL 109 TIME 10 CEND TITLE=S DOF, GROUND MOTION SUBTIE LARGE STIFENESS 工ABEL= X VIA FORCING SET 1=1 100 DLOAD=105 TSTEP=6 DISPLFALL ACCELSALL OLOAD=1
SOL 109 TIME 10 CEND TITLE=SDOF, GROUND MOTION SUBTI= LARGE MASS LABEL = X VIA X SET 1=1 DLOAD=105 TSTEP=6 DISPL=ALL ACCEL=ALL OLOAD=1 BEGIN BULK DAREA,105,1,1,1.0E8 CONM2,3,2,,1.0 CONM2,999,1,,1.0E8 CROD,1,1,1,2 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 +TAB1,0.,0.,.025,0.,.05,0.,.075,0.,+TAB2 +TAB2,5.075,12.5,10.075,25.0,ENDT TLOAD1,105,105,,1,8 TSTEP,6,199,.025,1 ENDDATA Large Stiffness Method SOL 109 TIME 10 CEND TITLE=SDOF, GROUND MOTION SUBTI= LARGE STIFFNESS LABEL = X VIA FORCING SET 1=1 SPC = 100 DLOAD=105 TSTEP=6 DISPL=ALL ACCEL=ALL OLOAD=1