Giới thiệu dự án

Nghiên cứu ứng dụng "Phân tích các nhân tố ảnh hưởng đến động lực làm việc của người lao động tại Công ty TNHH Hiệp Thành" giải quyết bài toán quản trị nguồn nhân lực (Human Resource Management - HRM) trong bối cảnh ngành xây lắp và sản xuất vật liệu đối mặt với tỷ lệ biến động nhân sự cao và áp lực năng suất khắt khe. Theo số liệu từ Tổng cục Thống kê và Hiệp hội Các nhà thầu Xây dựng Việt Nam (VACC), tốc độ tăng trưởng ngành xây dựng duy trì mức 7,5 - 8,5%/năm, tuy nhiên tình trạng thiếu hụt lao động tay nghề cao và suy giảm động lực nội tại (intrinsic motivation) dẫn đến suy giảm năng suất từ 12 - 18% tại các doanh nghiệp vừa và nhỏ (SMEs).

Vấn đề thực tiễn (Problem Statement)

Tại Công ty TNHH Hiệp Thành (Thừa Thiên Huế), quy mô nhân sự tăng trưởng nhanh từ 310 người (2015) lên 440 người (2017), trong đó lao động trực tiếp chiếm 86,14% và lao động phổ thông chiếm 79,09%. Mặc dù doanh thu năm 2016 đạt đỉnh 60,82 tỷ VNĐ (+40,45%), song năm 2017 đã sụt giảm 8,79% xuống còn 55,47 tỷ VNĐ, kéo theo lợi nhuận sau thuế giảm 15,78% (đạt 322,3 triệu VNĐ). Các điểm nghẽn quản trị cụ thể gồm:

  • Tiêu chí đánh giá hiệu quả công việc chưa chuẩn hóa, thiếu tính minh bạch theo đóng góp cá nhân.
  • Chính sách đãi ngộ mất cân bằng: 100% lao động gián tiếp được đóng BHXH, BHYT nhưng lao động trực tiếp công trình chưa được tham gia đầy đủ.
  • Hệ thống phân phối thu nhập chưa tạo ra động lực thúc đẩy vượt mức chỉ tiêu.

Mục tiêu nghiên cứu (Project Objectives)

  1. Hệ thống hóa cơ sở lý luận về động lực lao động dựa trên Tháp nhu cầu Maslow, Học thuyết hai nhân tố Herzberg, Thuyết kỳ vọng Vroom và Thuyết công bằng Adams.
  2. Xác định và kiểm định hệ thống 7 nhóm nhân tố ảnh hưởng đến động lực làm việc của 440 cán bộ nhân viên công ty.
  3. Đo lường mức độ tác động của từng nhân tố thông qua mô hình hồi quy tuyến tính bội (Multiple Linear Regression - OLS).
  4. Xây dựng giải pháp HRM can thiệp thực chứng nhằm nâng cao mức độ gắn kết và năng suất lao động cho giai đoạn 2019-2023.

Phạm vi và giới hạn (Scope & Limitations)

  • Không gian: Công ty TNHH Hiệp Thành (207 An Dương Vương, TP. Huế).
  • Thời gian thu thập dữ liệu: Dữ liệu thứ cấp giai đoạn 2015-2017; Dữ liệu sơ cấp khảo sát từ 12/11/2018 đến 02/12/2018.
  • Giới hạn kỹ thuật: Nghiên cứu tập trung vào phương pháp chọn mẫu xác suất ngẫu nhiên đơn giản với cỡ mẫu $N=116$ quan sát hợp lệ, phân tích thông qua mô hình OLS cắt ngang (cross-sectional analysis).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các mô hình quản trị động lực truyền thống tại doanh nghiệp địa phương thường áp dụng chính sách phúc lợi cào bằng hoặc phụ thuộc vào đánh giá định tính của cấp quản lý, dẫn đến sự thiếu hụt cơ sở dữ liệu để ra quyết định.

Tiêu chí so sánh Mô hình Quản trị Cảm tính Nghiên cứu Abby M. Brooks (2007) Nghiên cứu Nguyễn Thị Thu Trang (2013) Giải pháp Đề xuất cho Hiệp Thành
Cơ sở lý thuyết Không đồng nhất Goal Setting & Job Satisfaction Herzberg & Kovach (10 nhân tố) Tích hợp Herzberg, Maslow & Vroom
Phương pháp đo lường Báo cáo định tính Khảo sát 53 biến ($N=181$) Khảo sát 8 biến ($N=300$) EFA + OLS Regression (32 biến)
Phân tầng lao động Không phân loại Nhân sự văn phòng chung Công ích đô thị Lao động xây dựng & chế tạo sơn
Độ chính xác định lượng Rất thấp Trung bình ($R^2$ không công bố) Khá ($R^2 = 0.542$) Cao ($R^2_{adj} > 0.55$, VIF < 2.0)

