ALE/CEL中重力場波動是需要解決的關(guān)鍵問題。

以CPT貫入問題為例，實現(xiàn)ALE和CEL法中土體重力的施加。（附件中附CAE，inp文件）

1. ALE法，模型圖如圖1，CPT直徑0.036m，土體長寬均為1.5m，黏土容重8kN/m³。
2. 重力平衡后的豎向應(yīng)力S22云圖如圖2所示，土體沉降如圖3，可以看出S22成層分布，土體沉降較小。
3. 
   

     圖1 ALE法CPT貫入模型圖             圖2 ALE法重力平衡后S22云圖 

CPT貫入后的豎向應(yīng)力S22云圖如圖4所示，下部土體受CPT貫入影響很小，應(yīng)該保持初始應(yīng)力狀態(tài)，為檢驗重力平衡效果，取遠(yuǎn)離CPT的土體底部不同深度的單元，畫出S22在貫入過程中的變化，如圖5，可以發(fā)現(xiàn)，S22曲線基本保持水平，重力場穩(wěn)定，說明重力平衡效果良好。

      圖3 重力平衡后土體沉降           圖4 ALE法CPT貫入后S22云圖 

         圖5 CPT貫入過程中遠(yuǎn)離CPT的土單元的S22（驗證重力平衡效果）

2.CEL法（step1重力，step2保持CPT靜止來驗證重力分布是否正確），模型圖如圖6，重力平衡后的豎向應(yīng)力S22云圖如圖7所示。重力施加完后，保持CPT靜止10S，最終的S22云圖如8，土單元S22的變化如圖9.

S22云圖成層分布，在整個過程（STEP2/Dynamic explicit分析步）中S22曲線基本保持水平，重力場穩(wěn)定，說明重力平衡效果良好。

    圖6 CEL法CPT模型圖             圖7 CEL法重力平衡后S22云圖

      圖8 CEL法重力平衡后保持CPT靜止10S后的S22云圖

        圖9 CEL法CPT保持靜止過程中土單元的S22（驗證重力平衡效果）