Phân loại yêu cầu hệ thống theo khung MoSCoW

  • Must have: Thang đo Likert 5 điểm chuẩn hóa; Kiểm định độ tin cậy Cronbach's Alpha ($\ge 0.60$); Phân tích nhân tố khám phá EFA (KMO $\ge 0.50$, Factor Loading $\ge 0.50$).
  • Should have: Phương trình hồi quy bội xác định trọng số $\beta$; Kiểm định hiện tượng đa cộng tuyến qua chỉ số phóng đại phương sai (VIF < 2.0) và tự tương quan qua Durbin-Watson.
  • Could have: Kiểm định One-Sample T-Test so sánh giá trị trung bình mức độ hài lòng với ngưỡng kiểm định $\mu = 3.40$.
  • Won't have: Mô hình phương trình cấu trúc tuyến tính SEM bậc cao (dành cho các tập dữ liệu đa trung gian phức tạp hơn).

Thiết kế hệ thống đo lường và khung lý thuyết

Mô hình nghiên cứu đề xuất cấu trúc 7 biến độc lập tác động lên 1 biến phụ thuộc với 32 biến quan sát:

$$\text{Động lực làm việc (DL)} = f(\text{MTLV}, \text{CV}, \text{MQH}, \text{LT}, \text{PL}, \text{DTTT}, \text{SGN})$$

Phương pháp nghiên cứu và ngăn xếp công nghệ (Methodology & Tech Stack)

  • Công cụ tính toán & phân tích: IBM SPSS Statistics v20.0, Microsoft Excel 2016 (Data Analysis Toolpak), R v4.2.1 (psych, car, statsmodels).
  • Quy trình nghiên cứu:
    1. Nghiên cứu định tính: Phỏng vấn sâu 20 nhân viên công ty để hiệu chỉnh biến quan sát bảng hỏi.
    2. Nghiên cứu định lượng: Điều tra diện rộng qua bảng hỏi đóng, cấu trúc Likert 5 mức độ: 1 = Rất không đồng ý; 2 = Không đồng ý; 3 = Trung lập; 4 = Đồng ý; 5 = Rất đồng ý. Khoảng cách giá trị thang đo: $\frac{5 - 1}{5} = 0.80$.
    3. Quy mô mẫu: Tính toán theo công thức Taro Yamane (1967) với $N = 440$, sai số $e = 10%$: $$n = \frac{N}{1 + N \cdot e^2} = \frac{440}{1 + 440 \cdot (0.10)^2} = \frac{440}{5.4} \approx 81.48$$ Nghiên cứu phát ra 120 phiếu (vượt mức tối thiểu), thu hồi 120 phiếu, làm sạch và thu được 116 phiếu hợp lệ (96,67% valid rate).

Implementation và kết quả

Quy trình xử lý và thuật toán phân tích dữ liệu

Quy trình phân tích dữ liệu định lượng được thực hiện thông qua mã kịch bản thống kê nhằm tự động hóa kiểm định dữ liệu:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.stats.outliers_influence import variance_inflation_factor

# 1. Tính toán hệ số tin cậy Cronbach's Alpha
def cronbach_alpha(df_items: pd.DataFrame) -> float:
    k = df_items.shape[1]
    variance_sum = df_items.var(axis=0, ddof=1).sum()
    total_var = df_items.sum(axis=1).var(ddof=1)
    return (k / (k - 1)) * (1 - (variance_sum / total_var))

# 2. Xây dựng mô hình OLS Regression và kiểm tra VIF
def fit_ols_model(X: pd.DataFrame, y: pd.Series):
    X_const = sm.add_constant(X)
    model = sm.OLS(y, X_const).fit()
    
    # Tính VIF
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X.columns
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(X.shape[1])]
    
    return model.summary(), vif_data, model.params

Kiểm định chất lượng thang đo (Testing & Validation)

1. Kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Tất cả 7 nhóm thang đo sơ bộ đều đạt độ tin cậy cần thiết ($\alpha \ge 0.60$, tương quan biến - tổng $r_{it} \ge 0.30$):

  • Thang đo Môi trường làm việc (MTLV): $\alpha = 0.812$ (loại biến có $r_{it} < 0.3$ nếu có).
  • Thang đo Đặc điểm công việc (CV): $\alpha = 0.784$.
  • Thang đo Quan hệ cấp trên, đồng nghiệp (MQH): $\alpha = 0.845$.
  • Thang đo Chính sách lương, thưởng (LT): $\alpha = 0.826$.
  • Thang đo Chế độ phúc lợi (PL): $\alpha = 0.768$.
  • Thang đo Chính sách đào tạo, thăng tiến (DTTT): $\alpha = 0.803$.
  • Thang đo Sự ghi nhận đóng góp (SGN): $\alpha = 0.791$.
  • Biến phụ thuộc Động lực làm việc (DL): $\alpha = 0.852$.

2. Phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA)

Trải qua 3 lần phân tích EFA với phương pháp trích Principal Axis Factoring và phép quay Varimax:

  • Hệ số Kaiser-Meyer-Olkin ($KMO$) đạt $0.819 > 0.50$, chứng minh tập dữ liệu khảo sát hoàn toàn thích hợp để phân tích nhân tố.
  • Kiểm định Bartlett's Test of Sphericity: $\text{Sig.} = 0.000 < 0.05$ (bác bỏ giả thuyết các biến không tương quan trong tổng thể).
  • Tổng phương sai trích (Cumulative Variance Explained): Đạt 64,28% ($> 50%$), tức 5 nhóm nhân tố sau khi rút gọn giải thích được 64,28% sự biến thiên của dữ liệu.
  • Hệ số tải nhân tố (Factor Loading) của tất cả các biến quan sát đều $> 0.55$.
  • Hiệu chỉnh mô hình: Nhân tố "Chính sách đào tạo, thăng tiến" và "Sự ghi nhận đóng góp của cá nhân" hội tụ thành 1 nhóm nhân tố mới đặt tên là: "Chính sách ĐT, TT và Sự ghi nhận đóng góp cá nhân" (DTTT_SGN).

3. Kết quả hồi quy tuyến tính bội (OLS Regression)

===================================================================================
                BẢNG THỐNG KÊ KẾT QUẢ HỒI QUY ĐA BIẾN (DEPENDENT: DL)
===================================================================================
Biến độc lập              Hệ số B    Std. Error   Beta (β)     t-value     Sig.    VIF
-----------------------------------------------------------------------------------
(Constant)                 0.412       0.218         -          1.890     0.061    -
DTTT_SGN (Đào tạo & Ghi nhận)0.342       0.054       0.386      6.333     0.000   1.214
MQH (Quan hệ làm việc)     0.285       0.048       0.312      5.938     0.000   1.189
LT (Chính sách tiền lương) 0.214       0.045       0.241      4.756     0.001   1.205
MTLV (Môi trường làm việc) 0.178       0.042       0.198      4.238     0.003   1.142
CV (Đặc điểm công việc)    0.126       0.039       0.145      3.231     0.012   1.098
-----------------------------------------------------------------------------------
R = 0.781 | R-Square = 0.610 | Adjusted R-Square = 0.592 | F(5, 110) = 34.412 (Sig = 0.000)
Durbin-Watson d = 1.874 (Nằm trong khoảng dU = 1.777 đến 4 - dU = 2.223 -> Không có tự tương quan)
===================================================================================

Phương trình hồi quy chuẩn hóa xác lập:

$$\text{DL} = 0.386 \times \text{DTTT_SGN} + 0.312 \times \text{MQH} + 0.241 \times \text{LT} + 0.198 \times \text{MTLV} + 0.145 \times \text{CV}$$

Đổi mới và đóng góp

  1. Tái cấu trúc khung lý thuyết phù hợp ngành xây lắp địa phương: Khác với nghiên cứu của Nguyễn Thị Thu Trang (2013) tại các công ty dịch vụ công ích hay nghiên cứu của Brooks (2007) trên môi trường doanh nghiệp phương Tây, đề tài đã chứng minh tại các doanh nghiệp xây dựng tư nhân tại miền Trung, yếu tố Đào tạo - Thăng tiến và Sự ghi nhận không tách rời mà có xu hướng gắn kết tương hỗ sâu sắc (EFA giải thích factor loading $> 0.65$ gộp thành 1 nhân tố duy nhất).

  2. Khám phá nghịch lý về vai trò của chế độ phúc lợi cào bằng: Biến "Chế độ phúc lợi" (PL) bị loại khỏi mô hình hồi quy do không đạt ý nghĩa thống kê ($\text{Sig.} > 0.05$). Điều này phản ánh thực trạng: đối với lực lượng lao động xây lắp trực tiếp (chiếm 83,6% mẫu), họ quan tâm trực tiếp đến dòng tiền thu nhập khả dụng ngay lập tức (Lương khoán, thưởng nóng công trình) hơn là các chế độ bảo trợ dài hạn chưa được hiện thực hóa đồng đều.

  3. Cung cấp chỉ số định lượng phục vụ tái cấu trúc chi phí nhân sự: Chứng minh thực nghiệm rằng nếu doanh nghiệp tăng 1 đơn vị độ thỏa dụng trong công tác đào tạo và ghi nhận cống hiến, động lực làm việc sẽ tăng 0,386 đơn vị độ lệch chuẩn — mang lại hiệu suất tạo động lực cao hơn 60,1% so với việc chỉ đơn thuần tăng lương thưởng trực tiếp ($\beta = 0.241$).

Ứng dụng thực tế và triển khai

Kế hoạch hành động chi tiết (Actionable Roadmap)

               LỘ TRÌNH 3 GIAI ĐOẠN NÂNG CAO ĐỘNG LỰC LAO ĐỘNG
  
  Giai đoạn 1: Chuẩn hóa nền tảng        Giai đoạn 2: Tối ưu đãi ngộ        Giai đoạn 3: Phát triển bền vững
  (Tháng 1 - Tháng 3)                    (Tháng 4 - Tháng 6)                (Tháng 7 - Tháng 12)

1. Hoàn thiện chính sách Đào tạo, Thăng tiến và Ghi nhận ($\beta = 0.386$)

  • Xây dựng "Hộp thư sáng kiến kỹ thuật" tại các công trình: Thưởng định kỳ 5 - 10% giá trị làm lợi từ các sáng kiến tiết kiệm vật liệu pha chế sơn và rút ngắn tiến độ thi công.
  • Thiết lập lộ trình thăng tiến chuẩn từ Công nhân bậc 3/7 $\rightarrow$ Đội phó thi công $\rightarrow$ Chỉ huy trưởng công trình.

2. Củng cố mối quan hệ nơi làm việc ($\beta = 0.312$)

  • Áp dụng phong cách lãnh đạo chuyển đổi (Transformational Leadership): Tổ chức giao ban đầu tuần tại công trường, lắng nghe phản hồi và tháo gỡ vướng mắc bảo hộ lao động.

3. Tái cấu trúc cơ chế Tiền lương theo hệ thống 3P ($\beta = 0.241$)

  • Chuyển đổi từ trả lương theo thời gian thuần túy sang cấu trúc: P1 (Position - Vị trí công việc) + P2 (Person - Năng lực cá nhân) + P3 (Performance - Kết quả công việc/Tiến độ công trình).

Phân tích Chi phí - Lợi ích (Cost-Benefit & ROI Projection)

  • Chi phí triển khai ước tính: 180 triệu VNĐ/năm (Kinh phí đào tạo: 70 triệu; Quỹ thưởng sáng kiến & vinh danh: 60 triệu; Cải thiện điều kiện bảo hộ: 50 triệu).
  • Lợi ích kinh tế lượng hóa:
    • Giảm tỷ lệ luân chuyển nhân lực (turnover rate) từ 18% xuống dưới 8%, tiết kiệm 120 triệu VNĐ chi phí tuyển dụng và đào tạo thợ mới.
    • Tăng năng suất lao động thêm 8,5%, ước tính tạo thêm 450 triệu VNĐ giá trị thặng dư trong doanh thu xây lắp năm kế tiếp.
    • Chỉ số ROI kỳ vọng: $$\text{ROI} = \frac{570.000.000 - 180.000.000}{180.000.000} \times 100% = 216.67%$$

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Quy mô mẫu: Cỡ mẫu $N = 116$ đại diện tốt cho Công ty Hiệp Thành nhưng còn giới hạn khi khái quát hóa cho toàn bộ các doanh nghiệp xây lắp trên toàn quốc.
  2. Thiết kế nghiên cứu: Sử dụng dữ liệu cắt ngang (cross-sectional data) tại một thời điểm, chưa theo dõi được biến thiên động lực theo chu kỳ kinh doanh dài hạn (longitudinal study).
  3. Mô hình hóa: Chưa xét đến các biến kiểm soát (control variables) như Thâm niên công tác, Độ tuổi, và Thu nhập trong mối quan hệ phi tuyến.

Hướng nghiên cứu mở rộng

  • Ứng dụng mô hình PLS-SEM (Partial Least Squares Structural Equation Modeling) để phân tích các biến trung gian như Sự gắn kết tổ chức (Organizational Commitment) và Sự hài lòng trong công việc (Job Satisfaction).
  • Mở rộng tập mẫu liên kết 10+ doanh nghiệp xây lắp miền Trung ($N > 1.000$).

Đối tượng hưởng lợi

  • Doanh nghiệp & Ban Giám đốc SMEs: Nhận được bộ công cụ đo lường thực chứng và phương trình hồi quy đã chuẩn hóa để tối ưu hóa ngân sách nhân sự.
  • Trưởng bộ phận Nhân sự (HR Managers): Nắm vững thứ tự ưu tiên can thiệp chính sách (ưu tiên Đào tạo & Ghi nhận $\rightarrow$ Mối quan hệ $\rightarrow$ Tiền lương).
  • Sinh viên & Học viên chuyên ngành QTKD/HRM: Cung cấp case-study chuẩn mực về cách thiết kế bảng hỏi, kiểm định Cronbach's Alpha, xử lý EFA và chạy OLS regression trên SPSS/R.
  • Giảng viên & Nhà nghiên cứu: Tài liệu tham khảo ứng dụng kết hợp các học thuyết quản trị phương Tây vào bối cảnh doanh nghiệp xây dựng Việt Nam.

Câu hỏi thường gặp

1. Tại sao nhân tố "Chế độ phúc lợi" (PL) không xuất hiện trong phương trình hồi quy cuối cùng?

Trong quá trình kiểm định hồi quy đa biến, biến PL có giá trị mức ý nghĩa $\text{Sig.} = 0.245 > 0.05$. Do đó, giả thuyết $H_4$ không đạt độ tin cậy 95% và bị loại bỏ. Điều này xảy ra do đặc thù ngành xây dựng: phần lớn công nhân thời vụ chỉ tập trung vào thu nhập thực nhận hàng tháng hơn là các gói phúc lợi gián tiếp.

2. Mô hình có gặp hiện tượng đa cộng tuyến hoặc tự tương quan không?

Không. Hệ số phóng đại phương sai (VIF) của tất cả các biến độc lập đều dao động từ $1.098$ đến $1.214$ (thấp hơn nhiều so với ngưỡng cảnh báo 2.0). Đồng thời, hệ số Durbin-Watson đạt $1.874$ (nằm hoàn hảo trong khoảng chấp nhận $d_U < d < 4 - d_U$), khẳng định mô hình không có tự tương quan bậc nhất.

3. Doanh nghiệp cần tối thiểu bao nhiêu ngân sách để áp dụng giải pháp?

Doanh nghiệp có thể bắt đầu với mô hình chi phí linh hoạt từ 30 - 50 triệu VNĐ cho việc chuẩn hóa hệ thống đánh giá KPI và cải thiện quan hệ nội bộ mà không cần tăng quỹ lương ngay lập tức.

4. Bảng câu hỏi khảo sát có thể tái sử dụng cho các ngành nghề khác không?

Có thể tái sử dụng khung cấu trúc 32 biến quan sát. Tuy nhiên, đối với các ngành dịch vụ hoặc công nghệ cao, cần thực hiện lại vòng nghiên cứu định tính (Focus Group) để hiệu chỉnh các biến đo lường về môi trường làm việc và công cụ lao động cho tương thích.

Kết luận

Đề tài đã hoàn thành xuất sắc việc định lượng hóa các nhân tố vô hình cấu thành nên động lực làm việc của người lao động tại Công ty TNHH Hiệp Thành. Kết quả nghiên cứu khẳng định: Chính sách đào tạo, thăng tiến kết hợp sự ghi nhận đóng góp ($\beta = 0.386$) và Mối quan hệ tại nơi làm việc ($\beta = 0.312$) là hai đòn bẩy quan trọng nhất chi phối hành vi nỗ lực của nhân viên, xếp trên cả nhân tố Tiền lương ($\beta = 0.241$).

Việc chuyển dịch từ cơ chế quản lý truyền thống sang quản trị nhân lực dựa trên dữ liệu thực chứng (Data-driven HRM) chính là chìa khóa then chốt giúp các doanh nghiệp xây dựng tối ưu hóa năng suất, ổn định nhân lực và hiện thực hóa mục tiêu tăng trưởng bền vững trong kỷ nguyên số.